Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  868 1663 (N34E27 L=181)  2516    573        6 ( 655/1546-1058/1773)  -4.6 E48-E11 N29-N38
      2 1147 1545 (N33E04 L=204)   401     85        0 (1112/1513-1193/1568)  51.3 E06-E01 N29-N34
      3  480 1471 (N14E46 L=162)   884    183        3 ( 458/1381- 520/1545) -73.3 E49-E41 N08-N18
      4  612 1416 (N13E34 L=174)   308     77        1 ( 579/1385- 635/1453)  78.7 E35-E33 N10-N15
      5  361 1372 (N06E53 L=156)   516    135        1 ( 306/1357- 405/1388)  12.5 E57-E48 N03-N06
      6  257 1312 (N00E61 L=147)   479    103        2 ( 230/1289- 289/1335)   0.0 E64-E56 S01-N01
      7 1788  983 (N12W53 L=261)  1634    240        1 (1722/ 884-1822/1067) -36.9 W45-W58 N06-N15
      8  274 1002 (S18E51 L=157)   359    100        0 ( 239/ 984- 303/1017)  15.9 E54-E47 S19-S16
      9  739  986 (S10E16 L=192)   452    152        1 ( 661/ 964- 806/1007)  26.6 E20-E12 S11-S07
     10 1545  713 (S09W40 L=248)   311    107        0 (1532/ 675-1567/ 763) -38.7 W37-W42 S10-S06
     11  990  621 (S27W07 L=215)   369     56        0 ( 953/ 600-1025/ 643)  11.3 W04-W09 S27-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