vebyk

vebyk performs ordinary kriging and can be easily adapted to other kriging methods.

The program is designed to interpolate values on a regular two-dimensional grid using ordinary kriging. By slightly modifying the code it is also possible to use simple kriging. The grid of estimated or "kriged" values is rectangular and spans the range of coordinates in the dataset of sampled values. The spacing between estimated points for x- and y-axes can be specified individually and the sampled values do not have to follow a spatial order.
For the search neighborhood, the number of points used for interpolation can be specified. Arbitrarily oriented spatial anisotropy of the covariance function can be accounted for. The current implementation is based on the von Karman family of covariance function. The correlation length a and the exponent nu of the von Karman function can be specified.

Finally, "handles" are provided to switch on a wait bar or a cross validation mode for comparison of the interpolated value at a sampled location.

Please cite my diploma thesis when you publish based on this software (it also contains a description of the algorithm):
Sidler, R., 2003: Kriging and conditional geostatistical simulation based on scale-invariant covariance models. Diploma thesis, ETH Zurich.