Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  914 1766 (N37E33 L= 16)  2609    547        2 ( 880/1646- 976/1923) -35.9 E51-E22 N28-N50
      2 1366 1399 (N26W10 L= 59)   334    110        0 (1336/1338-1382/1444) -56.3 W06-W11 N21-N28
      3 1768 1398 (N41W43 L= 93)   594     82        1 (1736/1365-1794/1417)  15.9 W39-W46 N38-N43
      4 1844 1341 (N40W52 L=101)   874     53        0 (1823/1308-1868/1374) -36.9 W49-W54 N38-N42
      5  793 1334 (N06E21 L= 28)   400     51        0 ( 767/1300- 820/1364) -71.6 E22-E19 N04-N07
      6 1859 1256 (N35W53 L=103)   642    118        0 (1825/1231-1880/1279)  18.4 W49-W55 N32-N36
      7 1062 1147 (N02E01 L= 49)   201     90        0 (1015/1137-1114/1156)  26.6 E03-W01 N01-N03
      8 1147 1135 (N04W04 L= 54)   194     13        0 (1128/1114-1164/1156)  90.0 W02-W04 N02-N05
      9 1677 1031 (N14W39 L= 89)   774    215        0 (1641/ 956-1731/1079)  82.9 W36-W43 N09-N17
     10 1588 1024 (N11W33 L= 83)   489    115        1 (1556/ 973-1618/1076) -51.3 W29-W35 N07-N13
     11  851  940 (S14E07 L= 42)   254     54        2 ( 836/ 900- 871/ 974) -76.0 E08-E06 S16-S12
     12 1703  886 (N07W46 L= 95)   741    130        0 (1650/ 836-1740/ 942)  69.4 W42-W48 N02-N10
     13  390  891 (S29E37 L= 12)   706    185        1 ( 323/ 841- 451/ 938)  37.9 E42-E32 S33-S25
     14 1836  877 (N11W58 L=108)   781    104        1 (1811/ 828-1865/ 924)  74.1 W56-W60 N07-N14
     15  699  848 (S24E14 L= 36)   634    213        1 ( 602/ 816- 797/ 893)  39.3 E20-E09 S27-S18
     16  874  696 (S28W01 L= 50)   442     75        1 ( 833/ 672- 911/ 728)  59.0 E01-W02 S30-S25
     17  652  620 (S39E12 L= 38)   427     66        1 ( 617/ 603- 699/ 646)   8.1 E15-E08 S40-S37
     18 1191  293 (S42W40 L= 90)   355     60        0 (1159/ 274-1221/ 314)   9.5 W37-W43 S43-S41
  

    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