第四章空间形态与空间格局分解.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 变换函数插值 根据一个或多个空间参量的经验方程进行整体空间插值，也是经常使用的空间插值方法，这种经验方程称为变换函数。 冲积平原的土壤重金属污染与几个重要因子有关，其中距污染源（河流）的距离，和高程两个因子最重要。一般情况，携带重金属的粗粒泥沙沉积在河滩上，携带重金属的细粒泥沙沉淀在低洼的在洪水期容易被淹没的地方，而那些洪水频率低的地方，由于携带重金属污染泥沙颗粒比较少，受到污染轻。由于距河流的距离和高程是比较容易得到的空间变量，可以用各种重金属含量与它们的经验方程进行空间插值，以改进对重金属污染的预测。 举 例 本例回归方程的形式如下： 式中是z x 某种重金属含量（ppm），b0…bn是回归系数，p1…pn是独立空间变量，本例p1是距河流的距离因子，p2是高程因子。 局部插值方法只使用邻近的数据点来估计未知点的值，包括几个步骤： 1）定义一个邻域或搜索范围； 2）搜索落在此邻域范围的数据点； 3）选择表达这有限个点的空间变化的数学函数； 4）为落在规则格网单元上的数据点赋值。 重复这个步骤直到格网上的所有点赋值完毕。 使用局部插值方法需要注意的几个方面是： 所使用的插值函数； 邻域的大小、形状和方向； 数据点的个数； 数据点的分布方式是规则的还是不规则的。 4. 空间插值 4.3 空间插值方法 －局部插值方法 最近邻点法：泰森多边形方法 泰森多边形（Thiessen，又叫Dirichlet 或Voronoi多边形）采用了一种极端的边界内插方法，只用最近的单个点进行区域插值。泰森多边形按数据点位置将区域分割成子区域，每个子区域包含一个数据点，各子区域到其内数据点的距离小于任何到其它数据点的距离，并用其内数据点进行赋值。连接所有数据点的连线形成Delaunay三角形，与不规则三角网TIN具有相同的拓扑结构。 GIS和地理分析中经常采用泰森多边形进行快速的赋值，实际上泰森多边形的一个隐含的假设是任何地点的气象数据均使用距它最近的气象站的数据。而实际上，除非是有足够多的气象站，否则这个假设是不恰当的，因为降水、气压、温度等现象是连续变化的，用泰森多边形插值方法得到的结果图变化只发生在边界上，在边界内都是均质的和无变化的。 移动平均插值方法：距离倒数插值 距离倒数插值方法综合了泰森多边形的邻近点方法和趋势面分析的渐变方法的长处，它假设未知点x0处属性值是在局部邻域内中所有数据点的距离加权平均值。距离倒数插值方法是加权移动平均方法的一种。加权移动平均方法的计算公式如下： 加权移动平均公式最简单的形式称为线性插值，公式如下： 距离倒数插值方法是GIS软件根据点数据生成栅格图层的最常见方法。距离倒数法计算值易受数据点集群的影响，计算结果经常出现一种孤立点数据明显高于周围数据点的“鸭蛋”分布模式，可以在插值过程中通过动态修改搜索准则进行一定程度的改进。 样条函数插值方法 在计算机用于曲线与数据点拟合以前，绘图员是使用一种灵活的曲线规逐段的拟合出平滑的曲线。这种灵活的曲线规绘出的分段曲线称为样条。与样条匹配的那些数据点称为桩点，绘制曲线时桩点控制曲线的位置。曲线规绘出的曲线在数学上用分段的三次多项式函数来描述这种曲线，其连接处有连续的一阶和二阶连续导数。 样条函数是数学上与灵活曲线规对等的一个数学等式，是一个分段函数，进行一次拟合只有与少数点拟合，同时保证曲线段连接处连续。这就意味着样条函数可以修改少数数据点配准而不必重新计算整条曲线，趋势面分析方法做不到这一点，（图4.29）。 样条函数的局部特征 （a：当二次样条曲线的一个点位置变化时，只需要重新计算四段曲线； b：而一次样条曲线的一个点位置变化时，只需要重新计算两段曲线） 一般的分段多项式p x 定义为： p x pi x xi x xi+1 i 1, 2, 3…, k-1 j 0, 1, 2, ..., r-1; i 1, 2, ... , k-1 x1, ..., xk-1将区间x0，xk分成k个子区间，这些分割点称“断点”，曲线上具有这些x值的点称为“节”。 函数pi x 为小于等于m次的多项式。r项用来表示样条函数的约束条件： r 0，无约束； r 1，函数连续且对它的导数无任何约束； r m-1，区间[x0 xk]可用一个多项式表示； r m，约束条件最多。