Modeling the Echelle Spectra Continuum with Alpha Shapes and Local Regression Fitting

Abstract

Continuum normalization of echelle spectra is an important data analysis step that is difficult to automate. Polynomial fitting requires a reasonably high order model to follow the steep slope of the blaze function. However, in the presence of deep spectral lines, a high order polynomial fit can result in ripples in the normalized continuum that increase errors in spectral analysis. Here, we present two algorithms for flattening the spectrum continuum. The Alpha-shape Fitting to Spectrum algorithm (AFS) is completely data-driven, using an alpha shape to obtain an initial estimate of the blaze function. The Alpha-shape and Lab Source Fitting to Spectrum algorithm (ALSFS) incorporates a continuum constraint from a lab source reference spectrum for the blaze function estimation. These algorithms are tested on a simulated spectrum, where we demonstrate improved normalization compared to polynomial regression for continuum fitting. We show an additional application, using the algorithms for mitigation of spatially correlated quantum efficiency variations and fringing in the CCD detector of the EXtreme PREcision Spectrometer (EXPRES).

References in corpus (4)

• The Gaia FGK Benchmark Stars - High resolution spectral library
• Insights on the Spectral Signatures of Stellar Activity and Planets from PCA
• Impact of micro-telluric lines on precise radial velocities and its correction
• Towards a Self-calibrating, Empirical, Light-Weight Model for Tellurics in High-Resolution Spectra

Cited by in corpus (2)

• RASSINE: Interactive tool for normalising stellar spectra I. Description and performance of the code
• First detection of ozone in the mid-infrared at Mars: implications for methane detection