Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  506 1243 (N31E27 L= 37)  1697    359        3 ( 363/1188- 669/1279) -16.3 E40-E16 N24-N34
      2  430 1029 (N20E37 L= 27)  1356    332        1 ( 309/ 959- 576/1116) -39.5 E46-E27 N12-N27
      3  197  952 (N20E61 L=  3)  1103    173        2 ( 151/ 872- 257/1022) -47.7 E66-E55 N14-N25
      4 1835  907 (S25W52 L=116)   426     50        1 (1810/ 877-1855/ 936)  33.7 W49-W55 S26-S24
      5 1499  833 (S18W22 L= 86)   290    115        0 (1485/ 771-1510/ 889)  56.3 W20-W24 S21-S15
      6 1634  719 (S30W31 L= 94)   290     84        0 (1608/ 684-1663/ 755)  36.9 W28-W33 S31-S28
      7 1697  576 (S43W36 L=100)  1239    227        2 (1598/ 519-1789/ 634)  -2.5 W25-W48 S44-S40
      8  411  435 (S14E61 L=  3)  2707    362        9 ( 320/ 363- 541/ 511) -47.5 E70-E51 S20-S07
      9  936  286 (S37E31 L= 32)  1888    233        5 ( 849/ 201-1013/ 373)  19.2 E44-E21 S41-S31
     10  526  289 (S26E64 L=  0)  1094     76        1 ( 484/ 258- 554/ 320) -90.0 E68-E58 S28-S23
     11  786  200 (S39E50 L= 13)   406     73        0 ( 770/ 179- 813/ 237)  39.8 E53-E46 S41-S36
  

    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