Isotropic covariance functions for disparate spatial resolutions
Lucas da Cunha Godoy, Marcos Oliveira Prates, Fernando Andrés Quintana, Jun Yan
TL;DR
This work introduces the ball-Hausdorff distance, a distance between bounded spatial sets defined via the Hausdorff distance between their minimum enclosing balls. Under length-space assumptions, it has a closed-form expression in terms of Chebyshev centers and radii, and if the base metric is conditionally negative definite, bh is CND, enabling isometric Hilbert-space embeddings and a broad class of valid isotropic covariances for set-indexed random fields. The authors show that common covariance families, including Matérn and powered exponential, remain valid after appropriate transformations, while achieving substantial computational gains by reducing set relationships to centers and radii. Numerical results demonstrate improved covariance validity and drastic speedups over standard Hausdorff-based approaches, especially on non-Euclidean domains like the sphere. The framework supports mixed-resolution spatial data problems (change of support, data fusion, misalignment) and points to future avenues in cross-covariance modeling and differentiability analysis on manifolds.
Abstract
Distances between sets arise naturally in data fusion problems involving both point-referenced and areal observations, as well as in set-indexed stochastic processes more broadly. However, commonly used constructions of distances on sets, including those derived from the Hausdorff distance, generally fail to be conditionally negative definite, precluding their use in isotropic covariance models. We propose the ball-Hausdorff distance, defined as the Hausdorff distance between the minimum enclosing balls of bounded sets in a metric space. For length spaces, we derive an explicit representation of this distance in terms of the associated centers and radii. We show that the ball-Hausdorff distance is conditionally negative definite whenever the underlying metric is conditionally negative definite, which implies, via Schoenberg's theorem, an isometric embedding into a Hilbert space. As a consequence, broad classes of isotropic covariance functions, including the Matérn and powered exponential families, are valid for random fields indexed by sets. The resulting construction reduces set-to-set dependence to low-dimensional geometric summaries, leading to substantial computational simplifications in covariance evaluation.
