Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1031 1801 (N42E25 L=236)   350     61        0 (1008/1775-1064/1823)  90.0 E26-E23 N39-N45
      2  689 1762 (N30E46 L=214)  1579    183        5 ( 620/1706- 757/1817)  73.3 E50-E42 N25-N36
      3  837 1687 (N28E31 L=230)  1409    292        1 ( 754/1649- 946/1756)  22.6 E37-E23 N24-N32
      4  423 1572 (N12E55 L=205)   617    146        0 ( 387/1548- 481/1610)   5.7 E59-E49 N10-N14
      5 1547 1447 (N32W23 L=284)  1651    424        0 (1422/1300-1697/1602) -25.6 W11-W37 N26-N39
      6  500 1527 (N11E47 L=214)   451     88        0 ( 467/1502- 532/1548)  53.1 E49-E43 N08-N12
      7 1650 1466 (N37W32 L=293)   376     97        0 (1606/1431-1686/1500)  -7.1 W27-W35 N35-N38
      8  583 1408 (N05E36 L=225)   913    105        2 ( 557/1350- 611/1454)  90.0 E38-E33 N01-N08
      9 1808 1328 (N33W49 L=310)   516     48        0 (1777/1309-1834/1347)   0.0 W45-W53 N32-N35
     10  566 1147 (S10E30 L=231)   901    256        3 ( 475/1071- 671/1221)  50.2 E34-E24 S16-S04
     11 1218 1150 (N04W08 L=268)   401     89        0 (1172/1125-1261/1179)  45.0 W05-W09 N02-N07
     12  163 1000 (S25E66 L=195)  1965    170        4 ( 131/ 921- 189/1084)  78.7 E71-E61 S30-S20
     13  343  917 (S28E44 L=217)   739    213        3 ( 255/ 865- 419/ 972)  26.6 E53-E37 S31-S23
     14 1202  873 (S11W14 L=274)   415    147        0 (1114/ 850-1259/ 892)  18.4 W08-W17 S13-S09
     15  561  779 (S32E24 L=237)  2435    523        4 ( 366/ 650- 787/ 866)  30.5 E42-E08 S44-S24
     16 1601  501 (S20W55 L=316)   858    189        0 (1523/ 456-1681/ 559)  55.0 W51-W59 S24-S14
     17 1311  520 (S28W33 L=294)   547    182        0 (1230/ 502-1395/ 544)  14.0 W26-W38 S29-S25
     18 1489  353 (S31W58 L=319)  2534    287        4 (1352/ 325-1618/ 390)  23.5 W46-W69 S36-S25
  

    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