Blind deconvolution in astronomy with adaptive optics: the parametric marginal approach

Abstract

ABSTRACT One of the major limitations of using adaptive optics (AO) to correct image post-processing is the lack of knowledge about the system’s point spread function (PSF). The PSF is not always available as direct imaging on isolated point-like objects, such as stars. The use of AO telemetry to predict the PSF also suffers from serious limitations and requires complex and yet not fully operational algorithms. A very attractive solution is to estimate the PSF directly from the scientific images themselves, using blind or myopic post-processing approaches. We demonstrate that such approaches suffer from severe limitations when a joint restitution of object and PSF parameters is performed. As an alternative, here we propose a marginalized PSF identification that overcomes this limitation. In this case, the PSF is used for image post-processing. Here we focus on deconvolution, a post-processing technique to restore the object, given the image and the PSF. We show that the PSF estimated by marginalization provides good-quality deconvolution. The full process of marginalized PSF estimation and deconvolution constitutes a successful blind deconvolution technique. It is tested on simulated data to measure its performance. It is also tested on experimental AO images of the asteroid 4-Vesta taken by the Spectro-Polarimetric High-contrast Exoplanet Research (SPHERE)/Zurich Imaging Polarimeter (Zimpol) on the Very Large Telescope to demonstrate application to on-sky data.