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

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、析取克里金法　　析