Galaxy Image Restoration with Shape Constraint
arXiv:2101.10021 · doi:10.1007/s00041-021-09880-9
Abstract
Images acquired with a telescope are blurred and corrupted by noise. The blurring is usually modeled by a convolution with the Point Spread Function and the noise by Additive Gaussian Noise. Recovering the observed image is an ill-posed inverse problem. Sparse deconvolution is well known to be an efficient deconvolution technique, leading to optimized pixel Mean Square Errors, but without any guarantee that the shapes of objects (e.g. galaxy images) contained in the data will be preserved. In this paper, we introduce a new shape constraint and exhibit its properties. By combining it with a standard sparse regularization in the wavelet domain, we introduce the Shape COnstraint REstoration algorithm (SCORE), which performs a standard sparse deconvolution, while preserving galaxy shapes. We show through numerical experiments that this new approach leads to a reduction of galaxy ellipticity measurement errors by at least 44%.
arXiv comment22 pages, 6 figures, 1 table, accepted in Journal of Fourier Analysis and Applications
References in corpus (16)
- Weak Gravitational Lensing
- A Method for Weak Lensing Observations
- Shear calibration biases in weak lensing surveys
- Weak lensing for precision cosmology
- Bayesian Galaxy Shape Measurement for Weak Lensing Surveys -I. Methodology and a Fast Fitting Algorithm
- The effects of seeing on Sersic profiles. II. The Moffat PSF
- Bayesian Lensing Shear Measurement
- Precision simulation of ground-based lensing data using observations from space
- Deep learning dark matter map reconstructions from DES SV weak lensing data
- Deep Learning for space-variant deconvolution in galaxy surveys
- Biases in, and corrections to, KSB shear measurements
- Space variant deconvolution of galaxy survey images
- Radio Weak Lensing Shear Measurement in the Visibility Domain - I. Methodology
- The effects of environment on the intrinsic shape of galaxies
- Radio Galaxy Shape Measurement with Hamiltonian Monte Carlo in the Visibility Domain
- Weak Lensing Measurements in Simulations of Radio Images