Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  368 1449 (N06E54 L=177)   971    188        1 ( 310/1395- 431/1510)  56.3 E57-E50 N01-N10
      2  501 1480 (N09E44 L=187)   427    105        1 ( 456/1461- 552/1508)  -6.3 E48-E39 N08-N10
      3 1774 1346 (N32W46 L=276)   400    111        0 (1761/1290-1786/1402)  77.5 W44-W48 N27-N36
      4 1879 1278 (N31W61 L=291)  1367    263        0 (1809/1192-1937/1402) -25.3 W52-W71 N27-N37
      5 1544 1209 (N16W27 L=257)   330     87        1 (1501/1190-1594/1230)   9.5 W24-W30 N14-N16
      6 1638  626 (S14W50 L=281)   420     75        0 (1611/ 595-1675/ 645)  59.0 W48-W52 S16-S11
      7 1698  496 (S19W65 L=295)   619     88        1 (1658/ 480-1749/ 510)  24.0 W59-W70 S20-S16
      8 1480  376 (S32W54 L=285)   344     71        0 (1450/ 360-1521/ 392)  31.0 W50-W57 S33-S30
  

    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