Interpolation[Kriging]
Constructor
Calling Sequence
Parameters
Description
Examples
Compatibility
Kriging(points, data)
points
-
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.
data
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
Overview
Set the variogram for a Kriging object
with⁡Interpolation:
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:
ptp
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.
SetVariogram⁡k,Spherical⁡1,40,4
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.
k⁡3,0
11.4000000000000004
The value at other points is interpolated.
k⁡2,0
7.96736181702214985
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:
pk
plots:-display⁡ptp,pk
We can also find out what the modeled variance is at various points.
k⁡3,0,output=variance
0.
k⁡2,0,output=variance
14.8894610336107736
We can also display this variance, or use the variance to color the visualization of the interpolated surface.
plot3d⁡k⁡x,y,output=variance,x=0..3,y=0..3
pk2≔plot3d⁡k⁡x,y,x=0..3,y=0..3,color=k⁡x,y,output=variance25,0.8,0.8,colortype=HSV:
plots:-display⁡ptp,pk2
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:
k⁡0.2,0.3
−2.75173577049669937
The Interpolation[Kriging]/Constructor command was introduced in Maple 2018.
For more information on Maple 2018 changes, see Updates in Maple 2018.
See Also
Interpolating at a point
Kriging
SetVariogram
Download Help Document