空间数据插值

空间数据插值第一页，共二十八页，2022年，8月28日所谓空间数据插值,即通过探寻收集到的样点/样方数据的规律,外推/内插到整个研究区域为面数据的方法.即根据已知区域的数据求算待估区域值,影响插值精度的主要因素就是插值法的选取.第二页，共二十八页，2022年，8月28日空间数据插值方法的基本原理:任何一种空间数据插值法都是基于空间相关性的基础上进行的。即空间位置上越靠近,则事物或现象就越相似,空间位置越远,则越相异或者越不相关，体现了事物/现象对空间位置的依赖关系。（，南京师范大学地理科学学院地理信息系统专业网络课程教程）由于经典统计建模通常要求因变量是纯随机独立变量，而空间插值则要求插值变量具备某种程度的空间自相关性的具随机性和结构性的区域化变量。即区域内部是随机的，与位置无关的，而在整体的空间分布上又是有一定的规律可循的，这也是不宜用简单的统计分析方法进行插值预估的原因。从而空间统计学应用而生。无论用哪种插值方法，根据统计学假设可知，样本点越多越好，而样本的分布越均匀越好。第三页，共二十八页，2022年，8月28日常用的空间数据插值方法之一：趋势面分析趋势面分析（Trendanalyst）。严格来说趋势面分析并不是在一种空间数据插值法。它是根据采样点的地理坐标X，Y值与样点的属性Z值建立多元回归模型，前提假设是，Z值是独立变量且呈正态分布，其回归误差与位置无关。根据自行设置的参数可建立线性、二次…或n次多项式回归模型，从而得到不同的拟合平面，可以是平面，亦可以是曲面。精度以最小二乘法进行验证。第四页，共二十八页，2022年，8月28日常用的空间数据插值方法之一：趋势面分析趋势面分析中，将Z值分解成如下等式：

Z=++全局趋势局部变异全区的，规模较大的地理过程的反映变化缓慢规模较小的局部区域的地理过程变化较快，局部异常，由于空间数据不具备重复抽样条件，所以通常将后两项合并。趋势值即回归值，而后两项将合并到拟合残差中。在趋势面拟合中，空间位置以平面坐标为佳，即将经纬度坐标转换为以米为单位的平面大地坐标。通常趋势面分析用于分析趋势和异常而不追求高的拟合精度，一般达到60-80%，阶数在1-4之间即可。拟合精度按R^2系数和F值检验。随机干扰抽样误差或观测误差，不含系统误差第五页，共二十八页，2022年，8月28日常用的空间数据插值方法之一：趋势面分析由上述可知，趋势面分析是经典统计学在点数据进行空间展面上的应用，属于全局多项式插值，即对整个研究区域用一个多项式进行拟合。它的缺点在于：当研究区域范围较大，地形很复杂时，需要用高阶多项式拟合以提高精度，但高阶将增加其计算成本，因而需要进行改进。第六页，共二十八页，2022年，8月28日常用的空间数据插值方法之二：局部多项式插值局部多项式插值（LocalPolynomialInterpolation）：用多个多项式进行拟合。每个多项式都只在特定重叠的邻近区域内有效，通过设定搜索半径和方向的来定义邻近区域。

