Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  434 1658 (N50E48 L= 91)  1171    143        2 ( 394/1603- 476/1726) -23.4 E70-E39 N44-N57
      2 1154 1595 (N34W17 L=157)   486    101        1 (1111/1562-1199/1634) -51.3 W15-W19 N30-N36
      3  578 1311 (N22E26 L=114)   584    146        0 ( 508/1285- 643/1333)   5.7 E30-E20 N21-N23
      4  348 1259 (N22E44 L= 95)   264     42        2 ( 329/1240- 373/1278) -21.8 E46-E41 N20-N23
      5 1728  909 (S15W46 L=186)  1528    295        1 (1694/ 773-1753/1050)  90.0 W43-W48 S23-S06
      6  247  710 (S10E61 L= 79)   711    100        0 ( 221/ 660- 277/ 755)  35.5 E64-E56 S12-S07
      7 1756  574 (S38W57 L=196)  1632    212        2 (1678/ 509-1843/ 648)  24.0 W48-W67 S42-S33
  

    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