Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1516 1717 (N57W44 L=163)  1098    119        1 (1469/1673-1568/1748)  17.2 W33-W62 N51-N60
      2 1551 1613 (N47W41 L=159)   436     60        0 (1517/1594-1579/1630)  -8.1 W36-W43 N45-N48
      3  864 1404 (N23E15 L=103)   711    177        1 ( 830/1340- 920/1485) -60.9 E18-E11 N18-N27
      4  751 1333 (N17E21 L= 97)   525    145        1 ( 713/1262- 780/1405)  90.0 E22-E19 N12-N21
      5  665  785 (S18E20 L= 98)   574    133        0 ( 614/ 737- 720/ 813)  40.6 E23-E16 S21-S15
      6 1207  665 (S19W16 L=134)   224     62        2 (1171/ 645-1236/ 684) -11.3 W13-W18 S20-S17
      7 1326  405 (S36W32 L=151)   416     95        0 (1291/ 380-1372/ 440)  51.3 W29-W34 S37-S32
  

    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