Spectro-Perfectionism: An Algorithmic Framework for Photon Noise-Limited Extraction of Optical Fiber Spectroscopy

Abstract

We describe a new algorithm for the "perfect" extraction of one-dimensional spectra from two-dimensional (2D) digital images of optical fiber spectrographs, based on accurate 2D forward modeling of the raw pixel data. The algorithm is correct for arbitrarily complicated 2D point-spread functions (PSFs), as compared to the traditional optimal extraction algorithm, which is only correct for a limited class of separable PSFs. The algorithm results in statistically independent extracted samples in the 1D spectrum, and preserves the full native resolution of the 2D spectrograph without degradation. Both the statistical errors and the 1D resolution of the extracted spectrum are accurately determined, allowing a correct chi-squared comparison of any model spectrum with the data. Using a model PSF similar to that found in the red channel of the Sloan Digital Sky Survey spectrograph, we compare the performance of our algorithm to that of cross-section based optimal extraction, and also demonstrate that our method allows coaddition and foreground estimation to be carried out as an integral part of the extraction step. This work demonstrates the feasibility of current- and next-generation multi-fiber spectrographs for faint galaxy surveys even in the presence of strong night-sky foregrounds. We describe the handling of subtleties arising from fiber-to-fiber crosstalk, discuss some of the likely challenges in deploying our method to the analysis of a full-scale survey, and note that our algorithm could be generalized into an optimal method for the rectification and combination of astronomical imaging data.

9 pages, 4 figures, emulateapj; minor corrections and clarifications; to be published in the PASP

Cited by in corpus (19)

• Stellar Population Inference with Prospector
• Overview of the Instrumentation for the Dark Energy Spectroscopic Instrument
• The Early Data Release of the Dark Energy Spectroscopic Instrument
• DESI Bright Galaxy Survey: Final Target Selection, Design, and Validation
• The Spectroscopic Data Processing Pipeline for the Dark Energy Spectroscopic Instrument
• The wide-field, multiplexed, spectroscopic facility WEAVE: Survey design, overview, and simulated implementation
• The BOSS Lyman-alpha Forest Sample from SDSS Data Release 9
• Performance Verification of the EXtreme PREcision Spectrograph
• An Extreme Precision Radial Velocity Pipeline: First Radial Velocities from EXPRES
• No escaping helium from 55 Cnc e
• ATOCA: an algorithm to treat order contamination. Application to the NIRISS SOSS mode
• Optimal 1D Ly-$\alpha$ Forest Power Spectrum Estimation I: DESI-Lite Spectra
• Characterization of the ESPRESSO Line-Spread Function and Improvement of the Wavelength Calibration Accuracy
• Simulating image coaddition with the Nancy Grace Roman Space Telescope: I. Simulation methodology and general results
• A framework to measure the properties of intergalactic metal systems with two-point flux statistics
• Covariance-regularized reconstruction of data cubes in integral field spectroscopy and application to MaNGA data
• Two-dimensional Multi-fiber Spectrum Image Correction Based on Machine Learning Techniques
• Review of detector requirements: some challenges for the present
• DESI DR1 Lyα 1D power spectrum: The Fast Fourier Transform estimator measurement