Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  623 1757 (N32E49 L=286)  1633    278        2 ( 558/1676- 668/1811) -35.5 E56-E40 N26-N36
      2  505 1720 (N28E57 L=277)   645     97        2 ( 468/1678- 547/1750) -80.5 E60-E55 N24-N30
      3  554 1679 (N25E49 L=286)  1043    167        1 ( 505/1626- 609/1711)  69.4 E53-E45 N20-N28
      4 1396 1623 (N39W13 L=348)   279     45        0 (1368/1608-1424/1637)  11.3 W10-W15 N38-N40
      5  576 1461 (N11E37 L=298)  1623    315        1 ( 505/1343- 670/1597)  86.8 E39-E35 N03-N21
      6  735 1456 (N13E26 L=309)   469    110        1 ( 709/1406- 764/1498)  90.0 E27-E24 N10-N16
      7  933 1339 (N10E11 L=324)   501    141        2 ( 887/1278- 978/1408)  90.0 E13-E09 N05-N14
      8 1841 1222 (N24W55 L= 29)   547     71        1 (1817/1184-1861/1261) -53.1 W51-W57 N21-N26
      9  659 1226 (S02E25 L=309)   210     79        0 ( 628/1193- 679/1251) -36.9 E27-E23 S03-N00
     10 1549 1208 (N15W28 L=  3)   470    145        0 (1504/1135-1574/1247) -50.2 W24-W30 N10-N16
     11 1376 1197 (N10W17 L=352)   261     51        1 (1361/1161-1394/1228) -56.3 W15-W18 N08-N11
     12  204 1080 (S16E60 L=275)  1269    188        3 ( 167/1019- 265/1148)  41.6 E64-E54 S19-S11
     13 1780 1016 (N10W50 L= 25)   465    114        0 (1765/ 944-1796/1061) -74.1 W48-W50 N05-N12
     14  579  803 (S28E25 L=310)  1493    531        1 ( 388/ 684- 691/ 991)  47.5 E41-E18 S38-S14
     15 1824  886 (N04W57 L= 32)   333     70        1 (1804/ 850-1839/ 917) -76.0 W54-W59 N01-N05
     16 1786  663 (S09W63 L= 38)   674     76        1 (1742/ 639-1819/ 689)  12.5 W57-W66 S10-S07
     17 1597  511 (S23W53 L= 28)  1695    290        3 (1534/ 396-1650/ 643) -56.3 W47-W58 S30-S15
  

    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