Normalized to: Politsch, C.
[1]
oai:arXiv.org:1908.07151 [pdf] - 2034397
Trend Filtering -- I. A Modern Statistical Tool for Time-Domain
Astronomy and Astronomical Spectroscopy
Submitted: 2019-08-19, last modified: 2020-01-10
The problem of denoising a one-dimensional signal possessing varying degrees
of smoothness is ubiquitous in time-domain astronomy and astronomical
spectroscopy. For example, in the time domain, an astronomical object may
exhibit a smoothly varying intensity that is occasionally interrupted by abrupt
dips or spikes. Likewise, in the spectroscopic setting, a noiseless spectrum
typically contains intervals of relative smoothness mixed with localized higher
frequency components such as emission peaks and absorption lines. In this work,
we present trend filtering, a modern nonparametric statistical tool that yields
significant improvements in this broad problem space of denoising $spatially$
$heterogeneous$ signals. When the underlying signal is spatially heterogeneous,
trend filtering is superior to any statistical estimator that is a linear
combination of the observed data---including kernel smoothers, LOESS, smoothing
splines, Gaussian process regression, and many other popular methods.
Furthermore, the trend filtering estimate can be computed with practical and
scalable efficiency via a specialized convex optimization algorithm, e.g.
handling sample sizes of $n\gtrsim10^7$ within a few minutes. In a companion
paper, we explicitly demonstrate the broad utility of trend filtering to
observational astronomy by carrying out a diverse set of spectroscopic and
time-domain analyses.
[2]
oai:arXiv.org:2001.03552 [pdf] - 2034594
Trend Filtering -- II. Denoising Astronomical Signals with Varying
Degrees of Smoothness
Submitted: 2020-01-10
Trend filtering---first introduced into the astronomical literature in Paper
I of this series---is a state-of-the-art statistical tool for denoising
one-dimensional signals that possess varying degrees of smoothness. In this
work, we demonstrate the broad utility of trend filtering to observational
astronomy by discussing how it can contribute to a variety of spectroscopic and
time-domain studies. The observations we discuss are (1) the Lyman-$\alpha$
forest of quasar spectra; (2) more general spectroscopy of quasars, galaxies,
and stars; (3) stellar light curves with planetary transits; (4) eclipsing
binary light curves; and (5) supernova light curves. We study the
Lyman-$\alpha$ forest in the greatest detail---using trend filtering to map the
large-scale structure of the intergalactic medium along quasar-observer lines
of sight. The remaining studies share broad themes of: (1) estimating
observable parameters of light curves and spectra; and (2) constructing
observational spectral/light-curve templates. We also briefly discuss the
utility of trend filtering as a tool for one-dimensional data reduction and
compression.