ArcGIS教程：几种克里金法的概述.docx

地理国情监测云平台
ArcGIS教程：几种克里金法的概述
　　1、普通克里金法
　　普通克里金法假设模型为
　　Z(s) = µ + ε(s)
　　其中，µ 是一个未知常量。对于普通克里金法，我们所关心的主要问题之一就是对常量平均值的假设是否合理。有时有很充分的科学依据来拒绝该假设。不过，作为一种简单的预测方法，它具有显著的灵活性。下图所举的是处于某一空间维度中的示例：

　　从图上看，数据好像是从山谷或山体的线横断面中采集的高程值。而且，好像数据在左侧变化更显著，而在右侧则变得更平滑。事实上，该数据是在平均值 µ 为常量的情况下基于普通克里金法模型模拟得到的。虚线给出的是平均值，该平均值是是真值但是是未知的。因此，普通克里金法可用于似乎带有某种趋势的数据。单凭数据无法确定已观测到的模式是否是自相关(µ 为常量的情况下，在误差 ε(s) 之间)或趋势(µ(s) 随 s 变化)所造成的。
普通克里金法可以使用半变异函数或协方差(用于表达空间自相关的数学形式)，使用变换和移除趋势，还允许测量误差。
2、简单克里金法
简单克里金法假设模型为：
　　Z(s) = µ + ε(s)
　　其中，µ 是已知常量
　　例如，下图中使用的数据与普通克里金法和泛克里金法概念介绍中所使用的数据相同，观测数据以实心圆的形式给出：

　　虚线表示的已知常量为 µ。这点可以与普通克里金法进行比较。对于简单克里金法，因为假设确切已知 µ，那么也确切已知数据位置上的 ε(s)。对于普通克里金法，如果估算了 µ，那么也会估算 ε(s)。如果已知 ε(s)，可以比估算 ε(s)时更好地估算自相关。通常，已知确切平均值 µ 的假设是不现实的。但是，有时候，假定一个基于物理的模型能够给出已知趋势却是有意义的。由此可以使用模型和观测值的差值(称为残差)，并且假设残差中的趋势已知为零，可以在残差上使用简单克里金法。
简单克里金法可以使用半变异函数或协方差(用于表达自相关的数学形式)和变换，并且允许测量误差。
3、指示克里金法
　　指示克里金法假设模型为
　　I(s) = µ + ε(s),
　　其中，µ 是一个未知常量，I(s) 是一个二进制变量。二进制数据的创建可利用连续数据的阈值实现，或者观测数据可以为 0 或 1。例如，假设存在一个由某点是否为森林栖息地或非森林栖息地的相关信息组成的样本，则其中二进制变量用来指示这两种类别。使用二进制变量时，指示克里金法的处理过程与普通克里金法相同。
　　在下图中，已使用了解阈值中介绍的阈值将数据转换为二进制值。

　　用空心方块给出了观测的二进制数据。虚线表示所有指示变量的未知平均值，即 µ。可以将这一点与普通克里金法进行比较。在使用普通克里金法时，假设 ε(s) 是自相关的。请注意，因为指示变量为 0 或 1，所以插值结果将位于 0 和 1之间，而且基于指示克里金法的预测可解释为变量是 1 的概率或属于 1 所指示的类别的概率。如果创建指示变量时使用了阈值，则生成的插值地图会显示超出(或低于)阈值的概率。
　　通过选择多个阈值可以为同一数据集创建多个指示变量。在本例中，一个阈值创建主要指示变量，而另一个指示变量则用作协同克里金法中的次要变量。
指示克里金法可使用半变异函数或协方差，它们都是用于表达自相关的数学形式。
4、概率克里金法
　　概率克里金法假设模型为
　　I(s) = I(Z(s) > ct) = µ1 + ε1(s) Z(s) = µ2 + ε2(s),
　　其中 µ1 和 µ2 为未知常量，I(s) 是通过使用阈值指示 I(Z(s) > ct) 创建的二进制变量。请注意，现在有两种类型的随机误差：ε1(s) 和 ε2(s)，因此它们各自存在自相关，并且它们之间存在互相关。概率克里金法要实现指示克里金法相同的功能很吃力，而使用协同克里金法进行尝试则可更好地实现。
　　例如，在下图中普通克里金法、泛克里金法、简单克里金法和指示克里金法概念使用相同的数据，请注意标注为Z(u=9) 的基准的指示变量为 I(u) = 0，标注为 Z(s=10) 的基准的指示变量为 I(s) = 1。

