Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1159 1498 (N33W19 L= 60)   632    174        3 (1105/1418-1216/1536)  31.0 W13-W24 N28-N34
      2 1174 1223 (N16W13 L= 53)   319    130        1 (1135/1197-1224/1257) -31.0 W10-W15 N13-N18
      3 1348 1046 (N02W20 L= 60)   264     77        0 (1328/1005-1364/1080)  56.3 W18-W21 N00-N03
      4 1854  985 (S14W58 L= 99)  1130    240        2 (1802/ 930-1892/1062)  36.9 W52-W62 S17-S09
      5  230  926 (N12E59 L=342)   388     54        0 ( 211/ 898- 252/ 961) -76.0 E60-E56 N09-N13
      6 1849  896 (S19W57 L= 97)   461    121       -1 (1821/ 858-1870/ 939)  26.6 W53-W59 S21-S17
      7  244  752 (N01E62 L=339)  1048    302        3 ( 183/ 722- 324/ 785)   3.2 E70-E52 S01-N03
      8 1720  657 (S31W43 L= 84)   879    209        1 (1667/ 589-1761/ 758)  48.4 W38-W47 S35-S24
      9  287  578 (S11E64 L=336)  2030    310        8 ( 224/ 495- 334/ 657)  84.8 E71-E58 S16-S05
     10 1494  451 (S41W23 L= 64)  1231    257        1 (1386/ 368-1605/ 546)  26.6 W13-W33 S45-S35
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y