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