Area-to-point Kriging with inequality-type data

In practical applications of area-to-point spatial interpolation, inequality constraints, such as non-negativity or more general constraints on the maximum and/or minimum attribute value, should be taken into account. The geostatistical framework proposed in this paper deals with the spatial interpolation problem of downscaling areal data under such constraints, while: (1) explicitly accounting for support differences between sample data and unknown values, (2) guaranteeing coherent (mass-preserving) predictions, and (3) providing a measure of reliability (uncertainty) for the resulting predictions. The formal equivalence between Kriging and spline interpolation allows solving constrained area-to-point interpolation problems via quadratic programming (QP) algorithms, after accounting for the support differences between various constraints involved in the problem formulation. In addition, if inequality constraints are enforced on the entire set of points discretizing the study domain, the numerical algorithms for QP problems are applied only to selected locations where the corresponding predictions violate such constraints. The application of the proposed method of area-to-point spatial interpolation with inequality constraints in one and two dimension is demonstrated using realistically simulated data.