def visibilities_from_image(self, image):
real_visibilities = self.real_visibilities_from_image(image=image)
imag_visibilities = self.imag_visibilities_from_image(image=image)
return vis.Visibilities(visibilities_1d=np.stack(
(real_visibilities, imag_visibilities), axis=-1))
def mapped_reconstructed_visibilities(self):
visibilities = inversion_util.mapped_reconstructed_visibilities_from(
transformed_mapping_matrix=self.transformed_mapping_matrix,
reconstruction=self.reconstruction,
)
return vis.Visibilities(visibilities=visibilities)
def visibilities_from_image(self, image):
"""
...
"""
warnings.filterwarnings("ignore")
return vis.Visibilities(visibilities=self.forward(
image.binned.native[::-1, :]) # flip due to PyNUFFT internal flip
)
def mapped_reconstructed_visibilities(self):
real_visibilities = inversion_util.mapped_reconstructed_data_from_mapping_matrix_and_reconstruction(
mapping_matrix=self.transformed_mapping_matrices[0],
reconstruction=self.reconstruction,
)
imag_visibilities = inversion_util.mapped_reconstructed_data_from_mapping_matrix_and_reconstruction(
mapping_matrix=self.transformed_mapping_matrices[1],
reconstruction=self.reconstruction,
)
return vis.Visibilities(visibilities_1d=np.stack(
(real_visibilities, imag_visibilities), axis=-1))
def visibilities_from_image(self, image):
if self.preload_transform:
visibilities = transformer_util.visibilities_via_preload_jit_from(
image_1d=image.binned,
preloaded_reals=self.preload_real_transforms,
preloaded_imags=self.preload_imag_transforms,
)
else:
visibilities = transformer_util.visibilities_jit(
image_1d=image.binned,
grid_radians=self.grid,
uv_wavelengths=self.uv_wavelengths,
)
return vis.Visibilities(visibilities=visibilities)
def modified_visibilities_from_visibilities(self, visibilities):
interferometer = copy.deepcopy(self)
interferometer.data = vis.Visibilities(visibilities=visibilities)
return interferometer