Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  801 1862 (N44E54 L=321)   957     72        2 ( 766/1834- 835/1892)  -7.1 E64-E47 N41-N46
      2  516 1750 (N28E67 L=308)  1144    114        3 ( 475/1707- 559/1788) -74.1 E70-E63 N24-N31
      3  403 1610 (N16E64 L=311)  1190    148        3 ( 343/1553- 446/1647)  -9.9 E79-E56 N12-N17
      4  904 1575 (N26E23 L=352)   633    154        0 ( 882/1509- 940/1632) -77.5 E25-E20 N21-N30
      5  651  936 (S18E20 L=355)   392    107        1 ( 612/ 910- 695/ 979)  -9.5 E23-E17 S19-S16
      6  580  917 (S21E24 L=351)   239     53        1 ( 560/ 896- 609/ 943) -14.0 E26-E22 S21-S19
      7  784  711 (S28E06 L=  9)   823    309        1 ( 643/ 678- 921/ 744)  32.0 E15-W01 S33-S22
      8 1306  343 (S36W45 L= 60)  3448    462        7 (1130/ 272-1480/ 431)  49.9 W32-W52 S45-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