Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1151 1820 (N57E15 L= 36)  1664    234        1 (1051/1771-1248/1855)  49.1 E23-E07 N48-N64
      2 1769 1109 (N19W50 L=101)  2144    339        6 (1681/1015-1843/1184) -30.3 W41-W58 N12-N23
      3  236 1130 (S10E58 L=353)   849    159        1 ( 190/1061- 289/1181)  48.8 E62-E52 S15-S06
      4 1664 1110 (N17W40 L= 91)   900    203        1 (1598/1028-1714/1164) -24.4 W34-W45 N12-N19
      5 1907 1012 (N17W70 L=121)   812     84        1 (1888/ 983-1927/1066)  78.7 W66-W74 N14-N19
      6 1216  614 (S22W21 L= 72)   342    119        0 (1168/ 580-1276/ 654)  50.2 W18-W23 S24-S18
      7 1705  505 (S18W66 L=117)   702     95        0 (1673/ 463-1737/ 536)  68.2 W63-W68 S20-S15
      8 1481  441 (S27W48 L= 99)  1061    190        2 (1399/ 417-1559/ 476)  35.0 W42-W53 S29-S22
      9 1160  237 (S50W40 L= 91)   537     80        1 (1124/ 208-1197/ 260)   0.0 W33-W44 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