Kriging

Command line

The Kriging operation can be directly executed by typing one of the following expressions on the command line of the Main window:

  

OUTMAP

=

MapKrigingSimple(InputPointMap, Georef, SemiVarModel)

OUTMAP

=

MapKrigingSimple(InputPointMap, Georef, SemiVarModel, 1 | 0)

OUTMAP

=

MapKrigingSimple(InputPointMap, Georef, SemiVarModel, 1 | 0, No | Average | Firstval)

OUTMAP

=

MapKrigingSimple(InputPointMap, Georef, SemiVarModel, 1 | 0, Average | Firstval , Tolerance)

OUTMAP

=

MapKrigingSimple(InputPointMap, Georef, SemiVarModel, , No | Average | Firstval)

OUTMAP

=

MapKrigingSimple(InputPointMap, Georef, SemiVarModel, , Average | Firstval , Tolerance)

OUTMAP

=

MapKrigingSimple(InputPointMap, Georef, SemiVarModel, plane | sphere)

OUTMAP

=

MapKrigingSimple(InputPointMap, Georef, SemiVarModel, plane | sphere, 1 | 0)

OUTMAP

=

MapKrigingSimple(InputPointMap, Georef, SemiVarModel, plane | sphere, 1 | 0, No | Average | Firstval)

OUTMAP

=

MapKrigingSimple(InputPointMap, Georef, SemiVarModel, plane | sphere, 1 | 0, Average | Firstval , Tolerance)

OUTMAP

=

MapKrigingSimple(InputPointMap, Georef, SemiVarModel, plane | sphere, , No | Average | Firstval)

OUTMAP

=

MapKrigingSimple(InputPointMap, Georef, SemiVarModel, plane | sphere, , Average | Firstval , Tolerance)

  

OUTMAP

=

MapKrigingOrdinary(InputPointMap, Georef, SemiVarModel, LimDist)

OUTMAP

=

MapKrigingOrdinary(InputPointMap, Georef, SemiVarModel, LimDist, 1 | 0)

OUTMAP

=

MapKrigingOrdinary(InputPointMap, Georef, SemiVarModel, LimDist, 1 | 0, min, max)

OUTMAP

=

MapKrigingOrdinary(InputPointMap, Georef, SemiVarModel, LimDist, 1 | 0, min, max, No | Average | Firstval )

OUTMAP

=

MapKrigingOrdinary(InputPointMap, Georef, SemiVarModel, LimDist, 1 | 0, min, max, Average | Firstval, Tolerance)

OUTMAP

=

MapKrigingOrdinary(InputPointMap, Georef, SemiVarModel, plane | sphere, LimDist)

OUTMAP

=

MapKrigingOrdinary(InputPointMap, Georef, SemiVarModel, plane | sphere, LimDist, 1 | 0)

OUTMAP

=

MapKrigingOrdinary(InputPointMap, Georef, SemiVarModel, plane | sphere, LimDist, 1 | 0, min, max)

OUTMAP

=

MapKrigingOrdinary(InputPointMap, Georef, SemiVarModel, plane | sphere, LimDist, 1 | 0, min, max, No | Average | Firstval )

OUTMAP

=

MapKrigingOrdinary(InputPointMap, Georef, SemiVarModel, plane | sphere, LimDist, 1 | 0, min, max, Average | Firstval, Tolerance)

where:

OUTMAP

is the name of the output raster map containing the Kriging estimates.

MapKrigingSimple

is the command to start the Kriging operation, using the Simple Kriging method.

MapKrigingOrdinary

is the command to start the Kriging operation, using the Ordinary Kriging method.

InputPointMap

is the name of the input point map with a value domain. To use a value attribute column from an attribute table which is linked to an ID point map, use InputPointMap.Columnname.

Georef

is the name of an existing georeference for the output raster map.

SemiVarModel

Model(nugget, sill, range) | Power(nugget, slope, pow)

This expression defines the semi-variogram model that should be used and the expected parameters.

Model

Spherical | Exponential | Gaussian | Wave | RatQuad | Circular

nugget

value for the nugget, according to your semi-variogram (real value).

sill

value for the sill, according to your semi-variogram (real value).

range

value for the range, according to your semi-variogram (real value > 0).

slope

when using the Power model: value for the 'slope'. When pow is specified as 1 (i.e. thus using a linear model), then slope is the direction coefficient, i.e. Dg/Dh.

pow

when using the Power model: an exponent 0 < real value < 2. When you use value 1, the Power model will become linear and the slope will be constant. If the power exponent is 2 the assumed stochastic model (�randomness�) is not always justifiable and the interpolation can become pathological.

plane | sphere

optional parameter to calculate distances in a plane (plane), or on the sphere by using the projection of the coordinate system as used by the georeference of the output raster map (sphere). If not specified, distances will be calculated in a plane.

LimDist

For Ordinary Kriging: a value for the limiting distance: points that are farther away from an output pixel than the limiting distance will not be used in the Kriging equations.

1 | 0

Specify whether or not an error map has to be created. The error map will contain the square root of the Kriging error variance values, i.e. standard deviations per pixel. When this parameter is not specified, no error map will be created.

1

Create an error map. The error map will obtain the same name as specified for the output Kriging map followed by the additional string _Error.

0

Do not create an error map.

min, max

For Ordinary Kriging: Optional parameters to specify the minimum and maximum number of points that should be taken into account during the calculation of each Kriging estimate/prediction, i.e. the minimum and maximum number of points within the limiting distance of an output pixel that should be taken into account in the calculation of the output value for that pixel. 1 <= min. nr. of points (integer) <= max. nr. of points (integer).

  • When the minimum is not specified, value 1 will be used. When, for an output pixel, less points than the specified minimum are found within the specified limiting distance, no Kriging is performed: the output pixel will be assigned the undefined value.
  • When the maximum is not specified, value 16 will be used. When, for an output pixel, more points than the specified maximum are found within the specified limiting distance, only the specified maximum number of points that are nearest to the output pixel will be used.

Tip: To speed up the algorithm when using large numbers of points within each limiting distance, it is advised to first rasterize the point values and then to use the Kriging from Raster operation.

No | Average | Firstval

 

Choose how to handle possible coinciding points or duplicates. When no method is specified, Average will be used.

No

No removal of duplicates. Mind: when there are coinciding points, the Kriging system may become unstable or unsolvable.

Average

Values of duplicates are replaced by the average (arithmetic mean); the coordinates of the first point are taken as new position.

Firstval

Values of duplicates (and their coordinates) are replaced by those of the first point.

Tolerance

Optional parameter to specify the distance within which 2 points are considered to coincide. When not specified, a tolerance of 0.1 m will be used.

When the definition symbol = is used, a dependent output map is created; when the assignment symbol := is used, the dependency link is immediately broken after the output map has been calculated.

Example:

To perform Ordinary Kriging on the attribute values in column MyAttrib of an attribute table linked to point map MyPoints,

 

use the following expression:

 

OUTMAP

=

MapKrigingOrdinary(MyPoints.MyAttrib, MyGeoref, Spherical(10,70,800), 200, 1, 5, 35, Firstval, 0.10)

When you wish to use spherical distances, the expression will read:

 

OUTMAP

=

MapKrigingOrdinary(MyPoints.MyAttrib, MyGeoref, Spherical(10,70,800), sphere, 200, 1, 5, 35, Firstval, 0.10)

See also: