Ejemplo n.º 1
0
    def test__visualizes_interferometer_using_configs(
        self,
        masked_interferometer_7,
        general_config,
        include_all,
        plot_path,
        plot_patch,
    ):

        visualizer = vis.PhaseInterferometerVisualizer(
            masked_dataset=masked_interferometer_7, image_path=plot_path
        )

        visualizer = visualizer.new_visualizer_with_preloaded_critical_curves_and_caustics(
            preloaded_critical_curves=include_all.preloaded_critical_curves,
            preloaded_caustics=include_all.preloaded_caustics,
        )

        visualizer.visualize_interferometer()

        assert plot_path + "subplots/subplot_interferometer.png" in plot_patch.paths
        assert plot_path + "interferometer/visibilities.png" in plot_patch.paths
        assert plot_path + "interferometer/u_wavelengths.png" not in plot_patch.paths
        assert plot_path + "interferometer/v_wavelengths.png" not in plot_patch.paths
        assert plot_path + "interferometer/primary_beam.png" in plot_patch.paths
Ejemplo n.º 2
0
    def __init__(self, masked_interferometer, settings, cosmology, results=None):

        super(Analysis, self).__init__(
            masked_dataset=masked_interferometer,
            settings=settings,
            cosmology=cosmology,
            results=results,
        )

        self.visualizer = visualizer.PhaseInterferometerVisualizer(
            masked_dataset=masked_interferometer
        )

        result = ag_analysis.last_result_with_use_as_hyper_dataset(results=results)

        if result is not None:

            self.hyper_galaxy_visibilities_path_dict = (
                result.hyper_galaxy_visibilities_path_dict
            )

            self.hyper_model_visibilities = result.hyper_model_visibilities

        else:

            self.hyper_galaxy_visibilities_path_dict = None
            self.hyper_model_visibilities = None
Ejemplo n.º 3
0
    def __init__(self,
                 masked_interferometer,
                 cosmology,
                 image_path=None,
                 results=None):

        super(Analysis, self).__init__(cosmology=cosmology, results=results)

        self.visualizer = visualizer.PhaseInterferometerVisualizer(
            masked_dataset=masked_interferometer, image_path=image_path)

        self.visualizer.visualize_hyper_images(
            hyper_galaxy_image_path_dict=self.hyper_galaxy_image_path_dict,
            hyper_model_image=self.hyper_model_image,
        )

        self.masked_dataset = masked_interferometer

        result = analysis_data.last_result_with_use_as_hyper_dataset(
            results=results)

        if result is not None:

            self.hyper_galaxy_visibilities_path_dict = (
                result.hyper_galaxy_visibilities_path_dict)

            self.hyper_model_visibilities = result.hyper_model_visibilities

        else:

            self.hyper_galaxy_visibilities_path_dict = None
            self.hyper_model_visibilities = None
Ejemplo n.º 4
0
    def test__source_and_lens__visualizes_fit_using_configs(
        self,
        masked_interferometer_7,
        masked_interferometer_fit_x2_plane_inversion_7x7,
        include_all,
        plot_path,
        plot_patch,
    ):

        visualizer = vis.PhaseInterferometerVisualizer(
            masked_dataset=masked_interferometer_7, image_path=plot_path)

        visualizer = visualizer.new_visualizer_with_preloaded_critical_curves_and_caustics(
            preloaded_critical_curves=include_all.preloaded_critical_curves,
            preloaded_caustics=include_all.preloaded_caustics,
        )

        visualizer.visualize_fit(
            fit=masked_interferometer_fit_x2_plane_inversion_7x7,
            during_analysis=True)

        assert plot_path + "subplots/subplot_fit_interferometer.png" in plot_patch.paths
        assert plot_path + "fit_interferometer/visibilities.png" in plot_patch.paths
        assert plot_path + "fit_interferometer/noise_map.png" not in plot_patch.paths
        assert (plot_path + "fit_interferometer/signal_to_noise_map.png"
                not in plot_patch.paths)
        assert (plot_path + "fit_interferometer/model_visibilities.png"
                in plot_patch.paths)
        assert (plot_path +
                "fit_interferometer/residual_map_vs_uv_distances_real.png"
                not in plot_patch.paths)
        assert (
            plot_path +
            "fit_interferometer/normalized_residual_map_vs_uv_distances_real.png"
            in plot_patch.paths)
        assert (plot_path +
                "fit_interferometer/chi_squared_map_vs_uv_distances_real.png"
                in plot_patch.paths)

        #    assert plot_path + "subplots/subplot_inversion.png" in plot_patch.paths
        assert plot_path + "inversion/reconstructed_image.png" in plot_patch.paths
        assert plot_path + "inversion/reconstruction.png" in plot_patch.paths
        assert plot_path + "inversion/errors.png" not in plot_patch.paths
        #  assert plot_path + "inversion/residual_map.png" not in plot_patch.paths
        #  assert plot_path + "inversion/normalized_residual_map.png" not in plot_patch.paths
        #  assert plot_path + "inversion/chi_squared_map.png" in plot_patch.paths
        assert (plot_path + "inversion/regularization_weight_map.png"
                not in plot_patch.paths)
        assert (plot_path + "inversion/interpolated_reconstruction.png"
                not in plot_patch.paths)
        assert plot_path + "inversion/interpolated_errors.png" in plot_patch.paths
