source: trunk/templates/doublefourier.html@ 17

<%inherit file="base.html"/>

<!-- do the plotting -->
<%
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()

%>

<!-- link to the plot from html -->

% for wn in context['data']['doublefourier']['frequency axis'][:1]:
    <img src='dbsky1${wn}.png' alt="The sky model should appear here<br>">
    <img src='dbskybeams1${wn}.png' alt="The sky*beams should appear here<br>">
    <img src='skyfft${wn}.png' alt="The skyfft should appear here<br>">
    <img src='skyfftimag${wn}.png' alt="The skyfft.imag should appear here<br>">
% 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
%>
    <img src='${plotfile}' alt="The interferogram should appear here<br>">
% endfor