def plot(bvis):
plt.clf()
plot_uvcoverage(
[bvis], title='UV Coverage {source:s}'.format(source=bvis.source))
plt.tight_layout()
plt.show(block=False)
plt.clf()
plot_visibility(
[bvis],
title='Visibility amplitude {source:s}'.format(source=bvis.source))
plt.tight_layout()
plt.show(block=False)
return bvis
guard_band_image=1.0,
delA=0.02)
advice = rsexecute.compute(advice_list, sync=True)
pb_npixel = 1024
d2r = numpy.pi / 180.0
pb_cellsize = d2r * FOV_deg / pb_npixel
cellsize = advice['cellsize']
if show:
bvis_list = rsexecute.compute(bvis_graph, sync=True)
plot_azel(bvis_list, title=basename)
plt.savefig('azel.png')
plt.show(block=False)
vis_list = rsexecute.compute(bvis_graph, sync=True)
plot_uvcoverage(vis_list, title=basename)
plt.savefig('uvcoverage.png')
plt.show(block=False)
# Now construct the components
print("Constructing components")
original_components, offset_direction = create_simulation_components(
context,
phasecentre,
frequency,
pbtype, [0.0, 0.0],
flux_limit,
pbradius * HWHM,
pb_npixel,
pb_cellsize,
fov=10)