Search

This study introduces a novel frequency-based data-driven controller for adaptive optics, using power spectral density for optimization while ensuring stability criteria. It addresses disturbance rejection, command amplitude constraints and system transfer functions through convex optimization to obtain an optimal control in an infinite input response filter form. Evaluated within the SAXO+ project, it demonstrates efficacy under diverse atmospheric conditions and operational scenarios. The proposed controller is tested in both standard and disentangled adaptive optics schemes, showcasing its adaptability and performance. Experimental validation is conducted using the COMPASS simulation tool, affirming the controller's promise for enhancing adaptive optics systems in real-world applications.

We tackle covariance estimation in low-sample scenarios, employing a structured covariance matrix with shrinkage methods. These involve convexly combining a low-bias/high-variance empirical estimate with a biased regularization estimator, striking a bias-variance trade-off. Literature provides optimal settings of the regularization amount through risk minimization between the true covariance and its shrunk counterpart. Such estimators were derived for zero-mean statistics with i.i.d. diagonal regularization matrices accounting for the average sample variance solely. We extend these results to regularization matrices accounting for the sample variances both for centered and non-centered samples. In the latter case, the empirical estimate of the true mean is incorporated into our shrinkage estimators. Introducing confidence weights into the statistics also enhance estimator robustness against outliers. We compare our estimators to other shrinkage methods both on numerical simulations and on real data to solve a detection problem in astronomy., Comment: Accepted to EUSIPCO 2024

Exoplanet detections and characterizations via direct imaging require high contrast and high angular resolution. These requirements typically require (i) cutting-edge instrumental facilities, (ii) optimized differential imaging to introduce a diversity in the signals of the sought-for objects, and (iii) dedicated processing algorithms to further eliminate the residual stellar leakages. Substantial efforts have been undertaken on the design of more efficient post-processing algorithms but their performance remains upper-bounded at shorter angular separations due to the the lack of diversity induced by the processing of each epoch of observations individually. We propose a new algorithm that is able to combine several observations of the same star by accounting for the Keplerian orbital motion across epochs of the sought-for sources in order to constructively co-add their weak signals. The proposed algorithm, PACOME, integrates an exploration of the plausible orbits within a statistical detection and estimation formalism. It is extended to a multi-epoch combination of the maximum likelihood framework of PACO, which is a mono-epoch post-processing algorithm. We derive a reliable multi-epoch detection criterion, interpretable both in terms of probability of detection and of false alarm. We tested the proposed algorithm on several datasets obtained from the VLT/SPHERE instrument with IRDIS and IFS. By resorting to injections of synthetic exoplanets, we show that PACOME is able to detect sources remaining undetectable in mono-epoch frameworks. The gain in detection sensitivity scales as high as the square root of the number of epochs. We also applied PACOME on a set of observations from the HR 8799 star hosting four known exoplanets, which are detected with very high signal-to-noise ratios. In addition, its implementation is efficient, fast, and fully automatized., Comment: Accepted for publication in A&A

Context. Adaptive optics (AO) is now a tool commonly deployed in astronomy. The real time correction of the atmospheric turbulence that AO enables allows telescopes to perform close to the diffraction limit at the core of their point spread function (PSF). Among other factors, AO-corrected PSFs depend on the ability of the wavefront corrector (WFC), generally a deformable mirror, to fit the incident wavefront corrugations. Aims. In this work, we focus on this error introduced by the WFC, the so-called fitting error. To date, analytical models only depend on the WFC cut-off frequency, and Monte Carlo simulations are the only solution for studying the impact of the WFC influence function shape on the AO-corrected PSF. We aim to develop an analytical model accounting for the influence function shape. Methods. We first obtain a general analytical model of the fitting error structure function. With additional hypotheses, we then derive an analytical model of the AO-corrected power spectral density. These two analytical solutions are compared with Monte Carlo simulations on different ideal profiles (piston, pyramid, Gaussian) as well as with real hardware (DM192 from ALPAO). Results. Our analytical predictions show a very good agreement with the Monte Carlo simulations. We show that in the image plane, the depth of the correction as well as the transition profile between the AO-corrected area and the remaining turbulent halo depend on the influence functions of the WFC. We also show that the generally assumed hypothesis of stationarity of the AO correction is actually not met. Conclusions. As the fitting error is the intrinsic optimal limit of an AO system, our analytical model allows for the assessment of the theoretical limits of extreme AO systems limited by the WFC in high-contrast imaging through a context where other errors become comparable.

