Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1216 1568 (N30W24 L=138)   280     69        0 (1184/1539-1249/1594) -68.2 W23-W25 N27-N32
      2 1316 1465 (N22W27 L=141)   345     62        1 (1288/1435-1348/1486) -45.0 W25-W29 N19-N23
      3  432 1354 (N32E35 L= 79)   326    123        0 ( 386/1320- 476/1372)   7.1 E39-E31 N30-N33
      4 1884 1016 (S15W63 L=177)  1296    223        2 (1811/ 960-1926/1083)  11.3 W53-W69 S17-S10
      5 1339  954 (S08W17 L=131)   384    108        1 (1316/ 905-1357/1010) -81.9 W15-W18 S11-S04
      6  177  819 (N04E66 L= 48)  1184    123        2 ( 121/ 788- 218/ 869)  10.0 E76-E59 N02-N06
      7  862  748 (S12E15 L= 99)   890    257        0 ( 801/ 634- 927/ 844)  39.8 E20-E08 S17-S06
      8 1830  764 (S30W57 L=171)   564     58        2 (1801/ 738-1854/ 793)  26.6 W52-W60 S31-S28
      9  339  699 (S06E52 L= 62)   548     93        1 ( 310/ 662- 378/ 732)  33.7 E56-E48 S07-S03
     10 1747  646 (S36W48 L=162)  1582    250        1 (1674/ 580-1851/ 727) -23.6 W39-W71 S43-S28
     11  359  570 (S14E56 L= 58)  1219    138        1 ( 319/ 527- 401/ 638)  22.6 E62-E50 S16-S09
  

    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