Abstract
Context.Direct imaging is a method of choice for probing the close environment of young stars. Even with the coupling of adaptive optics and coronagraphy, the direct detection of off-axis sources such as circumstellar disks and exoplanets remains challenging due to the required high contrast and small angular resolution. Angular differential imaging (ADI) is an observational technique that introduces an angular diversity to help disentangle the signal of off-axis sources from the residual signal of the star in a post-processing step.Aims.While various detection algorithms have been proposed in the last decade to process ADI sequences and reach high contrast for the detection of point-like sources, very few methods are available to reconstruct meaningful images of extended features such as circumstellar disks. The purpose of this paper is to describe a new post-processing algorithm dedicated to the reconstruction of the spatial distribution of light (total intensity) received from off-axis sources, in particular from circumstellar disks.Methods.Built on the recentPACOalgorithm 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 problems 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 byREXPACOcompared to standard post-processing methods such as cADI. In particular, the modeling of the spatial covariances proves useful 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 ofREXPACOto produce images of the light distribution with reduced artifacts. Finally, we show howREXPACOcan be combined withPACOto disentangle the signal of circumstellar disks from the signal of candidate point-like sources.Conclusions.REXPACOis a novel post-processing algorithm for reconstructing images of the circumstellar environment from high contrast ADI sequences. It produces numerically deblurred images and exploits the spatial covariances of the stellar leakages and of the noise to efficiently eliminate this nuisance term. The processing is fully unsupervised, all tuning parameters being directly estimated from the data themselves.
Highlights
Circumstellar disks are at the heart of the planet formation processes
We develop a regularized image reconstruction algorithm (REXPACO) following an inverse problems approach based on a forward image formation model of the off-axis sources in the Angular differential imaging (ADI) sequences
We have introduced a new post-processing algorithm, named REXPACO, dedicated to the reconstruction of flux distributions from ADI data sets in the presence of extended features, like circumstellar disks
Summary
Circumstellar disks are at the heart of the planet formation processes. Direct imaging in the near infrared is a unique method for addressing this question because planets can be detected together with the disk environment. Recent protoplanetary and debris disk studies performed in intensity or in polarimetry have focused on morphological disk characteristics such as the presence of spirals (Benisty et al 2015; Ren et al 2018a; Muro-Arena et al 2020), asymmetries, warps (Dawson et al 2011; Kluska et al 2020) and gaps (van Boekel et al 2017) that could be signposts for the presence of exoplanets. In this context, protoplanetary disks allow unique studies of the exoplanet-disk interactions (Keppler et al 2018; Haffert et al 2019; Mesa et al 2019). In order to understand the physical processes governing these disks, it is essential to reconstruct their surface-brightness or scattering phase function and to disentangle, in a postprocessing step, the disks and the exoplanets
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.