Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  965 1609 (N39E11 L=127)  1026    182        1 ( 896/1553-1046/1650) -21.0 E17-E04 N35-N41
      2  481 1514 (N27E47 L= 91)  1696    224        5 ( 411/1434- 534/1590) -29.4 E56-E39 N21-N31
      3 1029 1392 (N24E03 L=135)   517    153        0 (1007/1329-1049/1466)  90.0 E04-E01 N20-N29
      4  709 1386 (N21E25 L=113)   240     76        0 ( 675/1365- 749/1400)  21.8 E26-E21 N19-N21
      5 1804 1214 (N20W57 L=195)   362    107        0 (1783/1168-1825/1266) -63.4 W55-W59 N16-N22
      6  624  778 (S17E24 L=114)  2960    626        4 ( 397/ 586- 803/ 990)  43.1 E44-E12 S32-S02
      7 1286  818 (S09W18 L=156)   307    121        1 (1234/ 777-1315/ 853)  51.3 W15-W19 S11-S06
      8 1859  710 (S11W71 L=209)  1288    128        5 (1823/ 657-1892/ 769)  81.9 W67-W77 S14-S07
      9 1612  580 (S21W48 L=186)  2324    491        5 (1493/ 428-1726/ 707)  62.2 W42-W55 S32-S12
     10  289  516 (S40E66 L= 72)   744     70        0 ( 269/ 481- 318/ 548)  16.7 E77-E57 S43-S37
     11 1387  429 (S34W35 L=173)   393     71        0 (1351/ 413-1425/ 447)  -9.5 W31-W37 S34-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