改进的粒子群优化克里金模型，普通的克里金模型MATLAB数据包与遗传算法的结合

克里金法时一种用于空间插值的地学统计方法。 克里金法用半变异测定空间要素，要素即自相关要素。 半变异公式为： 其中γ(h) 是已知点 xi 和 xj 的半变异，***h***表示这两个点之间的距离，z是属性值。 假设不存在漂移，普通克里金法重点考虑空间相关因素，并用拟合的半变异直接进行插值。 估算某测量点z值的通用方程为： 式中，z0是待估计值，zx是已知点x的值，Wx是每个已知点关联的权重，s是用于估计的已知点数目。 权重可以由一组矩阵方程得到。 此程序对半变异进行拟合时采用的时最简单的正比例函数拟合 数据为csv格式 保存格式如下： 第一行为第一个点以此类推

这个函数完成了我从函数 variogram 和 variogramfit 开始的。 它不是用于执行克里金法的高度优化的代码，但可能具有教育价值。 请注意，要使用克里金法，您应该下载最新版本的 variogramfit。 克里金法使用普通克里金法来插入在坐标 x 和 y 位于未采样位置 xi, yi 的位置处测量的变量 z。 该函数需要包含有关变异函数的所有必要信息的变量 vstruct。 vstruct 是函数 variogramfit 的第四个输出参数。 这是一个基本但易于使用的函数来执行普通的克里金插值。 我称之为基本的，因为它总是包含所有观察值来估计未采样位置的值。 当样本位置不在自相关范围内但需要像 ak 最近邻搜索算法或类似的东西时，这可能没有必要。 因此，这些算法最适用于相对较少数量的观察（100-500）。 对于更大量的观察，我建议使用 GSTAT。

将三维开源几何内核Open CASCADE和基于地质统计学的普通克里金(Ordinary Kriging)插值算法从底层相结合，设计了在CAD几何操作核心之上加入地质统计学插值算法的三维水文地质建模软件开发架构。利用Open CASCADE的三维图形渲染、可视化交互、编辑等功能和普通克里金的地质统计学插值功能，以Visual Studio为开发工具，以C++和Python为开发语言，SQLite作为地质数据库设计并开发了Hydrogeo3D矿井水文地质建模软件，从而实现了建模过程局部细节可编辑的功能。

二、普通克里金(OK) (Ordinary Kriging)

该程序执行在excel文件中获得的一组点的2D克里金法。

主要介绍了python 普通克里金（Kriging）法的实现，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧

在网上找了非常多资源没见过教怎么用MATLAB编程实现普通克里金法的，作为初学者感觉上手有一定难度。这个PDF详细介绍了普通克里金法的推导过程，并且对如何用MATLAB实现普通克里金法进行了详细讲解，如要自己改写克里金法可以进行参考学习。 该方法是对普通克里金法进行直接的实现的，如要学习更有扩展性的编程方法参考Engineering Design via Surrogate Modelling,A Practice Guide这本书。

普通克里金插值的详细步骤，主要是通过自己学习普通克里金插值对其计算过程的一个描述，对于刚接触克里金插值的人来说，是很好的一个了解资料，可以快速掌握，并根据资料编程实现，资料中省去了克里金插值中较为繁琐的公式推导。若要看公式的详细推导，文档中也给出了可参考的文献以及部分伪代码。