High-contrast long-slit spectrographs can be used to characterize exoplanets. High-contrast long-slit spectroscopic data are however corrupted by stellar leakages which largely dominate other signals and make the process of extracting the companion spectrum very challenging. This paper presents a complete method to calibrate the spectrograph and extract the signal of interest. The proposed method is based on a flexible direct model of the high-contrast long-slit spectroscopic data. This model explicitly accounts for the instrumental response and for the contributions of both the star and the companion. The contributions of these two components and the calibration parameters are jointly estimated by solving a regularized inverse problem. This problem having no closed-form solution, we propose an alternating minimization strategy to effectively find the solution. We have tested our method on empirical long-slit spectroscopic data and by injecting synthetic companion signals in these data. The proposed initialization and the alternating strategy effectively avoid the self-subtraction bias, even for companions observed very close to the coronagraphic mask. Careful modeling and calibration of the angular and spectral dispersion laws of the instrument clearly reduce the contamination by the stellar leakages. In practice, the outputs of the method are mostly driven by a single hyper-parameter which tunes the level of regularization of the companion SED., Comment: Paper under review by Astronomy & Astrophysics

We present the results obtained with an end-to-end simulator of an Extreme Adaptive Optics (XAO) system control loop. It is used to predict its on-sky performances and to optimise the AO loop algorithms. It was first used to validate a novel analytical model of the fitting error, a limit due to the Deformable Mirror (DM) shape. Standard analytical models assume a sharp correction under the DM cutoff frequency, disregarding the transition between the AO corrected and turbulence dominated domains. Our model account for the influence function shape in this smooth transition. Then, it is well-known that Shack-Hartmann wavefront sensors (SH-WFS) have a limited spatial bandwidth, the high frequencies of the wavefront being seen as low frequencies. We show that this aliasing error can be partially compensated (both in terms of Strehl ratio and contrast) by adding priors on the turbulence statistics in the framework of an inverse problem approach. This represents an alternative to the standard additional optical filter used in XAO systems. In parallel to this numerical work, a bench was aligned to experimentally test the AO system and these new algorithms comprising a DM192 ALPAO deformable mirror and a 15x15 SH-WFS. We present the predicted performances of the AO loop based on end-to-end simulations.

The Evanescent Wave Coronagraph (EvWaCo) is an achromatic coronagraph mask with adjustable size over the spectral domain [600nm, 900nm] that will be installed at the Thai National Observatory. We present in this work the development of a bench to characterise its Extreme Adaptive Optics system (XAO) comprising a DM192 ALPAO deformable mirror (DM) and a 15x15 Shack-Hartmann wavefront sensor (SH-WFS). In this bench, the turbulence is simulated using a rotating phase plate in a pupil plane. In general, such components are designed using a randomly generated phase screen. Such single realisation does not necessarily provide the wanted structure function. We present a solution to design the printed pattern to ensure that the beam sees a strict and controlled Kolmogorov statistics with the correct 2D structure function. This is essential to control the experimental conditions in order to compare the bench results with the numerical simulations and predictions. This bench is further used to deeply characterise the full 27 mm pupil of the ALPAO DM using a 54x54 ALPAO SH-WFS. We measure the average shape of its influence functions as well as the influence function of each single actuator to study their dispersion. We study the linearity of the actuator amplitude with the command as well as the linearity of the influence function profile. We also study the actuator offsets as well as the membrane shape at 0-command. This knowledge is critical to get a forward model of the DM for the XAO control loop.

In interferometry, the quality of the reconstructed image depends on the algorithm used and its parameters, and users often need to compare the results of several algorithms to disentangle artifacts from actual features of the astrophysical object. Such comparisons can rapidly become cumbersome, as these software packages are very different. OImaging is a graphical interface intended to be a common frontend to image reconstruction software packages. With OImaging, the user can now perform multiple reconstructions within a single interface. From a given dataset, OImaging allows benchmarking of different image reconstruction algorithms and assessment of the reliability of the image reconstruction process. To that end, OImaging uses the IMAGE-OI OIFITS extension proposed to standardize communication with image reconstruction algorithms.

