Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  669 1828 (N40E59 L= 60)  1467    121        2 ( 636/1773- 708/1882) -12.3 E80-E48 N35-N43
      2  891 1665 (N33E27 L= 92)   346     86        1 ( 868/1625- 920/1700)  90.0 E28-E24 N29-N36
      3  359 1557 (N14E63 L= 56)  3467    386        5 ( 252/1466- 460/1642)  51.8 E70-E54 N07-N21
      4  996 1091 (N00E03 L=116)  1102    305        1 ( 954/ 939-1047/1216) -52.4 E07-W03 S07-N07
      5 1797 1082 (N19W51 L=170)  1692    283        1 (1712/ 970-1865/1186) -26.6 W41-W61 N13-N23
      6  476 1109 (S10E35 L= 83)   430    134        0 ( 418/1093- 539/1136)  26.6 E39-E31 S11-S06
      7 1588  608 (S13W47 L=166)  1531    305        1 (1460/ 544-1723/ 666)  40.9 W39-W56 S19-S06
      8 1581  410 (S25W60 L=178)  1807    175        4 (1496/ 385-1675/ 436)  19.4 W51-W70 S27-S20
      9  935  375 (S45W13 L=132)   577     80        1 ( 891/ 353- 978/ 407)   7.1 W09-W17 S46-S42
  

    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