- Timestamp:
- Sep 25, 2014 12:53:16 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/dirtyimage.py
r68 r71 19 19 20 20 # must be a better way of doing this 21 argx = b_x * (wn * 100.0) * spatial_axis / 206265.022 argy = b_y * (wn * 100.0) * spatial_axis / 206265.023 valx = numpy.exp(2 j * math.pi * argx)24 valy = numpy.exp(2 j * math.pi * argy)21 argx = numpy.radians(b_x * (wn * 100.0) * spatial_axis / 3600.0) 22 argy = numpy.radians(b_y * (wn * 100.0) * spatial_axis / 3600.0) 23 valx = numpy.exp(2.0j * math.pi * argx) 24 valy = numpy.exp(2.0j * math.pi * argy) 25 25 26 26 vis = spectrum[wn_spectrum==wn] … … 36 36 37 37 # normalise 38 image /= numpy.sum(beam) 39 beam /= numpy.sum(beam) 38 # The sum of the dirty beam should equal the 0 baseline vis measurement: 0 39 # Normalising the volume would preserve flux under convolution but 40 # is problematic as the volume is 0, so normalise the peak to 1. 41 beam /= 2.0 * len(b_x_list) 42 43 # likewise for dirty map 44 image /= 2.0 * len(b_x_list) 40 45 41 46 return image, beam … … 88 93 for iwn,wn in enumerate(wavenumber): 89 94 # submit jobs 90 print 'starting dirty for', wn91 95 indata = (b_x_list, b_y_list, spectra, wn_spectra, iwn, wn, 92 96 spatial_axis, npix,)
Note:
See TracChangeset
for help on using the changeset viewer.