Sampling and mapping of heterogeneous surfaces: multi-resolution tiling adjusted to spatial variability

Abstract

Mapping by sampling and prediction of local and regional values of two-dimensional surfaces is a frequent, complex task in geographical information systems. This article describes a method for the approximation of two-dimensional surfaces by optimizing sample size, arrangement and prediction accuracy simultaneously. First, a grid of an ancillary data set is approximated by a quadtree to determine a predefined number of homogeneous mapping units. This approximation is optimal in the sense of minimizing Kullback-divergence between the quadtree and the grid of ancillary data. Then, samples are taken from each mapping unit. The performance of this sampling has been tested against other sampling strategies (regular and random) and found to be superior in reconstructing the grid using three interpolation techniques (inverse squared Euclidean distance, kriging, and Thiessen-polygonization). Finally, the discrepancy between the ancillary grid and the surface to be mapped is modelled by different levels and spatial structures of noise. Conceptually this method is advantageous in cases when sampling strata cannot be well defined a priori and the spatial structure of the phenomenon to be mapped is not known, but ancillary information (e.g., remotely-sensed data), corresponding to its spatial pattern, is available.