《ArcGIS教程：几种克里金法的概述》由会员分享，可在线阅读，更多相关《ArcGIS教程：几种克里金法的概述（6页珍藏版）》请在金锄头文库上搜索。
1、ArcGIS 教程：几种克里金法的概述1、普通克里金法普通克里金法假设模型为Z(s)=卩 + (s)其中，4是一个未知常量。对于普通克里金法，我们所关心的主要问题之一就是对常量平均值的假设 是否合理。有时有很充分的科学依据来拒绝该假设。不过，作为一种简单的预测方法，它具有显著的灵活 性。下图所举的是处于某一空间维度中的示例：从图上看，数据好像是从山谷或山体的线横断面中采集的高程值。而且，好像数据在左侧变化更显著， 而在右侧则变得更平滑。事实上，该数据是在平均值 4 为常量的情况下基于普通克里金法模型模拟得到的 虚线给出的是平均值，该平均值是是真值但是是未知的。因此，普通克里金法可用于似乎带有某
2、种趋势的 数据。单凭数据无法确定已观测到的模式是否是自相关(4 为常量的情况下，在误差 (s) 之间)或趋势 (4(s) 随 s 变化)所造成的。普通克里金法可以使用半变异函数或协方差(用于表达空间自相关的数学形式)，使用变换和移除趋势 还允许测量误差。2、简单克里金法简单克里金法假设模型为：Z(s)=卩 + (s)其中，4是已知常量例如，下图中使用的数据与普通克里金法和泛克里金法概念介绍中所使用的数据相同，观测数据以实 心圆的形式给出：虚线表示的已知常量为 4。这点可以与普通克里金法进行比较。对于简单克里金法，因为假设确切已 知4,那么也确切已知数据位置上的 (s)。对于普通克里金法，如果估
3、算了 4,那么也会估算 (s)。 如果已知 (s),可以比估算 (s)时更好地估算自相关。通常，已知确切平均值4的假设是不现实的。 但是，有时候，假定一个基于物理的模型能够给出已知趋势却是有意义的。由此可以使用模型和观测值的 差值(称为残差)，并且假设残差中的趋势已知为零，可以在残差上使用简单克里金法。简单克里金法可以使用半变异函数或协方差(用于表达自相关的数学形式)和变换，并且允许测量误差。3、指示克里金法指示克里金法假设模型为I(s) = 4 + (s),其中，4是一个未知常量，I(s)是一个二进制变量。二进制数据的创建可利用连续数据的阈值实现， 或者观测数据可以为 0 或 1。例如，假设
4、存在一个由某点是否为森林栖息地或非森林栖息地的相关信息组 成的样本，则其中二进制变量用来指示这两种类别。使用二进制变量时，指示克里金法的处理过程与普通 克里金法相同。在下图中，已使用了解阈值中介绍的阈值将数据转换为二进制值。用空心方块给出了观测的二进制数据。虚线表示所有指示变量的未知平均值，即山可以将这一点与 普通克里金法进行比较。在使用普通克里金法时，假设 (s)是自相关的。请注意，因为指示变量为0或 1，所以插值结果将位于 0 和 1之间，而且基于指示克里金法的预测可解释为变量是 1 的概率或属于 1 所 指示的类别的概率。如果创建指示变量时使用了阈值，则生成的插值地图会显示超出(或低于)
5、阈值的概率。通过选择多个阈值可以为同一数据集创建多个指示变量。在本例中，一个阈值创建主要指示变量，而 另一个指示变量则用作协同克里金法中的次要变量。指示克里金法可使用半变异函数或协方差，它们都是用于表达自相关的数学形式。4、概率克里金法概率克里金法假设模型为I(s) = I(Z(s) ct) = P-1 + l(s) Z(s) = P2 + 2(s),其中M和42为未知常量，I(s)是通过使用阈值指示I(Z(s) ct)创建的二进制变量。请注意， 现在有两种类型的随机误差： 1(s)和 2(s)，因此它们各自存在自相关，并且它们之间存在互相关。 概率克里金法要实现指示克里金法相同的功能很吃力，
6、而使用协同克里金法进行尝试则可更好地实现。例如，在下图中普通克里金法、泛克里金法、简单克里金法和指示克里金法概念使用相同的数据，请 注意标注为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。概率克里金法尝试利用原始数据中除二进制变量之外的其他信息。但是，这也存在一些 代价。必须要进行更多的估算，包括估算每个变量的自相关
7、和互相关。然而，每次估算未知的自相关参数 时，都会引入更多的不确定性，因此概率克里金法可能不值得付出额外努力。概率克里金法可以使用半变异函数或协方差(用于表达自相关的数学形式)、交叉协方差(用于表达互相 关的数学形式)和变换，但是不允许测量误差。5、析取克里金法析取克里金法假设的模型为f(z(s)= m + s (s),其中，M是一个未知常量，f(Z(s)是Z(s)的一个任意函数。请注意，您可以写成f(Z(s) =I(Z(s) ct),这样，指示克里金法就成为析取克里金法的一种特殊情况。在Geostatistical Analyst中使用析取 克里金法时，您既可预测值本身，也可预测指示器。在
8、Geostatistical Analyst 中，提供的 g(Z(s0)函数其实就是 Z(s0)本身和 I(Z(sO) ct)。一 般来说，相比普通克里金法，析取克里金法可以做更多事情。尽管回报更丰厚，但成本也更高。析取克里 金法要求接受二元正态分布假设和fi(Z(si)函数的近似值;但是很难对假设进行验证，而且从数学和计 算角度来看，解决方案都很复杂。析取克里金法可使用半变异函数或协方差(用于表达自相关的数学公式)以及变换，但是它不允许出现 测量误差。6、协同克里金法协同克里金法使用多种变量类型的信息。主要的感兴趣变量是Z1,可利用Z1的自相关性和Z1与所 有其他变量类型间的互相关性进行更好
9、的预测。使用其他变量的信息进行预测固然很好，但也有其弊端。 协同克里金法需要进行更多的估计，包括估计每一变量的自相关性以及所有的互相关性。从理论上讲，克 里金法最保险，因为即使没有互相关性，也可以对Z1的自相关性进行估计。但是，每一次对未知自相关 参数进行估计都会产生更多的变化，从而使得为获得精确性作出的努力得不偿失。普通协同克里金法对模型做了假设Zl(s) = m + l(s) Z2(s)=卩2 + 2(s),其中，M和42为未知常量。请注意，此时存在两种类型的随机误差， 1(s)和 2(s),因此，它 们各自具有自相关性且两者之间存在互相关性。同普通克里金法一样，普通协同克里金法尝试对 Z
10、1(s0) 进 行预测，但是为了使预测更精确，还使用了协变量 Z2(s) 中的信息。例如，下图中的数据与普通克里金法 使用的数据相同，只是在此处额外添加了另外一个变量。请注意，数据Z1和Z2具有自相关性。还应注意，当Z1低于其平均值41时，Z2往往会高于其平 均值42,反之亦然。因此，Z1和Z2具有负的互相关性。在本例中，每一位置s都有对应的Z1(s)和 Z2(s);但这是不必要的，每一变量类型都可以具有其自身独特的一组位置。主要的感兴趣变量是Z1,自相 关性和互相关性都用来进行更好的预测。其他协同克里金法 - 泛协同克里金法、简单协同克里金法、指示协同克里金法、概率协同克里金法和 析取协同克里金法 - 为上述方法的扩展方法，适用于有多个数据集的情况。例如，可通过以下方式实现指 示协同克里金法：为数据使用多个阀值，而后使用每个阀值的二进制数据来预测主要感兴趣的阀值。因此， 指示协同克里金法与概率克里金法类似，但前者对异常值和其他不稳定数据敏感度较低。协同克里金法可以使用半变异函数或协方差（用于表示自相关性的数学形式）、互协方差（用于表示互相 关性的数学形式）和变换以及趋势移除;并在执行普通协同克里金法、简单协同克里金法或泛协同克里金法 时允许存在测量误差。
