High-order aberration compensation with Multi-frame Blind Deconvolution and Phase Diversity image restoration techniques

Context. For accurately measuring intensities and determining magnetic field strengths of small-scale solar (magnetic) structure, knowledge of and compensation for the point spread function is crucial. For images recorded with the Swedish 1-meter Solar Telescope, restoration with Multi-Frame Blind Deconvolution and Joint Phase Diverse Speckle methods lead to remarkable improvements in image quality but granulation contrasts that are too low, indicating additional stray light. Aims. We propose a method to compensate for stray light from high-order atmospheric aberrations not included in MFBD and JPDS processing. Methods. To compensate for uncorrected aberrations, a reformulation of the image restoration process is proposed that allows the average effect of hundreds of high-order modes to be compensated for by relying on Kolmogorov statistics for these modes. The applicability of the method requires simultaneous measurements of Fried's parameter r0. The method is tested with simulations as well as real data and extended to include compensation for conventional stray light. Results. We find that only part of the reduction of granulation contrast in SST images is due to uncompensated high-order aberrations. The remainder is still unaccounted for and attributed to stray light from the atmosphere, the telescope with its re-imaging system and to various high-altitude seeing effects. Conclusions. We conclude that statistical compensation of high-order modes is a viable method to reduce the loss of contrast occurring when a limited number of aberrations is explicitly compensated for with MFBD and JPDS processing. We show that good such compensation is possible with only 10 recorded frames. The main limitation of the method is that already MFBD and JPDS processing introduces high-order compensation that, if not taken into account, can lead to over-compensation.Comment: in press in Astronomy & Astrophysic

Image Reconstruction with Analytical Point Spread Functions

The image degradation produced by atmospheric turbulence and optical aberrations is usually alleviated using post-facto image reconstruction techniques, even when observing with adaptive optics systems. These techniques rely on the development of the wavefront using Zernike functions and the non-linear optimization of a certain metric. The resulting optimization procedure is computationally heavy. Our aim is to alleviate this computationally burden. To this aim, we generalize the recently developed extended Zernike-Nijboer theory to carry out the analytical integration of the Fresnel integral and present a natural basis set for the development of the point spread function in case the wavefront is described using Zernike functions. We present a linear expansion of the point spread function in terms of analytic functions which, additionally, takes defocusing into account in a natural way. This expansion is used to develop a very fast phase-diversity reconstruction technique which is demonstrated through some applications. This suggest that the linear expansion of the point spread function can be applied to accelerate other reconstruction techniques in use presently and based on blind deconvolution.Comment: 10 pages, 4 figures, accepted for publication in Astronomy & Astrophysic

A Demonstration of Wavefront Sensing and Mirror Phasing from the Image Domain

In astronomy and microscopy, distortions in the wavefront affect the dynamic range of a high contrast imaging system. These aberrations are either imposed by a turbulent medium such as the atmosphere, by static or thermal aberrations in the optical path, or by imperfectly phased subapertures in a segmented mirror. Active and adaptive optics (AO), consisting of a wavefront sensor and a deformable mirror, are employed to address this problem. Nevertheless, the non-common-path between the wavefront sensor and the science camera leads to persistent quasi-static speckles that are difficult to calibrate and which impose a floor on the image contrast. In this paper we present the first experimental demonstration of a novel wavefront sensor requiring only a minor asymmetric obscuration of the pupil, using the science camera itself to detect high order wavefront errors from the speckle pattern produced. We apply this to correct errors imposed on a deformable microelectromechanical (MEMS) segmented mirror in a closed loop, restoring a high quality point spread function (PSF) and residual wavefront errors of order $\sim 10$ nm using 1600 nm light, from a starting point of $\sim 300$ nm in piston and $\sim 0.3$ mrad in tip-tilt. We recommend this as a method for measuring the non-common-path error in AO-equipped ground based telescopes, as well as as an approach to phasing difficult segmented mirrors such as on the \emph{James Webb Space Telescope} primary and as a future direction for extreme adaptive optics.Comment: 9 pages, 6 figure

Time-sequential Pipelined Imaging with Wavefront Coding and Super Resolution

Wavefront coding has long offered the prospect of mitigating optical aberrations and extended depth of field, but image quality and noise performance are inevitably reduced. We report on progress in the use of agile encoding and pipelined fusion of image sequences to recover image quality

Stokes imaging polarimetry using image restoration at the Swedish 1-m Solar Telescope

Aims: We aim to achieve high spatial resolution as well as high polarimetric sensitivity, using an earth-based 1m-class solar telescope, for the study of magnetic fine structure on the Sun. Methods: We use a setup with 3 high-speed, low-noise cameras to construct datasets with interleaved polarimetric states, particularly suitable for Multi-Object Multi-Frame Blind Deconvolution image restorations. We discuss the polarimetric calibration routine as well as various potential sources of error in the results. Results: We obtained near diffraction limited images, with a noise level of approximately 10^(-3) I(cont). We confirm that dark-cores have a weaker magnetic field and at a lower inclination angle with respect to the solar surface than the edges of the penumbral filament. We show that the magnetic field strength in faculae-striations is significantly lower than in other nearby parts of the faculae.Comment: Accepted for publication in Astronomy & Astrophysics, 12 pages, 11 figure

Learning Wavefront Coding for Extended Depth of Field Imaging

Depth of field is an important factor of imaging systems that highly affects the quality of the acquired spatial information. Extended depth of field (EDoF) imaging is a challenging ill-posed problem and has been extensively addressed in the literature. We propose a computational imaging approach for EDoF, where we employ wavefront coding via a diffractive optical element (DOE) and we achieve deblurring through a convolutional neural network. Thanks to the end-to-end differentiable modeling of optical image formation and computational post-processing, we jointly optimize the optical design, i.e., DOE, and the deblurring through standard gradient descent methods. Based on the properties of the underlying refractive lens and the desired EDoF range, we provide an analytical expression for the search space of the DOE, which is instrumental in the convergence of the end-to-end network. We achieve superior EDoF imaging performance compared to the state of the art, where we demonstrate results with minimal artifacts in various scenarios, including deep 3D scenes and broadband imaging

Speckle Control with a remapped-pupil PIAA-coronagraph

The PIAA is a now well demonstrated high contrast technique that uses an intermediate remapping of the pupil for high contrast coronagraphy (apodization), before restoring it to recover classical imaging capabilities. This paper presents the first demonstration of complete speckle control loop with one such PIAA coronagraph. We show the presence of a complete set of remapping optics (the so-called PIAA and matching inverse PIAA) is transparent to the wavefront control algorithm. Simple focal plane based wavefront control algorithms can thus be employed, without the need to model remapping effects. Using the Subaru Coronagraphic Extreme AO (SCExAO) instrument built for the Subaru Telescope, we show that a complete PIAA-coronagraph is compatible with a simple implementation of a speckle nulling technique, and demonstrate the benefit of the PIAA for high contrast imaging at small angular separation.Comment: 6 figures, submitted to PAS