The direct detection and characterization of planetary and substellar companions at small angular separations is a rapidly advancing field. Dedicated high-contrast imaging instruments deliver unprecedented sensitivity, enabling detailed insights into the atmospheres of young low-mass companions. In addition, improvements in data reduction and PSF subtraction algorithms are equally relevant for maximizing the scientific yield, both from new and archival data sets. We aim at developing a generic and modular data reduction pipeline for processing and analysis of high-contrast imaging data obtained with pupil-stabilized observations. The package should be scalable and robust for future implementations and in particular well suitable for the 3-5 micron wavelength range where typically (ten) thousands of frames have to be processed and an accurate subtraction of the thermal background emission is critical. PynPoint is written in Python 2.7 and applies various image processing techniques, as well as statistical tools for analyzing the data, building on open-source Python packages. The current version of PynPoint has evolved from an earlier version that was developed as a PSF subtraction tool based on PCA. The architecture of PynPoint has been redesigned with the core functionalities decoupled from the pipeline modules. Modules have been implemented for dedicated processing and analysis steps, including background subtraction, frame registration, PSF subtraction, photometric and astrometric measurements, and estimation of detection limits. The pipeline package enables end-to-end data reduction of pupil-stabilized data and supports classical dithering and coronagraphic data sets. As an example, we processed archival VLT/NACO L' and M' data of beta Pic b and reassessed the planet's brightness and position with an MCMC analysis, and we provide a derivation of the photometric error budget.

Context. High-contrast exoplanet imaging is a rapidly growing field as can be seen through the significant resources invested. In fact, the detection and characterization of exoplanets through direct imaging is featured at all major ground-based observatories. Aims. We aim to improve the signal-to-noise ratio (SNR) achievable for ground-based, adaptive-optics assisted exoplanet imaging by applying sophisticated post-processing algorithms. In particular, we investigate the benefits of including time domain information. Methods. We introduce a new speckle-suppression technique in data post-processing based on wavelet transformation. This technique explicitly considers the time domain in a given data set (specifically the frequencies of speckle variations and their time dependence) and allows us to filter-out speckle noise. We combine our wavelet-based algorithm with state-of-the-art principal component analysis (PCA) based PSF subtraction routines and apply it to archival data sets of known directly imaged exoplanets. The data sets were obtained in the L filter where the short integration times allow for a sufficiently high temporal sampling of the speckle variations. Results. We demonstrate that improvements in the peak SNR of up to forty to sixty percent can be achieved. We also show that, when combined with wavelet-denoising, the PCA PSF model requires systematically smaller numbers of components for the fit to achieve the highest SNR. The improvement potential is, however, data set dependent or, more specifically, closely linked to the field rotation available in a given data set: larger amounts of rotation allow for a better suppression of the speckle noise. Conclusions. We have demonstrated that by applying advanced data post-processing techniques, the contrast performance in archival high-contrast imaging data sets can be improved.