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()