Nonparametric spectral density estimation from irregularly sampled data

We introduce a nonparametric spectral density estimator for continuous-time and continuous-space processes measured at fully irregular locations. Our estimator is constructed using a weighted nonuniform Fourier sum whose weights yield a high-accuracy quadrature rule with respect to a user-specified window function. The resulting estimator significantly reduces the aliasing seen in periodogram approaches and least squares spectral analysis, sidesteps the dangers of ill-conditioning of the nonuniform Fourier inverse problem, and can be adapted to a wide variety of irregular sampling settings. After a discussion of methods for computing the necessary weights and a theoretical analysis of sources of bias, we close with demonstrations of the method's efficacy, including for processes that exhibit very slow spectral decay and for processes in multiple dimensions.