An Alternating Minimization Algorithm with Trajectory for Direct Exoplanet Detection -- The AMAT Algorithm

Abstract

Effective image post-processing algorithms are vital for the successful direct imaging of exoplanets. Standard PSF subtraction methods use techniques based on a low-rank approximation to separate the rotating planet signal from the quasi-static speckles, and rely on signal-to-noise ratio maps to detect the planet. These steps do not interact or feed each other, leading to potential limitations in the accuracy and efficiency of exoplanet detection. We aim to develop a novel approach that iteratively finds the flux of the planet and the low-rank approximation of quasi-static signals, in an attempt to improve upon current PSF subtraction techniques. In this study, we extend the standard L2 norm minimization paradigm to an L1 norm minimization framework to better account for noise statistics in the high contrast images. Then, we propose a new method, referred to as Alternating Minimization Algorithm with Trajectory, that makes a more advanced use of estimating the low-rank approximation of the speckle field and the planet flux by alternating between them and utilizing both L1 and L2 norms. For the L1 norm minimization, we propose using L1 norm low-rank approximation, a low-rank approximation computed using an exact block-cyclic coordinate descent method, while we use randomized singular value decomposition for the L2 norm minimization. Additionally, we enhance the visibility of the planet signal using a likelihood ratio as a postprocessing step. Numerical experiments performed on a VLT/SPHERE-IRDIS dataset show the potential of AMAT to improve upon the existing approaches in terms of higher S/N, sensitivity limits, and ROC curves. Moreover, for a systematic comparison, we used datasets from the exoplanet data challenge to compare our algorithm to other algorithms in the challenge, and AMAT with likelihood ratio map performs better than most algorithms tested on the exoplanet data challenge.

Submitted to A&A, 20 pages, 20 figures

References in corpus (29)

• Angular Differential Imaging: a Powerful High-Contrast Imaging Technique
• A giant planet imaged in the disk of the young star Beta Pictoris
• Detection and Characterization of Exoplanets and Disks using Projections on Karhunen-Loeve Eigenimages
• The Gemini Planet Imager: First Light
• SPHERE: the exoplanet imager for the Very Large Telescope
• 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
• Imaging Extrasolar Giant Planets
• PynPoint: An Image Processing Package for Finding Exoplanets
• VIP: Vortex Image Processing package for high-contrast direct imaging
• Spectral and atmospheric characterization of 51 Eridani b using VLT/SPHERE
• Direct exoplanet detection and characterization using the ANDROMEDA method: Performance on VLT/NaCo data
• Non-negative Matrix Factorization: Robust Extraction of Extended Structures
• Improving and Assessing Planet Sensitivity of the GPI Exoplanet Survey with a Forward Model Matched Filter
• The Orbit and Transit Prospects for $\beta$ Pictoris b constrained with One Milliarcsecond Astrometry
• Searching for companions down to 2 AU from beta Pictoris using the L'-band AGPM coronagraph on VLT/NACO
• VLT/SPHERE robust astrometry of the HR8799 planets at milliarcsecond-level accuracy Orbital architecture analysis with PyAstrOFit
• Low-rank plus sparse decomposition for exoplanet detection in direct-imaging ADI sequences. The LLSG algorithm
• Improving Signal to Noise in the Direct Imaging of Exoplanets and Circumstellar Disks
• MAYONNAISE: a morphological components analysis pipeline for circumstellar disks and exoplanets imaging in the near infrared
• Supervised detection of exoplanets in high-contrast imaging sequences
• STIM map: detection map for exoplanets imaging beyond asymptotic Gaussian residual speckle noise
• On the Complexity of Robust PCA and ${\ell }_1$-norm Low-Rank Matrix Approximation
• Comparing Apples with Apples: Robust Detection Limits for Exoplanet High-Contrast Imaging in the Presence of non-Gaussian Noise
• RSM detection map for direct exoplanet detection in ADI sequences
• Iterative Angular Differential Imaging (IADI): An exploration of recovering disk structures in scattered light with an iterative ADI approach
• Improved Contrast in Images of Exoplanets using Direct SNR Optimization
• Inverse-problem versus principal component analysis methods for angular differential imaging of circumstellar disks. The mustard algorithm
• NA-SODINN: a deep learning algorithm for exoplanet image detection based on residual noise regimes