第四节 空间数据处理-5.pptVIP

第五节 空间数据的内插方法 空间数据内插概念 设已知一组空间数据，它们可以是离散点的形式，也可以是分区数据的形式，空间数据的内插就是从这些数据中找到一个函数关系式，使该关系最好地逼近这些已知的空间数据，并能根据该函数关系式推求出区域范围内其他任意点或任意分区的值。 空间数据内插方法 根据已知点和已知分区数据的不同，将空间数据内插分为 点的内插 区域的内插 第五节 空间数据的内插方法 一、点内插 插值过程中，由于取样的数据点呈离散分布形式，或者数据点虽然按格网排列，但格网的密度不能满足使用的要求，这样就需要以数据点为基础进行插值运算。 插值运算的手段是选择一个合理的数学模型，使用已知点上的信息求出函数的待定系数。通常采用： 逐点内插 局部函数内插 第五节 空间数据的内插方法 （1）局部函数内插 局部函数法通常以格网小块为加密区，采用低次项函数拟合地表面。 线性内插 双线性多项式内插 双三次多项式（样条函数） 第五节 空间数据的内插方法 ①线性内插 当分块插值区的地表面为一平面，按直线比例内插待定点的高程 。 数学模型 使用最靠近待定点的3个数据点来定义函数。设待定点高程为zp，则线性内插函数为： zp = a0 + a1 x + a2 y  将3个数据点的坐标值代入，联立求解出系数a0 ， a1 ， a2 ，待定点在给予平面坐标xp，yp之后，即可求出内插高程zp . 第五节 空间数据的内插方法 ②双线性多项式内插 当分块插值区内待定点的高程在轴x（和y ）平行的方向上与坐标y（和x）成直线比例关系。 数学模型 设待定点的高程为zp，则按双线性多项式内插的函数式为： zp = a0 + a1 x + a2 y + a3 x y 此时取最靠近的4个数据点定义函数的待定参数。 第五节 空间数据的内插方法 ③双三次多项式（样条函数） 当分块插值区用双三次多项式即样条函数模拟地表面时。 数学模型 设待定点的高程为zp，函数式为： zp= f(x，y)=a1x3y3 + a2x2y3 + a3xy3 + a4y3 + a5x3y2 + a6x2y2 + a7xy2 + a8y2 +a9x3y + a10x2y + a11xy + a12y + a13x3 + a14x2 + a15x + a16 需要4个格网节点的数据，以及x和y方向的斜率、曲面的扭曲来定义函数的待定值。双三次多项式样条函数是正方形格网布点时内插加密常用方法之一。 第五节 空间数据的内插方法 （2）逐点内插 移动拟合法是典型的逐点内插法。移动拟合法是指对每一个待定点取用一个多项式曲面拟合该点附近的地表面。 此时取待定点作平面坐标的原点，并用待定点为圆心，以 R 为半径的圆内取数据点来定义函数的待定系数。 数学模型 设取二次多项式来拟合，则待求点的高程可写成一般式为： zp = Ax2 + Bxy + Cy2 + Dx + Ey + F 上式中有6个待定系数，至少取用6个点的数据代入方程，求解待定系数。 第五节 空间数据的内插方法 当取用的数据点多于6个时，则以数据点高程z作为观测值，可列出误差方程式 Q = A x2 + B x y + C y2 + D x + E y + F - Z 并以数据点到待求点的距离给予适当的权重，按最小二乘法平差原理求解。权值应与距离成反比，间距愈近，对待求点测定值的影响应愈大。 如取 W = 1 / d2或( (R-d) /d )2 式中：d为待定点到数据点间的水平距离， R为定义函数待定参数时所求的圆半径。 该方法称为按距离加权最小二乘内插算法。 第五节 空间数据的内插方法 点内插实现过程 ①数据取样： 数据整理过程，包括简单的数据插值运算、数据检查等； ②数据内插方法选择： 根据数据特点，选择数学模型进行数据内插运算； ③数据记录： 将内插结果按一定的数据结构形式，存储起来，供系统使用。 第五节 空间数据的内插方法 主要方法： 叠置法 比重法 区域内插主要解决离散属性数据问题。研究的目标是从已知分区的数据(如社会经济数据)中推出同一地区的另一组分区数据。 如： 已知某地区各县中历年的人口分布数据，但因行政区划分使该地区中某些县的边界线发生了变化，现在需要推算新行政区中历年的人口分布数据，就可用这种插值方法。 二、区域内插 人口 面积 A 35 7 B 30 6 C 10 3 A