Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  907 1666 (N46E12 L=129)  2972    588       14 ( 773/1584-1072/1724)  12.3 E23-W02 N39-N51
      2 1440 1655 (N46W36 L=177)   917    197        1 (1372/1631-1520/1679)   3.8 W28-W43 N44-N48
      3 1257 1628 (N44W18 L=159)   482    103        2 (1224/1587-1296/1665)  50.2 W15-W22 N40-N47
      4 1635 1541 (N37W50 L=190)   444     95        0 (1607/1525-1661/1560)   9.5 W46-W52 N35-N38
      5  422 1459 (N29E47 L= 94)   314    123        1 ( 410/1403- 437/1510) -82.9 E48-E45 N24-N32
      6 1281 1423 (N29W17 L=158)   216     58        2 (1254/1397-1303/1450) -33.7 W15-W19 N27-N30
      7  949 1305 (N21E06 L=135)   438     95        2 ( 907/1279-1010/1334) -15.9 E08-E01 N18-N22
      8  899 1228 (N16E08 L=132)   689    249        5 ( 799/1201-1008/1263)  13.0 E14-E01 N13-N17
      9  463 1190 (N12E38 L=103)   399    177        3 ( 426/1115- 509/1260) -48.8 E41-E33 N07-N15
     10  823 1141 (N10E13 L=128)   478    194        1 ( 767/1088- 925/1203) -31.0 E16-E06 N06-N13
     11 1546 1064 (N06W34 L=175)   361    112        0 (1506/1023-1596/1103) -39.8 W31-W37 N03-N08
     12 1237  903 (S04W13 L=154)   261     87        0 (1220/ 851-1256/ 947)  80.5 W12-W14 S07-S01
     13  249  866 (S10E57 L= 84)   497    165        2 ( 208/ 836- 302/ 909) -18.4 E61-E51 S11-S07
     14 1548  781 (S12W35 L=176)   371    107        2 (1517/ 734-1594/ 808)  45.0 W33-W38 S14-S09
     15  311  773 (S15E51 L= 90)   500     70        2 ( 280/ 744- 338/ 798)  -9.5 E54-E48 S16-S13
     16  464  403 (S41E47 L= 94)  1197    296        0 ( 406/ 364- 503/ 463)  32.0 E56-E40 S45-S35
     17 1290  378 (S40W23 L=164)  1307    370        3 (1143/ 328-1450/ 426)  11.8 W11-W36 S44-S36
  

    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