显然，局部多项式插值是对全局多项式，即趋势面拟合的一大改进。这里涉及到一个搜索邻域的概念。第七页，共二十八页，2022年，8月28日从空间自相关性的概念可知，空间上越靠近，属性就越相似，相关性也越高。那么，两个样点间在多远的距离内所具备相关性可以不考虑，或者其相关将消失呢？可以根据经验或专业背景找出这么一个阈值，作为邻近区域的半径。同时，如果其自相关性在不同的方向上消失的距离值也不同的话，将还需要设置一个方向值以及长短两个半径值，此时的邻近区域将呈椭圆。（如当属性值受风向影响较大时，应当将风向角度设置为搜索方向，即长半径所在的方向）空间数据插值之邻近区域：第八页，共二十八页，2022年，8月28日空间数据插值之邻近区域：通过半径和方向可以定义出一个以待估点为中心的区域（圆或者椭圆）。此外，还可以通过限制参与某待估点值进行预测的样点数来定义邻近区域。即参与某点预测的最多样点数和最少样点数。在由半径和方向决定的区域内包含到的样点数为0时，则扩大搜索区域使其达到最小样点数值。第九页，共二十八页，2022年，8月28日空间数据插值之各向异性：在设定邻近区域时，提到了一个方向参数。即当空间相关性沿各个方向上的消失距离都一致时，其邻近区域应该是一个圆，如图a，叫各向同性。否则，如图b，在西南-东北方向上的消失距离明显小于东南-西北方向，则其邻近区域应当是一个平行于东南-西北方向的椭圆，其方向角度（AngleDirection）设为长轴与X轴的角度值。图b的现象即各向异性（Anisotropy)。(图片来源：ArcgisDesjktopHelp文件)图中的Range（变程）参数，即自相关消失或不予考虑的半径值。图b中的MinorRange,最小变程，即相关性消失得最快的方向上的半径值，而MajorRange，最大变程即相关性消失最慢的方向上的半径值。图a图b第十页，共二十八页，2022年，8月28日常用的空间数据插值方法之三：移动平均插值法（MovingAverage)移动平均插值法，通过设定邻近区域，取该区域内样点的平均值作为待估点的值。移动平均插值适用于样点分布均匀、密集，而且变化缓慢的情况下，对缺失值进行填补。移动平均主要用于消除随机干扰，即局部降噪功能。移动平均插值的优势在于计算简便快速，但适用范围较窄。第十一页，共二十八页，2022年，8月28日常用的空间数据插值方法之四：线性三角网法（TriangulaionwithLinearInterpolation）线性三角网法是最佳的Delaunay三角形，连续样点数据间的连线形成三角形，覆盖整个研究区域，所有三角形的边都不相交。（即与构建TIN文件的原理一致）线性三角网法将在整个研究区域内均匀分配数据，地图上的稀疏区域会形成截然不同的三角面。第十二页，共二十八页，2022年，8月28日常用的空间数据插值方法之五：最近邻点插值法（NearestNeighbor)最近邻点插值法，又称泰森多边形（Thiessen或Voronoi多边形）分析法。即在每个样点数据周边生成一个邻近区域，即Thiessen多边形，使得每个多边形内的任意一点离其内部的样点最近，在多边形内插值时只有其中心样点参与运算，如图：最近邻点插值法同样只适用于样点分布均匀、紧密完整，且只有少数缺失值时，对缺失值进行填补第十三页，共二十八页，2022年，8月28日常用的空间数据插值方法之六：自然邻近插值法（NaturalNeighbor)自然邻近插值法是对泰森多边形插值法的改进。它对研究区域内各点都赋予一个权重系数，插值时使用邻点的权重平均值决定待估点的权重。每完成一次估值就将新值纳入原样点数据集重新计算泰松多边形并重新赋权重，再对下一待估点进行估值运算。对于由样点数据展面生成栅格数据而言，通过设置栅格大小（cellsize)来决定自然邻近插值中的泰森多边形的运行次数n，即，设整个研究区域的面积area，则有：n=area/cellsize可设置各向异性参数（半径和方向）来辅助权重系数的计算。

第十四页，共二十八页，2022年，8月28日常用的空间数据插值方法之七：反距离权重插值法(InverseDistanceWeighting,IDW)反距离权重插值综合了泰森多边形的自然邻近法和多元回归渐变方法的长处，在插值时为待估点Z值为邻近区域内所有数据点都的距离加权平均值，当有各向异性时，还要考虑方向权重。权重函数与待估点到样点间的距离的U次幂成反比，即随着距离增大，权重呈幂函数递减。且对某待估点而言，其所有邻域的样点数的权重和为1。决定反距离权重插值法结果的参数包括距离的U次幂值的确定，同时还取决于确定邻近区域的所使用的方法。此外，为消除样点数据的不均匀分布的影响，还可设置引入一个平滑参数，以保证没有哪个样点被赋予全部的权重，即使得插值运算时尽可能不只有一个样点参与运算。第十五页，共二十八页，2022年，8月28日常用的空间数据插值方法之七：反距离权重插值法(InverseDistanceWeighting,IDW)IDW是一种全局插值法，即全部样点都参与某一待估点的Z值的估算；IDW的适用于呈均匀分布且密集程度足以反映局部差异的样点数据集；IDW与之前介绍的插值法的不同之处在于，它是一种精确的插值法，即插值生成的表面中预测的样点值与实测样点值完全相等。第十六页，共二十八页，2022年，8月28日常用的空间数据插值方法之八：最小曲率法(MinimumCurvature)

最小曲率插值法，非精确插值法。其插值基准是生成一个具有最小曲率（即弯曲度最小），且到各样点的Z值的距离最小的曲面。

