Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1012 1875 (N60W33 L=225)  1758    244        2 ( 933/1809-1120/1939)   5.7 W19-W59 N55-N62
      2 1228 1742 (N42W37 L=228)   296     58        0 (1206/1716-1255/1763)  15.9 W33-W40 N39-N42
      3 1325 1735 (N38W44 L=235)   334     49        0 (1296/1718-1350/1754) -26.6 W40-W46 N37-N39
      4  255 1170 (N25E52 L=139)   732    162        1 ( 240/1095- 276/1252) -68.2 E54-E49 N20-N30
      5 1677 1025 (S11W41 L=232)   320     84        0 (1651/ 994-1704/1067) -59.0 W39-W42 S13-S08
      6 1136  670 (S21E01 L=191)   548    259        0 (1038/ 599-1244/ 735)  11.3 E08-W07 S23-S17
      7  252  661 (S04E63 L=128)   454     60        0 ( 235/ 628- 269/ 690)  31.0 E65-E59 S06-S02
      8  890  356 (S36E26 L=165)  1167    317        2 ( 749/ 338-1050/ 387) -21.8 E35-E14 S39-S30
  

    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