Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1202 1761 (N55W04 L=197)   537    129        2 (1149/1735-1244/1787)  29.1 E01-W08 N51-N58
      2 1367 1707 (N53W21 L=214)   556    155        1 (1303/1691-1424/1722)   0.0 W14-W27 N50-N54
      3  463 1582 (N29E51 L=142)   577    137        1 ( 438/1523- 493/1634) -69.4 E52-E49 N24-N32
      4  339 1460 (N19E57 L=136)   442     96        5 ( 307/1438- 373/1482)  18.4 E60-E53 N17-N20
      5  741 1362 (N18E23 L=171)   462    207        1 ( 715/1267- 767/1461) -67.4 E25-E19 N12-N24
      6 1632 1405 (N32W40 L=234)   402     89        2 (1609/1373-1655/1433) -45.0 W38-W42 N29-N34
      7 1715 1382 (N32W48 L=241)   618     70        2 (1691/1352-1746/1411) -36.9 W45-W50 N29-N33
      8 1553 1057 (N09W34 L=227)   270     65        0 (1527/1028-1578/1084) -26.6 W31-W35 N07-N10
      9 1808 1028 (N10W56 L=250)   671    152        4 (1778/ 946-1832/1076)  90.0 W53-W59 N04-N12
     10  295  871 (S18E50 L=144)  1237    301        5 ( 220/ 782- 352/ 973) -45.0 E57-E44 S22-S12
     11  182  877 (S19E62 L=131)   541    111        2 ( 166/ 842- 196/ 911) -76.0 E64-E60 S20-S16
     12  235  800 (S23E56 L=137)   403    120        3 ( 221/ 739- 246/ 846) -90.0 E57-E54 S26-S19
     13  418  763 (S23E39 L=155)   405    113        2 ( 394/ 717- 442/ 828) -90.0 E40-E36 S25-S18
     14 1758  680 (S12W59 L=252)  3990    647       11 (1670/ 568-1847/ 812)  76.0 W52-W67 S19-S03
     15 1108  671 (S20W10 L=203)   330    123        1 (1041/ 654-1175/ 689)  14.0 W05-W13 S20-S18
     16 1424  562 (S23W34 L=227)   857    171        2 (1357/ 500-1499/ 621)  55.0 W30-W37 S28-S18
     17  802  488 (S37E09 L=184)   430     54        1 ( 772/ 462- 835/ 516)  38.7 E11-E06 S38-S34
  

    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