A Blind Deconvolution Algorithm using Zernike-polynomial-based Phase Fitting

The main inherent property of medical ultrasound imaging is speckle noise which generally obscures and reduces the diagnostic image resolution and contrast. Consequently, the substantial improvement of ultrasound images is an important prerequisite, whenever ultrasound is used as one of the most utilized diagnostic modalities. The problem of reconstruction of ultrasound images by means of blind deconvolution has long been recognized as a prerequisite in ultrasound imaging processing in the recent decades. Recently, the Total Variation (TV) regularization method has become extremely effective approach for image reconstruction, especially for restoring edges of the blurring image. In this paper, we present a new blind iterative TV deconvolution algorithm for reconstructing ultrasound images from blurry and noisy observations. First, it proposes the initial estimation of the point-spread function (PSF) based on a generalization of two-dimensional homomorphic filtering in cepstrum domain. It is demonstrated that the initial PSF can be effectively estimated by applying a proper smoothing low-pass filtering in cepstrum domain. Second, it introduces a novel blind iterative TV deconvolution which is derived from an alternating minimization algorithm. Fast Fourier Transform (FFT) is used in the pre-iteration computation. The innovative blind deconvolution is based on either concurrent or successive estimation of the PSF function and the image of interest. The iterative scheme is devised to recover the image and simultaneously identify the PSF function. The estimated PSF and restored image will be close to real values in the subsequent iterative deconvolution. Final, tests on phantom and clinical images have proven our novel blind iterative TV deconvolution gives more positive results of higher spatial resolution and better defined tissue structures than other deconvolution methods.

Aim: The objective of this project was to evaluate high resolution images from an adaptive optics retinal imager through comparisons with standard film-based and standard digital fundus imagers. Methods: A clinical prototype adaptive optics fundus imager (AOFI) was used to collect retinal images from subjects with various forms of retinopathy to determine whether improved visibility into the disease could be provided to the clinician. The AOFI achieves low-order correction of aberrations through a closed-loop wavefront sensor and an adaptive optics system. The remaining high-order aberrations are removed by direct deconvolution using the point spread function (PSF) or by blind deconvolution when the PSF is not available. An ophthalmologist compared the AOFI images with standard fundus images and provided a clinical evaluation of all the modalities and processing techniques. All images were also analyzed using a quantitative image quality index. Results: This system has been tested on three human subjects (one normal and two with retinopathy). In the diabetic patient vascular abnormalities were detected with the AOFI that cannot be resolved with the standard fundus camera. Very small features, such as the fine vascular structures on the optic disc and the individual nerve fiber bundles are easily resolved by the AOFI. Conclusion: This project demonstrated that adaptive optic images have great potential in providing clinically significant detail of anatomical and pathological structures to the ophthalmologist.

Context. Modeling the optical point spread function (PSF) is particularly challenging for adaptive optics (AO)-assisted observations owing to the its complex shape and spatial variations. Aims. We aim to (i) exhaustively demonstrate the accuracy of a recent analytical model from comparison with a large sample of imaged PSFs, (ii) assess the conditions for which the model is optimal, and (iii) unleash the strength of this framework to enable the joint estimation of atmospheric parameters, AO performance and static aberrations. Methods. We gathered 4812 on-sky PSFs obtained from seven AO systems and used the same fitting algorithm to test the model on various AO PSFs and diagnose AO performance from the model outputs. Finally, we highlight how this framework enables the characterization of the so-called low wind effect on the Spectro-Polarimetic High contrast imager for Exoplanets REsearch (LWE; SPHERE) instrument and piston cophasing errors on the Keck II telescope. Results. Over 4812 PSFs, the model reaches down to 4% of error on both the Strehl-ratio (SR) and full width at half maximum (FWHM). We particularly illustrate that the estimation of the Fried’s parameter, which is one of the model parameters, is consistent with known seeing statistics and follows expected trends in wavelength using the Multi Unit Spectroscopic Explorer instrument (λ6/5) and field (no variations) from Gemini South Adaptive Optics Imager images with a standard deviation of 0.4 cm. Finally, we show that we can retrieve a combination of differential piston, tip, and tilt modes introduced by the LWE that compares to ZELDA measurements, as well as segment piston errors from the Keck II telescope and particularly the stair mode that has already been revealed from previous studies. Conclusions. This model matches all types of AO PSFs at the level of 4% error and can be used for AO diagnosis, post-processing, and wavefront sensing purposes.

