Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1393 1738 (N53W41 L=207)   548     64        0 (1357/1714-1423/1765)  14.0 W35-W47 N50-N55
      2 1037 1736 (N54W02 L=167)   359     49        0 (1015/1709-1062/1759)  45.0 W00-W04 N51-N55
      3  933 1730 (N53E09 L=156)   407     46        0 ( 904/1714- 960/1747)  -9.5 E12-E06 N51-N54
      4  963 1491 (N34E04 L=161)   253     30        0 ( 943/1471- 991/1507)  18.4 E05-E02 N32-N35
      5  473 1259 (N18E38 L=127)   496    205        0 ( 412/1199- 545/1291)  29.1 E42-E32 N13-N19
      6  213 1065 (N05E61 L=104)   463     55        1 ( 187/1039- 243/1096)  23.2 E64-E57 N02-N06
      7  176  971 (S02E66 L= 99)   819     78        2 ( 143/ 948- 210/ 997)  11.3 E71-E61 S02-N00
      8  513  564 (S27E38 L=127)  1661    429        1 ( 399/ 414- 652/ 717)  33.3 E59-E24 S39-S15
      9 1385  682 (S18W24 L=189)   287    125        0 (1316/ 667-1442/ 699)   6.3 W19-W28 S19-S17
     10  801  489 (S31E17 L=149)  3029    591        6 ( 599/ 346- 996/ 676)  34.8 E37-E01 S43-S18
     11 1606  512 (S31W47 L=212)   475     65        0 (1574/ 487-1633/ 542)  36.9 W43-W49 S32-S28
     12 1527  321 (S47W53 L=218)  1460     93        4 (1489/ 284-1567/ 365) -19.1 W44-W70 S52-S42
  

    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