Changeset 28 for trunk/templates/doublefourier.html
- Timestamp:
- May 20, 2014, 3:20:46 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/templates/doublefourier.html
r25 r28 3 3 <!-- do the plotting --> 4 4 <% 5 import matplotlib.pyplot as plt6 5 import os.path 7 import string8 9 _valid_chars = "_.%s%s" % (string.ascii_letters, string.digits)10 def _char_replacer(s):11 '''A small utility function that echoes the argument or returns '_' if the12 argument is in a list of forbidden characters.13 '''14 if s not in _valid_chars:15 return '_'16 return s17 18 def sanitize(text):19 filename = ''.join(_char_replacer(c) for c in text)20 return filename21 22 6 23 7 frequency_axis = context['data']['doublefourier']['frequency axis'] 24 spatial_axis = context['data']['doublefourier']['spatial axis']25 spatial_freq_axis = context['data']['doublefourier']['spatial frequency axis']26 8 27 #for i,wn in enumerate(frequency_axis): 28 for i,wn in enumerate(frequency_axis[:1]): 29 plt.figure() 9 # construct html pages with sky model plane plots 10 for iwn,wn in enumerate(frequency_axis): 11 link = os.path.join(context['dirname'], 'skymodelt0plane%s.html' % (wn)) 12 with open(link, 'w') as f: 13 template = mylookup.get_template('skymodelplane.html') 14 context_copy = context.kwargs 15 context_copy['wn'] = wn 16 context_copy['iwn'] = iwn 17 context_copy['stagename'] = 'doublefourier' 18 context_copy['resultname'] = 'sky at time 0' 19 f.write(template.render(**context_copy)) 30 20 31 sky_image = context['data']['doublefourier']['sky at time 0'][:,:,i] 32 plt.imshow(sky_image, interpolation='nearest', origin='lower', 33 aspect='equal', extent=[spatial_axis[0], spatial_axis[-1], 34 spatial_axis[0], spatial_axis[-1]]) 35 plt.axis('image') 21 for iwn,wn in enumerate(frequency_axis): 22 link = os.path.join(context['dirname'], 23 'skymodelpbt0plane%s.html' % (wn)) 24 with open(link, 'w') as f: 25 template = mylookup.get_template('skymodelplane.html') 26 context_copy = context.kwargs 27 context_copy['wn'] = wn 28 context_copy['iwn'] = iwn 29 context_copy['stagename'] = 'doublefourier' 30 context_copy['resultname'] = 'sky*beams at time 0' 31 f.write(template.render(**context_copy)) 36 32 37 plt.colorbar(orientation='vertical') 38 ## plt.xlabel(beam_image.axes[0].title) 39 ## plt.ylabel(beam_image.axes[1].title) 40 ## plt.title(beam_image.title) 33 for iwn,wn in enumerate(frequency_axis): 34 link = os.path.join(context['dirname'], 35 'skymodelfftt0plane%s.html' % (wn)) 36 with open(link, 'w') as f: 37 template = mylookup.get_template('skymodelplane.html') 38 context_copy = context.kwargs 39 context_copy['wn'] = wn 40 context_copy['iwn'] = iwn 41 context_copy['stagename'] = 'doublefourier' 42 context_copy['resultname'] = 'skyfft at time 0' 43 f.write(template.render(**context_copy)) 41 44 42 ## plt.grid() 43 plt.savefig(os.path.join(context['dirname'], 'dbsky1%s.png' % wn)) 44 plt.close() 45 link = os.path.join(context['dirname'], 'sourcespectrumbaseline.html') 46 with open(link, 'w') as f: 47 template = mylookup.get_template('sourcespectrum.html') 48 context_copy = context.kwargs 49 context_copy['sourceid'] = 'baseline' 50 context_copy['spectrum'] = context['data']['doublefourier']\ 51 ['skyfft spectrum at time 0'] 52 print 'before' 53 f.write(template.render(**context_copy)) 54 print 'after' 45 55 46 47 #for i,wn in enumerate(frequency_axis): 48 for i,wn in enumerate(frequency_axis[:1]): 49 plt.figure() 50 51 sky_image = context['data']['doublefourier']['sky*beams at time 0'][:,:,i] 52 plt.imshow(sky_image, interpolation='nearest', origin='lower', 53 aspect='equal', extent=[spatial_axis[0], spatial_axis[-1], 54 spatial_axis[0], spatial_axis[-1]]) 55 plt.axis('image') 56 57 plt.colorbar(orientation='vertical') 58 ## plt.xlabel(beam_image.axes[0].title) 59 ## plt.ylabel(beam_image.axes[1].title) 60 ## plt.title(beam_image.title) 61 62 ## plt.grid() 63 plt.savefig(os.path.join(context['dirname'], 'dbskybeams1%s.png' % wn)) 64 plt.close() 65 66 #for i,wn in enumerate(frequency_axis): 67 for i,wn in enumerate(frequency_axis[:1]): 68 plt.figure() 69 70 sky_image = context['data']['doublefourier']['skyfft at time 0'][:,:,i].real 71 plt.imshow(sky_image, interpolation='nearest', origin='lower', 72 aspect='equal', extent=[spatial_freq_axis[0], spatial_freq_axis[-1], 73 spatial_freq_axis[0], spatial_freq_axis[-1]]) 74 plt.axis('image') 75 76 plt.colorbar(orientation='vertical') 77 ## plt.xlabel(beam_image.axes[0].title) 78 ## plt.ylabel(beam_image.axes[1].title) 79 ## plt.title(beam_image.title) 80 81 ## plt.grid() 82 plt.savefig(os.path.join(context['dirname'], 'skyfft%s.png' % wn)) 83 plt.close() 84 85 #for i,wn in enumerate(frequency_axis): 86 for i,wn in enumerate(frequency_axis[:1]): 87 plt.figure() 88 89 sky_image = context['data']['doublefourier']['skyfft at time 0'][:,:,i].imag 90 plt.imshow(sky_image, interpolation='nearest', origin='lower', 91 aspect='equal', extent=[spatial_axis[0], spatial_axis[-1], 92 spatial_axis[0], spatial_axis[-1]]) 93 plt.axis('image') 94 95 plt.colorbar(orientation='vertical') 96 ## plt.xlabel(beam_image.axes[0].title) 97 ## plt.ylabel(beam_image.axes[1].title) 98 ## plt.title(beam_image.title) 99 100 ## plt.grid() 101 plt.savefig(os.path.join(context['dirname'], 'skyfftimag%s.png' % wn)) 102 plt.close() 103 104 # interferograms 105 baseline_interferograms = context['data']['doublefourier']['baseline interferograms'] 106 for baseline,interferogram in baseline_interferograms.items(): 107 plt.figure() 108 plt.plot(interferogram) 109 plotfile = 'interferogram%s.png' % list(baseline) 110 plotfile = sanitize(plotfile) 111 plt.savefig(os.path.join(context['dirname'], plotfile)) 112 plt.close() 113 56 link = os.path.join(context['dirname'], 'baselineffttest.html') 57 with open(link, 'w') as f: 58 template = mylookup.get_template('sourcespectrum.html') 59 context_copy = context.kwargs 60 context_copy['sourceid'] = 'testfft' 61 context_copy['spectrum'] = context['data']['doublefourier']\ 62 ['test FTS at time 0'] 63 print 'before' 64 f.write(template.render(**context_copy)) 65 print 'after' 114 66 %> 115 67 116 <!-- link to the plot from html-->68 <!-- link to the plot pages from this page --> 117 69 118 % for wn in context['data']['doublefourier']['frequency axis'][:1]: 119 <img src='dbsky1${wn}.png' alt="The sky model should appear here<br>"> 120 <img src='dbskybeams1${wn}.png' alt="The sky*beams should appear here<br>"> 121 <img src='skyfft${wn}.png' alt="The skyfft should appear here<br>"> 122 <img src='skyfftimag${wn}.png' alt="The skyfft.imag should appear here<br>"> 70 <h3>Sky Model Planes at t=0</h3> 71 <div class="row"> 72 % for wn in frequency_axis: 73 <% 74 wnformatted = '%6.4gcm-1' % wn 75 %> 76 <div class="col-xs-6 col-md-3"> 77 <a href="skymodelt0plane${wn}.html" class="thumbnail"> 78 <img src="sky_at_time_0plane${wn}.png" alt="..."> 79 <div> 80 <h4>${wnformatted}</h4> 81 </div> 82 </a> 83 </div> 123 84 % endfor 85 </div> 124 86 125 % for baseline in baseline_interferograms.keys(): 126 <% 127 plotfile = 'interferogram%s.png' % list(baseline) 128 plotfile = sanitize(plotfile) 129 %> 130 <img src='${plotfile}' alt="The interferogram should appear here<br>"> 87 <h3>Sky Model Planes*Pb at t=0</h3> 88 <div class="row"> 89 % for wn in frequency_axis: 90 <% 91 wnformatted = '%6.4gcm-1' % wn 92 %> 93 <div class="col-xs-6 col-md-3"> 94 <a href="skymodelpbt0plane${wn}.html" class="thumbnail"> 95 <img src="sky_beams_at_time_0plane${wn}.png" alt="..."> 96 <div> 97 <h4>${wnformatted}</h4> 98 </div> 99 </a> 100 </div> 131 101 % endfor 102 </div> 103 104 <h3>Sky Model FFT at t=0</h3> 105 <div class="row"> 106 % for wn in frequency_axis: 107 <% 108 wnformatted = '%6.4gcm-1' % wn 109 %> 110 <div class="col-xs-6 col-md-3"> 111 <a href="skymodelfftt0plane${wn}.html" class="thumbnail"> 112 <img src="skyfft_at_time_0plane${wn}.png" alt="..."> 113 <div> 114 <h4>${wnformatted}</h4> 115 </div> 116 </a> 117 </div> 118 % endfor 119 </div> 120 121 <h3>Measured spectrum at t=0</h3> 122 <div class="col-xs-6 col-md-3"> 123 <a href="sourcespectrumbaseline.html" class="thumbnail"> 124 <img src="sourcespectrumbaseline.png" alt="..."> 125 </a> 126 </div> 127 128 <h3>Measured interferogram at t=0</h3> 129 <div class="col-xs-6 col-md-3"> 130 <a href="baselineffttest.html" class="thumbnail"> 131 <img src="sourcespectrumtestfft.png" alt="..."> 132 </a> 133 </div>
Note:
See TracChangeset
for help on using the changeset viewer.