We propose to use the blind deconvolution and its modification to extract the point spread function (PSF) of layered metamaterials from a SNOM measurement. The measurement results are processed using the blind deconvolution algorithm to reconstruct the real-valued non-coherent PSF, or using the modified blind deconvolution introduced in this paper to reconstruct the complex-valued coherent PSF. The two algorithms are tested in simulations with a layered metamaterial deposited on a thin metallic mask with test apertures. We show that the modified algorithm is capable of recovering the approximate shape of complex PSF with a sub-wavelength full width at half maximum from a measurement in which the apertures are larger than the wavelength.

Accurate astrometry and photometry of saturated and coronagraphic point spread functions (PSFs) are fundamental to both ground- and space-based high contrast imaging projects. For ground-based adaptive optics imaging, differential atmospheric refraction and flexure introduce a small drift of the PSF with time, and seeing and sky transmission variations modify the PSF flux distribution. For space-based imaging, vibrations, thermal fluctuations and pointing jitters can modify the PSF core position and flux. These effects need to be corrected to properly combine the images and obtain optimal signal-to-noise ratios, accurate relative astrometry and photometry of detected objects as well as precise detection limits. Usually, one can easily correct for these effects by using the PSF core, but this is impossible when high dynamic range observing techniques are used, like coronagrahy with a non-transmissive occulting mask, or if the stellar PSF core is saturated. We present a new technique that can solve these issues by using off-axis satellite PSFs produced by a periodic amplitude or phase mask conjugated to a pupil plane. It will be shown that these satellite PSFs track precisely the PSF position, its Strehl ratio and its intensity and can thus be used to register and to flux normalize the PSF. A laboratory experiment is also presented to validate the theory. This approach can be easily implemented in existing adaptive optics instruments and should be considered for future extreme adaptive optics coronagraph instruments and in high-contrast imaging space observatories.

The Integrated Optical System (IOS) is an extreme adaptive optics system designed for NASA’s Laser Com- munication Relay Demonstration mission. There is a great deal of overlap between the requirements for laser communication AO and high-contrast exoplanet imaging AO systems. Both require very high Strehl ratios with narrow fields of view. This overlap allows the IOS to serve as a testbed and technology demonstrator for astronomical extreme adaptive optics systems. There are several example technologies from the IOS that are already making the transition to astronomical AO systems. The first is that the real time controller based on Direct Memory Access transfer between the WFS camera link frame-grabber and a DSP board is being reused on the upgrade to PALM-3000 AO system at Palomar Observatory. This enables the system to minimize latency by bypassing the CPU and its inherent timing jitter. Technologies like this will be crucial to enabling high contrast imaging on the next generation of extremely large telescopes. In addition, the IOS measures Fried’s parameter from wavefraont measures in near real time. This technology has already been deployed to PALM-3000. The main function of Laser Communication AO systems is to couple the incoming light into single mode fiber. This is the same configuration that will be used by AO coupled radial velocity spectrographs. The adaptive optics system is a woofer/tweeter design, with one deformable mirror correcting for low spatial frequencies with large amplitude and a second deformable mirror correcting for high spatial frequencies with small amplitude. The system uses a Shack-Hartmann wavefront sensor. The system has achieved first light and is undergoing commissioning. We will present an overview of the system design and initial performance.

