Approche inverse pour la reconstruction des environnements circumstellaires en polarimétrie avec l'instrument d'imagerie directe ESO / VLT SPHERE IRDIS
Institution:
LyonDisciplines:
Directors:
Abstract EN:
Circumstellars environments observation is a key for the comprehension of planet formation. If the very large telescopes allow the resolution of these environments, their observation is difficult due to the high contrast between the environment and their host stars. In fact the host stars are 1000 to 10 000 times brighter than the environment, even 10 000 000 times brighter for exoplanets. When images of these circumstellar environnements are acquired in direct imaging, the signal of the environnements mixed to star light residuals. Yet, the light of the environment is partially linearly polarized while the light od the star is unpolarized. The instrument Infrared Dual-band Imaging and Spectroscopy (IRDIS) of the European Southern Observatory’s (ESO) Spectro-Polarimeter High-contrast Expolanet REsearch (SPHERE) instrument, installed at one of the four Very Large Telescopes (VLT) in Atacama in Chile, acquires datasets where the polarization is modulated according to a known angles cycle. It is then possible, by combinations of the data, to extract the polarized signal of the environment from the unpolarized residual light of the stars and unpolarized light of the disks. The stat-of-the-art methods to extract such signal do not take optimally into account the photon noise statistics of the data, which dominate the signal of interest, nor the read out noise of the detector. Moreover, if any image from a rotation cycle is missing, the rest of the cycle is not used. Finally, any centering and rotation of the data or deconvolution by the PSF is generally performed in separated steps from the data reduction. The bad pixels and dead pixels are interpolated before the processing. The consequence of such approach is that the propagation of the errors in the data is not controlled.The « inverse problem » methods allow such processing while controlling the error propagation in the reconstructions. These approaches have never been developed, so far, for high contrast direct imaging in polarimetry. My goal in this thesis is to optimally reconstruct, from the polarimetric data of the instrument ESO/VLT-SPHERE IRDIS, maps of the circumstellar environments polarized light, the ascociated polarization angles and the unpolarized star light residuals and circumstellar environments light. First, I develop a nonlinear physical model of the data, pixelwise independent, parametric in these quantities of interest, or linear with respect to the Stokes parameters, from which they can be estimated. Throughout this thesis, I complete the model by adding centering, rotations and convolutions, making it pixelwise dependent. The parameters are then estimated by the minimization of an objective function, derived from the co-log-likelyhood of the data, under some constraint, such as positivity constraint or epigraphical constraint, and regularizations as smooth and non-smooth Total Variation and the Shatten norm on the Hessian. This methods are all applied on simulated datasets, created to reproduce typical astrophysical datasets obtained in circumstellar environment polarimetrical direct imaging. Depending of the properties of the functions considered in the objective function, the research of its minimum is done with different algorithms as the Variable Metric Limited Memory and Bound algorithm, Forward-Backward with backtracking and the preconditioned primal-dual Condat-Vu algorithm with backtracking. I also use the Stein Unbiased Risk Estimator to auto-tune the weights of the regularization. In the results, I show that the use of a complete direct model of the data, taking in account the recentering, the rotations and the convolution and the estimation of its parameters from a constraint problem, taking in account the measure precision and the missing data reduces the error on the estimation maps in such astrophysics context
Abstract FR:
L'étude des environnements circumstellaires permet d'en apprendre plus sur la formation des exoplanètes. Si les avancées instrumentales permettent la résolution de ces environnements, leur observation en imagerie directe est difficile du fait du grand contraste entre les environnements et leurs étoiles hôtes. En effet, celles-ci sont 1000 à 10 000 fois plus brillantes, ou encore 10 000 000 fois dans le cas des exoplanètes. Lors de l'acquisition directe d'image de ces environnements, leur signal est mélangé au résidu de lumière stellaire. Or, la lumière de l'environnements circumstellaires est partiellement polarisée linéairement, tandis que le résidu de lumière stellaire ne l'est pas. Le sous-instrument Infrared Dual-band Imaging and Spectroscopy (IRDIS) de l'instrument de l'European Southern Observatory (ESO) appelé Spectro-Polarimeter High-contrast Expolanet REsearch (SPHERE), situé sur l'un des quatre Very Large Telescopes (VLT) dans le désert d'Atacama au Chili, acquiert des jeux de données où la polarisation linéaire est modulée en rotation, selon plusieurs cycles d'angles connus. Ainsi, par combinaison de ces données, il est possible de démélanger la lumière diffusée par l'environnement circumstellaire et la lumière de l'étoile. Les méthodes de l'état-de-l'art permettant ce démélange, ne prennent en compte ni le bruit de photon, qui domine le signal d'intérêt, ni le bruit de lecture du détecteur. De plus, si une image du cycle de rotation est manquante, les autres images du cycle sont supprimées. Par ailleurs, le traitement par interpolation des pixels morts, le recentrage des images et toutes éventuelles rotations et déconvolutions des données, sont fait indépendamment du démélange. De ce fait, la propagation des erreurs n'est pas contrôlée.Les méthodes de types « problèmes inverses » permettent, à partir d'un modèle direct des données, de procéder au démélange tout en ayant le contrôle sur la propagation des erreurs dans les reconstructions. Une telle approche n'a jamais été développée pour l'imagerie directe en polarimétrie. Le but de ma thèse, est de reconstruire de manière optimale, à partir des données polarimétriques de l’instrument ESO/VLT-SPHERE IRDIS, des cartes de la lumière polarisée des environnements circumstellaires, des angles de polarisation associés et des résidus lumineux de l'étoile et de la lumière non-polarisée de l'environnement. Je propose tout d’abord un modèle physique séparable des données, paramétré non-linéairement en ces quantité d’intérêt et linéairement en les paramètres de Stokes, à partir duquel ces quantités peuvent être estimées. Ensuite, je complexifie le modèle en y incluant recentrage, rotations et convolution, le rendant ainsi non-séparable. Les paramètres sont alors estimés par la minimisation de la co-log-vraissemblance des données, sous contraintes, de positivité dans le cas non-linéaire ou épigraphique dans le cas linéaire. Dans le cas non-séparable, un terme de régularisation est également ajouté, comme la variation totale et la norme de Shatten sur le Hessien. Selon les propriétés des fonctions considérées dans le critère, je procède à sa minimisation par différents algorithmes tels que l’algorithme Variable Metric Limited Memory and Bound, Forward-Backward avec backtracking et l’algorithme primal-dual préconditionné Condat-Vu avec backtracking. Je propose également une méthode de réglage automatique des poids des régularisations à partir de l’estimateur non-biaisé du risque de Stein (SURE). L’ensemble de ces méthodes d'estimation sont appliquées sur des données synthétiques et sur données astrophysiques. Je montre alors que l’utilisation d’un modèle directe complet des données, prenant en compte le recentrage, les rotations et la convolution, ainsi que l’estimation de ses paramètres à partir d’un critère régularisé sous contrainte, prenant en compte la précision des mesures et les données manquantes, permet de réduire l’erreur faite sur l’estimation des observables