Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  973 1608 (N39E11 L=127)   973    162        2 ( 908/1555-1051/1651) -18.4 E16-E03 N35-N41
      2  486 1514 (N27E46 L= 91)  1730    222        6 ( 416/1431- 539/1596) -29.1 E57-E38 N21-N31
      3 1036 1384 (N24E03 L=135)   443    161        1 (1017/1314-1055/1463)  90.0 E04-E01 N19-N29
      4 1805 1226 (N20W58 L=195)   370    112        0 (1784/1178-1822/1276) -71.6 W56-W59 N17-N23
      5  632  778 (S17E23 L=114)  3000    621        3 ( 406/ 590- 805/ 988)  46.0 E43-E12 S32-S02
      6 1292  810 (S09W19 L=156)   350    126        1 (1239/ 773-1321/ 845)  45.0 W15-W20 S11-S06
      7 1625  685 (S14W45 L=182)   423    103        0 (1600/ 660-1655/ 704) -11.3 W42-W47 S15-S12
      8 1599  524 (S25W49 L=186)  1625    307        3 (1493/ 430-1696/ 635)  60.6 W42-W53 S32-S16
      9  291  516 (S40E65 L= 72)   797     75        0 ( 270/ 477- 321/ 547)  16.7 E78-E57 S43-S37
     10 1395  424 (S34W36 L=173)   403     70        0 (1360/ 406-1429/ 440)  -9.5 W32-W38 S35-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