Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1150 1758 (N42E17 L=198)   378     56        0 (1129/1726-1166/1785) -51.3 E18-E14 N39-N44
      2  759 1714 (N27E39 L=176)   309     70        1 ( 729/1677- 783/1739) -23.2 E42-E35 N24-N28
      3 1657 1631 (N54W33 L=249)   633     58        0 (1631/1600-1687/1653)  39.3 W29-W40 N49-N58
      4  388  960 (S26E40 L=175)  1098    277        3 ( 353/ 840- 448/1086)  79.4 E41-E36 S33-S17
      5  251  834 (S37E54 L=161)   732    166        0 ( 196/ 772- 302/ 887) -21.0 E61-E47 S39-S33
      6 1706  654 (S06W55 L=270)   344     43        2 (1684/ 629-1727/ 678) -33.7 W53-W57 S07-S05
      7  317  603 (S51E50 L=165)   463     69        0 ( 297/ 581- 342/ 636) -31.0 E53-E46 S52-S48
      8 1089  249 (S51W41 L=256)   812     99        0 (1052/ 212-1116/ 306) -52.1 W35-W45 S54-S45
  

    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