Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1368 1754 (N55W02 L=156)   903     95        2 (1336/1703-1405/1789)  90.0 E01-W05 N48-N59
      2 1059 1710 (N39E20 L=134)   644    121        1 ( 996/1681-1109/1737)   5.7 E25-E15 N37-N40
      3  731 1542 (N19E33 L=122)  1168    223        3 ( 691/1472- 789/1655) -37.6 E39-E26 N14-N25
      4 1527 1612 (N49W21 L=175)   548     67        1 (1495/1592-1566/1636)  14.0 W17-W25 N46-N50
      5  306 1526 (N09E68 L= 86)   909    133        0 ( 263/1480- 334/1589) -82.9 E71-E65 N05-N13
      6  844 1365 (N11E19 L=135)   327     52        1 ( 812/1341- 870/1392)  56.3 E20-E17 N08-N12
      7 1780 1129 (N23W48 L=203)   577    160        1 (1724/1097-1823/1179) -18.4 W42-W53 N21-N25
      8  520  888 (S24E28 L=126)   319    104        0 ( 484/ 862- 563/ 924)  -8.1 E31-E24 S25-S22
      9 1283  744 (S13W23 L=177)   736    266        0 (1245/ 612-1320/ 850) -75.1 W20-W25 S21-S06
     10 1506  640 (S13W41 L=195)   390     58        2 (1474/ 616-1541/ 657)  53.1 W37-W43 S14-S10
     11  708  604 (S37E09 L=145)   308     71        1 ( 678/ 581- 739/ 629)  53.1 E11-E07 S38-S34
     12  843  346 (S51W11 L=165)   358     42        0 ( 819/ 323- 868/ 371)  78.7 W08-W12 S53-S48
  

    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