克里金插值法

克里金插值法

克里金插值法

克里金插值法又称空间局部插值法，是以变异函数理论和结构分析为基础，在有限区域内对区域化变量进行无偏最优估计的一种方法，是地统计学的主要内容之一，由南非矿产工程师 D. Matheron于1951年在寻找金矿时首次提出，法国著名统计学家G. Matheron随后将该方法理论化、系统化，并命名为Kriging，即克里金插值法。

1 克里金插值法原理

克里金插值法的适用范围为区域化变量存在空间相关性，即如果变异函数和结构分析的结果表明区域化变量存在空间相关性，则可以利用克里金插值法进行内插或外推。其实质是利用区域化变量的原始数据和变异函数的结构特点，对未知样点进行线性无偏、最优估计，无偏是指偏差的数学期望为0，最优是指估计值与实际值之差的平方和最小[1]。因此，克里金插值法是根据未知样点有限领域内的若干已知样本点数据，在考虑了样本点的形状、大小和空间方位，与未知样点的相互空间关系，以及变异函数提供的结构

信息之后，对未知样点进行的一种线性无偏最优估计。

假设研究区域a 上研究变量Z （x ），在点x i ∈A （i=1，2，……，n ）处属性值为Z （x i ），则待插点x 0∈A 处的属性值Z （x 0）的克里金插值结果Z*（x 0）是已知采样点属性值Z （x i ）（i=1，2，……，n ）的加权和，即：

)()(10*

i n i i x Z x Z ∑==λ

（1）

式中i

λ是待定权重系数。 其中Z(x i )之间存在一定的相关关系，这种相关性除与距离有关外，还与其相对方向变化有关，克里金插值方法将研究的对象称“区域化变量”

针对克里金方法无偏、最小方差条件可得到无偏条件可得待定权系数i

λ (i=1，2，……，n)满足关系式：

11=∑=n i i λ

（2）

以无偏为前提，kriging 方差为最小可得到求解待定权系数i

λ的方程组：

⎪⎪⎩

⎪⎪⎨⎧=⋯⋯==+∑∑==1)n ,2,1)(,(),(101n i i j j i n i i j x x C x x C λμλ， （3）

式中，C （x i ，x j ）是Z(x i )和Z(x j )的协方差

函数。

2 方法步骤

克里金插值法的应用步骤如下：

1、输入原始数据，即采样点，下面以输入

三个采样点求待估插值为例来进行说明。如图1所示：

图1 采样点图示

2、网格化，选择区域的范围和网格的大小，对区域进行网格化处理。

3、数据检验与分析，根据采样值是否合乎

实际情况，剔除明显差异点。

4、直方图的计算，直方图有助于掌握区域

变化的分布规律，以便决定是否对原始数据进行转换。

5、利用变异函数进行变异函数计算，了解

变量的空间结构。

6、克里金插值估计

（1）待估点权重系数估计

利用多边形估计的方法，首先确定离待估点

最近的采样点的权重，根据公式（4）进行采样点权重估计：

∑=

++=n i w w i i d c d c 111

λ （4）

（2）根据搜索策略选择合适的参估点，如

图2：

图2 参估点图示

（3）根据已经求出的变异函数以及采样点

数量，三个采样点列出三个等式，求出方程组的系数，公式为：

⎥⎥⎥⎦

⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡)3,0()2,0()1,0()3,3()2,3()1,3()3,2()2,2()1,2()3,1()2,1()1,1(321C C C C C C C C C C C C λλλ （5）

（4）分析在各向同性条件下改变块金值与

在块金值相同条件下改变各向异性对权重值的影响[2]。各向同性条件下改变块金值时对权重值的影响效果如图3（a ），在块金值相同条件下改变各向异性对权重值带来的影响如图3（b ）：

（

a ）

（b ） 图3 各向同性条件下改变块金值与在块金值相同条件下改变各向异性对权重值的影响

（5）根据求出的权重值，代入公式（1），即可求得评估领域内n 个采样值的线性组合[2]。 克里金插值法的方法路线图如下：

图4 方法路线图 3 克里金插值法分类及适用类型

克里金插值法主要有以下几种类型：普通克

里金（Ordinary Kriging ）、简单克里金（Simple Kriging ）、泛克里金（Universal Kriging ）、协同克里金（Co-Kriging ）、对数正态克里金（Logistic Normal Kriging ）、指示克里金（Indicator 导入数

数据分

是否

是否

根据数计算样计算样

按距离数据变泛克里按组统计

绘制

绘制拟合计算克是

否

是 否 进行预

Kriging ）、概率克里金（Probability Kriging ）和析取克里金（Disjunctive Kriging ）等[1]。

克里金插值法可以简单地表达为：

)()()(s s s Z εμ+=

(6)

式中，s 为不同位置的点，可以人为是用经纬度表示的空间坐标；Z （s ）为s 处的变量值，它可以分解为确定趋势值)(s μ和自相关随机误差)(s ε。通过对这个公式进行变化，可以生成克里金插值法的不同类型。

首先，对于趋势值)(s μ，可以简单地赋予一个常量，即在任何位置s 处)(s μ=μ，如果μ是未知的，这便是普通克里金基本模型；)(s μ也可表示为空间坐标的线性函数，如：

xy y x y x s 52423210)(ββββββμ+++++=

(7)

如果趋势面方程中的回归系数是未知的，则形成泛克里金模型；如果在任何时候趋势已知的（如所有系数和协方差均已知），无论趋势常量与否，都会形成简单克里金模型。

其次，无论趋势如何复杂，)(s μ仍无法获得很好的预测，在这种情况下需要对误差项)(s ε进