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',