An Adaptive Optic (AO) system is a fundamental requirement of 8m-class telescopes. We know that in order to obtain the maximum possible resolution allowed by these telescopes we need to correct the atmospheric turbulence. Thanks to adaptive optic systems we are able to use all the effective potential of these instruments, drawing all the information from the universe sources as best as possible. In an AO system there are two main components: the wavefront sensor (WFS) that is able to measure the aberrations on the incoming wavefront in the telescope, and the deformable mirror (DM) that is able to assume a shape opposite to the one measured by the sensor. The two subsystem are connected by the reconstructor (REC). In order to do this, the REC requires a “common language" between these two main AO components. It means that it needs a mapping between the sensor-space and the mirror-space, called an interaction matrix (IM). Therefore, in order to operate correctly, an AO system has a main requirement: the measure of an IM in order to obtain a calibration of the whole AO system. The IM measurement is a 'mile stone' for an AO system and must be done regardless of the telescope size or class. Usually, this calibration step is done adding to the telescope system an auxiliary artificial source of light (i.e a fiber) that illuminates both the deformable mirror and the sensor, permitting the calibration of the AO system. For large telescope (more than 8m, like Extremely Large Telescopes, ELTs) the fiber based IM measurement requires challenging optical setups that in some cases are also impractical to build. In these cases, new techniques to measure the IM are needed. In this PhD work we want to check the possibility of a different method of calibration that can be applied directly on sky, at the telescope, without any auxiliary source. Such a technique can be used to calibrate AO system on a telescope of any size. We want to test the new calibration technique, called “sinusoidal modulation technique”, on the Large Binocular Telescope (LBT) AO system, which is already a complete AO system with the two main components: a secondary deformable mirror with by 672 actuators, and a pyramid wavefront sensor. My first phase of PhD work was helping to implement the WFS board (containing the pyramid sensor and all the auxiliary optical components) working both optical alignments and tests of some optical components. Thanks to the “solar tower” facility of the Astrophysical Observatory of Arcetri (Firenze), we have been able to reproduce an environment very similar to the telescope one, testing the main LBT AO components: the pyramid sensor and the secondary deformable mirror. Thanks to this the second phase of my PhD thesis: the measure of IM applying the sinusoidal modulation technique. At first we have measured the IM using a fiber auxiliary source to calibrate the system, without any kind of disturbance injected. After that, we have tried to use this calibration technique in order to measure the IM directly “on sky”, so adding an atmospheric disturbance to the AO system. The results obtained in this PhD work measuring the IM directly in the Arcetri solar tower system are crucial for the future development: the possibility of the acquisition of IM directly on sky means that we are able to calibrate an AO system also for extremely large telescope class where classic IM measurements technique are problematic and, sometimes, impossible. Finally we have not to forget the reason why we need this: the main aim is to observe the universe. Thanks to these new big class of telescopes and only using their full capabilities, we will be able to increase our knowledge of the universe objects observed, because we will be able to resolve more detailed characteristics, discovering, analyzing and understanding the behavior of the universe components.

Context.The direct imaging of potentially habitable exoplanets is one prime science case for the next generation of high contrast imaging instruments on ground-based, extremely large telescopes. To reach this demanding science goal, the instruments are equipped with eXtreme Adaptive Optics (XAO) systems which will control thousands of actuators at a framerate of kilohertz to several kilohertz. Most of the habitable exoplanets are located at small angular separations from their host stars, where the current control laws of XAO systems leave strong residuals.Aims.Current AO control strategies such as static matrix-based wavefront reconstruction and integrator control suffer from a temporal delay error and are sensitive to mis-registration, that is, to dynamic variations of the control system geometry. We aim to produce control methods that cope with these limitations, provide a significantly improved AO correction, and, therefore, reduce the residual flux in the coronagraphic point spread function (PSF).Methods.We extend previous work in reinforcement learning for AO. The improved method, called the Policy Optimization for Adaptive Optics (PO4AO), learns a dynamics model and optimizes a control neural network, called a policy. We introduce the method and study it through numerical simulations of XAO with Pyramid wavefront sensor (PWFS) for the 8-m and 40-m telescope aperture cases. We further implemented PO4AO and carried out experiments in a laboratory environment using Magellan Adaptive Optics eXtreme system (MagAO-X) at the Steward laboratory.Results.PO4AO provides the desired performance by improving the coronagraphic contrast in numerical simulations by factors of 3–5 within the control region of deformable mirror and PWFS, both in simulation and in the laboratory. The presented method is also quick to train, that is, on timescales of typically 5–10 s, and the inference time is sufficiently small (<ms) to be used in real-time control for XAO with currently available hardware even for extremely large telescopes.

