def main():
m31image = create_test_image(frequency=frequency, cellsize=cellsize)
nchan, npol, ny, nx = m31image.data.shape
vt = load(3)
m31image.wcs.wcs.crval[0] = vt.phasecentre.ra.deg
m31image.wcs.wcs.crval[1] = vt.phasecentre.dec.deg
m31image.wcs.wcs.crpix[0] = float(nx // 2)
m31image.wcs.wcs.crpix[1] = float(ny // 2)
fig = show_image(m31image)
fig.savefig('%s_created_img.png' % sys.argv[2])
dump(1, m31image)
def main():
vt = load(1)
#TODO move cellsize and npixel to Drop parameters
model = create_image_from_visibility(vt, cellsize=0.001, npixel=256)
dirty, sumwt = invert_2d(vt, model)
psf, sumwt = invert_2d(vt, model, dopsf=True)
fig = show_image(dirty)
png = '%s_dirty_img.png' % sys.argv[4]
fig.savefig(png)
dump(2, psf)
dump(3, dirty)
def main():
dirty = load(2)
psf = load(1)
comp, residual = deconvolve_cube(dirty, psf, niter=1000, threshold=0.001,
fracthresh=0.01, window='quarter', gain=0.7, scales=[0, 3, 10, 30])
restored = restore_cube(comp, psf, residual)
fig = show_image(comp)
plt.title('Solution')
fig.savefig('%s_Sol.png' % sys.argv[5])
fig = show_image(residual)
plt.title('Residual')
fig.savefig('%s_Residual.png' % sys.argv[6])
fig = show_image(restored)
plt.title('Restored')
fig.savefig('%s_Restored.png' % sys.argv[7])
dump(3, comp)
dump(4, residual)
def test_show_image(self):
show_image(self.m31image)
def main():
img = load(1)
fig = show_image(img)
fig.savefig(sys.argv[2])
nchan=1,
channel_bandwidth=[numpy.sum(channel_bandwidth)],
cellsize=cellsize,
phasecentre=phasecentre)
# In[10]:
dirty, sumwt = invert_function(predicted_vis,
model,
vis_slices=vis_slices,
dopsf=False,
context='timeslice')
# In[11]:
show_image(dirty)
plt.show()
# In[12]:
deconvolved, residual, restored = ical(block_vis=blockvis,
model=model,
vis_slices=vis_slices,
timeslice='auto',
algorithm='hogbom',
niter=1000,
fractional_threshold=0.1,
threshold=0.1,
context='timeslice',
nmajor=5,
gain=0.1,