One of the major limitations of adaptive optics (AO) corrected image post-processing is the lack of knowledge on the system point spread function (PSF). The PSF is not always available as a direct imaging on isolated point like objects such as stars. Its prediction using AO telemetry also suffers from serious limitations and requires complex and yet not fully operational algorithms. A very attractive solution consists in a direct PSF estimation using the scientific images themselves thanks to 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 we propose here a marginalized PSF identification that overcomes this limitation. Then 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 marginalisation 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 adaptive optics images of the asteroid 4-Vesta by VLT/SPHERE/Zimpol to demonstrate application to on-sky data.

Context. Adaptive optics (AO) systems greatly increase the resolution of large telescopes, but produce complex point spread function (PSF) shapes, varying in time and across the field of view. This PSF must be accurately known since it provides crucial information about optical systems for design, characterisation, diagnostics and image post processing. Aims. We develop here a model of the AO long exposure PSF, adapted to various seeing conditions and any AO system. This model is made to match accurately both the core of the PSF and its turbulent halo. Methods. The PSF model we develop is based on a parsimonious parameterization of the phase power spectral density with only five parameters to describe circularly symmetric PSFs and seven parameters for asymmetrical ones. Moreover, one of the parameters is directly the Fried parameter r0 of the turbulence s strength. This physical parameter is an asset in the PSF model since it can be correlated with external measurements of the r0, such as phase slopes from the AO real time computer (RTC) or site seeing monitoring. Results. We fit our model against endtoend simulated PSFs using OOMAO tool, and against on sky PSFs from the SPHERE ZIMPOL imager and the MUSE integral field spectrometer working in AO narrowfield mode. Our model matches the shape of the AO PSF both in the core and the halo, with a sub 1 percent relative error for simulated and experimental data. We also show that we retrieve the r0 parameter with subcentimeter precision on simulated data. For ZIMPOL data, we show a correlation of 97 percent between our r0 estimation and the RTC estimation. Finally, MUSE allows us to test the spectral dependency of the fitted r0 parameter. It follows the theoretical $\lambda^{6/5}$ evolution with a standard deviation of 0.3 cm. Evolution of other PSF parameters, such as residual phase variance or aliasing, is also discussed.

Context. Over the past decades, several interplanetary missions have studied small bodies in situ, leading to major advances in our understanding of their geological and geophysical properties. These missions, however, have had a limited number of targets. Among them, the NASA Dawn mission has characterised in detail the topography and albedo variegation across the surface of asteroid (4) Vesta down to a spatial resolution of approximately 20 m/pixel scale. Aims. Here our aim was to determine how much topographic and albedo information can be retrieved from the ground with VLT/SPHERE in the case of Vesta, having a former space mission (Dawn) providing us with the ground truth that can be used as a benchmark. Methods. We observed Vesta with VLT/SPHERE/ZIMPOL as part of our ESO large programme (ID 199.C-0074) at six different epochs, and deconvolved the collected images with a parametric point spread function (PSF). We then compared our images with synthetic views of Vesta generated from the 3D shape model of the Dawn mission, on which we projected Vesta's albedo information. Results. We show that the deconvolution of the VLT/SPHERE images with a parametric PSF allows the retrieval of the main topographic and albedo features present across the surface of Vesta down to a spatial resolution of approximately 20-30 km. Contour extraction shows an accuracy of approximately 1 pixel (3.6 mas). The present study provides the very first quantitative estimate of the accuracy of ground-based adaptive-optics imaging observations of asteroid surfaces. Conclusions. In the case of Vesta, the upcoming generation of 30-40m telescopes (ELT, TMT, GMT) should in principle be able to resolve all of the main features present across its surface, including the troughs and the north-south crater dichotomy, provided that they operate at the diffraction limit.

Ground-based astronomy is severely limited by the atmospheric turbulence, resulting in a large Point-Spread-Function (PSF) and poor imaging resolution. Even imaging with Adaptive Optics (AO) cannot completely correct the aberrated wavefront, and a residual turbulence still corrupts the observation. Thus the consequences of the turbulence on the PSF is of first interest when building any ground-based telescope. The Power Spectral Density (PSD) of a spatially stationary turbulent phase carries all the information needed for describing the long-exposure PSF. We then develop an analytical description of the long-exposure PSF as a series expansion of the aberrated phase PSD. Our description of the PSF given the PSD of the phase is a simple theoretical way to describe the impact of turbulence on the PSF. We also show accordance with previous papers when restricting our model to its first expansion order. Finally we derive applications of our formula to some particular cases, such as Kolmogorov or von-K\'arm\'an models, or the AO correction impact on the PSF.