%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