def test_create_list_spectral(self):
if not self.casacore_available:
return
msfile = rascil_path("data/vis/ASKAP_example.ms")
vis_by_channel = list()
nchan_ave = 16
nchan = 192
for schan in range(0, nchan, nchan_ave):
max_chan = min(nchan, schan + nchan_ave)
v = create_visibility_from_ms(msfile, range(schan, max_chan))
vis_by_channel.append(v[0])
assert len(vis_by_channel) == 12
for v in vis_by_channel:
assert v.vis.data.shape[-1] == 4
assert v.polarisation_frame.type == "linear"
def test_create_list_slice_visibility(self):
if not self.casacore_available:
return
msfile = rascil_path("data/vis/ASKAP_example.ms")
vis_by_channel = list()
nchan_ave = 16
nchan = 192
for schan in range(0, nchan, nchan_ave):
max_chan = min(nchan, schan + nchan_ave)
v = create_visibility_from_ms(msfile, start_chan=schan, end_chan=max_chan - 1)
nchannels = len(numpy.unique(v[0].frequency))
assert nchannels == nchan_ave
vis_by_channel.append(v[0])
assert len(vis_by_channel) == 12
for v in vis_by_channel:
assert v.vis.data.shape[-1] == 4
assert v.polarisation_frame.type == "linear"
def test_create_list_average_slice_visibility(self):
if not self.casacore_available:
return
msfile = rascil_path("data/vis/ASKAP_example.ms")
vis_by_channel = list()
nchan_ave = 16
nchan = 192
for schan in range(0, nchan, nchan_ave):
max_chan = min(nchan, schan + nchan_ave)
v = create_visibility_from_ms(msfile, start_chan=schan, end_chan=max_chan - 1, average_channels=True)
nchannels = len(numpy.unique(v[0].frequency))
assert nchannels == 1
vis_by_channel.append(v[0])
assert len(vis_by_channel) == 12
for ivis, v in enumerate(vis_by_channel):
assert v.vis.data.shape[-1] == 4
assert v.polarisation_frame.type == "linear"
assert numpy.max(numpy.abs(v.vis)) > 0.0, ivis
assert numpy.max(numpy.abs(v.flagged_vis)) > 0.0, ivis
assert numpy.sum(v.weight) > 0.0, ivis
assert numpy.sum(v.flagged_weight) > 0.0, ivis