Angular and spectral differential imaging is an observational technique of choice to investigate the immediate vicinity of stars. The relative angular motion and spectral scaling between on-axis and off-axis sources are exploited by post-processing techniques to separate two components: the residual star light and the light coming from surrounding objects such as circumstellar disks or point-like objects. This paper introduces a new algorithm to jointly unmix these two components and deconvolve disk images. The proposed algorithm is based on a statistical modeling of the residual star light, accounting for its spatial and spectral correlations. While critical, these correlations are not modeled and compensated for by existing reconstruction algorithms. We leverage dedicated shrinkage techniques to estimate the large number of parameters of our model of correlations in a data-driven fashion. We show that the resulting separable model of the spatial and spectral covariances captures very accurately the star light, thus allowing its efficient suppression. We apply our method on several datasets from the VLT/SPHERE instrument and compare the performances against standard algorithms of the field (cADI, PCA). We show that accounting for the multiple correlations of the data significantly enhances the reconstruction quality, leading to better preservation of both the disk morphology and photometry. Thanks to its unique joint spectral modeling of the data, the proposed algorithm can reconstruct disks having a circular symmetry (e.g., rings, spirals) at intensities one million times fainter than the star, without resorting to additional reference datasets free from the off-axis objects., Comment: Accepted to Monthly Notices of the Royal Astronomical Society

Polarimetric imaging is one of the most effective techniques for high-contrast imaging and characterization of circumstellar environments. These environments can be characterized through direct-imaging polarimetry at near-infrared wavelengths. The SPHERE/IRDIS instrument installed on the Very Large Telescope in its dual-beam polarimetric imaging (DPI) mode, offers the capability to acquire polarimetric images at high contrast and high angular resolution. However dedicated image processing is needed to get rid of the contamination by the stellar light, of instrumental polarization effects, and of the blurring by the instrumental point spread function. We aim to reconstruct and deconvolve the near-infrared polarization signal from circumstellar environments. We use observations of these environments obtained with the high-contrast imaging infrared polarimeter SPHERE-IRDIS at the VLT. We developed a new method to extract the polarimetric signal using an inverse approach method that benefits from the added knowledge of the detected signal formation process. The method includes weighted data fidelity term, smooth penalization, and takes into account instrumental polarization. The method enables to accurately measure the polarized intensity and angle of linear polarization of circumstellar disks by taking into account the noise statistics and the convolution of the observed objects by the instrumental point spread function. It has the capability to use incomplete polarimetry cycles which enhance the sensitivity of the observations. The method improves the overall performances in particular for low SNR/small polarized flux compared to standard methods.

Aims. The purpose of this paper is to describe a new post-processing algorithm dedicated to the reconstruction of the spatial distribution of light received from off-axis sources, in particular from circumstellar disks. Methods. Built on the recent PACO algorithm dedicated to the detection of point-like sources, the proposed method is based on the local learning of patch covariances capturing the spatial fluctuations of the stellar leakages. From this statistical modeling, we develop a regularized image reconstruction algorithm (REXPACO) following an inverse problem approach based on a forward image formation model of the off-axis sources in the ADI sequences. Results. Injections of fake circumstellar disks in ADI sequences from the VLT/SPHERE-IRDIS instrument show that both the morphology and the photometry of the disks are better preserved by REXPACO compared to standard postprocessing methods like cADI. In particular, the modeling of the spatial covariances proves usefull in reducing typical ADI artifacts and in better disentangling the signal of these sources from the residual stellar contamination. The application to stars hosting circumstellar disks with various morphologies confirms the ability of REXPACO to produce images of the light distribution with reduced artifacts. Finally, we show how REXPACO can be combined with PACO to disentangle the signal of circumstellar disks from the signal of candidate point-like sources. Conclusions. REXPACO is a novel post-processing algorithm producing numerically deblurred images of the circumstellar environment. It exploits the spatial covariances of the stellar leakages and of the noise to efficiently eliminate this nuisance term.

An adaptive optics system with a single deformable mirror is being implemented on the THEMIS 90cm solar telescope. This system is designed to operate in the visible and is required to be as robust as possible in order to deliver the best possible correction in any atmospheric conditions, even if wavefronts are sensed on some low-contrast solar granulation. In extreme conditions, the images given by the subapertures of the Shack-Hartmann wavefront sensor get randomly blurred in space, in the set of subapertures, and the distribution of blurred images is rapidly changing in time, some of them possibly fading away. The algorithms we have developed for such harsh conditions rely on inverse problem approach. As an example, with the gradients of the wavefronts, the wavefront sensor also estimates their errors, including their covariance. This information allows the control loop to promptly optimize itself to the fast varying conditions, both in space (wavefront reconstruction) and in time. A major constraint is to fit the calculations in a low-cost multi-core CPU. An overview of the algorithms in charge of implementing this strategy is presented, focusing on wavefront sensing.

