Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1479 1558 (N32W44 L=319)  1594    348        2 (1324/1493-1627/1597) -27.6 W31-W54 N24-N36
      2  551 1453 (N38E31 L=244)   314     66        1 ( 518/1431- 585/1478) -29.7 E34-E27 N36-N40
      3  334 1375 (N35E52 L=223)  1533    281        2 ( 223/1314- 429/1427) -12.3 E72-E40 N29-N38
      4  961 1347 (N26W00 L=275)   457    215        1 ( 866/1334-1068/1357) -13.0 E06-W07 N24-N27
      5 1587 1290 (N14W42 L=317)   470    111        1 (1562/1238-1611/1346)  71.6 W39-W43 N10-N16
      6 1257 1159 (N11W16 L=291)   307    113        2 (1188/1137-1302/1179) -26.6 W12-W18 N09-N13
      7 1784  948 (S10W52 L=327)   828    197        1 (1717/ 913-1867/ 972) -17.4 W45-W62 S14-S07
      8 1681  870 (S13W42 L=317)   820    198        3 (1642/ 795-1712/ 956)  56.3 W38-W45 S17-S07
      9 1750  808 (S18W49 L=324)  1149    255        3 (1721/ 727-1797/ 911)  85.2 W45-W54 S23-S11
     10  385  739 (S05E47 L=228)  1114    231        5 ( 327/ 655- 432/ 824) -84.3 E51-E43 S10-N00
     11 1323  746 (S15W15 L=290)   324     72        0 (1290/ 732-1364/ 765) -36.9 W13-W17 S16-S13
     12  285  653 (S11E60 L=215)  2316    339        7 ( 206/ 546- 401/ 752) -40.9 E67-E51 S17-S04
     13 1234  617 (S22W08 L=283)  1571    396        4 (1190/ 516-1334/ 752)  86.4 W05-W14 S29-S13
     14  335  695 (S08E53 L=222)   355     36        2 ( 316/ 669- 353/ 717)  45.0 E54-E51 S09-S06
     15  307  591 (S15E61 L=214)   410     73        1 ( 276/ 571- 336/ 616)   7.1 E64-E56 S16-S13
     16  368  481 (S23E60 L=215)  1018     95        2 ( 339/ 440- 393/ 527)  24.8 E68-E54 S25-S19
     17  983  441 (S31E11 L=264)   322    127        1 ( 929/ 395-1022/ 486)  32.0 E16-E08 S33-S27
     18 1404  384 (S43W21 L=296)   863    101        1 (1352/ 355-1446/ 426)  20.6 W16-W24 S45-S40
  

    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