Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  654 1739 (N26E49 L= 31)   499     80        1 ( 637/1705- 674/1786) -26.6 E55-E45 N23-N29
      2 1051 1643 (N30E17 L= 63)   553    158        0 ( 991/1569-1100/1695) -22.6 E23-E11 N26-N32
      3  719 1626 (N20E37 L= 43)   907    215        2 ( 639/1555- 796/1680) -13.2 E45-E28 N16-N22
      4  155 1137 (S19E69 L= 11)  1097    100        1 ( 102/1098- 184/1171)  15.9 E79-E64 S22-S17
      5  152 1034 (S26E68 L= 12)   936    114        1 ( 108/1007- 197/1056)   0.0 E77-E60 S27-S24
      6 1722  655 (S06W57 L=137)   753     82        2 (1694/ 617-1756/ 699) -59.0 W54-W59 S08-S03
      7 1402  433 (S28W45 L=125)   708    269        0 (1280/ 399-1496/ 497)  58.4 W39-W49 S34-S21
      8 1643  474 (S17W62 L=142)   498     96        0 (1611/ 453-1676/ 490)  45.0 W59-W65 S19-S15
      9 1557  365 (S26W64 L=144)  3524    432        3 (1355/ 277-1682/ 430)  53.1 W53-W77 S37-S17
  

    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