Kriging in Geography

| |

Spatial kriging belongs to a family of spatial interpolation methods often used to estimate raster values. Typically, points such as elevation or other values not typically represented in a continuous surface can be estimated using known data from surrounding points. 

What is different about kriging, relative to other methods such as inverse distance weighting (IDW), is the method uses a statistical model rather a deterministic weighted estimation. Therefore, kriging provides both an estimate but also a measure that can be used to evaluate the certainty for the predicted value. 

Spatial kriging does have a lot of similarity to IDW, in that both methods using surrounding distances to points to derive a weighted estimate. One main difference, however, is kriging’s weighted values are based on both distance and spatial arrangement of points that surround the area of interest. 

Steps in Applying Kriging

There are two main steps in applying kriging, mainly creating variograms that have covariance functions that formulate a spatial autocorrelation and then determine a predicted value. For each pair of locations, that is the point of interest and all other points, the distance difference is squared in the variograph.

The user typically also selects a function to fit the surface model to, such as a spherical surface or among from many other possibilities, to allow for the estimated weighted distances to be fit to. This makes spatial kriging typically utilize a form of Gaussian surface regression models as part of its estimated surface values.

The end product of kriging is a prediction surface that covers a geographic area based on the value of sampled locations with in that study area.

Two maps showing chloride water sample levels.
On the left is a map of sample chloride levels in water samples. On the right is a surface prediction map created using kriging. Source: USGS.

Additionally, as spatial models often have some distance where their efficacy could be limited, users should be aware where, that is what value ranges, given models could have limitations based on an efficacy.

With Ordinary Kriging, a standard error surface map is also created that spatially visualizes the uncertainty related to the predicted values.

A map with mostly yellow and shades fo orange to show the uncertainty in a predicted value for each area of the map.
A surface error map showing the uncertainty of the surface prediction map for water chloride levels. Map: USGS.

Ordinary and Universal Kriging

There are two main methods used in kriging, called ordinary and universal kriging. Ordinary means you assume the constant mean is unknown, and is generally commonly used by users, while in universal methods a polynomial function could be used for kriging values that may have a given bias, such as direction and regions of data concentration, in order to estimate points.[1]

Using Kriging to Interpolate Missing Geospatial Data

As for applications of these methods, many problems involving incomplete raster data have typically applied a form of spatial kriging, where the goal involves estimating what missing data might be.

For instance, in a relatively recent research paper, scientists used kriging to estimate traffic pollution in given areas. This can be used by planners who are attempting to estimate pollution levels on streets based on varying traffic patterns and using existing known data collection points.[2] 

Spatial kriging has also been used to estimate house prices based on a surface consisting of known house prices, where results are not only available for planning and real estate purposes, but they can also be used to estimate if, for some reason, a given house price is either well-above or well-below an expected trend line.

This can help both consumers and others looking at a series of discrete housing price data where data gaps mean estimates for surrounding regions are often needed.[3]

Some scholars, while finding kriging useful, have also begun to compare spatial kriging methods to different machine learning algorithms, in particular deep learning neural networks and random forest models. In these cases, learning-based algorithms have shown relatively greater accuracy levels than standard kriging methods, particularly if training data are robust or sufficiently developed for given problems.[4] However, spatial kriging could still be of greater benefit where given data are difficult to obtain for training or there might be a limited understanding of a given phenomenon.

Overall, new applications of spatial kriging are always appearing. Spatial kriging has demonstrated to be a strong statistical method that allows for reasonable estimates, based on given fit models, for undetermined, typically raster values.

New methods have appeared in recent years to replace spatial kriging, such as deep learning techniques; however, spatial kriging is still considered a well suited method for many applications. 

Read next: Ten Reasons Why Spatial Data Quality Matters More Now Than Ever


[1]    For more on spatial kriging and how it can be used in estimating values, see: Montero, J. M., Fernández-Avilés, G., & Mateu, J. (2015). Spatial and spatio-temporal geostatistical modeling and kriging. 

[2]    For more on traffic pollution estimates using spatial kriging, see:  Pinto, J. A., Kumar, P., Alonso, M. F., Andreão, W. L., Pedruzzi, R., Espinosa, S. I., & de Almeida Albuquerque, T. T. (2020). Kriging method application and traffic behavior profiles from local radar network database: A proposal to support traffic solutions and air pollution control strategies. Sustainable Cities and Society56, 102062.

[3]    For more on spatial kriging and how it can be used to estimate house price data, see:  Chica-Olmo, J., & Cano-Guervos, R. (2020). Does my house have a premium or discount in relation to my neighbors? A regression-kriging approach. Socio-Economic Planning Sciences72, 100914.

[4]    For more on using machine learning techniques relative to kriging-based approaches, see:  Derdouri, A., & Murayama, Y. (2020). A comparative study of land price estimation and mapping using regression kriging and machine learning algorithms across Fukushima prefecture, Japan. Journal of Geographical Sciences30(5), 794–822.