def test_ver_images(self):
'''Test image production.'''
scanset = ScanSet('test.hdf5')
scanset.remove_column('Feed0_RCP-filt')
images = scanset.calculate_images(no_offsets=True, direction=1)
img = images['Feed0_RCP']
def test_rough_image_nooffsets_nofilt(self):
'''Test image production.'''
scanset = ScanSet('test.hdf5')
scanset.remove_column('Feed0_RCP-filt')
images = scanset.calculate_images(no_offsets=True)
img = images['Feed0_RCP']
if HAS_MPL:
fig = plt.figure('img')
plt.imshow(img, origin='lower')
plt.colorbar()
plt.savefig('img_nooff_nofilt.png')
plt.close(fig)
def test_find_scan_through_pixel1(self):
scanset = ScanSet('test.hdf5')
for i in scanset.chan_columns:
scanset.remove_column(i + '-filt')
images = scanset.calculate_images()
ysize, xsize = images['Feed0_RCP'].shape
_, _, _, _, _, _, _, coord = \
scanset.find_scans_through_pixel(xsize//2, 0, test=True)
dec_scan = os.path.join(self.obsdir_dec,
self.dec_scans[self.n_dec_scans // 2])
assert dec_scan in coord
assert coord[dec_scan] == 'dec'
ra_scan = os.path.join(self.obsdir_ra, 'Ra0.fits')
assert ra_scan in coord
assert coord[ra_scan] == 'ra'