The convergence of many proximal algorithms involving a gradient descent relies on its Lipschitz constant. To avoid computing it, backtracking rules can be used. While such a rule has already been designed for the forward-backward algorithm (FBwB), this scheme is not flexible enough when a non-differentiable penalization with a linear operator is added to a constraint. In this work we propose a backtracking rule for the primal-dual scheme (PDwB), and evaluate its performance for the epigraphical constrained high dynamical reconstruction in high contrast polarimetric imaging, under TV penalization., Comment: in Proceedings of iTWIST'20, Paper-ID: 12, Nantes, France, December, 2-4, 2020

During the last two decades, the first generation of beam combiners at the Very Large Telescope Interferometer has proved the importance of optical interferometry for high-angular resolution astrophysical studies in the near- and mid-infrared. With the advent of 4-beam combiners at the VLTI, the u-v coverage per pointing increases significantly, providing an opportunity to use reconstructed images as powerful scientific tools. Therefore, interferometric imaging is already a key feature of the new generation of VLTI instruments, as well as for other interferometric facilities like CHARA and JWST. It is thus imperative to account for the current image reconstruction capabilities and their expected evolutions in the coming years. Here, we present a general overview of the current situation of optical interferometric image reconstruction with a focus on new wavelength-dependent information, highlighting its main advantages and limitations. As an Appendix we include several cookbooks describing the usage and installation of several state-of-the art image reconstruction packages. To illustrate the current capabilities of the software available to the community, we recovered chromatic images, from simulated MATISSE data, using the MCMC software SQUEEZE. With these images, we aim at showing the importance of selecting good regularization functions and their impact on the reconstruction., Comment: Accepted for publication in Experimental Astronomy as part of the topical collection: Future of Optical-infrared Interferometry in Europe

We present a new formulation of a family of proximity operators that generalize the projector step for phase retrieval. These proximity operators for noisy intensity measurements can replace the classical "noise free" projection in any projection-based algorithm. They are derived from a maximum likelihood formulation and admit closed form solutions for both the Gaussian and the Poisson cases. In addition, we extend these proximity operators to undersampled intensity measurements. To assess their performance, these operators are exploited in a classical Gerchberg Saxton algorithm. We present numerical experiments showing that the reconstructed complex amplitudes with these proximity operators perform always better than using the classical intensity projector while their computational overhead is moderate.

This paper provides a general introduction to the problem of image reconstruction from interferometric data. A simple model of the interferometric observables is given and the issues arising from sparse Fourier data are discussed. The effects of various regularizations are described. In the proposed general framework, most existing algorithms can be understood. For an astronomer, such an understanding is crucial not only for selecting and using an algorithm but also to ensure correct interpretation of the resulting image., Comment: 20 pages, 5 figures

The detection of exoplanets in coronographic images is severely limited by residual starlight speckles. Dedicated post-processing can drastically reduce this "stellar leakage" and thereby increase the faintness of detectable exoplanets. Based on a multi-spectral series expansion of the diffraction pattern, we derive a multi-mode model of the residuals which can be exploited to estimate and thus remove the residual speckles in multi-spectral coronographic images. Compared to other multi-spectral processing methods, our model is physically grounded and is suitable for use in an (optimal) inverse approach. We demonstrate the ability of our model to correctly estimate the speckles in simulated data and demonstrate that very high contrasts can be achieved. We further apply our method to removing speckles from a real data cube obtained with the SPHERE IFS instrument., Comment: accepted for publication in MNRAS on 25th of August 2017, 17 pages, 15 figures

