source: trunk/templates/skymodelplane.html@ 47

Last change on this file since 47 was 47, checked in by JohnLightfoot, 10 years ago

almost working

File size: 2.2 KB
Line 
1<%inherit file="base.html"/>
2
3<!-- do the plotting -->
4<%
5import matplotlib.pyplot as plt
6import numpy as np
7import os.path
8import string
9
10_valid_chars = "_.%s%s" % (string.ascii_letters, string.digits)
11
12def _char_replacer(s):
13 '''A small utility function that echoes the argument or returns '_' if the
14 argument is in a list of forbidden characters.
15 '''
16 if s not in _valid_chars:
17 return '_'
18 return s
19
20def sanitize(text):
21 filename = ''.join(_char_replacer(c) for c in text)
22 return filename
23
24
25plt.figure()
26
27iwn = context['iwn']
28stagename = context['stagename']
29resultname = context['resultname']
30
31sky_image = context['data'][stagename][resultname][:,:,iwn]
32spatial_axis = context['data'][stagename]['spatial axis']
33wn = context['wn']
34
35if sky_image.dtype==np.complex:
36 plt.subplot(211)
37 plt.imshow(sky_image.real, interpolation='nearest', origin='lower',
38 aspect='equal', extent=[spatial_axis[0], spatial_axis[-1],
39 spatial_axis[0], spatial_axis[-1]])
40# vmax=np.max(sky_image.real)*1.1,
41# vmin=np.min(sky_image.real)*0.9)
42 plt.colorbar(orientation='vertical')
43 plt.axis('image')
44 plt.title('Sky model real %6.4gcm-1' % wn)
45
46 plt.subplot(212)
47 plt.imshow(sky_image.imag, interpolation='nearest', origin='lower',
48 aspect='equal', extent=[spatial_axis[0], spatial_axis[-1],
49 spatial_axis[0], spatial_axis[-1]])
50# vmax=np.max(sky_image.imag)*1.1,
51# vmin=np.min(sky_image.imag)*0.9)
52 plt.colorbar(orientation='vertical')
53 plt.axis('image')
54 plt.title('Sky model imag %6.4gcm-1' % wn)
55
56else:
57 plt.imshow(sky_image, interpolation='nearest', origin='lower',
58 aspect='equal', extent=[spatial_axis[0], spatial_axis[-1],
59 spatial_axis[0], spatial_axis[-1]])
60# vmax=np.max(sky_image)*1.1,
61# vmin=np.min(sky_image)*0.9)
62 plt.colorbar(orientation='vertical')
63 plt.axis('image')
64 plt.title('Sky model real %6.4gcm-1' % wn)
65
66filename = '%splane%s.png' % (resultname, wn)
67filename = sanitize(filename)
68filename = os.path.join(context['dirname'], filename)
69plt.savefig(filename)
70plt.close()
71%>
72
73<!-- link to the plot from html -->
74
75<img src='${os.path.basename(filename)}' alt="The sky model should appear here<br>">
Note: See TracBrowser for help on using the repository browser.