Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1238 1834 (N57E17 L=354)  1317    163        3 (1170/1790-1313/1870)  -3.0 E26-E07 N54-N59
      2  598 1712 (N26E53 L=318)  2652    314        9 ( 562/1567- 648/1831) -25.2 E74-E39 N16-N34
      3 1654 1629 (N57W38 L= 49)   760     73        2 (1633/1582-1672/1666)  42.5 W33-W48 N52-N63
      4 1891  889 (N13W68 L= 79)  1018    121        2 (1873/ 831-1907/ 939)  90.0 W65-W71 N09-N16
      5  492  854 (S28E30 L=341)   309    116        0 ( 481/ 796- 506/ 909) -63.4 E31-E27 S30-S24
      6  641  743 (S31E16 L=355)   220    121        0 ( 586/ 728- 695/ 765)  35.5 E20-E13 S32-S27
      7  786  672 (S31E05 L=  6)   243     67        0 ( 755/ 654- 823/ 693)  36.9 E06-E02 S32-S28
      8 1394  295 (S35W57 L= 68)  1181    151        2 (1326/ 273-1459/ 330)  10.6 W47-W63 S37-S32
  

    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