deconvolve#

aiapy.psf.deconvolve( smap, *, psf: Array | None = None, iterations: int = 25, clip_negative: bool = True, )[source]#

Deconvolve an AIA image with the point spread function.

Perform image deconvolution on an AIA image with the instrument point spread function using the Richardson-Lucy deconvolution algorithm [1].

Note

If the jax package is installed it will be used to accelerate the computation. jax can use CPUs or GPUs, see their documentation for instructions. For more information on PSF deconvolution on a GPU, see [2].

Parameters:

smap (Map) – An AIA image.
psf (array-like, optional) – The point spread function. Defaults to None and it will be calculated with aiapy.psf.calculate_psf.
iterations (int, optional) – Number of iterations in the Richardson-Lucy algorithm, defaults to 25.
clip_negative (bool, optional) – If the image has negative intensity values, set them to zero. Defaults to True.

Returns:

Deconvolved AIA image

Return type:

Map