a list, Vector, listlist, Array, or Matrix containing the coordinates of the points corresponding to data values. If this is a listlist, each inner list must contain n entries, where n is the dimensionality of the data. If this is an Array or Matrix, then there must be n columns, with each row corresponding to one set of coordinates. If it is a plain list or a Vector, then the input points are 1-dimensional.
a list, Array, or Matrix containing the data, which must be arranged in the same order as the corresponding coordinates appear in points
The Kriging command returns a Kriging object. See this help page for a general mathematical description of the Kriging process.
Input sample points must not contain duplicates. The presence of duplicate points can lead to unexpected results.
The following help pages describe the Kriging object and its methods further:
apply Kriging at a grid of values
display the empirical variogram
fit parameters of a variogram
Interpolate at a point
Set the variogram for a Kriging object
Specify some points and some data manually.
points ≔ 0,0,1,0,3,0,1,1,2,1,3,2,0,3,2,3
data ≔ 7.9,7.7,11.4,2.2,3.0,5.7,1.4,4.6
We can visualize these points and data values as follows.
ptp ≔ plotspointplot3d⁡seq⁡op⁡pointsi,datai,i=1..nops⁡points:
Create a Kriging object.
k ≔ Kriging⁡points,data
k≔Krⅈgⅈng ⅈntⅇrpolatⅈon obȷⅇct wⅈth 8 samplⅇ poⅈntsVarⅈogram: Sphⅇrⅈcal(4.35,32.49,2.236067977)
This uses an estimate to set the variogram used. This is often useful if we do not have a model for the variogram, but if we do, we can set the variogram manually.
Krⅈgⅈng ⅈntⅇrpolatⅈon obȷⅇct wⅈth 8 samplⅇ poⅈntsVarⅈogram: Sphⅇrⅈcal(1,40,4)
If we evaluate k at one of the input points, we get the corresponding value back.
The value at other points is interpolated.
We can visualize the interpolated surface as follows; or we can include the data points, too.
pk ≔ plot3d⁡k⁡x,y,x=0..3,y=0..3:
We can also find out what the modeled variance is at various points.
We can also display this variance, or use the variance to color the visualization of the interpolated surface.
pk2 ≔ plot3d⁡k⁡x,y,x=0..3,y=0..3,color=k⁡x,y,output=variance25,0.8,0.8,colortype=HSV:
Maple also contains functionality for generating data that is spatially correlated according to a given variogram.
points,data ≔ KrigingGenerateSpatialData⁡Spherical⁡1,10,1
Create a Kriging object:
k≔Krⅈgⅈng ⅈntⅇrpolatⅈon obȷⅇct wⅈth 30 samplⅇ poⅈntsVarⅈogram: Sphⅇrⅈcal(1.25259453854486,13.6487615617241,.5525536774)
Use the Kriging object to interpolate at a given point:
The Interpolation[Kriging]/Constructor command was introduced in Maple 2018.
For more information on Maple 2018 changes, see Updates in Maple 2018.
Interpolating at a point
Download Help Document