克里金法的工作原理
Resource Center ? 专业库 ? 地理处理 ? 地理处理工具参考 ? Spatial Analyst 工具箱 ? 插值工具集 ? 插值工具集概念
克里金法是通过一组具有 z 值的分散点生成估计表面的高级地统计过程。与插值工具集中的其他插值方法不同，选择用于生成输出表面的最佳估算方法之前，有效使用克里金法工具涉及 z 值表示的现象的空间行为的交互研究。
什么是克里金法？
IDW（反距离加权法）和样条函数法插值工具被称为确定性插值方法，因为这些方法直接基于周围的测量值或确定生成表面的平滑度的指定数学公式。第二类插值方法由地统计方法（如克里金法）组成，该方法基于包含自相关（即，测量点之间的统计关系）的统计模型。因此，地统计方法不仅具有产生预测表面的功能，而且能够对预测的确定性或准确性提供某种度量。
克里金法假定采样点之间的距离或方向可以反映可用于说明表面变化的空间相关性。克里金法工具可将数学函数与指定数量的点或指定半径内的所有点进行拟合以确定每个位置的输出值。克里金法是一个多步过程；它包括数据的探索性统计分析、变异函数建模和创建表面，还包括研究方差表面。当您了解数据中存在空间相关距离或方向偏差后，便会认为克里金法是最适合的方法。该方法通常用在土壤科学和地质中。
克里金法公式
由于克里金法可对周围的测量值进行加权以得出未测量位置的预测，因此它与反距离权重法类似。这两种插值器的常用公式均由数据的加权总和组成：
?
其中：
Z(si) = 第 i 个位置处的测量值 λi = 第 i 个位置处的测量值的未知权重
s0 = 预测位置 N = 测量值数
在反距离权重法中，权重λi仅取决于预测位置的距离。但是，使用克里金方法时，权重不仅取决于测量点之间的距离、预测位置，还取决于基于测量点的整体空间排列。要在权重中使用空间排列，必须量化空间自相关。因此，在普通克里金法中，权重 λi 取决于测量点、预测位置的距离和预测位置周围的测量值之间空间关系的拟合模型。以下部分将讨论如何使用常用克里金法公式创建预测表面地图和预测准确性地图。
使用克里金法创建预测表面地图
要使用克里金法插值方法进行预测，有两个任务是必需的： 找到依存规则。 ? 进行预测。
?
要实现这两个任务，克里金法需要经历一个两步过程：
1. 创建变异函数和协方差函数以估算取决于自相关模型（拟合模型）的统计相关性（称为空间自相关）值。 2. 预测未知值（进行预测）。
由于这两个任务是不同的，因此可以确定克里金法使用了两次数据：第一次是估算数据的空间自相关，第二次是进行预测。
变异分析
拟合模型或空间建模也称为结构分析或变异分析。在测量点结构的空间建模中，以经验半变异函数的图形开始，针对以距离 h 分隔的所有位置对，通过以下方程进行计算：
Semivariogram(distanceh) = 0.5 * average{(valuei – valuej}2] 该公式涉及到计算配对位置的差值平方。
下图显示了某个点（红色点）与所有其他测量位置的配对情况。会对每个测量点执行该过程。
计算配对位置的差值平方
通常，各位置对的距离都是唯一的，并且存在许多点对。快速绘制所有配对则变得难以处理。并不绘制每个配对，而是将配对分组为各个步长条柱单元。例如，计算距离大于 40 米但小于 50 米的所有点对的平均半方差。经验半变异函数是 y 轴上表示平均半变异函数值，x 轴上表示距离或步长的图（请参阅下图）。
经验半变异函数图示例
空间自相关量化时采用以下地理的基本原则：距离较近的事物要比距离较远的事物更相似。因此，位置对的距离越近（在半变异函数云的 x 轴上最左侧），具有的值就应该越相似（在半变异函数云的 y 轴上较低处）。位置对的距离变得越远（在半变异函数云的 x 轴上向右移动），就应该变得越不同，差值的平方就会更高（在半变异函数云的 y 轴上向上移动）。
根据经验半变异函数拟合模型
下一步是根据组成经验半变异函数的点拟合模型。半变异函数建模是空间描述和空间预测之间的关键步骤。克里金法的主要应用是预测未采样位置处的属性值。经验半变异函数可提供有关数据集的空间自相关的信息。但是，不提供所有可能的方向和距离的信息。因此，为确保克里金法预测的克里金法方差为正值，根据经验半变异函数拟合模型（即，连续函数或曲线）是很有必要的。该操作理论上类似于回归分析，在此回归分析中将根据数据点拟合连续线或曲线。
要根据经验半变异函数拟合模型，则选择用作模型的函数（例如，开始时上升并在距离变大而超过某一范围后呈现水平状态的球面类型）（请参阅下面的球面模型示例）。经验半变异函数上的点与模型有一些偏差；一些点在模型曲线上方，一些点在模型曲线下方。但是，如果添加一个相应的距离，每个点都会在线上方，或者如果添加另一个相应的距离，每个点都会在线下方，这两个距离值应该是相似的。有多种半变异函数模型可供选择。
半变异函数模型
克里金法工具提供了以下函数，可以从中选择用于经验半变异函数建模的函数：
? ? ? ? ?
圆
球面 指数 高斯 线性
所选模型会影响未知值的预测，尤其是当接近原点的曲线形状明显不同时。接近原点处的曲线越陡，最接近的相邻元素对预测的影响就越大。这样，输出曲面将更不平滑。每个模型都用于更准确地拟合不同种类的现象。 下图显示了两个常用模型并确定了函数的不同之处：
球面模型示例
该模型显示了空间自相关逐渐减小（等同于半方差的增加）到超出某个距离后自相关为零的过程。球面模型是最常用的模型之一。
球面模型示例
指数模型示例
该模型在空间自相关随距离的增加呈指数减小时应用。在这里，自相关仅会在无穷远处完全消失。指数模型也是常用模型。要选择使用哪个模型基于数据的空间自相关和数据现象的先验知识。
指数模型示例
有关更多数学模型的信息，请参见下面。
了解半变异函数 - 变程、基台和块金
正如前文所述，半变异函数显示了测量样本点的空间自相关。由于地理的基本原则（距离越近的事物就越相似），通常，接近的测量点的差值平方比距离很远的测量点的差值平方小。各位置对经调整后进行绘制，然后模型根据这些位置进行拟合。通常使用变程、基台和块金描述这些模型。
变程和基台