Ejemplo n.º 5
0
    def __init__(self, masked_interferometer, cosmology, image_path=None, results=None):

        super(Analysis, self).__init__(cosmology=cosmology, results=results)

        self.visualizer = visualizer.PhaseInterferometerVisualizer(
            masked_dataset=masked_interferometer, image_path=image_path
        )

        self.masked_dataset = masked_interferometer

        if results is not None and results.last is not None:

            self.hyper_galaxy_visibilities_path_dict = (
                results.last.hyper_galaxy_visibilities_path_dict
            )

            self.hyper_model_visibilities = results.last.hyper_model_visibilities
Ejemplo n.º 6
0
    def test__visualizes_interferometer_using_configs(self,
                                                      masked_interferometer_7,
                                                      include_all, plot_path,
                                                      plot_patch):
        visualizer = vis.PhaseInterferometerVisualizer(
            masked_dataset=masked_interferometer_7)

        visualizer = visualizer.new_visualizer_with_preloaded_critical_curves_and_caustics(
            preloaded_critical_curves=include_all.preloaded_critical_curves,
            preloaded_caustics=include_all.preloaded_caustics,
        )

        visualizer.visualize_interferometer(paths=af.Paths())

        assert (path.join(plot_path, "image", "interferometer",
                          "subplot_interferometer.png") in plot_patch.paths)
        assert (path.join(plot_path, "image", "interferometer",
                          "visibilities.png") in plot_patch.paths)
        assert (path.join(plot_path, "image", "interferometer",
                          "u_wavelengths.png") not in plot_patch.paths)
        assert (path.join(plot_path, "image", "interferometer",
                          "v_wavelengths.png") not in plot_patch.paths)
Ejemplo n.º 7
0
    def test__source_and_lens__visualizes_fit_using_configs(
        self,
        masked_interferometer_7,
        masked_interferometer_fit_x2_plane_inversion_7x7,
        include_all,
        plot_path,
        plot_patch,
    ):
        visualizer = vis.PhaseInterferometerVisualizer(
            masked_dataset=masked_interferometer_7)

        visualizer = visualizer.new_visualizer_with_preloaded_critical_curves_and_caustics(
            preloaded_critical_curves=include_all.preloaded_critical_curves,
            preloaded_caustics=include_all.preloaded_caustics,
        )

        visualizer.visualize_fit(
            paths=af.Paths(),
            fit=masked_interferometer_fit_x2_plane_inversion_7x7,
            during_analysis=True,
        )

        assert (path.join(
            plot_path,
            "image",
            "fit_interferometer",
            "subplot_fit_interferometer.png",
        ) in plot_patch.paths)
        assert (path.join(plot_path, "image", "fit_interferometer",
                          "visibilities.png") in plot_patch.paths)
        assert (path.join(plot_path, "image", "fit_interferometer",
                          "noise_map.png") not in plot_patch.paths)
        assert (path.join(plot_path, "image", "fit_interferometer",
                          "signal_to_noise_map.png") not in plot_patch.paths)
        assert (path.join(plot_path, "image", "fit_interferometer",
                          "model_visibilities.png") in plot_patch.paths)
        assert (path.join(
            plot_path,
            "image",
            "fit_interferometer",
            "residual_map_vs_uv_distances_real.png",
        ) not in plot_patch.paths)
        assert (path.join(
            plot_path,
            "image",
            "fit_interferometer",
            "normalized_residual_map_vs_uv_distances_real.png",
        ) in plot_patch.paths)
        assert (path.join(
            plot_path,
            "image",
            "fit_interferometer",
            "chi_squared_map_vs_uv_distances_real.png",
        ) in plot_patch.paths)

        #    assert path.join(plot_path, "image","inversion","subplot_inversion.png") in plot_patch.paths
        assert (path.join(plot_path, "image", "inversion",
                          "reconstructed_image.png") in plot_patch.paths)
        assert (path.join(plot_path, "image", "inversion",
                          "reconstruction.png") in plot_patch.paths)
        assert (path.join(plot_path, "image", "inversion", "errors.png")
                not in plot_patch.paths)
        #  assert path.join(plot_path, "image","inversion","residual_map.png") not in plot_patch.paths
        #  assert path.join(plot_path, "image","inversion","normalized_residual_map.png") not in plot_patch.paths
        #  assert path.join(plot_path, "image","inversion","chi_squared_map.png") in plot_patch.paths
        assert (path.join(plot_path, "image", "inversion",
                          "regularization_weight_map.png")
                not in plot_patch.paths)
        assert (path.join(plot_path, "image", "inversion",
                          "interpolated_reconstruction.png")
                not in plot_patch.paths)
        assert (path.join(plot_path, "image", "inversion",
                          "interpolated_errors.png") in plot_patch.paths)