Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1350  828 (S16W19 L=310)  1016    283        0 (1258/ 710-1418/ 954) -59.0 W13-W24 S23-S07
      2 1761  856 (S16W51 L=341)   610    113        1 (1717/ 810-1810/ 898)  26.6 W47-W55 S18-S13
      3 1938  837 (S17W76 L=  7)   918    102        1 (1911/ 797-1961/ 878)  29.1 W71-W83 S19-S14
      4  438  827 (S10E40 L=251)   255     38        0 ( 409/ 810- 463/ 849) -36.9 E41-E37 S11-S08
      5  729  775 (S15E20 L=270)   285     57        1 ( 701/ 744- 754/ 802)  45.0 E21-E18 S16-S13
      6  232  703 (S16E63 L=228)   586     89        3 ( 205/ 662- 255/ 738)  59.0 E66-E59 S18-S13
      7  329  660 (S19E53 L=237)   519     97        1 ( 307/ 611- 353/ 691)  68.2 E55-E50 S22-S17
      8 1151  564 (S31W06 L=296)   507    119        1 (1107/ 538-1205/ 592)   0.0 W02-W09 S32-S29
      9  919  290 (S51E17 L=274)   442     65        1 ( 890/ 272- 958/ 312)  15.9 E19-E12 S52-S48
  

    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