Optimal extraction of echelle spectra: getting the most from observations

Abstract

The price of instruments and observing time on modern telescopes is quickly increasing with the size of the primary mirror. Therefore, it is worth revisiting the data reduction algorithms to extract every bit of scientific information from observations. Echelle spectrographs are typical instruments in high-resolution spectroscopy, but attempts to improve the wavelength coverage and versatility of these instruments results in a complicated and variable footprint of the entrance slit projection onto the science detector. Traditional spectral extraction methods fail to perform a truly optimal extraction, when the slit image is not aligned with the detector columns but instead is tilted or even curved. We here present the mathematical algorithms and examples of their application to the optimal extraction and the following reduction steps for echelle spectrometers equipped with an entrance slit, that is imaged with various distortions, such as variable tilt and curvature. The new method minimizes the loss of spectral resolution, maximizes the signal-to-noise ratio, and efficiently identifies local outliers. In addition to the new optimal extraction we present order splicing and a more robust continuum normalization algorithms. We have developed and implemented new algorithms that create a continuum-normalized spectrum. In the process we account for the (variable) tilt/curvature of the slit image on the detector and achieve optimal extraction without prior assumptions about the slit illumination. Thus the new method can handle arbitrary image slicers, slit scanning, and other observational techniques aimed at increasing the throughput or dynamic range. We compare our methods with other techniques for different instruments to illustrate superior performance of the new algorithms compared to commonly used procedures.

arXiv comment18 pages, 23 figures, PyReduce code available at https://github.com/AWehrhahn/PyReduce

References in corpus (9)

• Astropy: A Community Python Package for Astronomy
• Binary companions of evolved stars in APOGEE DR14: Search method and catalog of ~5,000 companions
• X-shooter, the new wide band intermediate resolution spectrograph at the ESO Very Large Telescope
• Flat-relative optimal extraction. A quick and efficient algorithm for stabilised spectrographs
• An Extreme Precision Radial Velocity Pipeline: First Radial Velocities from EXPRES
• Precision and consistency of astrocombs
• New wavelength calibration of the HARPS spectrograph
• A Full Implementation of Spectro-Perfectionism for Precise Radial Velocity Exoplanet Detection: A Test Case With the MINERVA Reduction Pipeline
• Robust Reconstruction from Chopped and Nodded Images

Cited by in corpus (11)

• APERO: A PipelinE to Reduce Observations -- Demonstration with SPIRou
• CRIRES${}^{+}$ on sky at the ESO Very Large Telescope
• ATOCA: an algorithm to treat order contamination. Application to the NIRISS SOSS mode
• Automated SpectroPhotometric Image REDuction (ASPIRED)
• The ESPRESSO Redshift Drift Experiment I -- High-resolution spectra of the Lyman-$\alpha$ forest of QSO J052915.80-435152.0
• Abundance analysis of benchmark M dwarfs
• Review of detector requirements: some challenges for the present
• A First Look at CRIRES+: Performance Assessment and Exoplanet Spectroscopy
• Stable magnetic fields and changing starspots on Vega -- An ultra-deep decadal survey at Pic du Midi and OHP
• Method for tilt correction of calibration lines in high resolution spectra
• MOMOS: The Multi-Object MKID Optical Spectrometer Simulator and Data Reduction Package