二维插值的普通克里金法-基于地质统计学的MATLAB实现

简介：
本文介绍了利用MATLAB软件进行二维插值的一种方法——普通克里金法，并探讨了其在地质统计学中的应用与实践。 这个函数是我从variogram 和 variogramfit 函数开始开发的。它并非为执行克里金法的高度优化代码而设计，但可能具有教育价值。请注意，若要使用克里金法，请确保下载最新版本的 variogramfit。该函数采用普通克里金方法，在坐标 x 和 y 上未采样位置 xi, yi 处插入变量 z 的测量值。它需要一个包含所有必要变异函数信息的 vstruct 变量，vstruct 是由 variogramfit 函数输出的第四个参数提供。这是一个基本但易于使用的函数来执行普通克里金插值。我称之为基础是因为它总是使用所有的观察点来估计未采样位置处的值。然而，在样本点不在自相关范围内时（例如需要类似 ak 最近邻搜索算法的情况），这可能没有必要了。因此，这些方法最适合于相对少量的观测数据（大约100-500个）。对于大量观测数据，则建议使用 GSTAT。