High spatial resolution is the key for the understanding various astrophysical phenomena. But even with the future E-ELT, single dish instruments are limited to a spatial resolution of about 4 mas in the visible. For the closest objects within our Galaxy most of the stellar photosphere remains smaller than 1 mas. With the success of long baseline interferometry these limitations were soom overcome. Today low and high resolution interferometric instruments on the VLTI and CHARA offer an immense range of astrophysical studies. Combining more telescopes and moving to visible wavelengths broadens the science cases even more. With the idea of developing strong science cases for a future visible interferometer, we organized a science group around the following topics: pre-main sequence and main sequence stars, fundamental parameters, asteroseismology and classical pulsating stars, evolved stars, massive stars, active galactic nuclei (AGNs) and imaging techniques. A meeting was organized on the 15th and 16th of January, 2015 in Nice with the support of the Action Specific in Haute Resolution Angulaire (ASHRA), the Programme National en Physique Stellaire (PNPS), the Lagrange Laboratory and the Observatoire de la Cote d'Azur, in order to present these cases and to discuss them further for future visible interferometers. This White Paper presents the outcome of the exchanges. This book is dedicated to the memory of our colleague Olivier Chesneau who passed away at the age of 41., Comment: White Paper prospective. This book is dedicated to the memory of our colleague Olivier Chesneau who passed away at the age of 41. v2 includes some corrections to text

Assessing the quality of aperture synthesis maps is relevant for benchmarking image reconstruction algorithms, for the scientific exploitation of data from optical long-baseline interferometers, and for the design/upgrade of new/existing interferometric imaging facilities. Although metrics have been proposed in these contexts, no systematic study has been conducted on the selection of a robust metric for quality assessment. This article addresses the question: what is the best metric to assess the quality of a reconstructed image? It starts by considering several metrics, and selecting a few based on general properties. Then, a variety of image reconstruction cases is considered. The observational scenarios are phase closure and phase referencing at the Very Large Telescope Interferometer (VLTI), for a combination of two, three, four and six telescopes. End-to-end image reconstruction is accomplished with the MiRA software, and several merit functions are put to test. It is found that convolution by an effective point spread function is required for proper image quality assessment. The effective angular resolution of the images is superior to naive expectation based on the maximum frequency sampled by the array. This is due to the prior information used in the aperture synthesis algorithm and to the nature of the objects considered. The l1 norm is the most robust of all considered metrics, because being linear it is less sensitive to image smoothing by high regularisation levels. For the cases considered, this metric allows the implementation of automatic quality assessment of reconstructed images, with a performance similar to human selection., Comment: Key words: instrumentation: high angular resolution - instrumentation: interferosmeters - methods: data analysis - techniques: high angular resolution - techniques: image processing - techniques: interferometric. 18 pages, 8 figures, 4 tables

(abridged) Extreme adaptive optics (XAO) encounters severe difficulties to cope with the high speed (>1kHz), high accuracy and high order requirements for future extremely large telescopes. An innovative high order adaptive optics system using a self-referenced Mach-Zehnder wavefront sensor (MZWFS) allows counteracting these limitations. This sensor estimates very accurately the wavefront phase at small spatial scale by measuring intensity differences between two outputs, with a $\lambda /4$ path length difference between its two legs, but is limited in dynamic range due to phase ambiguity. During the past few years, such an XAO system has been studied by our team in the framework of 8-meter class telescopes. In this work, we report on our latest results with the XAO testbed recently installed in our lab, and dedicated to high contrast imaging with 30m-class telescopes (such as the E-ELT or the TMT). After reminding the principle of a MZWFS and describing the optical layout of our experiment, we will show the results of the assessment of the woofer-tweeter phase correctors, i.e., a Boston Micromachine continuous membrane deformable mirror (DM) and a Boulder Nonlinear Systems liquid crystal spatial light modulator (SLM). In particular, we will detail the calibration of the DM using Zygo interferometer metrology. Our method consists in the precise measurement of the membrane deformation while applying a constant deformation to 9 out of 140 actuators at the same time. By varying the poke voltage across the DM operating range, we propose a simple but efficient way of modeling the DM influence function using a Gaussian model. Finally, we show the DM flattening on the MZWFS allowing to compensate for low order aberrations., Comment: To appear in SPIE proceedings vol. 9617

Astronomical optical interferometers (OI) sample the Fourier transform of the intensity distribution of a source at the observation wavelength. Because of rapid atmospheric perturbations, the phases of the complex Fourier samples (visibilities) cannot be directly exploited , and instead linear relationships between the phases are used (phase closures and differential phases). Consequently, specific image reconstruction methods have been devised in the last few decades. Modern polychromatic OI instruments are now paving the way to multiwavelength imaging. This paper presents the derivation of a spatio-spectral ("3D") image reconstruction algorithm called PAINTER (Polychromatic opticAl INTErferometric Reconstruction software). The algorithm is able to solve large scale problems. It relies on an iterative process, which alternates estimation of polychromatic images and of complex visibilities. The complex visibilities are not only estimated from squared moduli and closure phases, but also from differential phases, which help to better constrain the polychromatic reconstruction. Simulations on synthetic data illustrate the efficiency of the algorithm., Comment: EUSIPCO, Aug 2015, NICE, France

