raise ex
coosys, racen, deccen = determine_phasecenter(ms=selfcal_ms, field=field)
phasecenter = "{0} {1}deg {2}deg".format(coosys, racen, deccen)
(dra, ddec, pixscale) = list(
determine_imsize(ms=selfcal_ms,
field=field,
phasecenter=(racen, deccen),
exclude_7m=exclude_7m,
only_7m=only_7m,
spw='all',
pixfraction_of_fwhm=1 / 8. if only_7m else 1 / 4.))
imsize = [dra, ddec]
cellsize = ['{0:0.2f}arcsec'.format(pixscale)] * 2
for key, value in list(imaging_parameters.items()):
if 'cell' not in imaging_parameters[key]:
imaging_parameters[key]['cell'] = cellsize
if 'imsize' not in imaging_parameters[key]:
imaging_parameters[key]['imsize'] = imsize
contimagename = os.path.join(imaging_root, basename) + "_" + arrayname
if do_bsens:
# just a sanity check
assert 'bsens' in contimagename
if not os.path.exists(contimagename + ".uvwave_vs_amp.png"):
# make a diagnostic plot to show the UV distribution
plotms(
vis=selfcal_ms,
xaxis='uvwave',
origin='contim_selfcal')
coosys, racen, deccen = determine_phasecenter(ms=selfcal_ms, field=field)
phasecenter = "{0} {1}deg {2}deg".format(coosys, racen, deccen)
(dra, ddec, pixscale) = list(
determine_imsize(ms=selfcal_ms,
field=field,
phasecenter=(racen, deccen),
exclude_7m=exclude_7m,
only_7m=only_7m,
spw='all',
pixfraction_of_fwhm=1 / 8. if only_7m else 1 / 4.))
imsize = [dra, ddec]
cellsize = ['{0:0.2f}arcsec'.format(pixscale)] * 2
for key, value in imaging_parameters.items():
if 'cell' not in imaging_parameters[key]:
imaging_parameters[key]['cell'] = cellsize
if 'imsize' not in imaging_parameters[key]:
imaging_parameters[key]['imsize'] = imsize
contimagename = os.path.join(imaging_root, basename) + "_" + arrayname
if do_bsens:
# just a sanity check
assert 'bsens' in contimagename
if not os.path.exists(contimagename + ".uvwave_vs_amp.png"):
# make a diagnostic plot to show the UV distribution
plotms(
vis=selfcal_ms,
xaxis='uvwave',