局部插值方法只使用邻近的数据点来估计未知点的值，包括几个步骤： 1）定义一个邻域或搜索范围；
2）搜索落在此邻域范围的数据点；
4）为落在规则格网单元上的数据点赋值。重复这个步骤直到格网上的所有点赋值完毕。
使用局部插值方法需要注意的几个方面是：所使用的插值函数；邻域的大小、形状和方向；数据点的个数；数据点的分布方式是规则的还是不规则的。
1）最近邻点法：泰森多边形方法
泰森多边形（Thiessen，又叫Dirichlet 或Voronoi多边形）采用了一种极端的边界内插方法，只用最近的单个点进行区域插值。泰森多边形按数据点位置将区域分割成子区域，每个子区域包含一个数据点，各子区域到其内数据点的距离小于任何到其它数据点的距离，并用其内数据点进行赋值。连接所有数据点的连线形成Delaunay三角形，与不规则三角网TIN具有相同的拓扑结构。
GIS和地理分析中经常采用泰森多边形进行快速的赋值，实际上泰森多边形的一个隐含的假设是任何地点的气象数据均使用距它最近的气象站的数据。而实际上，除非是有足够多的气象站，否则这个假设是不恰当的，因为降水、气压、温度等现象是连续变化的，用泰森多边形插值方法得到的结果图变化只发生在边界上，在边界内都是均质的和无变化的。
2）移动平均插值方法：距离倒数插值
距离倒数插值方法综合了泰森多边形的邻近点方法和趋势面分析的渐变方
法的长处，它假设未知点x0处属性值是在局部邻域内中所有数据点的距离加权平均值。距离倒数插值方法是加权移动平均方法的一种。加权移动平均方法的计算公式如下：
式中，权重系数由函数计算，要求当时，一般取倒数或负指数形式。其中最常见的形式是距离倒数加权函数，形式如下：
式中，x j为未知点，x i为已知数据点。
加权移动平均公式最简单的形式称为线性插值，公式如下：
距离倒数插值方法是GIS软件根据点数据生成栅格图层的最常见方法。距离倒数法计算值易受数据点
集的影响，计算结果经常出现一种孤立点数据明显高于周围数据点的“鸭蛋”分布模式，可以在插值过程中通过动态修改搜索准则进行一定程度的改进。
3）样条函数插值方法
在计算机用于曲线与数据点拟合以前，绘图员是使用一种灵活的曲线规逐段的拟合出平滑的曲线。这种灵活的曲线规绘出的分段曲线称为样条。与样条匹配的那些数据点称为桩点，绘制曲线时桩点控制曲线的位置。曲线规绘出的曲线在数学上用分段的三次多项式函数来描述这种曲线，其连接处有连续的一阶和二阶连续导数。
样条函数是数学上与灵活曲线规对等的一个数学等式，是一个分段函数，进行一次拟合只有与少数点拟合，同时保证曲线段连接处连续。这就意味着样条函
数可以修改少数数据点配准而不必重新计算整条曲线，趋势面分析方法做不到这一点，（图16）。
图8-16：样条函数的局部特征
（a：当二次样条曲线的一个点位置变化时，只需要重新计算四段曲线；
b：而一次样条曲线的一个点位置变化时，只需要重新计算两段曲线）
一般的分段多项式p(x)定义为：
p(x)=p i (x) x i < x < x i+1( i=1, 2, 3…, k-1 )
( j=0, 1, 2, ..., r-1; i=1, 2, ... , k-1 )
x1, ..., x k-1将区间x0，x k分成k个子区间，这些分割点称“断点”，曲线上具有这些x值的点称为“节”。函数p i(x)为小于等于m次的多项式。r项用来表示样条函数的约束条件：
r=0，无约束；
r=1，函数连续且对它的导数无任何约束；
r=m-1，区间[x0x k]可用一个多项式表示；
r=m，约束条件最多。
m=1,2,3时的样条分别为一次、二次、三次样条函数，其导数分别是0阶、1阶、2阶导数，二次样条函数的每个节点处必须有一阶连续导数，三次样条函数的每个节点初必须有二阶连续导数。r=m的简单样条只有k+m个自由度，r=m=3有着特殊的意义，因为它是三次多项式，该函数首次被人们称为样条函数。术语“三次样条”用于三维情况，此时进行曲面内插，而不是曲线内插。
由于离散子区间的范围较宽，可能是一条数字化的曲线，在这个范围内计算简单样条会引起一定的数学问题，因此在实际应用中都用B样条—一种特殊的样条函数。B样条是感兴趣区间以外均为零的其它样条的和，因此可按简单的方法用低次多项式进行局部拟合。
B样条经常用于数字化的线划在显示之前进行平滑处理，例如土壤、地质图上的各种边界，传统的制图总希望绘出较平滑的曲线。但是用B样条做多边形边界平滑也存在一些问题，特别是多边形面积和周长的计算，结果会与平滑前的不同。
综上所述，样条函数是分段函数，每次只用少量数据点，故插值速度快。样条函数与趋势面分析和移动平均方法相比，它保留了局部的变化特征。线性和曲面样条函数都在视觉上上得到了令人满意的结果。样条函数的一些缺点是：样条内插的误差不能直接估算，同时在实践中要解决的问题是样条块的定义以及如何在三维空间中将这些“块”拼成复杂曲面，又不引入原始曲面中所没有的异常现象等问题。
4）空间自协方差最佳插值方法：克里金插值
前面介绍的几个插值方法对影响插值效果的一些敏感性问题仍没有得到很
好的解决，例如趋势面分析的控制参数和距离倒数插值方法的权重对结果影响很大，这些问题包括：
●需要计算平均值数据点的数目；
●搜索数据点的邻域大小、方向和形状如何确定；
●有没有比计算简单距离函数更好的估计权重系数的方法；
●与插值有关的误差问题。
为解决这些问题，法国地理数学学家Georges Matheron和南非矿山工程师D.G.Krige研究了一种优化插值方法，用于矿山勘探。这个方法被广泛地应用于
地下水模拟、土壤制图等领域，成为GIS软件地理统计插值的重要组成部分。这种方法充分吸收了地理统计的思想，认为任何在空间连续性变化的属性是非常不规则的，不能用简单的平滑数学函数进行模拟，可以用随机表面给予较恰当的描述。这种连续性变化的空间属性称为“区域性变量”，可以描述象气压、高程及其它连续性变化的描述指标变量。这种应用地理统计方法进行空间插值的方法，被称为克里金（Kriging）插值。地理统计方法为空间插值提供了一种优化策略，即在插值过程中根据某种优化准则函数动态的决定变量的数值。Matheron ，Krige等人研究的插值方法着重于权重系数的确定，从而使内插函数处于最佳状态，即对给定点上的变量值提供最好的线性无偏估计。
克里金插值方法的区域性变量理论假设任何变量的空间变化都可以表示为
下述三个主要成分的和（图8-17）：
1）与恒定均值或趋势有关的结构性成分；
2）与空间变化有关的随机变量，即区域性变量；
3）与空间无关的随机噪声项或剩余误差项。
图8-17：区域变量理论将复杂的空间变化分为三个部分
(i)地形的平均特性；(ii)空间相关的不规则变化；(iii)随机的、局部的变化
另x为一维、二维或三维空间中的某一个位置，变量z在x处的值可又下式计算：
式中，m(x)是描述z(x)的结构性成分的确定性函数；是与空间变化有关的随机变化项，即区域性变量；是剩余误差项，空间上具有零平均值、与空间无关的高斯噪声项。
克里金方法的第一步是确定适当的m(x)函数，最简单的情况是m(x)等于采样区的平均值，距离矢量h分离的两点x，x+h之间的数学期望等于零：
式中z(x), z(x+h) 是随机变量z在x, x+h处的值，同时还假设两点之间的方差只与距离h有关，于是：
式中，是一种函数，称为半方差函数。
区域性变量理论的两个内在假设条件是差异的稳定性和可变性，一旦结构性成分确定后，剩余的差异变化属于同质变化，不同位置之间的差异仅是距离的函数。这样，区域性变量计算公式可以写成下式的形式：
半方差的估算公式如下：
式中，n为距离为h的采样点对的数目（n对点），采样间隔h也叫延迟。对应于h的的图被称为“半方差图”。图8-18表示一个典型的半方差图。
半方差是定量描述区域性变化的第一步，它为空间插值、优化采样方案提供了有益的信息。为了求得半方差图，必须先得到拟合半方差的理论模型，在半方差理论模型中：
1）延迟h的值较大的部分曲线呈水平方向。曲线的水平部分成为“梁（Sill）”。说明在延迟的这个范围内数据点没有空间相关性，因为所有的方差不随距离增减而变化。
2）曲线从的低值升到梁为止的延迟范围，称为“变程（Range）”。变程是半方差图最重要的部分，因为它描述了与空间有关的差异怎样随距离变化的。在变程范围内距离越近的点具有更相近的特征。变程给移动加权平均方法提供了一个确定窗口大小的方法。很显然，数据点和未知点之间的距离大于变
程范围，表明该数据点与未知点距离太远，对插值没有作用。
3）图中的拟合模型没有通过原点，而是在的正方向与坐标轴相截。
图8-18：半方差图
按半方差计算公式，当h=0时，必须为零。模型中出现的正值是剩余误差的估计值，它是与空间无关的噪声。称为“核(Nugget)”方差，是观测误差的和距离间隔很小的情况下的空间变化的组合。
当存在明显的变程和梁，同时核方差也很重要但数值不太大的情况下，可用球面模型进行半方差拟合（图8-19-1）。公式是：
0 < h < a
h >= a
h=0
式中，a是变程，h是延迟，c0+c1为梁。一般情况下用球面模型拟合效果比较理想。
如果存在明显的核方差和梁，而没有渐变的变程，则可用指数模型（图8-19-2）进行拟合：
如果核方差相对于与空间变化有关的随机变化很小的情况下，最好使用比较弯曲的曲线，如高斯曲线：（图8-19-4）
如果空间变化随变程渐变，但没有梁，则可用线性模型（图8-19-3）进行拟合：
式中，b为线的斜率。当变程的大小远超过人们希望的插值范围时，也用线性模型。
前面的讨论都假设地表特征在各个方向都是相同的，然而许多情况下空间变化中的都具有明显的方向性，这是就要用不同参数的模型来拟合半方差图。
图8-19：各种不同的变方差图
拟合后的半方差图重要的用途是确定局部内插需要的权重因子。确定的过程与加权移动插值方法类似，但不是按一种固定的函数计算，而是按采样点数据的半方差图的统计分析原理计算。即：
权重的选择应使是无偏估计，且估计的方差小于观测值的其它线性组合产生的方差。
的最小方差为：
只有下式成立时，才可获得的最小方差：
式中，是z在采样点x i, x j之间的半方差；是采样点x j和未知点x0之间的半方差，这两个量均可从已拟合模型的半方差图上得到。量计算最小方差需要的拉格朗日算子。
这个方法叫常规克里金插值。它是一个精确插值模型，内插值或最佳局部均值与数据点上的值一致。制图中常用比采样间隔更细的规则格网进行插值，内插值又可用前边提到的方法转换成等值线图。与此类似，估计的误差，又叫克里
以科学发展观为统领
I3992
I4 654
I5 536
解法过程用矩阵的表示如下：
领导人卡通形象
其中，A为数据点之间的半方差矩阵，b是每个数据点与未知点之间的半方差向量，λ为要计算的权重系数向量，Φ为解方程的拉格郎日算子。液基细胞学
首先计算这5个数据点之间的距离矩阵：
I 1 2 3 4 5
I1 0.0 5.099 9.899 5.000 3.162
I2 5.099 0.0 6.325 3.606 4.472
厦门槟榔小学I3 9.899 6.325 0.0 5.0 7.211
I4 5.0 3.606 5.0 0.0 2.236
I5 3.162 4.472 7.211 2.236 0.0
和它们与未知点之间的距离向量：
I I0
I1 4.243
I2 2.808
I3 5.657
I4 1.0
I5 2.0
将上述数值带入球面模型，得到相应的半方差（矩阵A, b）：
A= I 1 2 3 4 5 6
1 2.500 7.739 9.999 7.656 5.939 1.000细胞膜的渗透性
2 7.739 2.500 8.677 6.381 7.196 1.000
3 9.999 8.677 2.500 7.656 9.206 1.000
4 7.656 6.381 7.656 2.500 4.936 1.000
5 5.939 7.19
6 9.206 4.936 2.500 1.000
6 1.000 1.000 1.000 1.000 1.000 0.000
b = I 0
1 7.151
2 5.597
3 8.815
4 3.621
5 4.720
6 1.000
注意额外的第6行和第6列，是为了保证权重之和为1。
计算A的逆矩阵，得：
A-1= I 1 2 3 4 5 6
1 -0.17
2 0.050 0.022 -0.026 0.126 0.273
2 0.050 -0.167 0.032 0.077 0.007 0.207
吕宋岛
3 0.022 0.032 -0.111 0.066 -0.010 0.357
4 -0.026 0.077 0.066 -0.307 0.190 0.030
5 0.12
6 0.00
7 -0.010 0.190 -0.313 0.134
6 0.273 0.20
7 0.357 0.003 0.134 -6.873
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议