Skip to content
Snippets Groups Projects
Commit 426d1c31 authored by Julian Rüstig's avatar Julian Rüstig
Browse files

multi: some changes

parent bcc66910
Branches
No related tags found
1 merge request!50Draft: Mosaic imaging
......@@ -163,27 +163,6 @@ for ii in range(1):
plt.imshow(tot_response_adjoint[0, 0, 0], origin='lower')
plt.show()
# ff = (SKY_BEAMER @ REDUCER)(f).val
# fig, axes = plt.subplots(2, 2)
# ((a00, a01), (a10, a11)) = axes
# i00 = a00.imshow(f.val[0, 0, 0], origin='lower',
# vmin=f.val.min(), vmax=f.val.max())
# i01 = a01.imshow(tot_response_adjoint[0, 0, 0], origin='lower')
# i10 = a10.imshow(f.val[0, 0, 0, xi:xe, yi:ye], origin='lower',
# vmin=f.val.min(), vmax=f.val.max())
# i11 = a11.imshow(ff['fld0'], origin='lower',
# vmin=f.val.min(), vmax=f.val.max())
# plt.colorbar(i00, ax=a00)
# plt.colorbar(i01, ax=a01)
# plt.colorbar(i10, ax=a10)
# plt.colorbar(i11, ax=a11)
# a00.scatter([xi, xi, xe, xe], [yi, ye, yi, ye], c='orange')
# plt.show()
# lh = rve.ImagingLikelihood(obs, sky, 1e-7, False, nthreads=1)
lh = reduce(lambda x, y: x+y, energies)
lh = lh @ sky
......
[sky]
freq mode = single
polarization=I
space npix x = 128
space npix y = 128
space npix x = 256
space npix y = 256
space fov x = 25as
space fov y = 25as
# space fov x = 30.72as
# space fov y = 30.72as
image center ra = 13h37m00.75s
image center dec = -29d51m58.6s
image center frame = icrs
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment