Search

We are interested in creating statistical methods to provide informative summaries of random fields through the geometry of their excursion sets.To this end, we introduce an estimator for the length of the perimeter of excursion sets of random fields on $\mathbb R^2$ observed over regular square tilings. The proposed estimator acts on the empirically accessible binary digital images of the excursion regions and computes the length of a piecewise linear approximation of the excursion boundary. The estimator is shown to be consistent as the pixel size decreases, without the need of any normalization constant, and with neither assumption of Gaussianity nor isotropy imposed on the underlying random field. In this general framework, even when the domain grows to cover $\mathbb R^2$, the estimation error is shown to be of smaller order than the side length of the domain. For affine, strongly mixing random fields, this translates to a multivariate Central Limit Theorem for our estimator when multiple levels are considered simultaneously. Finally, we conduct several numerical studies to investigate statistical properties of the proposed estimator in the finite-sample data setting.

The excursion set of a $C^2$ smooth random field carries relevant information in its various geometric measures. From a computational viewpoint, one never has access to the continuous observation of the excursion set, but rather to observations at discrete points in space. It has been reported that for specific regular lattices of points in dimensions 2 and 3, the usual estimate of the surface area of the excursions remains biased even when the lattice becomes dense in the domain of observation. In the present work, under the key assumptions of stationarity and isotropy, we demonstrate that this limiting bias is invariant to the locations of the observation points. Indeed, we identify an explicit formula for the bias, showing that it only depends on the spatial dimension $d$. This enables us to define an unbiased estimator for the surface area of excursion sets that are approximated by general tessellations of polytopes in $\mathbb{R}^d$, including Poisson-Voronoi tessellations. We also establish a joint central limit theorem for the surface area and volume estimates of excursion sets observed over hypercubic lattices., 38 pages, 6 figures