Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1047 1766 (N42E16 L= 45)   303     46        0 (1014/1755-1080/1778)  31.0 E18-E13 N40-N43
      2  178 1296 (S02E68 L=352)  2048    322        5 ( 125/1225- 217/1368)  45.0 E73-E62 S06-N02
      3 1823 1325 (N29W51 L=111)   476     66        1 (1801/1300-1848/1347)   0.0 W48-W53 N27-N30
      4  143 1164 (S10E68 L=352)   928    214        1 ( 113/1126- 188/1212) -22.6 E73-E61 S12-S06
      5 1838  993 (N09W56 L=116)  1876    302        4 (1805/ 847-1858/1134) -58.0 W50-W62 N00-N16
      6  655  993 (S15E22 L= 39)   192     77        0 ( 638/ 958- 682/1033) -53.1 E22-E19 S16-S12
      7 1442  599 (S23W36 L= 96)   778    240        1 (1314/ 565-1538/ 651)  39.8 W28-W40 S27-S17
      8 1219  575 (S29W22 L= 82)   273     50        0 (1193/ 561-1257/ 597)   0.0 W19-W24 S30-S28
      9  995  519 (S37W09 L= 69)   286     63        2 ( 960/ 500-1036/ 538)  36.9 W06-W11 S39-S35
     10 1308  273 (S46W44 L=105)   760    116        3 (1254/ 247-1352/ 303)  23.2 W40-W48 S47-S44
     11  860  212 (S66W12 L= 72)   585     94        1 ( 839/ 186- 892/ 241) -14.0 W08-W16 S68-S62
  

    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