第三章 空间数据的处理——内插.pptVIP

§5 空间数据的内插方法 空间数据内插的概念 是根据一组已知的离散点数据或分区数据，从现有这些数据中找到一个函数关系式，使该关系式最好 地逼近这些已知的空间数据，并能根据该函数关系 式推求出任意点或分区的值。这种通过已知点数据 推求出其他未知点或未知区域数据的方法。 空间数据内插算法的分类： 按具体内容分两大类：点的内插和区域内插； 按使用已知点范围分：整体拟合和局部拟合。 趋势面拟合法 是利用回归分析思想，运用最小二乘法对数据 点进行的拟合，通过选择一个二元函数来逼近采样 数据的整体变化趋势，二元函数的表达式：  P为二元函数阶数，通常≥1 趋势面法 趋势面的阶数 p=0,水平面 f (x，y)=b0 p=1,倾斜面 f(x,y)=b0+b1x+b2y p=2,二次曲面 f(x,y)=b0+ b1x+b2y +b3x2+b4xy+b5y2 p=3,复杂曲面常用三次 f(x,y)=b0+ b1x+b2y +b3x2+b4xy+b5y2 +b6x3+b7x2y+b8xy2+b9y3 趋势面法 函数方程的求解 该二元函数必须满足观测值与拟合值之差的平方和 最小，即 然后用最小二乘法求解各式的系数bi 例：趋势面分析法的运用 下图显示五个已知值的气象站点，用x,y坐标，格网单元大 小为2000m的行和列来表示，求0号站的未知值。 求解步骤： 建立线形趋势面方程：z x,y=b0+b1x+b2y 根据最小二乘法得到b0，b1，b2的三个法方程： 改成矩阵形势并代入数据 双线性多项式内插法 适用情况 局部待插区域中，当采样点的特征值在x、y 方向分别按线性规律变化 内插函数 f(x,y) = Ax + By + Cxy + D 采样点的选择 4个已知采样点 尽量以内插点为中心均匀分布 离内插点距离最近 例如： 数据按正方形格网（边长为L） 节点布置，球内点P的高程ZP 已知：A(i,j， ZA) B(i+1,j, ZB) C(i,j+1, ZC) D(i+1,j+1, ZD) 求解思路： 先用点A和B及C和D两对点的高程， 线性内插出点M和N的高程ZM ,ZN ； 然后再由ZM , ZN 直线内插待求 点P的高程ZP. 移动拟合内插法 原理 1)待内插的点P为中心原点，按一定半径作圆； 2)选定一多项式内插函数，用圆内的采样点 解出函数参数； 3)使用内插函数计算待插点的特征值； 移动拟合内插法 假设取二次多项式来拟合，则待求点的特征值可写为一般式： f(x,y) = Ax2 + Bxy + Cy2 + Dx + Ey + F 通常做法是将坐标原点放置在待定点上，而采用的 数据点应落在半径为R的圆内（图） f(x,y) = Ax′2 + Bx′y′+ Cy′2 +Dx′+Ey′+F 按距离加权插值法(DEM) 设平面上分布一系列离散点，已知其坐标和高程为xi,yi，zi(I=1,2,…n),p(x,y)为任一格网点，根据周围离散点的高程，通过距离加权插值求得p点高程。这时： 区域内插算法——叠置法 是将目标区叠置在源区上，首先 确定两者面积的交集，然后计算 出目标区各个分区的内插值； 例如： 已知源区：3个县的人口统计数据 求目标区：3个流域的人口数 区域内插算法——叠置法 步骤 （1）将目标区叠加在源区上，计算两者的交集的面积 Ats （2）计算内插后目标区中各分区的人口值 Vt = Ats ( As )-1 Us Vt目标区的分区人口矩阵 Ats目标区与源区交集的面积矩阵 As 源区的分区面积矩阵 Us 源区的分区人口矩阵 区域内插算法——叠置法 计算结果式 Vtj = 上式适用条件 已知源区的各分区的某项绝对统计值，如：人口总数 区域内插算法——比重法 原理 根据平滑密度函数的原理，将源区的统计数据从各分区的均匀分布（同质性）转变为分区内的非均匀分布（非同质性）。 应用 非同质性代表着一般社会经济与自然现象的特点，如：行政线交界处气象数据等； 比重法 算法步骤 1）对源区（us）上栅格化，尺寸选择; 2）根据原图各分区的面积和栅格尺寸，将各分区的统计数据按比例赋与相应分区的各个格网点，并计算全图的统计数据总和U; 3）按公式计算相邻四个格网点的平均值赋给该格网点;