Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1671 1630 (N36W54 L=321)   637     55        1 (1648/1602-1692/1658)  78.7 W51-W56 N33-N38
      2 1791 1162 (N05W53 L=321)   352     70        1 (1768/1134-1814/1184)  31.0 W50-W55 N03-N06
      3 1843 1140 (N04W59 L=326)   579     83        3 (1824/1105-1864/1175) -78.7 W56-W61 N01-N06
      4  691 1058 (S00E20 L=247)   366    113        0 ( 649/1030- 739/1110) -31.0 E22-E17 S01-N02
      5 1250  820 (S16W14 L=281)  1277    285        0 (1128/ 713-1371/ 926)  38.7 W06-W21 S21-S09
      6 1798  750 (S20W56 L=323)   355     71        0 (1782/ 705-1812/ 783)  78.7 W54-W57 S22-S17
      7 1847  762 (S19W61 L=328)   456    115        3 (1824/ 725-1868/ 783)  36.9 W58-W63 S20-S17
      8 1734  519 (S35W56 L=323)  1175    128        0 (1682/ 486-1785/ 564)  32.0 W51-W59 S36-S31
  

    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