Spectral estimation for point processes and random fields

point patterns
Random Fields
spectral analysis
Authors
Affiliations

Jake P Grainger

École Polytechnique Fédérale de Lausanne

Tuomas Rajala

Natural Resources Institute Finland

David J Murrell

University College London

Sofia Olhede

École Polytechnique Fédérale de Lausanne

Published

December 18, 2025

Doi

Abstract

Spatial variables can be observed in many different forms, such as regularly sampled random fields (lattice data), point processes and randomly sampled spatial processes. Joint analysis of such collections of observations is clearly desirable, but complicated by the lack of an easily implementable analysis framework. We fill this gap by providing a multitaper analysis framework using coupled discrete and continuous data tapers, combined with the discrete Fourier transform for inference. Using this set of tools is important, as it forms the backbone for practical spectral analysis. In higher dimensions it is especially important not to be constrained to Cartesian product domains, and so we develop the methodology for spectral analysis using irregular domain data tapers and the tapered discrete Fourier transform. We discuss its fast implementation, as well as the asymptotic and large finite-domain properties. Estimators of partial association between different spatial processes are provided, as are principled methods to determine their significance, and we demonstrate their practical utility using a large-scale ecological dataset.

Fig. 4. The first column corresponds to the three point processes B. tovarensis, P. armata, U. pittieri and the gradient of the terrain. The 4 by 4 panel subplot on the right shows the magnitude coherence (upper triangle) and magnitude partial coherence (lower triangle) between the four processes. The spectral plots range from to 0.05 on each axis, and the data are shown with axes ranging from 0 to 1000 and 0 to 500.