Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  561 1697 (N24E54 L=222)  2610    195        6 ( 448/1595- 692/1777)  64.8 E57-E48 N15-N32
      2 1155 1491 (N26E05 L=271)  1771    404        3 (1056/1310-1207/1635) -46.7 E15-W02 N16-N34
      3 1859 1186 (N29W57 L=332)  1664    249        1 (1800/1076-1907/1287) -37.6 W49-W64 N23-N33
      4  358 1139 (S13E45 L=230)   338    100        0 ( 343/1091- 377/1191) -45.0 E47-E42 S15-S10
      5  182 1072 (S20E63 L=213)   991    170        0 ( 141/ 998- 248/1121)  36.9 E68-E55 S25-S16
      6  497  983 (S19E31 L=244)   528    161        0 ( 472/ 924- 533/1052) -82.9 E33-E28 S22-S14
      7  190  950 (S28E60 L=215)   754     72        1 ( 149/ 926- 234/ 967)   9.5 E66-E54 S28-S25
      8 1664  873 (N04W43 L=318)   936    223        1 (1589/ 807-1712/ 964)  70.0 W39-W45 S01-N10
      9  513  502 (S50E25 L=251)  1701    260        3 ( 417/ 429- 620/ 592) -12.8 E35-E12 S53-S45
     10 1676  458 (S17W66 L=341)  2083    291        1 (1558/ 389-1791/ 526)  43.2 W58-W77 S24-S09
  

    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