source_galaxy = al.Galaxy( redshift=1.0, light=al.lp.EllipticalSersic( centre=(0.0, 0.0), elliptical_comps=(0.096225, -0.055555), intensity=0.4, effective_radius=0.5, sersic_index=1.0, ), ) for instrument in ["vro", "euclid", "hst", "hst_up", "ao"]: imaging = instrument_util.load_test_imaging( data_name="lens_sie__source_smooth", instrument=instrument, psf_shape_2d=psf_shape_2d, ) mask = al.Mask.circular( shape_2d=imaging.shape_2d, pixel_scales=imaging.pixel_scales, sub_size=sub_size, radius=radius, ) masked_imaging = al.MaskedImaging(imaging=imaging, mask=mask) print("Light profile fit run times for image type " + instrument + "\n") print("Number of points = " + str(masked_imaging.grid.sub_shape_1d) + "\n")
from autolens.plot import plotters from test_autolens.simulate.imaging import instrument_util # In this tutorial, we'll introduce a new pixelization, called an adaptive-pixelization. This pixelization doesn't use # uniform grid of rectangular pixels, but instead uses ir'Voronoi' pixels. So, why would we want to do that? # Lets take another look at the rectangular grid, and think about its weakness. # Lets quickly remind ourselves of the image, and the 3.0" circular mask we'll use to mask it. imaging = instrument_util.load_test_imaging( data_name="lens_light_dev_vaucouleurs", instrument="vro") array = imaging.image mask = al.Mask.elliptical( shape=imaging.shape_2d, pixel_scales=imaging.pixel_scales, major_axis_radius=6.0, axis_ratio=0.5, phi=0.0, centre=(0.0, 0.0), ) aplt.Array(array=array, mask=mask, positions=[[(1.0, 1.0)]], centres=[[(0.0, 0.0)]]) imaging = instrument_util.load_test_imaging( data_name="lens_sis__source_smooth__offset_centre", instrument="vro") array = imaging.image mask = al.Mask.elliptical( shape=imaging.shape_2d,
import autolens as al import autolens.plot as aplt from test_autolens.simulate.imaging import instrument_util # In this tutorial, we'll introduce a new pixelization, called an adaptive-pixelization. This pixelization doesn't use # uniform grid of rectangular pixels, but instead uses ir'Voronoi' pixels. So, why would we want to do that? # Lets take another look at the rectangular grid, and think about its weakness. # Lets quickly remind ourselves of the image, and the 3.0" circular mask we'll use to mask it. imaging = instrument_util.load_test_imaging( data_name="lens_light_dev_vaucouleurs", instrument="vro" ) array = imaging.image array_overlay = al.Array.manual_2d(array=[[1.0, 2.0], [3.0, 4.0]], pixel_scales=5.05) aplt.Array(array=array, array_overlay=array_overlay)
import autolens as al from test_autolens.simulate.imaging import instrument_util imaging = instrument_util.load_test_imaging( data_name="lens_sie__source_smooth__offset_centre", instrument="vro") def fit_with_offset_centre(centre): mask = al.Mask.elliptical( shape_2d=imaging.shape_2d, pixel_scales=imaging.pixel_scales, major_axis_radius=3.0, axis_ratio=1.0, phi=0.0, centre=centre, ) # The lines of code below do everything we're used to, that is, setup an image and its grid, mask it, trace it # via a tracer, setup the rectangular mapper, etc. lens_galaxy = al.Galaxy( redshift=0.5, mass=al.mp.EllipticalIsothermal(centre=(2.0, 2.0), einstein_radius=1.2, elliptical_comps=(0.17647, 0.0)), ) source_galaxy = al.Galaxy( redshift=1.0, pixelization=al.pix.VoronoiMagnification(shape=(20, 20)), regularization=al.reg.Constant(coefficient=1.0),