Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  916 1884 (N64W33 L=358)  1901    223        4 ( 840/1812- 994/1959)   6.6 W18-W63 N58-N67
      2 1184 1737 (N41W36 L=  1)  2175    395        3 (1004/1696-1368/1779) -25.6 W21-W47 N34-N46
      3 1648 1620 (N21W63 L= 29)   597     40        0 (1618/1598-1672/1645)  -7.1 W59-W67 N19-N21
      4  788 1590 (N43E00 L=326)   348     95        0 ( 765/1546- 815/1627)  59.0 E02-W03 N40-N45
      5 1742 1489 (N11W63 L= 28)   616    102        0 (1709/1444-1767/1531) -90.0 W61-W64 N07-N13
      6 1557 1097 (S05W33 L=358)  1991    449        3 (1504/ 889-1587/1296)  60.3 W25-W38 S15-N06
      7 1875 1137 (S12W61 L= 27)   598    105        1 (1847/1078-1902/1171) -80.5 W58-W65 S14-S08
      8  723 1081 (N14E16 L=310)   871    226        1 ( 641/1027- 820/1160) -51.3 E19-E10 N08-N19
      9  130 1106 (N27E68 L=257)  1656    124        4 (  81/1051- 170/1156) -20.2 E82-E61 N23-N30
     10  246 1098 (N25E52 L=273)   498     72        2 ( 220/1062- 271/1138) -59.0 E55-E49 N22-N27
     11  334 1011 (N18E45 L=281)  1176    250        3 ( 237/ 922- 390/1095)  23.6 E55-E39 N14-N22
     12  325  864 (N10E48 L=277)   356     37        2 ( 297/ 842- 346/ 885) -33.7 E50-E46 N08-N10
     13  218  743 (N05E63 L=263)   438     55        1 ( 191/ 721- 245/ 765) -36.9 E65-E59 N03-N06
     14 1878  681 (S40W63 L= 29)   671     73        0 (1846/ 650-1906/ 698) -18.4 W56-W77 S44-S37
     15 1242  592 (S26W02 L=328)   538    195        0 (1184/ 528-1314/ 663)  26.6 E02-W08 S28-S23
     16 1745  568 (S44W45 L= 10)  1044    146        1 (1681/ 522-1809/ 614)   7.1 W36-W53 S45-S40
     17  345  524 (S10E60 L=266)   382     35        0 ( 319/ 510- 377/ 540)  -8.1 E63-E56 S10-S08
     18  726  249 (S36E48 L=278)   822    122        0 ( 677/ 213- 761/ 301)  12.1 E54-E40 S38-S32
  

    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