We present the Vortex Image Processing (VIP) library, a python package dedicated to astronomical high-contrast imaging. Our package relies on the extensive python stack of scientific libraries and aims to provide a flexible framework for high-contrast data and image processing. In this paper, we describe the capabilities of VIP related to processing image sequences acquired using the angular differential imaging (ADI) observing technique. VIP implements functionalities for building high-contrast data processing pipelines, encompass- ing pre- and post-processing algorithms, potential sources position and flux estimation, and sensitivity curves generation. Among the reference point-spread function subtraction techniques for ADI post-processing, VIP includes several flavors of principal component analysis (PCA) based algorithms, such as annular PCA and incremental PCA algorithm capable of processing big datacubes (of several gigabytes) on a computer with limited memory. Also, we present a novel ADI algorithm based on non-negative matrix factorization (NMF), which comes from the same family of low-rank matrix approximations as PCA and provides fairly similar results. We showcase the ADI capabilities of the VIP library using a deep sequence on HR8799 taken with the LBTI/LMIRCam and its recently commissioned L-band vortex coronagraph. Using VIP we investigated the presence of additional companions around HR8799 and did not find any significant additional point source beyond the four known planets. VIP is available at http://github.com/vortex-exoplanet/VIP and is accompanied with Jupyter notebook tutorials illustrating the main functionalities of the library.

Data processing constitutes a critical component of high-contrast exoplanet imaging. Its role is almost as important as the choice of a coronagraph or a wavefront control system, and it is intertwined with the chosen observing strategy. Among the data processing techniques for angular differential imaging (ADI), the most recent is the family of principal component analysis (PCA) based algorithms. PCA serves, in this case, as a subspace projection technique for constructing a reference point spread function (PSF) that can be subtracted from the science data for boosting the detectability of potential companions present in the data. Unfortunately, when building this reference PSF from the science data itself, PCA comes with certain limitations such as the sensitivity of the lower dimensional orthogonal subspace to non-Gaussian noise. Inspired by recent advances in machine learning algorithms such as robust PCA, we aim to propose a localized subspace projection technique that surpasses current PCA-based post-processing algorithms in terms of the detectability of companions at near real-time speed, a quality that will be useful for future direct imaging surveys. We used randomized low-rank approximation methods recently proposed in the machine learning literature, coupled with entry-wise thresholding to decompose an ADI image sequence locally into low-rank, sparse, and Gaussian noise components (LLSG). This local three-term decomposition separates the starlight and the associated speckle noise from the planetary signal, which mostly remains in the sparse term. We tested the performance of our new algorithm on a long ADI sequence obtained on beta Pictoris with VLT/NACO. Compared to a standard PCA approach, LLSG decomposition reaches a higher signal-to-noise ratio and has an overall better performance in the receiver operating characteristic space. (abridged).

(abridged) Vortex coronagraphs are among the most promising solutions to perform high contrast imaging at small angular separations. They feature a very small inner working angle, a clear 360 degree discovery space, have demonstrated very high contrast capabilities, are easy to implement on high-contrast imaging instruments, and have already been extensively tested on the sky. Since 2005, we have been designing, developing and testing an implementation of the charge-2 vector vortex phase mask based on concentric subwavelength gratings, referred to as the Annular Groove Phase Mask (AGPM). Science-grade mid-infrared AGPMs were produced in 2012 for the first time, using plasma etching on synthetic diamond substrates. They have been validated on a coronagraphic test bench, showing broadband peak rejection up to 500:1 in the L band, which translates into a raw contrast of about $6\times 10^{-5}$ at $2 \lambda/D$. Three of them have now been installed on world-leading diffraction-limited infrared cameras (VLT/NACO, VLT/VISIR and LBT/LMIRCam). During the science verification observations with our L-band AGPM on NACO, we observed the beta Pictoris system and obtained unprecedented sensitivity limits to planetary companions down to the diffraction limit ($0.1''$). More recently, we obtained new images of the HR 8799 system at L band during the AGPM first light on LMIRCam. After reviewing these first results obtained with mid-infrared AGPMs, we will discuss the short- and mid-term goals of the on-going VORTEX project, which aims to improve the performance of our vortex phase masks for future applications on second-generation high-contrast imagers and on future extremely large telescopes (ELTs).