Issue |
A&A
Volume 643, November 2020
|
|
---|---|---|
Article Number | A58 | |
Number of page(s) | 15 | |
Section | Astronomical instrumentation | |
DOI | https://doi.org/10.1051/0004-6361/202038679 | |
Published online | 03 November 2020 |
Joint estimation of atmospheric and instrumental defects using a parsimonious point spread function model
On-sky validation using state of the art worldwide adaptive-optics assisted instruments
1
Aix Marseille Univ, CNRS, CNES, LAM, Marseille, France
e-mail: olivier.martin@lam.fr
2
DOTA, ONERA, Université Paris Saclay, 91123 Palaiseau, France
Received:
17
June
2020
Accepted:
31
July
2020
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.
Key words: instrumentation: adaptive optics / methods: data analysis / techniques: image processing / atmospheric effects / methods: analytical
© O. Beltramo-Martin et al. 2020
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Adaptive optics (AO, Roddier 1999) is a game changer in the quest for high-angular resolution, especially for ground-based astronomical observations that face the presence of wavefront aberrations introduced by the atmosphere (Roddier 1981). Thanks to AO, the point spread function (PSF) delivered by an optical instrument is much narrower, by a factor up to 50 on the full width at half maximum (FWHM), than the seeing-limited scenario (Roddier 1981). Still, some correction residuals persist and render the AO PSF shape complex to model. Consequently, standard parametric models that reliably reproduce the seeing-limited PSFs, such as a Moffat function (Trujillo et al. 2001; Moffat 1969), become inefficient at describing the AO PSF. Moreover, contrary to seeing-limited observations, AO-corrected images suffer from the anisoplanatism effect (Fried 1982) that strengthens the spatial variations of the PSF on top of instrument defects. Determining the AO PSF is necessary for two major reasons. Firstly, understanding and accurately modeling the PSF morphology is key to diagnosing AO performance. From the PSF, we can identify the major contributors to the AO residual error (Beltramo-Martin et al. 2019; Ferreira et al. 2018; Martin et al. 2017). Secondly, the image delivered by an optical instrument depends on both the science object we want to characterize and the PSF. In order to estimate the interesting astrophysical quantities, one can either use a deconvolution technique (Fétick et al. 2020, 2019a; Benfenati et al. 2016; Flicker & Rigaut 2005; Mugnier et al. 2004; Fusco et al. 2003; Drummond 1998) or include the PSF as part of a model, as is performed in PSF-fitting astrometry/photometry retrieval techniques (Beltramo-Martin et al. 2019; Witzel et al. 2016; Schreiber et al. 2012; Diolaiti et al. 2000; Bertin & Arnouts 1996; Stetson 1987) and galaxy kinematics estimation (Puech et al. 2018; Bouché et al. 2015; Epinat et al. 2010) for instance.
Nevertheless, the PSF is not systematically straightforward to determine from the focal-plane image, particularly for observations made with a small field of view (FOV), or of extended objects, crowded populations or with a low signal-to-noise ratio (S/N). Alternative methods exist, such as PSF reconstruction (Beltramo-Martin et al. 2019, 2020; Wagner et al. 2018; Gilles et al. 2018; Ragland et al. 2018; Martin et al. 2016a; Ragland et al. 2016; Jolissaint et al. 2015; Exposito 2014; Clénet et al. 2008; Gendron et al. 2006; Veran et al. 1997), which is potentially very accurate (1% error level). However, this process lacks science verification and requires years to be fully implemented and operational as a stand-alone pipeline. The inherent complexity has inspired several efforts to produce simpler but still efficient reconstruction techniques (Fusco et al. 2020; Fétick et al. 2019b) that have been fully validated and are now in operation for Multi Unit Spectroscopic Explorer (MUSE) (Bacon et al. 2010) at the Very Large Telescope (VLT). The analytic AO PSF model described by Fétick et al. (2019b) has shown spectacular accuracy in reproducing the MUSE Narrow Field Mode (NFM) PSF as well as the PSF of the Zurich Imaging Polarimeter (ZIMPOL) (Schmid et al. 2018) that equips the Spectro-Polarimetic High contrast imager for Exoplanets REsearch (SPHERE) instrument (Beuzit et al. 2019) PSF and is an excellent candidate to rethink the way we perform PSF reconstruction, with a single flexible algorithm that complies with all kinds of AO correction. However, several questions arise as to its capabilities. Does this model really match any type of AO PSF, regardless of the observing conditions and even at very high Strehl-ratio (SR)? How much accuracy can we expect in a statistical sense and what are the fundamental limits? Can we improve the description of the PSF in the presence of static aberrations? The goal of this paper is to address these questions by (i) providing an exhaustive demonstration that this model complies with any type of AO system, telescope, and in optical and near infrared (NIR) wavelengths, (ii) assessing conditions for which the model is optimal for representing an AO PSF, and (iii) showing that this framework is robust enough to enable the joint estimation of AO residual and instrumental aberrations. In Sect. 2, we present the theoretical background and the model description as well as some upgrades we have implemented so as to treat new types of data and include static-aberration fitting. In Sect. 3, we present a statistical analysis of the model accuracy over 4812 PSFs with a discussion on the parameter estimates. Finally, in Sect. 4, we illustrate how the present model allows joint estimation of atmospheric conditions, AO performance and static aberrations, caused by the low wind effect (LWE) on SPHERE or segment cophasing error at Keck.
2. PSF model
This model was originally introduced by Fétick et al. (2019b) and validated on SPHERE/ZIMPOL and MUSE NFM data. In order to increase the range of applicability of such a model, we upgraded it by including (i) the model of the Apodized Lyot Coronagraph (APLC, Soummer 2005) used on SPHERE (focal-plane mask not considered) (ii) additional degrees of freedom to adjust static aberrations over a specific modal basis, which can be Zernike modes described on the pupil, tip-tilt and piston modes for each pupil area delimited by the spiders so as to describe the so-called LWE or piston modes for each pupil segment (Laginja et al. 2019; Leboulleux et al. 2018; Ragland 2018).
First of all, the model assumes the stationarity of the phase of the electric field in the pupil plane (Roddier 1981), which allows the system optical transfer function (OTF) to be split into a static part (telescope + internal aberrations) and an AO residual spatial filter
as follows,
where ρ is the separation vector within the pupil, λ is the imaging wavelength, and ρ/λ is the angular frequencies vector.
The static OTF derives from the instrument pupil mask 𝒫(r) and the optical path difference (OPD) map ΔStatic in the pupil plane, which gives for a monochromatic beam at wavelength λ
The SPHERE/IRDIS data we have treated (see Sect. 3.1 for a more complete description) were obtained during PSF calibration procedure using an off-axis stars. Therefore, the incoming beam was not getting through the focal-plane mask and the pupil mask model results from the multiplication of the apodizer (amplitude) function 𝒜 and the Lyot stop ℒ (Soummer 2005) as follows
The 2D functions 𝒜(r) and ℒ(r) are illustrated in Fig. 1.
![]() |
Fig. 1. From left to right: pupil plane apodizer and Lyot stop used during SPHERE/IRDIS observations with the N_ALC_YJH_S APLC (Vigan et al. 2010). |
In this paper, we split the static aberration contribution in two terms: the Non Common Path Aberrations (NCPA) calibration (non-coronagraphic) noted ΔNCPA (Jia et al. 2020; Vigan et al. 2019; Lamb et al. 2018; Vassallo et al. 2018; N’Diaye et al. 2016; Sitarski et al. 2014; Jolissaint et al. 2012; Sauvage et al. 2011; Robert et al. 2008; Mugnier et al. 2008; Blanc et al. 2003; Fusco et al. 2003), to which is added another contribution decomposed over a modal basis M to be specified. This latter could be simply Zernike modes over the whole pupil or per pupil area delimited by the spiders, so as to include potential LWE that may introduce severe asymmetries in the PSF structure. Analyses of SPHERE images (Milli et al. 2018; Sauvage et al. 2016, 2015) showed that the LWE mainly introduces piston, tip, and tilt differential aberrations between pupil areas separated by the spiders, which represents 12 parameters to be adjusted. Also, for segmented pupil telescopes like Keck telescopes, this basis can be segment piston or petal modes (Ragland 2018). We analyze the capacity of this model to identify such aberrations in Sect. 4. The static OPD is calculated following
where μStat(k) and Mk are respectively the coefficient and 2D shape of the kth mode over nm modes. To include segment piston errors, Mk must be defined as the pattern formed by the kth segment (36 in total for the Keck pupil) and ak is the piston value in meters. Eventually, one may estimate both the atmospheric parameters and static coefficients and we illustrate such a joint estimation on Keck data in Sect. 4.1.
Finally, the description of the AO residual spatial filter in Eq. (1) relies on two assumptions, which are, (i) the exposure is infinitely long, and (ii) the residual phase is a Gaussian statistical process. Thus,
is fully described by the residual phase power spectrum density (PSD) as follows (Roddier 1981):
where 𝒟AO, ℬAO, and 𝒲AO are the residual phase structure function, autocovariance function, and PSD respectively, k is the spatial frequencies vector and ℱ[x] is the 2D Fourier transform of x. In Fétick et al. (2019b), the PSD is described as a split function to separate the AO-corrected spatial frequencies from the uncorrected frequencies that follow a Kolmogorov’s −11/3 power law (Kolmogorov 1941)
where k = |k|, r0 is the Fried parameter (Fried 1966), kAO the equivalent AO cut-off frequency beyond which the AO correction no longer occurd, C a constant value and ℳ an asymmetric Moffat function (Moffat 1969) that depends on five shape parameters A, αx, αy, θ, and β and on the spatial frequency vector k = (kx, ky) as follows,
where ψ is a normalization factor that ensures that the integral of the Moffat function over the AO-corrected area is given by A: (Fétick et al. 2019b)
Contrary to the classical use of a Moffat function in astronomy (Trujillo et al. 2001; Moffat 1969), we stress that the Moffat function used in the present model serves in the description of the PSD and not the focal-plane PSF. This latter is deduced by the inverse Fourier transform of the OTF given in Eq. (1). Consequently, the PSF FWHM is actually driven by parameters A and r0 rather than αx, αy, and β. In order to really capture the physical meaning of this model, below we present a description of each of the seven parameters that the model relies on and their impact on the PSF:
– r0 constrains the uncorrected PSD that corresponds to the PSF wings, that is, the part of the PSF that remains untouched by any AO correction. For a system with nact × nact actuators to compensate for the incoming wavefront, this breaking occurs in the focal plane at approximately nact × λ/(2D), with D the primary mirror diameter and for a Nyquist-sampled PSF (Roddier 1999).
– C is a constant value that changes the gap between the AO-corrected part and the uncorrected high-spatial frequencies. This parameter plays a major role near the AO cutoff in both the PSD and PSF planes, and compares to the aliasing error that contaminates wavefront measurements due to the WFS discrete sampling (Bond et al. 2018a; Correia & Teixeira 2014; Jolissaint 2010; Rigaut et al. 1998).
– A is the total energy in nm2 contained in the Moffat PSD model (constant C not included) and thus characterizes the AO residual wavefront error that connects to the PSF SR, that is, the intensity of the PSF peak compared to the diffraction-limit scenario.
– αx, αy, and θ govern the elongation and skewness of the PSD as well as the direction of the elongation thanks to the angle θ. The PSD FWHM is proportional to αx, αy parameters; for the same amount of energy, a PSD with a larger FWHM indicates that the PSD flattens and the correction homogenizes across spatial frequencies.
– β represents the asymptotic slopes of the PSD at large spatial frequencies within the AO correction radius. From Fourier analysis of AO performance (Bond et al. 2018a; Correia & Teixeira 2014; Jolissaint 2010; Rigaut et al. 1998), we know that the PSD pattern introduced by the wavefront measurement noise follows a k−2 power law for WFSs sensitive to the first order derivative of the wavefront such as the Shack-Hartmann. Asymptotically, we would have β = 1 for a Shack-Hartmann based AO system if the measurement noise is the only error in the wavefront reconstruction process. Moreover, in the situation of an AO system correcting all atmospheric aberrations with the same relative level (which does not occur as the sensitivity of the wavefront sensor is aberration-dependent Fauvarque et al. 2016), we would obtain β = 11/6 ≃ 1.83 according to the von Kármánn expression of the atmospheric PSD (von Karman 1948). Nevertheless, larger values of β can be observed. Indeed, for very large values of β, the Moffat distribution converges towards a Gaussian shape (Trujillo et al. 2001), as does the PSF consequently. We expect this behavior with poor tip-tilt correction for instance or in the presence of strong telescope wind-shake or telescope/dome vibrations. As a summary, β should usually range between 1 and 1.83 for nominal AO correction but can reach higher values in the presence of non-atmospheric aberrations.
We illustrate in Fig. 2 azimuthal profiles of PSDs and corresponding PSFs for various sets of model parameters. This figure shows clearly how r0 drives the PSF halo and how the parameter A increases the area underneath the PSD curve, which corresponds to the residual wavefront error. We also observe that a larger value of β sharpens the PSD. Knowing that the total energy remains constant, this situation corresponds to less efficient low-order modes compensation. Finally, a larger value of α increases the PSD FWHM and distributes more energy over the highest AO-correct modes. The shape of the PSF is mainly conditioned by parameters A and r0, which are the two most significant parameters of this model. Parameters C, α, θ, and β are secondary parameters that shape the PSD in order to more accurately reproduce the PSF structure in various observing conditions and especially in the presence of sub-optimal AO control or non-atmospheric aberrations. Henceforth, this model is a combination of a PSF determination tool and an AO diagnosis facility gathered up together into a single and parsimonious analytical framework.
![]() |
Fig. 2. Azimuthal average in log-log scale of the PSD (left) and the PSF (right) obtained with four different settings for the model parameters, assuming that the PSD is symmetric, i.e. αx = αy = α and θ = 0. No static aberrations were added here and the parameter C was set to C = 10−4 rad2m2. The PSF maximum is normalized to the SR value. |
3. Model versus reality: statistical validation
3.1. Overview of data
We aim to provide the most exhaustive on-sky validation of the PSF analytical and parsimonious model proposed by Fétick et al. (2019b). To achieve this goal, we collected 4812 PSFs obtained from different observatories, AO systems, AO modes and spectral bands covering optical and NIR wavelengths. The GALACSI/MUSE and GeMS/GSAOI instruments deliver simultaneous PSFs at different wavelengths and field positions respectively. If we consider a single PSF for each observation they produce, we obtain a total of 1880 PSFs. However, although the realizations of atmospheric residual wavefront are not independent from a PSF to another within the same observation, these data allow to test the model in the presence of chromatic and field-dependent instrumental aberrations. In order to conserve the same fitting process for all data, we have fitted each PSF with the same starting point and regardless the results obtained on another PSF acquired during the same observation. A description of each data set is given below
SPHERE/ZIMPOL@VLT. SPHERE is a facility at the VLT that was installed in 2015 and relies on an extreme AO (EXAO) system (Sphere Ao for eXoplanets Observations; Fusco et al. 2016; Sauvage et al. 2015) to deliver a very efficient atmospheric correction level to three science instruments, which are ZIMPOL (Schmid et al. 2018), IRDIS (Dohlen et al. 2008), and IFS (Claudi et al. 2014). This set of data is composed of two subsets of PSFs from SPHERE/ZIMPOL. The first (26 PSFs) was obtained from observations of NGC 6121 in 2019 (Massari et al. 2020; Beltramo-Martin et al. 2020) with the ZIMPOL V-filter (central wavelength 554 nm, width 80.6 nm) and a pixel scale of 7.2 mas pixel−1 in the context of technical calibrations1 granted after the 2017 ESO calibration workshop2. The second set of 18 PSFs was acquired in 2018 with the N_R-filter (central wavelength 645.9 nm, width 56.7 nm) and a pixel scale of 3.6 mas pixel−1 during the ESO Large Program (ID 199.C-0074, PI P. Vernazza) (Vernazza et al. 2018). The data were reduced using the SPHERE Data and Reduction Handling pipeline (DRH) to extract the intensity image, subtract a bias frame, and correct for the flat-field. These data are particularly useful for testing the model in optical wavelengths and under strong atmospheric residual regime.
SPHERE/IRDIS@VLT. A total of 237 PSFs were obtained from the SPHERE Data Centre client (Delorme et al. 2017) using the Keyword Frame type set to IRD_SCIENCE_PSF_MASTER_CUBE. These PSFs were acquired over the last five years with the N_ALC_YJH_S APLC (Vigan et al. 2010) during PSF calibration and with a pixel scale of 12.5 mas pixel−1. These data were collected using the dual band filters DB_H23, and DB_K12, and they are useful for testing the model under very high SR regime and for validating the LWE retrieval presented in Sect. 4.2.
Keck AO/NIRC2@Keck II. We obtained 355 PSFs using the narrow field mode of NIRC2 with a pixel scale of 9.94 mas pixel−1 and using the Fe II and K cont filters. The Keck AO system on the Keck II telescope was operated in single conjugated AO (SCAO) mode using a natural (Wizinowich et al. 2000) or an on-axis laser (Wizinowich et al. 2006) guide star. These PSFs were obtained during PSF reconstruction engineering nights in 2013 and 2017 (Ragland et al. 2018, 2016) and such data are especially useful to validate the model under the influence of remaining piston cophasing errors, as we discuss in Sect. 4.1.
SOUL/LUCI@LBT. These data were delivered in 2020 from the to two LUCI NIR spectro-imagers assisted with the pyramid SCAO (PSCAO) SOUL AO system (Pinna et al. 2016) driven by a pyramid WFS (Ragazzoni & Farinato 1999). These PSFs were acquired using the H filter (1.653 μm) and with a sampling of 15 mas pixel−1. Although few (11 for our purpose) data sets have been obtained so far in comparison to others systems, these PSFs pave the way to a demonstration of PSF determination strategies in the presence of a pyramid WFS, which will be the baseline for the SCAO mode of HARMONI (Thatte 2017), MICADO (Davies et al. 2016) and METIS (Hippler et al. 2019) on the future 39 m Extremely Large Telescope (ELT).
CANARY/CAMICAZ@WHT. CANARY (Gendron et al. 2010; Myers et al. 2008) was designed as a pathfinder for demonstrating the reliability and robustness of the multi-object AO (MOAO) concept proposed for assisting very large field (> 1′) multi-object spectrographs, such as MOSAIC for the ELT (Hammer et al. 2016). Thanks to 26 nights of commissioning, tests, and validation, we collected 1268 PSFs in H-band and with a sampling of 30 mas pixel−1 using the NIR detector CAMICAZ (Gratadour et al. 2014). These data were obtained during phase B (Martin et al. 2017; Morris et al. 2014, 2013; Martin et al. 2013), which was dedicated to the demonstration of MOAO relying jointly on four (Rayleigh) laser guide star (LGS) and up to three natural guide star (NGS) in a 1’ FOV to perform the tomography using the Learn & Apply technique (Laidlaw et al. 2019; Martin et al. 2016b; Vidal et al. 2010). Among this quite large set of data, we have 522 PSFs in SCAO mode, 128 PSFs in ground layer AO (GLAO) mode and 618 MOAO PSFs. Such an archive is useful for testing the model on a 4.2 m class telescope and on tomographic and laser-assisted AO-corrected-PSFs.
GALACSI/MUSE NFM@VLT. MUSE is the ESO VLT second-generation wide-field integral field spectrograph operating in the visible (Bacon et al. 2010), covering a simultaneous spectral range of 465–930 nm and assisted by the ESO Adaptive Optics Facility (AOF, Oberti et al. 2018; Arsenault et al. 2008) including the GALACSI module (Ströbele et al. 2012). We focused our analysis on the narrow field mode (NFM) of MUSE, which delivers a laser tomography AO (LTAO)-corrected field covering a 7.5″ × 7.5″ FOV, providing near-diffraction-limited images with a sampling of 25 mas pixel−1. These PSFs were obtained during the commissioning phase in 2018 and are particularly useful for analyzing the model outputs with respect to the wavelength and demonstrating that the model also complies with under-sampled PSFs. Also, we performed a spectral binning to reach a spectral width of 5 nm (91 PSFs per cube if we remove the notch filter wavelengths), resulting in a total of 1986 PSFs.
GEMS/GSAOI@GEMINI. The Gemini South Adaptive Optics Imager (GSAOI) is a NIRd camera that benefits the correction provided by the Gemini Multi-conjugate Adaptive Optics (MCAO) System (GeMS, Neichel et al. 2014; Rigaut et al. 2014) on Gemini South. It delivers near-diffraction-limited images in the 0.9 – 2.4 μm wavelength range in a large FOV of 85″ × 85″. From the Gemini archive (Hirst & Cardenes 2017), we assembled a catalog of 911 isolated and nonsaturated PSFs extracted from 60 images of Trumpler 14 acquired in 2019 and using the J, H, K, and Brγ filters with a pixel scale of 20 mas pixel−1 (P.I. M. Andersen, observation program: GS-2019A-DD107). These data are particularly useful for investigating model parameters variations across a MCAO-corrected FOV.
In total, we created a dictionary of 4812 PSFs covering several observatories, AO correction types, optical and NIR wavelengths, as summarized in Table 1. Having this diversity of data is important for spanning the full range of possible AO correction levels and assessing which conditions must be met to achieve an accurate PSF representation.
Summary of PSFs obtained and processed for the analysis presented in this paper.
3.2. PSF fitting
In order to fit the model over the image and retrieve the associated parameters, we followed the same strict process for each of the 4812 PSFs in our dictionary, which is described as follows
– Define a model of the image including the PSD degrees of freedom (no static aberrations retrieval) μAO as well as four additional scalar factors δx, δy, γ, and ν in order to finely adjust the PSF position, flux, and constant background level,
The image model is calculated over a given number of pixels that is 10% larger than the on-sky images so as to mitigate aliasing effects due to the Fast Fourier Transform algorithm. The size of the on-sky image support is instrument-dependent and truncated in order to mitigate the noise contamination. When possible, we crop the on-sky image to conserve a FOV up to twice the AO cutoff.
– Define a criterion to minimize
where dij is the (i, j) pixel of the 2D image and Wij the weight matrix defined by
The weight matrix accounts for the noise variance, that is, both photon noise and read-out noise, and allows us to maximize the robustness of the fitting process (Mugnier et al. 2004).
– Perform the minimization using a nonlinear and iterative recipe based on the trust-reflective-region algorithm (Conn et al. 2000). We did not use specific regularization techniques on top of the weighting matrix as we have taken care of selecting good S/N images.
3.3. PSF model accuracy
In Fig. 3, we present a visual comparison of on-sky PSFs, best-match model and residual map for the seven instruments considered in this analysis. The model adapts to any kind of AO correction; the maps of residuals are visually similar to first order among all systems. On the Keck AO/NIRC2 image, we see some structures that are static speckles, probably introduced by a remaining cophasing error, and residual NCPA as illustrated in Sect. 4.1. On SOUL/LUCI and CANARY/CAMICAZ images, we also see a persistent pattern that can be explained by the presence of static aberrations not included in the model and the exposure time that was not sufficiently long (a few seconds of exposure) to average the atmospheric speckles.
![]() |
Fig. 3. Two dimensional comparison in log scale and over 64 × 64 pixels of (top) observed PSFs, (middle) fitted model, and (bottom) residual map for multiple types of AO systems and instruments working in either visible or NIR wavelengths. From left to right: SPHERE/ZIMPOL, SPHERE/IRDIS with the apodizer and the Lyot stop (no focal plane mask), Keck AO/NIRC2, SOUL/LUCI, CANARY/CAMICAZ in MOAO mode, GALACSI/MUSE NFM, and GEMS/GSAOI. All PSFs have been normalized to the sum of pixels. |
In Fig. 4, we illustrate the SR and FWHM obtained from the fitted model versus sky image-based estimates. The same algorithms (OTF integral for the SR, interpolation+contour for the FWHM) were used to calculate these metrics regardless of the nature of the data, either sky image or model. For all systems, we observe a remarkable correlation and a similar dispersion on both SR and FWHM and for all observing conditions, AO modes and AO correction levels. For very high SR values obtained with SPHERE/IRDIS, we start observing an overestimation of the SR suggesting that the model cannot fully retrieve some patterns in the image. At this level of correction, the instrumental contribution (residual NCPA, LWE for instance) of the PSF may dominate the PSF morphology while they are not included in the model, as for all 4812 handled data sets. Consequently, in order to guarantee an SR and FWHM accuracy at a few percent, it is not necessary to include a precise model of these instrumental defects for AO systems delivering SR up to ∼80%. We confirm this assumption in Sect. 4. By calculating the relative difference ((xmodel − xsky)/xsky) on SR and FWHM values over all the 4812 data sets, we measure a bias and a standard deviation (std) value of 0.7% and 4.0% on the SR estimation and −0.8% and 4.6% on the FWHM estimation. These numbers indicate that there is a marginal performance overestimation of 1% from the model (larger SR, lower FWHM), which fits the measurements uncertainty envelopes. As a result, this model achieve a PSF recovery at a 4% level.
![]() |
Fig. 4. Image SR/FWHM versus the same metrics retrieved on the fitted image using the same estimation process and for the 4812 PSFs treated for this analysis. Error bars on the SR are obtained from calculations presented in Martin et al. (2016a). Error bars on the FWHM are given from the contour estimation on interpolated images that are oversampled by up to a factor four to quantify the FWHM more accurately. |
Table 2 provides statistical results of the estimated SR and FWHM PSFs for all instruments, including the median values for all observations and the Pearson correlation factor. As suggested by Fig. 4, there is no specific bias, except for SPHERE/IRDIS for reasons mentioned above, as well as for SOUL/LUCI owing to the small amount of data we have access to so far. Overall, we conclude that (i) the model becomes biased for very high SR observations (SR > 80%), calling for the introduction of instrumental defects to improve the model accuracy, (ii) the SR is estimated with a 1-σ precision of 1%, and (iii) the FWHM is estimated with a 1-σ precision of 3 mas, which correspond to approximately to one-fifth down to one-tenth the pixel scale depending on the instrument.
Individual statistics per system and imaging wavelength on SR and FWHM estimates.
3.4. Exploitation of the model outcomes for diagnosing observing conditions and AO performance
As discussed above, fitting the shape of the residual phase PSD allows to retrieve atmospheric parameters and AO performance. Therefore, the goals in this section are to (i) confirm that the output parameters r0 and A follow expected trends and give confidence in the retrieval process, and (ii) provide statistics on α, β parameters to asses which values they should reach for a nominal AO correction and thus discriminate situations of sub-optimal AO correction. We have excluded the parameter θ from this analysis as it corresponds to a PSF orientation only.
3.4.1. Primary parameters estimates
The seeing is estimated from the PSF wings fitting relying on a Kolmogorov expression of the PSD. This measurement technique has proven to be robust (Fétick et al. 2019b) as it consists in determining a single parameter from a significant number of pixels. Thanks to the large redundancy across the pixels of the spatial signature that the algorithms is seeking out, this approach still gives meaningful results with moderate S/N in comparison to external profilers (Fétick et al. 2019b). However, contrary to these latter, this focal-plane-based seeing determination includes all turbulence effects that contribute to impact the PSF wings, such as the dome seeing (Lai et al. 2019; Conan et al. 2019), which the external profilers are not sensitive to as they are apart from the dome. Consequently, estimating the r0 from the focal-plane image allows to diagnose more accurately the AO performance in comparison to an external profiler.
The target here is to verify that the retrieved seeing values are consistent with what we know about the observing sites. The first verification we made is to analyze statistics on the seeing estimates presented in Table 3. Given that the Keck data were acquired in February, August, and September seasons, the median seeing at Mauna Kea is consistent with the literature (Ono et al. 2016; Miyashita et al. 2004). Similar observations can be drawn for La Palma by comparing to either CANARY telemetry data (Martin et al. 2016b) or the RoboDIMM (O’Mahony 2003). We also find consistency with analysis by Masciadri et al. (2014) for the Paranal site and by Tokovinin & Travouillon (2006) for Cerro Pachón.
Median and 1-σ standard-deviation of seeing values retrieved from the PSF-fitting process.
In addition, we highlight in Fig. 5 the r0 estimates at zenith and 500 nm obtained on a single image of GeMS/GSAOI and a single cube of GALACSI/MUSE. For the GEMS/GSAOI case, we have r0 measurements obtained from several PSFs distributed over the field. Given that r0 is assessed from the PSF wings that contain the non-AO-corrected spatial frequencies, we expect r0 estimates to be independent of the PSF position in the FOV. Similarly, we multiplied r0 estimates obtained on GALACSI/MUSE images by a factor (500/λc)6/5, where λc is the central wavelength of the considered image in the data cube, in order to disable the theoretical dependency of r0 with respect to λ.
![]() |
Fig. 5. r0 estimates at zenith and 500 nm with respect to the GEMS/GSAOI PSF field position (left) and with respect to the PSF GALACSI/MUSE wavelength (right). |
Using our model, we retrieve r0 as 14.2 cm ± 0.4 and 10.4 cm ± 0.2 cm for GEMS/GSAOI and GALACI/MUSE, respectively. The achieved 3% 1-σ error of our estimates shows that this method is robust and precise compared to telemetry-based techniques that typically reach 10% precision (Jolissaint et al. 2018).
For completeness, Fig. 6 shows histograms for the A parameter (the wavefront error) and the function A = f(r0) for the specific case of CANARY working in SCAO (i.e., the largest SCAO PSF dictionary we have). We clearly see that the wavefront standard-deviation error follows a law as we expect from the von Kármán PSD (von Karman 1948), which is proportional to
for a SCAO system whose on-axis PSF is not influenced by the atmospheric turbulence profile. This illustrates the agreement between the different retrieved parameters.
![]() |
Fig. 6. Left: histograms for the A parameter for each AO system. Right: wavefront error obtained on CANARY SCAO mode with respect to the retrieved r0. The solid line gives the trend in |
Moreover, the histogram of A estimates given in Fig. 6 shows that the retrieved wavefront errors vary within a meaningful range: the SPHERE AO system performs better than others, as expected from such an extreme AO system, with a median of 140 nm, which includes all the aberrations that blur the AO-corrected part of the PSF. The present SPHERE results gather PSFs obtained with SPHERE/ZIMPOL, acquired with a AO loop running at low frequency (300 Hz, mV ≃ 11), and SPHERE/IRDIS fir which 25% were contaminated by a strong LWE. Therefore, this wavefront error of 140 nm seems consistent with the literature (Sauvage et al. 2016). The Keck AO system reaches 230 nm, which also complies with Wizinowich (2012). From GALACSI/MUSE NFM images, we retrieve 285 nm of wavefront error, which agrees with the Oberti et al. (2018) analysis Moreover, the GALACSI system was not yet fully optimized to operate in LTAO and a recent acquisition in 2019 already showed improvements on SR. For CANARY/CAMICAZ, we obtain a wavefront errors of 320 nm, 396 nm and 417 nm in SCAO, MOAO, and GLAO mode respectively, which compares well with Martin et al. (2017) and Vidal et al. (2014). Finally, GeMS/GSAOI data unveil a median wavefront error of 570 nm from PSFs extracted from the field at 20″ up to 60″ off-axis, which corresponds to a SR of 7.8% at the edges of the field and complies with analysis of the performance of GeMS (Neichel et al. 2014; Rigaut et al. 2012).
We stress that this value of wavefront error is not deduced from the image SR but directly from the integral of the PSD, which can be determined from the model parameters (Fétick et al. 2019b). This corresponds to the real mathematical definition of the wavefront error and is not influenced by the Maréchal approximation (Parenti & Sasiela 1994) that is biased at low SR. Consequently, this model is also a robust focal plane-based wavefront error estimator.
3.4.2. Secondary parameters estimates
We have emphasized in Sect. 2 that parameters α and β govern the PSD shape and assist in the AO correction diagnosis.
Figure 7 presents histograms for each system; the β parameter has a relatively narrow histogram with a peak located at β = 1.82 and a 1-σ standard deviation of 0.6, while the α histograms seem particularly instrument-dependent with values from 0.001 rad.m for efficient AO correction to 1 rad.m for marginal AO correction (e.g., in the bluest visible wavelengths of MUSE). Our first conclusion here is that an optimized AO system should provide a PSD with a β parameter comprised between 1 and 1.9, as discussed in Sect. 2. For larger values of β, there is an excess of residual error into low-order spatial frequencies. For instance, we see that on Keck AO/NIRC2 images the median β increases from 1.7 up to 2.1 in NGS and LGS modes, respectively, indicating the presence of additional low-order modes introduced here mainly from the focal anisoplanatism (Wizinowich 2012; van Dam et al. 2006). In addition, we have observed cases with β > 2.5 on Keck data due to a strong wind-shake effect that was enlarging the PSF FWHM by a factor three compared to nominal performance. Both α and β parameters grow up in the presence of strong wind shake and are wavelength dependent as we see on GALACSI/MUSE histograms. However, from GEMS/GSAOI data analysis, we do not find clear trend with respect to the field position, which may be owing to the uniform correction across the field provided by the MCAO mode of GEMS.
![]() |
Fig. 7. Histograms for the β and α parameters for each system, as well as the averaged distribution. |
We illustrate here that those two parameters carry additional and relevant information on the AO correction. However, the exact connection with the AO status is not straightforward to identify. To enable this identification, we are currently developing a convolutionnal neural network (Herbel et al. 2018) that we train to estimate the model parameters from the AO control loop data, such as wavefront measurements. As we can collect a very large amount of data for the purpose of estimating a few tens of parameters, solving this problem is definitely achievable with dedicated simulations and on-sky data that we will continue to be collected regularly among observatories.
3.5. Discussion on the influence of exposure time and spectral bandwidth
One of the major assumptions in the model proposed by Fétick et al. (2019b) concerns the infinitely long exposure hypothesis. This model is therefore not capable of reproducing atmospheric speckles that average when taking a sufficiently long exposure. As the method relies on the second-order statistical moment of the residual phase, the time necessary to achieve a convergence of the PSD shape is highly dependent on atmospheric parameters but may be achieved in few seconds (Martin et al. 2012). Thanks to SOUL/LUCI data, we analyzed the PSF accuracy when fitted on short exposure images with integration times from 0.157 s to 60 s. We find that the model matches the PSF down to 1–2 s of exposure and for a PSF acquired at 1.6 μm, below which the parameters estimation begins to degenerate as illustrated in Fig. 8 (left) through the average of absolute parameters variations. We have also noticed that the PSF shape remains stable from few seconds exposure, which explains the stability of retrieved parameters.
Moreover, using GALACSI/MUSE data, we tested the model on a polychromatic image, from 2.5 nm up to 470 nm of spectral width by binning monochromatic PSF observed simultaneously with MUSE. When compensating for the beam dispersion (recentering PSFs and then stack), the parameter estimation does not deviate by more than 4% over the whole spectral width span as presented in Fig. 8 (right). Consequently, the presence of chromatic static aberrations and atmosphere chromaticity do not prevent the model from very accurate characterization of the PSF. When including the beam dispersion (no recentering before stacking), which is 350 mas (14 pixels) for the largest spectral width, the parameter estimation degrades up to 10% and goes beyond the 4% threshold after 300 nm of spectral width. As a result, the model remains robust and reliable, even in the presence of such a strong beam dispersion.
![]() |
Fig. 8. Left: average of absolute retrieved parameters with respect to the exposure time. The parameters are normalized by the value obtained from the long-exposure PSF (60 frames). Error bars are assessed by averaging results over the SOUL/LUCI data sets. Right: average retrieved parameters normalized by the median value over the whole span with respect to the spectral width. The envelope shows the ±4% of variations around the median value and the spectral bandwidth is systematically centered around 700 nm. Error bars are assessed by averaging results over the GALACSI/MUSE data sets. |
4. Joint atmospheric, AO performance and static aberrations retrieval
The previous section illustrates that this model accurately and reliably reproduces the AO residual component of the PSF. Considering that the model is parsimonious (7 AO parameters and 4 photometric and astrometric parameters) and the large amount of pixels (10 000 for a 100 × 100 image) that the model-fitting may rely on, one can attempt to retrieve more degrees of freedom, such as static aberrations.
In order to estimate the static coefficients, we replace in the criterion given in Eq. (10) by
, which becomes a parametric OTF term. Besides, determining the wavefront from a single long-exposure image suffer from strong degeneracies by essence and this approach does not overcome this problem by using a non-linear minimization algorithm to estimate the problem solution. We distinguish two sorts of degeneracy
– Sign ambiguity : this occurs when trying to fit pair modes for which the PSF is not sensitive to the sign (Mugnier et al. 2008). Therefore, we concentrate the static aberrations retrieval on piston, tip and tilt modes only. Moreover, we bound the solution domain between −λ/2 and λ/2 in order to overcome phase wrapping issues with piston modes.
– Local minima: two different wavefront patterns may produce a similar but slightly different PSF and make the PSF-fitting algorithm converge towards a local minimum of the criterion. To assess how much this problem affect our algorithm, we have compared the PSF model with various conditions of piston, tip and tilt levels in order to highlight possible combinations (10 000) of modes that could produce a similar PSF. From a vector of aberrations, we have tested different permutations of the elements of this vector and in 99% of cases, the relative PSF variation is larger than 1%, which is large enough to change the structure of the PSF and retrieve the correct wavefront map, as long as the S/N is sufficient (> 50). We also rely on an analysis from Gerwe et al. (2008) that shows that for a fit of 35 Zernike polynomials on a segmented pupil, there are no local minima as long as the initial guess for the static coefficient remains within ±0.2λ = 330 nm rms in H-band from the optimal solution. We are in a situation where the static aberrations we attempt to retrieve are already mitigated from the Keck telescope active control (Chanan et al. 1988) and the VLT spiders coating (Milli et al. 2018). Consequently, the aberrations level we must retrieve remains sufficiently weak to mitigate the presence of local minima. At a level greater than 300 nm (Keck AO residual is 280 nm) would impact drastically the PSF and degrade the telescope science exploitation so much that this aberration would be necessarily visible and mitigated as much as possible.
Henceforth, we are in a good situation to estimate piston, tip and tilt static modes on Keck and VLT images obtained on bright star.
4.1. Keck cophasing error retrieval
The analysis presented in this section focuses on a smaller sample of 129 Keck AO/NIRC2 data acquired in NGS mode and with a high SR (> 40%) in 2013 (Ragland et al. 2016). The residual NCPA map was calibrated at the beginning of the night. We compared model-fitting performance in three situations: (i) a fit of PSD and photometric and astrometric parameters (11 values), (ii) a fit of these parameters when including, following Eq. (2), the static aberration map calibrated at the beginning of the night and, (iii) a joint adjustment of the PSD, photometry/astrometry and the 36 piston values corresponding to each Keck pupil segment (47 parameters in total), including the calibration static aberration map.
Figure 9 shows the PSFs and the residual maps produced from the three model-fitting strategies considered in this section and obtained for three cases acquired in February, August, and September, 2013. We firstly conclude that accounting for the residual NCPA static aberrations improves the model fitting as revealed by the residual map. This is also confirmed by the results shown in Table 4 which gives the median and 1σ standard deviation of SR/FWHM estimates over the 129 data samples as well as the mean square error (MSE) determined from the residual map. We firstly observe that biases and dispersion values on the SR are larger than in Table 2, owing to the data sample we consider in the present analysis and that gathers only high SR data for which the residual atmospheric aberrations are weaker. This advocates for including the calibrated static aberrations map into the model in order to mitigate this bias in the model. From Fig. 9, we also observe that these static aberrations manifest mainly as static speckles in the PSF and consequently contribute marginally to shape the PSF core, which explain why the statistics on the FWHM estimates are not influenced by the reduction of the data sample. Moreover, the results shown in both Fig. 9 and Table 4 highlight that the model accuracy is slightly improved by including the NCPA residual map in the model, with still 50% improvement on the MSE. Furthermore, the cophasing map retrieval allows to drastically enhance the model-fitting results, for example by up to a factor three on the MSE compared to the first scenario, as also illustrated in Fig. 9 and Table 4. The fact that the model efficiently reproduces the speckles observed on sky images suggests that these speckles are mainly introduced by cophasing errors due to the telescope.
![]() |
Fig. 9. Comparison in log scale of on-sky PSFs and fitted model using the three implemented strategies. Each column correspond to a single data sample acquired at different epochs. For each model configuration, we report both the adjusted PSF and the residual map. This figure illustrates that the cophasing map retrieval improves drastically and consistently the residual on top of the gain brought by including the NCPA residual map. |
Median and 1-σ dispersion of SR and FWHM errors and mean square error obtained over the 129 H-band and high SR (> 40%) Keck/NIRC2 images treated in this analysis for the three implemented strategies: (i) PSD/stellar parameters retrieval, (ii) including the residual NCPA map in the model, and (iii) also including the estimation of Keck II cophasing errors.
In addition, in order to further demonstrate the strength of this framework, we present in Fig. 10 the retrieved cophasing maps for three specific cases. The corresponding NCPA map used in the model is presented in Ragland et al. (2016). We unmistakably observe the presence of a stair mode, that introduces a constant step between adjacent segment rows or columns and that has been already highlighted from data acquired in 2017 using an alternative technique (Ragland 2018).
![]() |
Fig. 10. Retrieved cophasing maps in nm from PSFs acquired during three different epochs in 2013. The telescope elevation and azimuth were respectively (79,240), (85,−82) and (43,240) from left to right. |
Moreover, each of those three images was obtained with different elevation-azimuth telescope positions that are reported in Table 5, which shows that the wavefront standard deviation value is very close to observations reported by Ragland (2018) for a 45° telescope elevation. Also, the stair case peak-to-valley (PV) energy increases with low telescope elevation and rotates with respect to the telescope azimuth. This observation reveals a connection between the presence of this stair mode and global flexure over the primary mirror that is controlled to compensate for the gravity effect. Deeper analyses will confirm the presence of such a residual flexure on the Keck telescoped and whether or not this framework now to allow us to correct for this stair mode to generally improve the scientific exploitation of the Keck, but also future segmented telescopes.
Summary of peak-to-valley (PV) and 1-σ standard deviation of retrieved cophasing error map regarding the corresponding telescope elevation/azimuth.
Measuring the piston map from the focal-plane image is feasible from results obtained with the present analysis, but necessitates to have a star bright enough in the field to calibrate these aberrations. We must ensure that the shape of the primary mirror is controlled without the aid of a focal-plane-based technique in order to achieve the best image quality regardless of the observed field. We are again in a situation where we have to retrieve a few parameters (47) from a potentially large amount of data delivered by segment position sensors, temperature, pressure and humidity sensors, atmospheric parameters near the dome and even AO telemetry, which can deliver information about the piston error when using a pyramid WFS (Bond et al. 2018b; Schwartz et al. 2018). It is not easy to investigate the connection between this ensemble of data and the piston map, and the use of neural networks for this purpose will be explored.
4.2. SPHERE low wind effect retrieval
The analysis presented in this section focuses on 176 SPHERE/IRDIS data obtained with good SR conditions (> 40%) so as to test two different fitting strategies, namely (i) fitting the PSD and photometry/astrometry parameters (11 parameters) and (ii) fitting these latter 11, and 12 additional parameters corresponding to piston, tip, and tilt of each VLT pupil area delimited by the spiders in order to account for the LWE. According to Sauvage et al. (2016), this description of the LWE allows the main impact of this effect to be reproduced, i.e., a strong PSF asymmetry.
Figure 11 presents the results of the two strategies implemented in our study over three particular cases for which a strong LWE is observed. We observe, especially in Fig. 11, that the sole parametrization of the PSD is not rich enough to precisely reproduce the strong asymmetry, which becomes a solved problem thanks to the aberration parametrization we propose. Regarding the estimated PSF shape and values given in Table 6, we show clear evidence that (i) describing the LWE as a combination of differential piston, tip, and tilt allows to accurately characterize the PSF asymmetry, and (ii) accounting for this description in the PSF model ensures that biases and dispersion on SR and FWHM estimates are drastically mitigated. Moreover, the average MSE value over the 176 data samples is significantly diminished, that is, by a factor two. On the three cases illustrated in Fig. 11, Table 7 also reports the SR and FWHM errors as well as the MSE obtained from the PSF-fitting. These results show clearly a gain on accuracy for these particular cases for which the LWE impact on the PSF is significant comparatively to the statistical trend observed on the 176 data sets. This gain is partially mitigated on the statistical analysis owing to the fact that only 25% of the data were noticeably contaminated by the LWE. Moreover, the atmospheric part of the PSF model can mimic PSF asymmetries by modifying the ratio αx/αy in Eq. (7). Despite it leads to an inaccurate representation of the LWE impact on the PSF, it permits to mitigate the SR and FWHM estimation error comparatively to a symmetric atmospheric PSF model. As a summary, fitting the twelve extra parameters to represent the LWE allows can enhance the PSF metrics estimation by a factor up to five.
![]() |
Fig. 11. Comparison in log scale of on-sky PSFs and fitted model considering the LWE or not. Each column correspond to a single data sample acquired at different epoch. For each model configuration, we report both the adjusted PSF and the residual map. This figure illustrates that the LWE retrieval, modelled by a differential piston, tip and tilt improves drastically and consistently the residual map. |
Median and 1-σ dispersion of SR and FWHM errors obtained over the 176 H-band and high SR (> 40%) SPHERE/IRDIS images treated in this analysis for the two implemented strategies: (i) PSD/stellar parameters retrieval only (11 parameters in total) and (ii) including piston, tip, and tilt retrieval for the four pupil segments (23 parameters in total).
In order to confirm the robustness of the LWE retrieval, we have compared the estimated static aberration map with ZELDA measurements (Vigan et al. 2019) taken during SPHERE commissioning nights in 2014. The PSF-fitting was performed using the differential tip-tilt sensor (DTTS) that delivers 32 × 32 pixels images (Baudoz et al. 2010; Sauvage et al. 2015) acquired simultaneously with SPHERE/IRDIS using the ZELDA focal-plane mask to measure optical aberrations within the pupil-plane. Moreover, in order to calibrate ZELDA measurements (that are in ADU) to reconstruct the phase, we followed the process described by Sauvage et al. (2015) and removed the mean ADU value over the pupil and adjusted a multiplicative factor to obtain the closest PSF possible from the DTTS observation as reported in Fig. 12. According to this calibration, we obtained standard deviations of 173 nm rms and 146 nm rms on the ZELDA map and the DTTS image-based map,respectively, which leads to a quadratic difference of 92 nm rms. This residual includes the internal aberrations in the IRDIS science path that do not impact the DTTS images and reach 50 nm and standard deviation (Beuzit et al. 2019). Also, this residual is likely mostly due to higher order aberrations not included in the differential piston, tip, and tilt of the static map model, which suggests that further improvements can be pursued to characterize the LWE. In conclusion, the framework presented in this paper is a powerful tool to assist in the estimation of internal aberrations. Furthermore, we are able to assess these aberrations from the DTTS image that delivers a post-AO and non-coronagraphic image regardless the coronagraphic mask employed during SPHERE/IRDIS observations. Consequently, this technique allows the joint estimation of atmospheric and instrumental defects using tip-tilt sensors measurements. Future work will address the extension of this strategy to ELT instruments that will rely on a 2 × 2 Shack-Hartmann WFS to measure low-order modes and that will provide PSFs from which we will be able to assess the telescope aberrations and calibrate a PSF model for science exploitation.
![]() |
Fig. 12. Top: from left to right: ZELDA measurements (114 nm std), retrieved static map (140 nm) from the PSF-fitting and residual (70 nm std) in nm. This comparison emphasizes the meaningfulness of the PSF-fitting outputs, which compares very well with dedicated measurements of the LWE. Bottom: from left to right: on-sky DTTS image, best-fitted PSF, ZELDA-based PSF. |
5. Conclusions
This paper revisits and improves the analytical framework proposed by Fétick et al. (2019b), which now includes a parametrization of static aberrations for a joint retrieval of atmospheric parameters, AO performance and static aberrations.
We demonstrate in an exhaustive manner, using 4812 PSFs obtained from four different observatories and seven optical or NIR instruments, that the proposed model matches the PSF of any AO flavor within 4% error, even for high-SR observations. We also illustrate that the retrieved parameters carry relevant information about the AO performance and the atmospheric conditions, especially seeing and wavefront error, that shows agreement with the literature.
Finally, we illustrate that this model, upgraded with additional degrees of freedom to estimate static aberrations, allows the atmospheric parameters, AO performance and static aberrations, to be retrieved simultaneously. Particularly, the framework presented in this paper allows us to assess (i) the Keck pupil segment piston errors and especially the presence of a stair mode that was already pointed out (Ragland 2018) and, (ii) the LWE on SPHERE/IRDIS images as a combination of differential piston, tip and tilt over the four VLT pupil segments delimited by the spiders.
This model is a unique tool that gathers an AO diagnosis and a PSF estimation facility in the simplest and the most parsimonious way possible. However, we have handled it as a parametric model so far and the next step of this work is to enable a forward estimation of its parameters from contextual data expect the focal-plane image. We emphasized that the connection of these parameters with the observing conditions is not easily made; nevertheless, thanks to the parsimony of this model, this regression problem consists in assessing a few parameters (up to a few tens) from a large amount of data, which is provided by the AO telemetry, all the sensors within the telescope and the dome and the external meteorological profilers, which can be achieved with the use of neural networks. We are currently developing convolutional neural networks capable of directly estimating the model outputs from either an imaged PSF or a subsample of AO telemetry and we will present this work in a dedicated publication.
Acknowledgments
This work has been partially funded by the French National Research Agency (ANR) program APPLY – ANR-19-CE31-0011. This work also benefited from the support of the WOLF project ANR-18-CE31-0018 of the and the OPTICON H2020 (2017-2020) Work Package 1. This work has made use of the SPHERE Data Centre, jointly operated by OSUG/IPAG (Grenoble), PYTHEAS/LAM/CeSAM (Marseille), OCA/Lagrange (Nice) and Observatoire de Paris/LESIA (Paris) and supported by a grant from Labex OSUG@2020 (Investissements d’avenir – ANR10 LABX56). Authors thank G. Fiorentino (INAF), F. Kerber (ESO), D. Massari (U. Bologna, INAF, Kapteyn Astronomical Institute), J. Milli (IPAG), and E. Tolstoy (Kapteyn Astronomical Institute) to have conducted the SPHERE/ZIMPOL observations of NGC 6121 in 2018 and provided the data. Authors thank P. Vernazza (LAM) to have delivered SPHERE/ZIMPOL PSF calibrators obtained during his large programme. Authors acknowledge the contribution of J. Milli (IPAG) to utilize the SPHERE Data centre pipeline. Authors thank F. Cantalloube (MPIA) and M. N’Diaye (Lagrange) to provide the APLC mask model. Authors thank Jean-François Sauvage (ONERA/LAM) for providing SPHERE DTTS and ZELDA measurements that served in the Sect. 4.2. Authors thank Sam Ragland (W.M. Keck Observatory) for giving access to NIRC2 images obtained during PSF reconstruction engineering nights on Keck II. Author thank G. Agapito (INAF) and E. Pinna (INAF) to communicate SOUL/LUCI commissioning data. Authors thank J. Vernet (ESO) and S. Oberti (ESO) for providing MUSE NFM commissioning data. Authors thank M. Andersen (Gemini), G. Sivo (Gemini) and A. Shugart (Gemini) for communicating GeMS/GSAOI observations of Trumpler 14 and supporting the data processing. Author express their gratitude to F. Pedreros-Bustos (LAM) for a fruitful revision of the manuscript.
References
- Arsenault, R., Madec, P. Y., Hubin, N., et al. 2008, Proc. SPIE, 7015, 701524 [CrossRef] [Google Scholar]
- Bacon, R., Accardo, M., Adjali, L., et al. 2010, Proc. SPIE, 7735, 773508 [Google Scholar]
- Baudoz, P., Dorn, R. J., Lizon, J.-L., et al. 2010, Proc. SPIE, 7735, 77355B [CrossRef] [Google Scholar]
- Beltramo-Martin, O., Correia, C. M., Ragland, S., et al. 2019, MNRAS, 487, 5450 [NASA ADS] [CrossRef] [Google Scholar]
- Beltramo-Martin, O., Marasco, A., Fusco, T., et al. 2020, MNRAS, 494, 775 [CrossRef] [Google Scholar]
- Benfenati, A., La Camera, A., & Carbillet, M. 2016, A&A, 586, A16 [CrossRef] [EDP Sciences] [Google Scholar]
- Bertin, E., & Arnouts, S. 1996, A&A, 117, 393 [Google Scholar]
- Beuzit, J. L., Vigan, A., Mouillet, D., et al. 2019, A&A, 631, A155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Blanc, A., Fusco, T., Hartung, M., Mugnier, L. M., & Rousset, G. 2003, A&A, 399, 373 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bond, C. Z., Correia, C. M., Sauvage, J.-F., et al. 2018a, Proc. SPIE, 10703, 107034M [Google Scholar]
- Bond, C. Z., Wizinowich, P., Chun, M., et al. 2018b, Proc. SPIE, 10703, 107031Z [Google Scholar]
- Bouché, N., Carfantan, H., Schroetter, I., Michel-Dansac, L., & Contini, T. 2015, AJ, 150, 92 [NASA ADS] [CrossRef] [Google Scholar]
- Chanan, G. A., Mast, T. S., & Nelson, J. E. 1988, Very Large Telescopes and their Instrumentation, 1, 421 [Google Scholar]
- Claudi, R., Giro, E., Turatto, M., et al. 2014, Proc. SPIE, 9147, 91471L [CrossRef] [Google Scholar]
- Clénet, Y., Lidman, C., Gendron, E., et al. 2008, Proc. SPIE, 7015, 701529 [CrossRef] [Google Scholar]
- Conan, R., van Dam, M., Vogiatzis, K., Das, K., & Bouchez, A. 2019, Sixth International Conference on Adaptive Optics for Extremely Large Telescopes, 70 [Google Scholar]
- Conn, A., Gould, N., & Toint, P. 2000, Trust Region Methods, MPS-SIAM Serieson Optimization (Society for Industrial and Applied Mathematics) [Google Scholar]
- Correia, C. M., & Teixeira, J. 2014, J. Opt. Soc. Am. A, 31, 2763 [CrossRef] [Google Scholar]
- Davies, R., Schubert, J., Hartl, M., et al. 2016, Proc. SPIE, 9908, 99081Z [CrossRef] [Google Scholar]
- Delorme, P., Meunier, N., Albert, D., et al. 2017, in SF2A-2017: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics, eds. C. Reylé, P. Di Matteo, F. Herpin, et al. [Google Scholar]
- Diolaiti, E., Bendinelli, O., Bonaccini, D., et al. 2000, Astrophysics Source Code Library [record ascl: 0011.001] [Google Scholar]
- Dohlen, K., Langlois, M., Saisse, M., et al. 2008, Proc. SPIE, 7014, 70143L [Google Scholar]
- Drummond, J. D. 1998, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, eds. D. Bonaccini, & R. K. Tyson, Proc. SPIE, 3353, 1030 [Google Scholar]
- Epinat, B., Amram, P., Balkowski, C., & Marcelin, M. 2010, MNRAS, 401, 2113 [NASA ADS] [CrossRef] [Google Scholar]
- Exposito, J. 2014, PhD Thesis, Université Paris Diderot, France [Google Scholar]
- Fauvarque, O., Neichel, B., Fusco, T., Sauvage, J.-F., & Girault, O. 2016, Optica, 3, 1440 [NASA ADS] [CrossRef] [Google Scholar]
- Ferreira, F., Gendron, E., Rousset, G., & Gratadour, D. 2018, A&A, 616, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Fétick, R. J., Jorda, L., Vernazza, P., et al. 2019a, A&A, 623, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Fétick, R. J. L., Fusco, T., Neichel, B., et al. 2019b, A&A, 628, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Fétick, R. J. L., Mugnier, L. M., Fusco, T., & Neichel, B. 2020, MNRAS, 496, 4209 [CrossRef] [Google Scholar]
- Flicker, R. C., & Rigaut, F. J. 2005, J. Opt. Soc. Am. A, 22, 504 [CrossRef] [Google Scholar]
- Fried, D. L. 1966, J. Opt. Soc. Am. (1917–1983), 56, 1380 [NASA ADS] [CrossRef] [Google Scholar]
- Fried, D. L. 1982, J. Opt. Soc. Am. (1917–1983), 72, 52 [NASA ADS] [CrossRef] [Google Scholar]
- Fusco, T., Mugnier, L. M., Conan, J. M., et al. 2003, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, eds. P. L. Wizinowich, D. Bonaccini, et al., Proc. SPIE, 4839, 1065 [Google Scholar]
- Fusco, T., Sauvage, J. F., Mouillet, D., et al. 2016, Proc. SPIE, 9909, 99090U [CrossRef] [Google Scholar]
- Fusco, T., Bacon, R., Kamann, S., et al. 2020, A&A, 635, A208 [CrossRef] [EDP Sciences] [Google Scholar]
- Gendron, E., Clénet, Y., Fusco, T., & Rousset, G. 2006, A&A, 457, 359 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gendron, E., Morris, T., Hubert, Z., et al. 2010, in Adaptive Optics Systems II, eds. B. Ellerbroek, M. Hart, N. Hubin, P. Wizinowich, Proc. SPIE, 7736 [Google Scholar]
- Gerwe, D., Johnson, M., & Calef, B. 2008, Advanced Maui Optical and Space Surveillance Technologies Conference, E28 [Google Scholar]
- Gilles, L., Wang, L., & Boyer, C. 2018, Proc. SPIE, 10703, 1070349 [Google Scholar]
- Gratadour, D., Puech, M., Vérinaud, C., et al. 2014, Proc. SPIE, 9148, 91486O [CrossRef] [Google Scholar]
- Hammer, F., Morris, S., Kaper, L., et al. 2016, in Ground-based and Airborne Instrumentation for Astronomy VI, Proc. SPIE, 9908, 990824 [Google Scholar]
- Herbel, J., Kacprzak, T., Amara, A., Refregier, A., & Lucchi, A. 2018, JCAP, 2018, 054 [Google Scholar]
- Hippler, S., Feldt, M., Bertram, T., et al. 2019, Exp. Astron., 47, 65 [CrossRef] [Google Scholar]
- Hirst, P., & Cardenes, R. 2017, in Astronomical Data Analysis Software and Systems XXV, eds. N. P. F. Lorente, K. Shortridge, & R. Wayth, ASP Conf. Ser., 512, 53 [Google Scholar]
- Jia, P., Wu, X., Yi, H., Cai, B., & Cai, D. 2020, AJ, 159, 183 [CrossRef] [Google Scholar]
- Jolissaint, L. 2010, J. Eur. Opt. Soc. Rapid Publ., 5, 10055 [CrossRef] [Google Scholar]
- Jolissaint, L., Mugnier, L. M., Neyman, C., Christou, J., & Wizinowich, P. 2012, Proc. SPIE, 8447, 844716 [CrossRef] [Google Scholar]
- Jolissaint, L., Ragland, S., & Wizinowich, P. 2015, Adaptive Optics for Extremely Large Telescopes IV (AO4ELT4), E93 [Google Scholar]
- Jolissaint, L., Ragland, S., Christou, J., & Wizinowich, P. 2018, Appl. Opt., 57, 7837 [NASA ADS] [CrossRef] [Google Scholar]
- Kolmogorov, A. N. 1941, Akademiia Nauk SSSR Doklady, 32, 16 [NASA ADS] [Google Scholar]
- Laginja, I., Leboulleux, L., Pueyo, L., et al. 2019, Proc. SPIE, 11117, 1111717 [Google Scholar]
- Lai, O., Withington, J. K., Laugier, R., & Chun, M. 2019, MNRAS, 484, 5568 [Google Scholar]
- Laidlaw, D. J., Osborn, J., Morris, T. J., et al. 2019, MNRAS, 483, 4341 [CrossRef] [Google Scholar]
- Lamb, M., Norton, A., Macintosh, B., et al. 2018, Proc. SPIE, 10703, 107035M [Google Scholar]
- Leboulleux, L., Sauvage, J.-F., Pueyo, L. A., et al. 2018, J. Astron. Telescopes Instrum. Syst., 4, 035002 [Google Scholar]
- Martin, O., Gendron, E., Rousset, G., & Vidal, F. 2012, in Adaptive Optics Systems III, Proc. SPIE, 8447, 84472A [CrossRef] [Google Scholar]
- Martin, O., Gendron, E., Rousset, G., et al. 2013, in Adaptive Optics for Extremely Large Telescopes, Florence, Italy, ed. C. Kidon, 3rd AO4ELT conf. [Google Scholar]
- Martin, O. A., Correia, C. M., Gendron, E., et al. 2016a, J. Astron. Telescopes Instrum. Syst., 2, 048001 [NASA ADS] [CrossRef] [Google Scholar]
- Martin, O. A., Correia, C. M., Gendron, E., et al. 2016b, Proc. SPIE, 9909, 99093P [CrossRef] [Google Scholar]
- Martin, O. A., Gendron, É., Rousset, G., et al. 2017, A&A, 598, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Masciadri, E., Lombardi, G., & Lascaux, F. 2014, MNRAS, 438, 983 [CrossRef] [Google Scholar]
- Massari, D., Marasco, A., Beltramo-Martin, O., et al. 2020, A&A, 634, L5 [CrossRef] [EDP Sciences] [Google Scholar]
- Milli, J., Kasper, M., Bourget, P., et al. 2018, Proc. SPIE, 10703, 107032A [Google Scholar]
- Miyashita, A., Takato, N., Usuda, T., Uraguchi, F., & Ogasawara, R. 2004, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, eds. J. Oschmann, & M. Jacobus, Proc. SPIE, 5489, 207 [Google Scholar]
- Moffat, A. F. J. 1969, A&A, 3, 455 [NASA ADS] [Google Scholar]
- Morris, T., Gendron, E., Basden, A., et al. 2013, in Proceedings of the Third AO4ELT Conference, eds. S. Esposito, L. Fini [Google Scholar]
- Morris, T., Gendron, E., Basden, A., et al. 2014, in Adaptive Optics Systems IV, Proc. SPIE, 9148, 1 [Google Scholar]
- Mugnier, L. M., Fusco, T., & Conan, J.-M. 2004, J. Opt. Soc. Am. A, 21, 1841 [NASA ADS] [CrossRef] [Google Scholar]
- Mugnier, L. M., Sauvage, J.-F., Fusco, T., Cornia, A., & Dandy, S. 2008, Opt. Express, 16, 18406 [NASA ADS] [CrossRef] [Google Scholar]
- Myers, R. M., Hubert, Z., Morris, T. J., et al. 2008, in Adaptive Optics Systems, eds. P. W. N. Hubin, C. Max, Proc. SPIE, 7015 [Google Scholar]
- N’Diaye, M., Vigan, A., Dohlen, K., et al. 2016, A&A, 592, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Neichel, B., Rigaut, F., Vidal, F., et al. 2014, MNRAS, 440, 1002 [NASA ADS] [CrossRef] [Google Scholar]
- Oberti, S., Kolb, J., Madec, P.-Y., et al. 2018, Proc. SPIE, 10703, 107031G [Google Scholar]
- O’Mahony, N. 2003, The Newsletter of the Isaac Newton Group of Telescopes, 7, 22 [Google Scholar]
- Ono, Y. H., Akiyama, M., Oya, S., et al. 2016, J. Opt. Soc. Am. A, 33, 726 [CrossRef] [Google Scholar]
- Parenti, R. R., & Sasiela, R. J. 1994, J. Opt. Soc. Am. A, 11, 288 [NASA ADS] [CrossRef] [Google Scholar]
- Pinna, E., Esposito, S., Hinz, P., et al. 2016, Proc. SPIE, 9909, 99093V [CrossRef] [Google Scholar]
- Puech, M., Evans, C. J., Disseau, K., et al. 2018, Proc. SPIE, 10702, 107028R [Google Scholar]
- Ragazzoni, R., & Farinato, J. 1999, A&A, 350, L23 [Google Scholar]
- Ragland, S. 2018, in Ground-based and Airborne Telescopes VII, SPIE Conf. Ser., 10700, 107001D [Google Scholar]
- Ragland, S., Jolissaint, L., Wizinowich, P., et al. 2016, in Adaptive Optics Systems V, Proc. SPIE, 9909, 99091P [CrossRef] [Google Scholar]
- Ragland, S., Dupuy, T. J., Jolissaint, L., et al. 2018, in Adaptive Optics Systems VI, Proc. SPIE, 10703, 107031J [Google Scholar]
- Rigaut, F. J., Véran, J. P., & Lai, O. 1998, in Adaptive Optical System Technologies, eds. D. Bonaccini, & R. K. Tyson, Proc. SPIE, 3353, 1038 [NASA ADS] [CrossRef] [Google Scholar]
- Rigaut, F., Neichel, B., Boccas, M., et al. 2012, Proc. SPIE, 8447, 84470I [CrossRef] [Google Scholar]
- Rigaut, F., Neichel, B., Boccas, M., et al. 2014, MNRAS, 437, 2361 [NASA ADS] [CrossRef] [Google Scholar]
- Robert, C., Fusco, T., Sauvage, J.-F., & Mugnier, L. 2008, Proc. SPIE, 7015, 70156A [CrossRef] [Google Scholar]
- Roddier, F. 1981, Progress in Optics (Amsterdam: North-Holland Publishing Co.), 19, 281 [NASA ADS] [CrossRef] [Google Scholar]
- Roddier, F. 1999, Adaptive Optics in Astronomy (Cambridge University Press) [CrossRef] [Google Scholar]
- Sauvage, J. F., Fusco, T., LeMignant, D., et al. 2011, Second International Conference on Adaptive Optics for Extremely Large Telescopes, 48, Online at http://ao4elt2.lesia.obspm.fr [Google Scholar]
- Sauvage, J. F., Fusco, T., Guesalaga, A., et al. 2015, Adaptive Optics for Extremely Large Telescopes IV (AO4ELT4), E9 [Google Scholar]
- Sauvage, J.-F., Fusco, T., Lamb, M., et al. 2016, Proc. SPIE, 9909, 990916 [Google Scholar]
- Schmid, H. M., Bazzon, A., Roelfsema, R., et al. 2018, A&A, 619, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schreiber, L., Diolaiti, E., Sollima, A., et al. 2012, Proc. SPIE, 8447, 84475V [CrossRef] [Google Scholar]
- Schwartz, N., Sauvage, J.-F., Correia, C., et al. 2018, Proc. SPIE, 10703, 1070322 [Google Scholar]
- Sitarski, B. N., Witzel, G., Fitzgerald, M., et al. 2014, in Adaptive Optics Systems IV, Proc. SPIE, 9148 [Google Scholar]
- Soummer, R. 2005, ApJ, 618, L161 [NASA ADS] [CrossRef] [Google Scholar]
- Stetson, P. B. 1987, PASP, 99, 191 [Google Scholar]
- Ströbele, S., La Penna, P., Arsenault, R., et al. 2012, Proc. SPIE, 8447, 844737 [CrossRef] [Google Scholar]
- Thatte, N. 2017, Adaptative Optics for Extremely Large Telescopes V [Google Scholar]
- Tokovinin, A., & Travouillon, T. 2006, MNRAS, 365, 1235 [NASA ADS] [CrossRef] [Google Scholar]
- Trujillo, I., Aguerri, J. A. L., Cepa, J., & Gutiérrez, C. M. 2001, MNRAS, 328, 977 [NASA ADS] [CrossRef] [Google Scholar]
- van Dam, M. A., Sasiela, R. J., Bouchez, A. H., et al. 2006, Proc. SPIE, 6272, 627231 [CrossRef] [Google Scholar]
- Vassallo, D., Farinato, J., Sauvage, J. F., et al. 2018, Proc. SPIE, 10705, 1070516 [Google Scholar]
- Veran, J.-P., Rigaut, F., Maitre, H., & Rouan, D. 1997, J. Opt. Soc. Am. A, 14, 3057 [NASA ADS] [CrossRef] [Google Scholar]
- Vernazza, P., Brozv, M., Drouard, A., et al. 2018, VizieR Online Data Catalog: J/A+A/618/A154 [Google Scholar]
- Vidal, F., Gendron, E., & Rousset, G. 2010, JOSA A, 27, A253 [NASA ADS] [CrossRef] [Google Scholar]
- Vidal, F., Gendron, É., Rousset, G., et al. 2014, A&A, 569, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Vigan, A., Moutou, C., Langlois, M., et al. 2010, MNRAS, 407, 71 [NASA ADS] [CrossRef] [Google Scholar]
- Vigan, A., N’Diaye, M., Dohlen, K., et al. 2019, A&A, 629, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- von Karman, T. 1948, Proc. Nat. Acad. Sci., 34, 530 [NASA ADS] [CrossRef] [Google Scholar]
- Wagner, R., Hofer, C., & Ramlau, R. 2018, J. Astron. Telescopes Instrum. Syst., 4, 049003 [Google Scholar]
- Witzel, G., Lu, J. R., Ghez, A. M., et al. 2016, in Proc. SPIE, Adaptive Optics Systems V, 9909, 99091O [CrossRef] [Google Scholar]
- Wizinowich, P. 2012, in Adaptive Optics Systems III, Proc. SPIE, 8447, 84470D [CrossRef] [Google Scholar]
- Wizinowich, P. L., Acton, D. S., Lai, O., et al. 2000, in Adaptive Optical Systems Technology, ed. P. L. Wizinowich, Proc. SPIE, 4007, 2 [NASA ADS] [CrossRef] [Google Scholar]
- Wizinowich, P. L., Le Mignant, D., Bouchez, A. H., et al. 2006, PASP, 118, 297 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Summary of PSFs obtained and processed for the analysis presented in this paper.
Individual statistics per system and imaging wavelength on SR and FWHM estimates.
Median and 1-σ standard-deviation of seeing values retrieved from the PSF-fitting process.
Median and 1-σ dispersion of SR and FWHM errors and mean square error obtained over the 129 H-band and high SR (> 40%) Keck/NIRC2 images treated in this analysis for the three implemented strategies: (i) PSD/stellar parameters retrieval, (ii) including the residual NCPA map in the model, and (iii) also including the estimation of Keck II cophasing errors.
Summary of peak-to-valley (PV) and 1-σ standard deviation of retrieved cophasing error map regarding the corresponding telescope elevation/azimuth.
Median and 1-σ dispersion of SR and FWHM errors obtained over the 176 H-band and high SR (> 40%) SPHERE/IRDIS images treated in this analysis for the two implemented strategies: (i) PSD/stellar parameters retrieval only (11 parameters in total) and (ii) including piston, tip, and tilt retrieval for the four pupil segments (23 parameters in total).
All Figures
![]() |
Fig. 1. From left to right: pupil plane apodizer and Lyot stop used during SPHERE/IRDIS observations with the N_ALC_YJH_S APLC (Vigan et al. 2010). |
In the text |
![]() |
Fig. 2. Azimuthal average in log-log scale of the PSD (left) and the PSF (right) obtained with four different settings for the model parameters, assuming that the PSD is symmetric, i.e. αx = αy = α and θ = 0. No static aberrations were added here and the parameter C was set to C = 10−4 rad2m2. The PSF maximum is normalized to the SR value. |
In the text |
![]() |
Fig. 3. Two dimensional comparison in log scale and over 64 × 64 pixels of (top) observed PSFs, (middle) fitted model, and (bottom) residual map for multiple types of AO systems and instruments working in either visible or NIR wavelengths. From left to right: SPHERE/ZIMPOL, SPHERE/IRDIS with the apodizer and the Lyot stop (no focal plane mask), Keck AO/NIRC2, SOUL/LUCI, CANARY/CAMICAZ in MOAO mode, GALACSI/MUSE NFM, and GEMS/GSAOI. All PSFs have been normalized to the sum of pixels. |
In the text |
![]() |
Fig. 4. Image SR/FWHM versus the same metrics retrieved on the fitted image using the same estimation process and for the 4812 PSFs treated for this analysis. Error bars on the SR are obtained from calculations presented in Martin et al. (2016a). Error bars on the FWHM are given from the contour estimation on interpolated images that are oversampled by up to a factor four to quantify the FWHM more accurately. |
In the text |
![]() |
Fig. 5. r0 estimates at zenith and 500 nm with respect to the GEMS/GSAOI PSF field position (left) and with respect to the PSF GALACSI/MUSE wavelength (right). |
In the text |
![]() |
Fig. 6. Left: histograms for the A parameter for each AO system. Right: wavefront error obtained on CANARY SCAO mode with respect to the retrieved r0. The solid line gives the trend in |
In the text |
![]() |
Fig. 7. Histograms for the β and α parameters for each system, as well as the averaged distribution. |
In the text |
![]() |
Fig. 8. Left: average of absolute retrieved parameters with respect to the exposure time. The parameters are normalized by the value obtained from the long-exposure PSF (60 frames). Error bars are assessed by averaging results over the SOUL/LUCI data sets. Right: average retrieved parameters normalized by the median value over the whole span with respect to the spectral width. The envelope shows the ±4% of variations around the median value and the spectral bandwidth is systematically centered around 700 nm. Error bars are assessed by averaging results over the GALACSI/MUSE data sets. |
In the text |
![]() |
Fig. 9. Comparison in log scale of on-sky PSFs and fitted model using the three implemented strategies. Each column correspond to a single data sample acquired at different epochs. For each model configuration, we report both the adjusted PSF and the residual map. This figure illustrates that the cophasing map retrieval improves drastically and consistently the residual on top of the gain brought by including the NCPA residual map. |
In the text |
![]() |
Fig. 10. Retrieved cophasing maps in nm from PSFs acquired during three different epochs in 2013. The telescope elevation and azimuth were respectively (79,240), (85,−82) and (43,240) from left to right. |
In the text |
![]() |
Fig. 11. Comparison in log scale of on-sky PSFs and fitted model considering the LWE or not. Each column correspond to a single data sample acquired at different epoch. For each model configuration, we report both the adjusted PSF and the residual map. This figure illustrates that the LWE retrieval, modelled by a differential piston, tip and tilt improves drastically and consistently the residual map. |
In the text |
![]() |
Fig. 12. Top: from left to right: ZELDA measurements (114 nm std), retrieved static map (140 nm) from the PSF-fitting and residual (70 nm std) in nm. This comparison emphasizes the meaningfulness of the PSF-fitting outputs, which compares very well with dedicated measurements of the LWE. Bottom: from left to right: on-sky DTTS image, best-fitted PSF, ZELDA-based PSF. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.