Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1373 1770 (N36W56 L=148)  1466    143        4 (1304/1730-1440/1799) -11.3 W48-W65 N34-N38
      2  695 1724 (N60W01 L= 92)   308     69        0 ( 657/1711- 722/1739)   0.0 E04-W04 N58-N61
      3 1706 1221 (S03W47 L=139)   291     69        1 (1687/1188-1718/1253) -90.0 W46-W48 S05-S01
      4  247 1109 (N29E53 L= 38)  1788    229        6 ( 140/1047- 343/1149) -13.5 E68-E43 N24-N31
      5 1888 1062 (S18W62 L=153)  1144    130        2 (1856/1006-1934/1109)  11.3 W56-W71 S21-S16
      6 1743  755 (S31W41 L=132)  2330    322        5 (1703/ 633-1775/ 915)  77.5 W36-W44 S39-S20
      7  221  874 (N16E59 L= 32)   466    128        1 ( 174/ 836- 277/ 902) -31.0 E65-E54 N12-N18
      8  358  827 (N10E47 L= 44)   646    151        1 ( 334/ 762- 390/ 891)  66.8 E48-E43 N06-N13
      9  235  769 (N09E61 L= 30)  1070    114        2 ( 209/ 708- 257/ 821)  45.0 E64-E57 N05-N12
     10 1517  481 (S42W17 L=108)   416     68        0 (1484/ 462-1550/ 505)  -9.5 W14-W20 S43-S40
     11  317  473 (S10E69 L= 22)   836     59        1 ( 282/ 451- 340/ 500)   0.0 E76-E64 S11-S08
     12  697  214 (S35E56 L= 35)  1554    127        0 ( 675/ 142- 724/ 274)  25.6 E74-E49 S42-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