Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  365 1279 (N21E47 L=356)   962    183        2 ( 281/1251- 441/1316)   0.0 E56-E40 N18-N23
      2 1820 1268 (N14W62 L=106)   925     91        2 (1782/1231-1854/1306) -20.6 W58-W66 N12-N16
      3 1337 1038 (N04W19 L= 63)   303    137        0 (1288/ 997-1396/1079)  35.5 W16-W23 N01-N06
      4  447  972 (N02E38 L=  6)   611    189        0 ( 376/ 920- 528/1023) -24.4 E43-E32 S01-N04
      5 1423  947 (S02W25 L= 68)   291     81        0 (1407/ 906-1444/ 989) -63.4 W23-W26 S04-N00
      6  996  913 (S03E02 L= 42)   331     84        1 ( 964/ 879-1025/ 959)  53.1 E04-W00 S04-N00
      7 1835  930 (S06W59 L=103)   832     97        1 (1780/ 913-1884/ 947)  -8.7 W53-W66 S07-S04
      8  633  845 (S06E25 L= 18)   542    178        0 ( 566/ 813- 677/ 907)  32.0 E30-E22 S07-S02
      9 1326  781 (S12W18 L= 62)  1048    422        0 (1153/ 683-1452/ 896)  32.3 W07-W26 S17-S05
     10  317  710 (S15E53 L=351)   742    201        1 ( 270/ 624- 357/ 790) -74.7 E55-E49 S20-S09
     11  401  649 (S19E46 L=358)   301     60        1 ( 385/ 619- 426/ 672)  31.0 E48-E43 S20-S16
  

    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