Astronomical optical interferometers sample the Fourier transform of the intensity distribution of a source at the observation wavelength. Because of rapid perturbations caused by atmospheric turbulence, the phases of the complex Fourier samples (visibilities) cannot be directly exploited. Consequently, specific image reconstruction methods have been devised in the last few decades. Modern polychromatic optical interferometric instruments are now paving the way to multiwavelength imaging. This paper is devoted to the derivation of a spatio-spectral (3D) image reconstruction algorithm, coined PAINTER (Polychromatic opticAl INTErferometric Reconstruction software). The algorithm relies on an iterative process, which alternates estimation of polychromatic images and of complex visibilities. The complex visibilities are not only estimated from squared moduli and closure phases, but also differential phases, which helps to better constrain the polychromatic reconstruction. Simulations on synthetic data illustrate the efficiency of the algorithm and in particular the relevance of injecting a differential phases model in the reconstruction., Comment: 12 pages, 10 figures, http://www.opticsinfobase.org/submit/review/copyright_permissions.cfm

Optical interferometers provide multiple wavelength measurements. In order to fully exploit the spectral and spatial resolution of these instruments, new algorithms for image reconstruction have to be developed. Early attempts to deal with multi-chromatic interferometric data have consisted in recovering a gray image of the object or independent monochromatic images in some spectral bandwidths. The main challenge is now to recover the full 3-D (spatio-spectral) brightness distribution of the astronomical target given all the available data. We describe a new approach to implement multi-wavelength image reconstruction in the case where the observed scene is a collection of point-like sources. We show the gain in image quality (both spatially and spectrally) achieved by globally taking into account all the data instead of dealing with independent spectral slices. This is achieved thanks to a regularization which favors spatial sparsity and spectral grouping of the sources. Since the objective function is not differentiable, we had to develop a specialized optimization algorithm which also accounts for non-negativity of the brightness distribution., Comment: This version has been accepted for publication in J. Opt. Soc. Am. A

We present the results of the fifth Interferometric Imaging Beauty Contest. The contest consists in blind imaging of test data sets derived from model sources and distributed in the OIFITS format. Two scenarios of imaging with CHARA/MIRC-6T were offered for reconstruction: imaging a T Tauri disc and imaging a spotted red supergiant. There were eight different teams competing this time: Monnier with the software package MACIM; Hofmann, Schertl and Weigelt with IRS; Thi\'ebaut and Soulez with MiRA ; Young with BSMEM; Mary and Vannier with MIROIRS; Millour and Vannier with independent BSMEM and MiRA entries; Rengaswamy with an original method; and Elias with the radio-astronomy package CASA. The contest model images, the data delivered to the contestants and the rules are described as well as the results of the image reconstruction obtained by each method. These results are discussed as well as the strengths and limitations of each algorithm.

With the advent of infrared long-baseline interferometers with more than two telescopes, both the size and the completeness of interferometric data sets have significantly increased, allowing images based on models with no a priori assumptions to be reconstructed. Our main objective is to analyze the multiple parameters of the image reconstruction process with particular attention to the regularization term and the study of their behavior in different situations. The secondary goal is to derive practical rules for the users. Using the Multi-aperture image Reconstruction Algorithm (MiRA), we performed multiple systematic tests, analyzing 11 regularization terms commonly used. The tests are made on different astrophysical objects, different (u,v) plane coverages and several signal-to-noise ratios to determine the minimal configuration needed to reconstruct an image. We establish a methodology and we introduce the mean-square errors (MSE) to discuss the results. From the ~24000 simulations performed for the benchmarking of image reconstruction with MiRA, we are able to classify the different regularizations in the context of the observations. We find typical values of the regularization weight. A minimal (u,v) coverage is required to reconstruct an acceptable image, whereas no limits are found for the studied values of the signal-to-noise ratio. We also show that super-resolution can be achieved with increasing performance with the (u,v) coverage filling. Using image reconstruction with a sufficient (u,v) coverage is shown to be reliable. The choice of the main parameters of the reconstruction is tightly constrained. We recommend that efforts to develop interferometric infrastructures should first concentrate on the number of telescopes to combine, and secondly on improving the accuracy and sensitivity of the arrays., Comment: 15 pages, 16 figures; accepted in A&A

