Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  938 1821 (N44E27 L= 36)   323    104        0 ( 911/1789- 960/1847) -18.4 E30-E24 N41-N46
      2 1023 1772 (N42E18 L= 45)   341     59        1 ( 989/1752-1058/1786)  26.6 E19-E15 N39-N43
      3 1805 1332 (N29W49 L=112)   547    123        1 (1770/1306-1843/1350)   7.1 W45-W53 N27-N30
      4  180 1287 (S02E67 L=356)   623    149        3 ( 157/1240- 197/1337) -59.0 E71-E64 S05-N00
      5 1812 1009 (N09W53 L=116)  1892    325        4 (1779/ 845-1837/1179) -63.4 W47-W57 N00-N18
      6 1434  617 (S22W35 L= 98)   621    140        0 (1364/ 585-1506/ 662)  53.1 W30-W38 S25-S17
      7 1257  285 (S47W40 L=103)  1166    213        3 (1158/ 253-1342/ 322)  14.9 W30-W47 S49-S44
      8  863  208 (S66W13 L= 76)   667     80        0 ( 838/ 181- 889/ 238) -20.0 W07-W18 S68-S63
  

    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