Space variant deconvolution of galaxy survey images

Abstract

Removing the aberrations introduced by the Point Spread Function (PSF) is a fundamental aspect of astronomical image processing. The presence of noise in observed images makes deconvolution a nontrivial task that necessitates the use of regularisation. This task is particularly difficult when the PSF varies spatially as is the case for the Euclid telescope. New surveys will provide images containing thousand of galaxies and the deconvolution regularisation problem can be considered from a completely new perspective. In fact, one can assume that galaxies belong to a low-rank dimensional space. This work introduces the use of the low-rank matrix approximation as a regularisation prior for galaxy image deconvolution and compares its performance with a standard sparse regularisation technique. This new approach leads to a natural way to handle a space variant PSF. Deconvolution is performed using a Python code that implements a primal-dual splitting algorithm. The data set considered is a sample of 10 000 space-based galaxy images convolved with a known spatially varying Euclid-like PSF and including various levels of Gaussian additive noise. Performance is assessed by examining the deconvolved galaxy image pixels and shapes. The results demonstrate that for small samples of galaxies sparsity performs better in terms of pixel and shape recovery, while for larger samples of galaxies it is possible to obtain more accurate estimates of the galaxy shapes using the low-rank approximation.

12 pages and 8 figures. To be published in A&A

References in corpus (8)

• The Cosmic Evolution Survey (COSMOS) -- Overview
• Bayesian Galaxy Shape Measurement for Weak Lensing Surveys -II. Application to Simulations
• LOFAR Sparse Image Reconstruction
• Low-rank plus sparse decomposition for exoplanet detection in direct-imaging ADI sequences. The LLSG algorithm
• Joint Planck and WMAP CMB Map Reconstruction
• High Resolution Weak Lensing Mass-Mapping Combining Shear and Flexion
• Stellar classification from single-band imaging using machine learning
• Super-resolution method using sparse regularization for point-spread function recovery

Cited by in corpus (8)

• Deep Learning for space-variant deconvolution in galaxy surveys
• PySAP: Python Sparse Data Analysis Package for Multidisciplinary Image Processing
• Sparse Lens Inversion Technique (SLIT): lens and source separability from linear inversion of the source reconstruction problem
• Deep Learning-based galaxy image deconvolution
• Learning to do multiframe wavefront sensing unsupervisedly: applications to blind deconvolution
• Galaxy Image Deconvolution for Weak Gravitational Lensing with Unrolled Plug-and-Play ADMM
• Galaxy Image Restoration with Shape Constraint
• A Distributed Learning Architecture for Scientific Imaging Problems