Half-sibling regression meets exoplanet imaging: PSF modeling and subtraction using a flexible, domain knowledge-driven, causal framework
arXiv:2204.03439 · doi:10.1051/0004-6361/202142529
Abstract
High-contrast imaging of exoplanets hinges on powerful post-processing methods to denoise the data and separate the signal of a companion from its host star, which is typically orders of magnitude brighter. Existing post-processing algorithms do not use all prior domain knowledge that is available about the problem. We propose a new method that builds on our understanding of the systematic noise and the causal structure of the data-generating process. Our algorithm is based on a modified version of half-sibling regression (HSR), a flexible denoising framework that combines ideas from the fields of machine learning and causality. We adapt the method to address the specific requirements of high-contrast exoplanet imaging data obtained in pupil tracking mode. The key idea is to estimate the systematic noise in a pixel by regressing the time series of this pixel onto a set of causally independent, signal-free predictor pixels. We use regularized linear models in this work; however, other (non-linear) models are also possible. In a second step, we demonstrate how the HSR framework allows us to incorporate observing conditions such as wind speed or air temperature as additional predictors. When we apply our method to four data sets from the VLT/NACO instrument, our algorithm provides a better false-positive fraction than PCA-based PSF subtraction, a popular baseline method in the field. Additionally, we find that the HSR-based method provides direct and accurate estimates for the contrast of the exoplanets without the need to insert artificial companions for calibration in the data sets. Finally, we present first evidence that using the observing conditions as additional predictors can improve the results. Our HSR-based method provides an alternative, flexible and promising approach to the challenge of modeling and subtracting the stellar PSF and systematic noise in exoplanet imaging data.
arXiv commentAccepted for publication in Astronomy & Astrophysics
References in corpus (37)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- Array Programming with NumPy
- Astropy: A Community Python Package for Astronomy
- scikit-image: Image processing in Python
- Binary companions of evolved stars in APOGEE DR14: Search method and catalog of ~5,000 companions
- Angular Differential Imaging: a Powerful High-Contrast Imaging Technique
- astroquery: An Astronomical Web-Querying Package in Python
- Detection and Characterization of Exoplanets and Disks using Projections on Karhunen-Loeve Eigenimages
- A Giant Planet Candidate near a Young Brown Dwarf
- A new algorithm for point spread function subtraction in high-contrast imaging: a demonstration with angular differential imaging
- Fundamental limitations of high contrast imaging set by small sample statistics
- PynPoint: An Image Processing Package for Finding Exoplanets
- VIP: Vortex Image Processing package for high-contrast direct imaging
- Detection and Characterization of Exoplanets using Projections on Karhunen-Loeve Eigenimages: Forward Modeling
- The Structure of High Strehl Ratio Point-Spread Functions
- Direct exoplanet detection and characterization using the ANDROMEDA method: Performance on VLT/NaCo data
- The near-infrared spectral energy distribution of β Pictoris b
- Improving and Assessing Planet Sensitivity of the GPI Exoplanet Survey with a Forward Model Matched Filter
- Searching for companions down to 2 AU from beta Pictoris using the L'-band AGPM coronagraph on VLT/NACO
- Scientific Domain Knowledge Improves Exoplanet Transit Classification with Deep Learning
- Low-rank plus sparse decomposition for exoplanet detection in direct-imaging ADI sequences. The LLSG algorithm
- PynPoint: a modular pipeline architecture for processing and analysis of high-contrast imaging data
- Characterizing the performance of the NIRC2 vortex coronagraph at W.M. Keck Observatory
- Improving Signal to Noise in the Direct Imaging of Exoplanets and Circumstellar Disks
- Supervised detection of exoplanets in high-contrast imaging sequences
- A New Standard for Assessing the Performance of High Contrast Imaging Systems
- TRAP: A temporal systematics model for improved direct detection of exoplanets at small angular separations
- STIM map: detection map for exoplanets imaging beyond asymptotic Gaussian residual speckle noise
- A Causal, Data-Driven Approach to Modeling the Kepler Data
- Spatial Context Awareness for Unsupervised Change Detection in Optical Satellite Images
- RSM detection map for direct exoplanet detection in ADI sequences
- A High-contrast Imaging Algorithm: Optimized Image Rotation and Subtraction
- Improved Contrast in Images of Exoplanets using Direct SNR Optimization
- Fitting very flexible models: Linear regression with large numbers of parameters
- Fainter and closer: finding planets by symmetry breaking
- ISPY - the NaCo Imaging Survey for Planets around Young stars: A young companion candidate embedded in the R CrA cloud
- The Mysterious Lives Of Speckles. I. Residual atmospheric speckle lifetimes in ground-based coronagraphs