Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1219 1842 (N58W47 L=119)   401     55        0 (1203/1813-1233/1870)  29.1 W41-W52 N55-N60
      2  885 1740 (N57W02 L= 74)   297     51        0 ( 867/1719- 915/1758)  21.8 W00-W05 N55-N58
      3  727 1723 (N58E16 L= 56)   743    186        1 ( 638/1697- 811/1746)   0.0 E27-E06 N55-N60
      4 1626 1633 (N32W65 L=137)   730     93        1 (1596/1591-1650/1681)  14.0 W58-W74 N30-N34
      5 1112 1511 (N35W14 L= 86)   221     62        0 (1082/1495-1154/1525)  -9.5 W11-W17 N34-N35
      6  895 1374 (N29E04 L= 68)   484    179        1 ( 835/1301- 929/1425) -60.9 E07-E01 N24-N33
      7  678 1338 (N30E19 L= 53)   320    129        0 ( 642/1302- 734/1386) -45.0 E21-E15 N26-N33
      8  153 1200 (N24E72 L=  0)  1028    108        0 ( 139/1154- 172/1262) -37.9 E77-E66 N21-N28
      9  278 1158 (N22E53 L= 19)   605     86        2 ( 261/1117- 295/1203) -78.7 E55-E51 N19-N24
     10 1841  971 (S11W58 L=130)   822    143        2 (1805/ 914-1868/1019) -66.8 W54-W61 S14-S07
     11 1799  922 (S13W53 L=125)   525    122        2 (1777/ 869-1817/ 980)  63.4 W50-W55 S15-S09
     12  656  862 (N01E25 L= 47)   357    100        0 ( 617/ 828- 703/ 893) -45.0 E26-E22 S01-N03
     13 1733  640 (S30W48 L=120)  2608    469        7 (1534/ 492-1873/ 749)  12.5 W29-W67 S36-S25
     14 1298  386 (S40W10 L= 83)   526    111        0 (1248/ 363-1348/ 414)  -7.1 W06-W14 S41-S38
     15 1135  330 (S42E03 L= 69)   438     71        0 (1105/ 312-1176/ 353)  11.3 E05-W00 S43-S39
     16  957  245 (S47E21 L= 51)  1202    151        2 ( 885/ 210-1025/ 282)  11.3 E29-E14 S49-S43
     17  800  182 (S51E42 L= 30)   717     57        1 ( 782/ 149- 818/ 214)  29.7 E51-E37 S55-S47
  

    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