We propose a new method of blind deconvolution (BD) suitable for the recovery of truncated blurred images. In practical cases we often encounter truncated images, but they are difficult to treat properly with conventional BD algorithms, since nontruncated blurred images are assumed. The proposed method is based on simulated annealing (SA), to which we introduce an adaptive masking process. After a truncated segment is cut out from the blurred image, the sharp edges surrounding the segment are removed by multiplying the initial window, whose shape is determined with the initial estimate of the point spread function (PSF). A more appropriate estimate for the PSF is obtained from the masked segment based on the SA algorithm, and according to the new estimate of the PSF, the shape of the window is varied. These procedures are iterated until we finally obtain the best possible estimate of the PSF. In the last step, the entire blurred image is deconvolved with the final estimate of the PSF to obtain the recovered image. The effectiveness of the proposed method is confirmed by computer simulation.

It is well known that using two or more deformable mirrors (DMs) can improve the compensation ability of an adaptive optics (AO) system. However, to keep the stability of an AO system, the correlation between the multiple DMs must be suppressed during the correction. In this paper, we proposed a slope-based decoupling algorithm to simultaneous control the multiple DMs. In order to examine the validity and practicality of this algorithm, a typical woofer-tweeter (W-T) AO system was set up. For the W-T system, a theory model was simulated and the results indicated in theory that the algorithm we presented can selectively make woofer and tweeter correct different spatial frequency aberration and suppress the cross coupling between the dual DMs. At the same time, the experimental results for the W-T AO system were consistent with the results of the simulation, which demonstrated in practice that this algorithm is practical for the AO system with dual DMs.

Adaptive optics (AO) enhances astronomical images by correcting perturbations in the optical system, resulting in complex point spread function (PSF) shapes, which provide critical information for optical design, calibration, and diagnostics. This paper presents a parametric model for an AO-corrected long-exposure PSF that adapts to complex shapes and various seeing conditions. A complementary estimation method is also introduced, addressing the inverse problem of estimating the atmospheric PSF or power spectral density (PSD) from observed PSF data. The model, based on Gaussian sum approximation (GSA), was tested using simulated PSF data from the OOMAO toolbox for Fried parameter (r0) values ranging from 0.10 to 0.20 m and on-sky data from the Very Large Telescope’s MUSE instrument. Two approaches were analyzed: (i) direct PSF estimation using GSA, and (ii) PSF estimation derived from atmospheric PSD estimation. Both methods were evaluated using root mean square error (RMSE). While approach (i) provided the best performance, approach (ii) produced more accurate atmospheric PSF estimates compared to the Moffat approximation and the Fetick method. Additionally, approach (ii) enabled the estimation of key physical parameters, such as the Fried parameter (r0) and residual AO variance ($\sigma _{AO}^2$). The proposed methods effectively estimate atmospheric PSF and PSD, even for complex PSF shapes influenced by spiders, central obstructions, or static aberrations.

In this paper we discuss the performance of low order adaptive optics (AO) systems. We present improved calculations of Strehl ratio achievable at various D/r<SUB>o</SUB> ratios for AO systems up to 7<SUP>th</SUP> order. Additionally we present calculations concerning the throughput of low order AO systems used in conjunction with spectrometers. We then show the results of a detailed simulation of the expected interaction of the University of Hawaii prototype AO system with spectrometers through the V,I,J,H,K optical bands. This simulation shows that AO spectroscopy performance should always be acceptable in the near IR (J,H,K), and often in the visible (V,I). We show the results of our initial efforts at deconvolving the AO PSF. These results indicate that deconvolution is likely to be quite widely applicable to AO images. Finally we discuss how it might be possible to extend the applicability and accuracy of deconvolution using the AO drive and wavefront error signals. The following results are for on-axis system, the effects of atmospheric anisoplanitism have not been addressed in this paper.

