Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1004 1687 (N42E16 L=161)   314     74        0 ( 966/1670-1037/1701)   0.0 E19-E12 N41-N43
      2  268 1446 (N13E66 L=111)  3493    403        6 ( 161/1307- 380/1565)  46.5 E76-E57 N03-N22
      3  293 1524 (N18E67 L=110)   663     83        1 ( 270/1496- 311/1551)  63.4 E69-E63 N16-N20
      4  390 1360 (N10E49 L=127)   641    125        1 ( 365/1319- 435/1397)  59.0 E51-E45 N07-N12
      5  875 1301 (N14E14 L=163)   197     29        1 ( 854/1281- 896/1320)  56.3 E14-E12 N12-N15
      6  389 1024 (S10E42 L=135)  1664    317        6 ( 353/ 901- 453/1166) -65.0 E45-E35 S16-S01
      7  202 1120 (S07E62 L=115)   402    130        1 ( 171/1086- 235/1164)  35.5 E65-E58 S09-S04
      8  826 1111 (N02E13 L=163)   313     73        0 ( 788/1089- 864/1138) -11.3 E15-E10 N00-N02
      9 1037 1029 (N00W01 L=177)   439    140        0 ( 978/1004-1104/1066)  33.7 E02-W04 S02-N03
     10 1798  732 (S06W62 L=238)   970    164        2 (1767/ 661-1819/ 802) -53.1 W57-W65 S09-S01
     11  250  737 (S31E54 L=122)   545    123        1 ( 214/ 714- 282/ 764)  20.6 E58-E50 S32-S28
     12 1449  598 (S19W36 L=212)   338    145        1 (1377/ 566-1512/ 649)  45.0 W31-W38 S22-S15
     13 1385  414 (S32W39 L=215)  1052    281        0 (1276/ 384-1497/ 462)   3.0 W28-W47 S34-S29
  

    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