We present a new algorithm, FRiM (FRactal Iterative Method), aiming at the reconstruction of the optical wavefront from measurements provided by a wavefront sensor. As our application is adaptive optics on extremely large telescopes, our algorithm was designed with speed and best quality in mind. The latter is achieved thanks to a regularization which enforces prior statistics. To solve the regularized problem, we use the conjugate gradient method which takes advantage of the sparsity of the wavefront sensor model matrix and avoids the storage and inversion of a huge matrix. The prior covariance matrix is however non-sparse and we derive a fractal approximation to the Karhunen-Loeve basis thanks to which the regularization by Kolmogorov statistics can be computed in O(N) operations, N being the number of phase samples to estimate. Finally, we propose an effective preconditioning which also scales as O(N) and yields the solution in 5-10 conjugate gradient iterations for any N. The resulting algorithm is therefore O(N). As an example, for a 128 x 128 Shack-Hartmann wavefront sensor, FRiM appears to be more than 100 times faster than the classical vector-matrix multiplication method., Comment: to appear in the Journal of the Optical Society of America A

This tutorial paper describes the problem of image reconstruction from interferometric data with a particular focus on the specific problems encountered at optical (visible/IR) wavelengths. The challenging issues in image reconstruction from interferometric data are introduced in the general framework of inverse problem approach. This framework is then used to describe existing image reconstruction algorithms in radio interferometry and the new methods specifically developed for optical interferometry., Comment: accepted for publication in IEEE Signal Processing Magazine

We investigate the problem of probing the local spatial structure of the magnetic field of the interstellar medium using multi-frequency polarized maps of the synchrotron emission at radio wavelengths. We focus in this paper on the three-dimensional reconstruction of the largest scales of the magnetic field, relying on the internal depolarization (due to differential Faraday rotation) of the emitting medium as a function of electromagnetic frequency. We argue that multi-band spectroscopy in the radio wavelengths, developed in the context of high-redshift extragalactic HI lines, can be a very useful probe of the 3D magnetic field structure of our Galaxy when combined with a Maximum A Posteriori reconstruction technique. When starting from a fair approximation of the magnetic field, we are able to recover the true one by using a linearized version of the corresponding inverse problem. The spectral analysis of this problem allows us to specify the best sampling strategy in electromagnetic frequency and predicts a spatially anisotropic distribution of posterior errors. The reconstruction method is illustrated for reference fields extracted from realistic magneto-hydrodynamical simulations.

Diffuse envelopes around Mira variables are among the most important sources influencing the chemical evolution of galaxies. However they represent an observational challenge because of their complex spectral features and their rapid temporal variability. We constrained the exact brightness distribution of the Mira star TLep with a model-independent analysis. We obtained single-epoch interferometric observations with a dataset continuous in the spectral domain (1.5-2.4mum) and in the spatial domain (baselines ranging from 11 to 96m). We performed a model independent image reconstruction for each spectral bin using the MIRA software. We completed the analysis by modeling the data with a simple star+layer model inspired from the images. Reconstructed images confirm the general picture of a central star partially obscured by the surrounding molecular shell of changing opacity. At 1.7mum, the shell becomes optically thin, with corresponding emission appearing as a ring circling the star. This is the first direct evidence of the spherical morphology of the molecular shell. Model fitting confirmed a spherical layer of constant size and changing opacity over the wavelengths. Rough modeling points to a continuum opacity within the shell, in addition to the CO and H2O features. Accordingly, it appeared impossible to model the data by a photosphere alone in any of the spectral bins., Comment: Accepted in A&A

One of the aims of next generation optical interferometric instrumentation is to be able to make use of information contained in the visibility phase to construct high dynamic range images. Radio and optical interferometry are at the two extremes of phase corruption by the atmosphere. While in radio it is possible to obtain calibrated phases for the science objects, in the optical this is currently not possible. Instead, optical interferometry has relied on closure phase techniques to produce images. Such techniques allow only to achieve modest dynamic ranges. However, with high contrast objects, for faint targets or when structure detail is needed, phase referencing techniques as used in radio interferometry, should theoretically achieve higher dynamic ranges for the same number of telescopes. Our approach is not to provide evidence either for or against the hypothesis that phase referenced imaging gives better dynamic range than closure phase imaging. Instead we wish to explore the potential of this technique for future optical interferometry and also because image reconstruction in the optical using phase referencing techniques has only been performed with limited success. We have generated simulated, noisy, complex visibility data, analogous to the signal produced in radio interferometers, using the VLTI as a template. We proceeded with image reconstruction using the radio image reconstruction algorithms contained in AIPS IMAGR (CLEAN algorithm). Our results show that image reconstruction is successful in most of our science cases, yielding images with a 4 milliarcsecond resolution in K band. (abridged), Comment: 11 pages, 36 figures