This paper presents the design of an Adaptive Optics (AO) system for retinal imaging applications. The development of retinal imaging systems allows for early diagnosis of eye diseases. Such systems can increase the quality of life of patients as well as curtail increasing health care costs through early eye disease detection and treatment. Until recently, AO systems have been prohibitively expensive and cumbersome. This has been mainly due to the size and cost of flexible membrane mirrors normally used as the aberration correction device. Recent developments in the technology of Microelectromechanical System (MEMS) based actuators allow the implementation of AO systems which would have been difficult to implement a few years ago due to exorbitant costs. The aim of this paper is to present the design of a compact and flexible low cost AO system using off the shelf components to measure and compensate for the aberrations of the eye. The design is based around the system’s main components which include a 52 channel magnetically actuated deformable membrane mirror, a Shack Hartmann wavefront sensor and a control system which runs on a single processor personal computer. All the components are commercially available. The use of the MEMS-based magnetically actuated mirror allows for increased resolution and force compared to conventional membrane mirrors designed mainly for use in astronomical applications. The performance of the closed-loop system is evaluated through experiments. Although designed as a diagnostic tool for eye diseases, such a system will find a number of applications in basic research in the visual sciences, including the study of microscopic structures in the living retina that could not be seen before. Optometrists, retinal surgeons, and ophthalmologists will also benefit from using such a system, through potential improvements on commonly used instruments such as phoropters and fundus cameras.

Solar Adaptive Optics (AO) shares many issues with night-time AO, but it also has its own particularities. The wavefront sensing is performed using correlations to efficiently work on the solar granulation as a reference. The field of view for that measurement usually is around 10". A sensor collecting such a wide field of view averages wavefront information from different sky directions, and the anisoplanatism thus has a peculiar impact on the performance of solar AO and MCAO systems. Since we are entering the era of large solar telescopes (European Solar Telescope, Advanced Technology Solar Telescope) understanding this issue is crucial to evaluate its impact on the performance of future AO systems. In this paper we model the correlating wide field sensor and the way it senses the high altitude turbulence. Thanks to this improved modelling, we present an analysis of the influence of this sensing on the performance of each AO configuration, conventional AO and MCAO. In addition to the analytical study, simulations similar to the case of the EST AO systems with FRiM-3D (the Fractal Iterative Method for Atmospheric tomography) are used in order to highlight the relative influence of design parameters. In particular, results show the performance evolution when increasing the telescope diameter. We analyse the effect of high altitude turbulence correlation showing that increasing the diameter of the telescope does not degrade the performance when correcting on the same spatial and temporal scales.

METIS is the Mid-infrared Extremely large Telescope Imager and Spectrograph, one of the first generation instruments of ESO’s 39m ELT. All scientific observing modes of METIS require adaptive optics (AO) correction close to the diffraction limit. Demanding constraints are introduced by the foreseen coronagraphy modes, which require highest angular resolution and PSF stability. Further design drivers for METIS and its AO system are imposed by the wavelength regime: observations in the thermal infrared require an elaborate thermal, baffling and masking concept. METIS will be equipped with a Single-Conjugate Adaptive Optics (SCAO) system. An integral part of the instrument is the SCAO module. It will host a pyramid type wavefront sensor, operating in the near-IR and located inside the cryogenic environment of the METIS instrument. The wavefront control loop as well as secondary control tasks will be realized within the AO Control System, as part of the instrument. Its main actuators will be the adaptive quaternary mirror and the field stabilization mirror of the ELT. In this paper we report on the phase B design work for the METIS SCAO system; the opto-mechanical design of the SCAO module as well as the control loop concepts and analyses. Simulations were carried out to address a number of important aspects, such as the impact of the fragmented pupil of the ELT on wavefront reconstruction. The trade-off that led to the decision for a pyramid wavefront sensor will be explained, as well as the additional control tasks such as pupil stabilization and compensation of non-common path aberrations.