<%inherit file="base.html"/> <% import matplotlib.pyplot as plt import os.path import string _valid_chars = "_.%s%s" % (string.ascii_letters, string.digits) def _char_replacer(s): '''A small utility function that echoes the argument or returns '_' if the argument is in a list of forbidden characters. ''' if s not in _valid_chars: return '_' return s def sanitize(text): filename = ''.join(_char_replacer(c) for c in text) return filename frequency_axis = context['data']['doublefourier']['frequency axis'] spatial_axis = context['data']['doublefourier']['spatial axis'] spatial_freq_axis = context['data']['doublefourier']['spatial frequency axis'] print 'spatial frequency axis', spatial_freq_axis #for i,wn in enumerate(frequency_axis): print 'frequency axis', frequency_axis for i,wn in enumerate(frequency_axis[:1]): plt.figure() sky_image = context['data']['doublefourier']['sky at time 0'][:,:,i] plt.imshow(sky_image, interpolation='nearest', origin='lower', aspect='equal', extent=[spatial_axis[0], spatial_axis[-1], spatial_axis[0], spatial_axis[-1]]) plt.axis('image') plt.colorbar(orientation='vertical') ## plt.xlabel(beam_image.axes[0].title) ## plt.ylabel(beam_image.axes[1].title) ## plt.title(beam_image.title) ## plt.grid() plt.savefig(os.path.join(context['dirname'], 'dbsky1%s.png' % wn)) plt.close() #for i,wn in enumerate(frequency_axis): for i,wn in enumerate(frequency_axis[:1]): plt.figure() sky_image = context['data']['doublefourier']['sky*beams at time 0'][:,:,i] plt.imshow(sky_image, interpolation='nearest', origin='lower', aspect='equal', extent=[spatial_axis[0], spatial_axis[-1], spatial_axis[0], spatial_axis[-1]]) plt.axis('image') plt.colorbar(orientation='vertical') ## plt.xlabel(beam_image.axes[0].title) ## plt.ylabel(beam_image.axes[1].title) ## plt.title(beam_image.title) ## plt.grid() plt.savefig(os.path.join(context['dirname'], 'dbskybeams1%s.png' % wn)) plt.close() #for i,wn in enumerate(frequency_axis): for i,wn in enumerate(frequency_axis[:1]): plt.figure() sky_image = context['data']['doublefourier']['skyfft at time 0'][:,:,i].real plt.imshow(sky_image, interpolation='nearest', origin='lower', aspect='equal', extent=[spatial_freq_axis[0], spatial_freq_axis[-1], spatial_freq_axis[0], spatial_freq_axis[-1]]) plt.axis('image') plt.colorbar(orientation='vertical') ## plt.xlabel(beam_image.axes[0].title) ## plt.ylabel(beam_image.axes[1].title) ## plt.title(beam_image.title) ## plt.grid() plt.savefig(os.path.join(context['dirname'], 'skyfft%s.png' % wn)) plt.close() #for i,wn in enumerate(frequency_axis): for i,wn in enumerate(frequency_axis[:1]): plt.figure() sky_image = context['data']['doublefourier']['skyfft at time 0'][:,:,i].imag plt.imshow(sky_image, interpolation='nearest', origin='lower', aspect='equal', extent=[spatial_axis[0], spatial_axis[-1], spatial_axis[0], spatial_axis[-1]]) plt.axis('image') plt.colorbar(orientation='vertical') ## plt.xlabel(beam_image.axes[0].title) ## plt.ylabel(beam_image.axes[1].title) ## plt.title(beam_image.title) ## plt.grid() plt.savefig(os.path.join(context['dirname'], 'skyfftimag%s.png' % wn)) plt.close() # interferograms baseline_interferograms = context['data']['doublefourier']['baseline interferograms'] for baseline,interferogram in baseline_interferograms.items(): plt.figure() plt.plot(interferogram) print 'base', baseline plotfile = 'interferogram%s.png' % list(baseline) print 'plotfile', plotfile plotfile = sanitize(plotfile) print 'plotfile', plotfile plt.savefig(os.path.join(context['dirname'], plotfile)) plt.close() %> % for wn in context['data']['doublefourier']['frequency axis'][:1]: % endfor % for baseline in baseline_interferograms.keys(): <% print 'baseline', baseline, type(baseline) plotfile = 'interferogram%s.png' % list(baseline) print 'plotfile', plotfile plotfile = sanitize(plotfile) print 'plotfile', plotfile %> % endfor