Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1083 1861 (N50E25 L= 82)   733    100        1 (1033/1842-1134/1878)  39.8 E28-E20 N47-N53
      2  664 1770 (N31E49 L= 58)   589    100        0 ( 641/1721- 684/1815) -32.5 E54-E43 N27-N34
      3  637 1656 (N22E42 L= 64)   464     82        1 ( 600/1631- 681/1676)  36.9 E44-E39 N20-N24
      4  383  650 (S43E41 L= 65)   287     58        0 ( 353/ 638- 407/ 673)  26.6 E44-E38 S44-S41
      5 1471  439 (S28W49 L=155)   810    226        0 (1412/ 365-1560/ 503)  79.7 W46-W52 S33-S22
      6  898  273 (S58W19 L=125)   346     55        0 ( 875/ 249- 925/ 302) -14.0 W15-W23 S59-S55
      7 1118  160 (S58W55 L=161)   885    114        2 (1091/ 128-1154/ 199) -35.0 W45-W63 S61-S54
  

    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