Extending AMIRAL's blind deconvolution of adaptive optics corrected images with Markov chain Monte Carlo methods

Adaptive optics (AO) corrected image restoration is particularly difficult, as it suffers from the lack of knowledge on the point spread function (PSF) in addition to usual difficulties. An efficient approach is to marginalize the object out of the problem and to estimate the PSF and (object and noise) hyperparameters only, before deconvolving the object using these estimates. Recent works have applied this marginal blind deconvolution method, combined to a parametric model of the PSF, to a series of AO corrected astronomical and satellite images. In this communication, we propose a new restoration method, which consists in choosing the Minimum Mean Square Error (MMSE) estimator and computing the latter thanks to a Markov chain Monte Carlo (MCMC) algorithm. We validate our method by means of realistic simulations, in two very different contexts: an astronomical observation on VLT/SPHERE and a ground-based LEO satellite observation on a 1.52m telescope. Finally, we present results on experimental images for both applications.