Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  857 1688 (N51E12 L= 55)   390     61        1 ( 829/1662- 890/1713)  31.0 E14-E09 N48-N52
      2 1211 1472 (N32W16 L= 83)   536    213        0 (1115/1450-1310/1499) -12.1 W08-W22 N30-N34
      3  237 1429 (N31E72 L=356)   875    116        1 ( 222/1388- 256/1491) -37.9 E78-E66 N27-N34
      4 1706 1364 (N21W54 L=121)   630    115        0 (1674/1310-1729/1412) -90.0 W51-W55 N17-N24
      5 1789 1303 (N17W61 L=128)  1915    254        4 (1742/1223-1847/1412) -71.6 W57-W66 N11-N23
      6  815 1363 (N26E13 L= 54)   225     90        0 ( 778/1335- 855/1380)  21.8 E15-E10 N24-N27
      7  852 1306 (N23E10 L= 57)   208     65        2 ( 823/1277- 880/1331) -26.6 E12-E08 N20-N24
      8  252  812 (S08E58 L=  9)   412     79        0 ( 240/ 773- 268/ 849)  68.2 E59-E56 S10-S05
      9 1518  744 (S16W32 L= 99)  1328    293        2 (1435/ 629-1587/ 846)  52.4 W26-W36 S23-S09
     10 1813  794 (S15W59 L=127)   393     38        1 (1793/ 770-1834/ 820)  56.3 W56-W61 S16-S13
     11  783  737 (S13E17 L= 51)   361     97        0 ( 748/ 695- 822/ 769)  33.7 E19-E13 S15-S10
     12  609  665 (S17E29 L= 38)   358     71        1 ( 587/ 641- 647/ 698) -33.7 E31-E26 S18-S14
     13 1688  535 (S32W53 L=121)  2536    369        1 (1598/ 431-1773/ 665)  60.6 W47-W60 S39-S22
     14  545  535 (S26E37 L= 30)   623    120        1 ( 528/ 477- 573/ 597)  60.9 E40-E34 S30-S21
     15  957  541 (S26E07 L= 61)   312    113        0 ( 906/ 528-1010/ 555)   0.0 E10-E02 S27-S25
     16 1501  389 (S42W39 L=107)   599     52        0 (1483/ 354-1521/ 423) -71.6 W36-W42 S45-S39
     17 1085  162 (S64W01 L= 68)   521     66        0 (1063/ 144-1114/ 182)  45.0 E01-W04 S66-S60
  

    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