Cleaning Images with Gaussian Process Regression

Abstract

Many approaches to astronomical data reduction and analysis cannot tolerate missing data: corrupted pixels must first have their values imputed. This paper presents astrofix, a robust and flexible image imputation algorithm based on Gaussian Process Regression (GPR). Through an optimization process, astrofix chooses and applies a different interpolation kernel to each image, using a training set extracted automatically from that image. It naturally handles clusters of bad pixels and image edges and adapts to various instruments and image types. For bright pixels, the mean absolute error of astrofix is several times smaller than that of median replacement and interpolation by a Gaussian kernel. We demonstrate good performance on both imaging and spectroscopic data, including the SBIG 6303 0.4m telescope and the FLOYDS spectrograph of Las Cumbres Observatory and the CHARIS integral-field spectrograph on the Subaru Telescope.

arXiv commentPublished by the Astronomical Journal. 16 pages, 12 figures. The algorithm is implemented in the Python package astrofix, which is available at https://github.com/HengyueZ/astrofix

References in corpus (14)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• Astropy: A Community Python Package for Astronomy
• Binary companions of evolved stars in APOGEE DR14: Search method and catalog of ~5,000 companions
• Cosmic-Ray Rejection by Laplacian Edge Detection
• The Gemini Planet Imager: First Light
• Overview of the Kepler Science Processing Pipeline
• SPHERE: the exoplanet imager for the Very Large Telescope
• A Fast Algorithm for Cosmic Rays Removal from Single Images
• Gemini Planet Imager Observational Calibrations I: Overview of the GPI Data Reduction Pipeline
• Cosmic Ray Rejection by Linear Filtering of Single Images
• deepCR: Cosmic Ray Rejection with Deep Learning
• Conceptual Design of the Coronagraphic High Angular Resolution Imaging Spectrograph (CHARIS) for the Subaru Telescope
• Evaluation of cosmic ray rejection algorithms on single-shot exposures

Cited by in corpus (5)

• Likelihood-Based Jump Detection and Cosmic Ray Rejection for Detectors Read Out Up-the-Ramp
• Gaussian Process Classification for Galaxy Blend Identification in LSST
• Discovery of the Binarity of Gliese 229B, and Constraints on the System's Properties
• Can AI Dream of Unseen Galaxies? Conditional Diffusion Model for Galaxy Morphology Augmentation
• Photometry on Structured Backgrounds: Local Pixelwise Infilling by Regression