　　如果要预测它们中间的位于 x 坐标 9.5 处的值，单独使用指示克里金法将给出接近 0.5 的预测值。但是，可以看出Z(s) 刚好高于阈值，而 Z(u) 却远低于阈值。因此，有理由相信位置 9.5 处的指示预测值应该小于 0.5。概率克里金法尝试利用原始数据中除二进制变量之外的其他信息。但是，这也存在一些代价。必须要进行更多的估算，包括估算每个变量的自相关和互相关。然而，每次估算未知的自相关参数时，都会引入更多的不确定性，因此概率克里金法可能不值得付出额外努力。
概率克里金法可以使用半变异函数或协方差(用于表达自相关的数学形式)、交叉协方差(用于表达互相关的数学形式)和变换，但是不允许测量误差。
5、析取克里金法
　　析取克里金法假设的模型为
　　f(Z(s)) = µ1 + ε(s),
　　其中，µ1 是一个未知常量，f(Z(s)) 是 Z(s) 的一个任意函数。请注意，您可以写成 f(Z(s)) = I(Z(s) > ct)，这样，指示克里金法就成为析取克里金法的一种特殊情况。在 Geostatistical Analyst 中使用析取克里金法时，您既可预测值本身，也可预测指示器。
　　在 Geostatistical Analyst 中，提供的 g(Z(s0)) 函数其实就是 Z(s0) 本身和 I(Z(s0) > ct)。一般来说，相比普通克里金法，析取克里金法可以做更多事情。尽管回报更丰厚，但成本也更高。析取克里金法要求接受二元正态分布假设和fi(Z(si)) 函数的近似值;但是很难对假设进行验证，而且从数学和计算角度来看，解决方案都很复杂。
析取克里金法可使用半变异函数或协方差(用于表达自相关的数学公式)以及变换，但是它不允许出现测量误差。
6、协同克里金法
　　协同克里金法使用多种变量类型的信息。主要的感兴趣变量是 Z1，可利用 Z1 的自相关性和 Z1 与所有其他变量类型间的互相关性进行更好的预测。使用其他变量的信息进行预测固然很好，但也有其弊端。协同克里金法需要进行更多的估计，包括估计每一变量的自相关性以及所有的互相关性。从理论上讲，克里金法最保险，因为即使没有互相关性，也可以对Z1 的自相关性进行估计。但是，每一次对未知自相关参数进行估计都会产生更多的变化，从而使得为获得精确性作出的努力得不偿失。
　　普通协同克里金法对模型做了假设
　　Z1(s) = µ1 + ε1(s) Z2(s) = µ2 + ε2(s),
　　其中，µ1 和 µ2 为未知常量。请注意，此时存在两种类型的随机误差，ε1(s) 和 ε2(s)，因此，它们各自具有自相关性且两者之间存在互相关性。同普通克里金法一样，普通协同克里金法尝试对 Z1(s0) 进行预测，但是为了使预测更精确，还使用了协变量 Z2(s) 中的信息。例如，下图中的数据与普通克里金法使用的数据相同，只是在此处额外添加了另外一个变量。
　　请注意，数据 Z1 和 Z2 具有自相关性。还应注意，当 Z1 低于其平均值 µ1 时，Z2 往往会高于其平均值 µ2，反之亦然。因此，Z1 和 Z2 具有负的互相关性。在本例中，每一位置 s 都有对应的 Z1(s) 和 Z2(s);但这是不必要的，每一变量类型都可以具有其自身独特的一组位置。主要的感兴趣变量是 Z1，自相关性和互相关性都用来进行更好的预测。

　　其他协同克里金法 - 泛协同克里金法、简单协同克里金法、指示协同克里金法、概率协同克里金法和析取协同克里金法 - 为上述方法的扩展方法，适用于有多个数据集的情况。例如，可通过以下方式实现指示协同克里金法：为数据使用多个阀值，而后使用每个阀值的二进制数据来预测主要感兴趣的阀值。因此，指示协同克里金法与概率克里金法类似，但前者对异常值和其他不稳定数据敏感度较低。
　　协同克里金法可以使用半变异函数或协方差(用于表示自相关性的数学形式)、互协方差(用于表示互相关性的数学形式)和变换以及趋势移除;并在执行普通协同克里金法、简单协同克里金法或泛协同克里金法时允许存在测量误差。
北京数字空间科技有限公司