The widely applied geostatistical interpolation methods of ordinary kriging (OK) or external drift kriging (EDK) interpolate the variable of interest to the unknown location  providing a linear estimator and an estimation variance as measure of uncertainty