影响最小曲率插值法精度的参数有：最大残差，通常允许残差在10%-1%之间最大循环次数，与栅格大小（cellsize）有关,通常设置为生成的栅格数量的一到两倍。第十七页，共二十八页，2022年，8月28日常用的空间数据插值方法之九：径向基函数插值法（RadialBasisFunction)所谓径向基函数即基函数是由单个变量的函数构成的，是一系列精确插值法的统称。该插值法中的单个变量是指待估点到样点间的距离H，其中每一插值法都是距离H的基函数。径向基函数是对最小曲率插值的改进，即属于精确的最小曲率插值法。第十八页，共二十八页，2022年，8月28日常用的空间数据插值方法之九：径向基函数插值法（RadialBasisFunction)径向基函数包括的多种函数有：倒转复二次函数（InverseMultiquadric),复对数（Multilog),复二次函数（Multiquadratic),自然三次样条函数（NaturalCubicSpline),薄板样条法函数（ThinPlateSpline);上述的每一函数式中都带有一个平滑因子R，即使得生成的曲面不至于太粗糙。在实际应用中，许多人都发现复二次函数的效果最佳。径向基函数比同为精确插值法IDW的优点在于，它可以计算出高于或低于样点Z值的预测值。通常俗称的样条插值法即径向基函数插值法。此后在实际应用中又发展出了多种样条插值法，包括GRASS软件的RST，RegulationSplinewithTension,ANUSPLINE样条插值软件自带的多种样条插值法；大大提升了样条插值的精度。第十九页，共二十八页，2022年，8月28日常用的空间数据插值方法之九：径向基函数插值法（RadialBasisFunction)径向基函数适用于样点数据集大、表面变化平缓的情况；当局部变异性大，且无法确定样点数据的准确性，或样点数据具很大不确定性时，不适用该技术。第二十页，共二十八页，2022年，8月28日常用的空间数据插值方法之十：地统计插值法（GeostatisticalAnalyst)前面提到的多种基于空间统计学的插值方法都属于确定性插值法。而另一类插值法就是地统计插值法，它是空间统计分析的一个分支。地统计与确定性插值的最大区别在于，地统计插值引入了概率模型，即地统计插值认为从一个统计模型不可能完全精确地得出预测值，所以在进行预测时，应该给出预测值的误差，即预测值在一定概率内合理。通常所说的地统计插值是指克里格插值法（Kriging)Z（s）=μ（s）+ε（s）S表示不同的位置点，可以是用经纬度表示的空间坐标。Z（s）是该位置点的属性值。μ（s）为确定趋势值，ε（s）为自相关随机误差。第二十一页，共二十八页，2022年，8月28日常用的空间数据插值方法之十：地统计插值法（GeostatisticalAnalyst)当要考虑多个协同变量的情况下，可采用协克里格插值法（Co-Kriging).则其计算公式将变为：Zj（s）=μj（s）+εj（s）表示的是第j个变量的情况。在协克里格中，只对主变量进行预估，但将在插值预估时引入不同变量间的随机误差项εj（s）的交叉相关性值，从而构建协同克里格模型。第二十二页，共二十八页，2022年，8月28日常用的空间数据插值方法之十：地统计插值法（GeostatisticalAnalyst）克里格插值是一个最优的无偏估计法。获得预测图并不要求数据呈正态分布。但当数据呈正态分布时，克里格插值法将是无偏估计法中效果最好的一种方法。因此，在进行克里格插值前，可先对非正态分布的数据进行转换，包括Log对数转换，Box-Cox转换，使之呈正态分布，然后再进行插值。第二十三页，共二十八页，2022年，8月28日常用的空间数据插值方法之十：地统计插值法（GeostatisticalAnalyst）根据样点数据统计特征的不同可将克里格分成多种不同的插值法：当样点数据是二进制值时，用指示克里格插值法进行概率预测；对样点数据进行了未知函数变换后，可用该变换函数进行析取克里格插值；当样点数据的趋势值μ（s）是一个未知常量时，用普通克里格；当样点数据的趋势可用一个多项式进行拟合，但回归系数未知时，用泛克里格插值法；当样点数据的趋势已知时，用简单克里格插值法；其中最常用的是普通克里格与泛克里格插值法；当加入了协变量进行插值时，则叫作协同普通克里格插值法和协同泛克里格插值法。第二十四页，共二十八页，2022年，8月28日常用的空间数据插值方法之十：地统计插值法（GeostatisticalAnalyst）同反距离权重插值法IDW一样,克里格插值法同样可以表示为:Z(x0)=∑λiZ(xi)Z(x0)为待估点的值,Z(xi)为待估点周围的已知样点值,λi为第i个已知点的权重所不同之处在于,IDW的权重为待估点与已知样点间距离的u次幂的倒数,而克里格的权重值不仅考虑待估点与已知样点、已知样点之间的距离,还考虑了其空间分布的方位。通过半变异函数来赋权重值。第二十五页，共二十八页，2022年，8月28日克里格插值之半变异函数和协方差函数（Semivariogram/covariance）：半变异函数和协方差函数都是空间自相关性的定量化表达函数。半变异函数的定义为：

r(h)=1/2var(Z(si)-Z(sj))

hij为样点si和

sj间的距离,Z(si)和Z(sj)分别代表样点的属性值；若把si和