Tutorial 5: Kriging¶
gstools.krige provides routines for Gaussian process regression, also known as kriging.
Kriging is a method of data interpolation based on predefined covariance models.
The aim of kriging is to derive the value of a field at some point , when there are fixed observed values at given points .
The resluting value at is calculated as a weighted mean:
The weights depent on the given covariance model and the location of the target point.
The different kriging approaches provide different ways of calculating .
The routines for kriging are almost identical to the routines for spatial random fields. First you define a covariance model, as described in Tutorial 2: The Covariance Model, then you initialize the kriging class with this model:
import gstools as gs # condtions cond_pos = [...] cond_val = [...] model = gs.Gaussian(dim=1, var=0.5, len_scale=2) krig = gs.krige.Simple(model, cond_pos=cond_pos, cond_val=cond_val, mean=1)
The resulting field instance
krig has the same methods as the
You can call it to evaluate the kriged field at different points,
you can plot the latest field or you can export the field and so on.
Provided Kriging Methods¶
The following kriging methods are provided within the
||External drift kriging (EDK).|
||Detrended simple kriging.|