pylab.imshow(pow,
origin='lower',
interpolation='nearest',
cmap=pylab.cm.binary,
extent=[0,nmodes,ell[0],ell[-1]],
aspect='auto')
pylab.xlabel('KL mode number')
pylab.ylabel(r'$\ell$')
return ell,pow
if __name__ == '__main__':
from shear_KL_source import params
params.load('../run/base_params.dat')
plot_bandpower( get_basis_filename() )
draw_lines = False
if draw_lines:
xlim = pylab.xlim()
ylim = pylab.ylim()
pylab.plot([0,900],[6144,6144],'--g')
pylab.plot([900,900],[0,6144],'--g')
pylab.savefig('figs/bandpower_marked.eps')
pylab.xlim(xlim)
pylab.ylim(ylim)
pylab.savefig('figs/fig02_bandpower.eps')
pylab.savefig('figs/fig02_bandpower.pdf')
if logplot:
pylab.semilogx(xvals,SN,label='Signal+Noise')
pylab.semilogx(xvals,S,label='Signal')
pylab.semilogx(xvals,N,'--k',label='Noise')
ax.xaxis.set_major_formatter(FuncFormatter(lambda x,*args: '%i'%x))
else:
pylab.plot(xvals,SN,label='Signal+Noise')
pylab.plot(xvals,S,label='Signal')
pylab.plot(xvals,N,'--k',label='Noise')
pylab.legend(loc=4)
pylab.xlim(1,4096)
pylab.grid(True,c='#AAAAAA')
pylab.xlabel('mode number')
pylab.ylabel('cumulative value')
if __name__ == '__main__':
import os
import sys
sys.path.append(os.path.abspath('../'))
from shear_KL_source import params
params.load('../run/base_params.dat')
from shear_KL_source.DES_tile.tools import get_basis_filename
plot_evals( get_basis_filename() )
pylab.savefig('figs/fig03_eigenvalues.eps')
pylab.savefig('figs/fig03_eigenvalues.pdf')
pylab.show()
origin='lower',
cmap=pylab.cm.RdGy,
extent=(0,60,0,60) )
pylab.title('n=%i' % (modes[i]+1))
#pylab.colorbar()
cmax = numpy.max(abs(evec))
pylab.clim(-cmax,cmax)
if i in (0,3,6):
#pylab.gca().yaxis.set_major_formatter(NullFormatter())
pylab.ylabel('arcmin')
if i in (6,7,8):
pylab.xlabel('arcmin')
if __name__ == '__main__':
import os
import sys
sys.path.append(os.path.abspath('../'))
from shear_KL_source import params
from shear_KL_source.DES_tile.tools import get_basis_filename
params.load('../run/base_params.dat')
plot_modes_9( get_basis_filename() )
pylab.savefig('figs/fig01_eigenmodes.eps')
pylab.savefig('figs/fig01_eigenmodes.pdf')
pylab.show()