Classically, optical and near-infrared interferometry have relied on closure phase techniques to produce images. Such techniques allow us to achieve modest dynamic ranges. In order to test the feasibility of next generation optical interferometers in the context of the VLTI-spectro-imager (VSI), we have embarked on a study of image reconstruction and analysis. Our main aim was to test the influence of the number of telescopes, observing nights and distribution of the visibility points on the quality of the reconstructed images. Our results show that observations using six Auxiliary Telescopes (ATs) during one complete night yield the best results in general and is critical in most science cases; the number of telescopes is the determining factor in the image reconstruction outcome. In terms of imaging capabilities, an optical, six telescope VLTI-type configuration and ~200 meter baseline will achieve 4 mas spatial resolution, which is comparable to ALMA and almost 50 times better than JWST will achieve at 2.2 microns. Our results show that such an instrument will be capable of imaging, with unprecedented detail, a plethora of sources, ranging from complex stellar surfaces to microlensing events., Comment: 11 pages, 26 figures

Instrumental projects that will improve the direct optical finding and characterisation of exoplanets have advanced sufficiently to trigger organized investigation and development of corresponding signal processing algorithms. The first step is the availability of field-of-view (FOV) models. These can then be submitted to various instrumental models, which in turn produce simulated data, enabling the testing of processing algorithms. We aim to set the specifications of a physical model for typical FOVs of these instruments. The dynamic in resolution and flux between the various sources present in such a FOV imposes a multiscale, independent layer approach. From review of current literature and through extrapolations from currently available data and models, we derive the features of each source-type in the field of view likely to pass the instrumental filter at exo-Earth level. Stellar limb darkening is shown to cause bias in leakage calibration if unaccounted for. Occurrence of perturbing background stars or galaxies in the typical FOV is unlikely. We extract galactic interstellar medium background emissions for current target lists. Galactic background can be considered uniform over the FOV, and it should show no significant drift with parallax. Our model specifications have been embedded into a Java simulator, soon to be made open-source. We have also designed an associated FITS input/output format standard that we present here., Comment: 9 pages (+5 of appendices), 7 figures, accepted for publication in Astronomy & Astrophysics

VITRUV is a second generation spectro-imager for the PRIMA enabled Very Large Telescope Interferometer. By combining simultaneously up to 8 telescopes VITRUV makes the VLTI up to 6 times more efficient. This operational gain allows two novel scientific methodologies: 1) massive surveys of sizes; 2) routine interferometric imaging. The science cases presented concentrate on the qualitatively new routine interferometric imaging methodology. The science cases are not exhaustive but complementary to the PRIMA reference mission. The focus is on: a) the close environment of young stars probing for the initial conditions of planet formation and disk evolution; b) the surfaces of stars tackling dynamos, activity, pulsation, mass-loss and evolution; c) revealing the origin of the extraordinary morphologies of Planetary Nebulae and related stars; d) studying the accretion-ejection structures of stellar black-holes (microquasars) in our galaxy; e) unveiling the different interacting components (torus, jets, BLRs) of Active Galactic Nuclei; and f) probing the environment of nearby supermassive black-holes and relativistic effects in the Galactic Center black-hole., Comment: 15 pages. The Power of Optical/IR Interferometry: Recent Scientific Results and 2nd Generation VLTI Instrumentation, Allemagne (2005) in press

The VITRUV project has the objective to deliver milli-arcsecond spectro-images of the environment of compact sources like young stars, active galaxies and evolved stars to the community. This instrument of the VLTI second generation based on the integrated optics technology is able to combine from 4 to 8 beams from the VLT telescopes. Working primarily in the near infrared, it will provide intermediate to high spectral resolutions and eventually polarization analysis. This paper summarizes the result from the concept study led within the Joint Research Activity advanced instruments of the OPTICON program., Comment: In "The Power of Optical/IR Interferometry: Recent Scientific Results and 2nd Generation VLTI Instrumentation", Allemagne (2005) in press