第五章 空间数据的处理.ppt

6.简述多边形数据栅格化的夹角之和法基本原理。7.什么是区域的内插? 叠置法和比重法有什么不同?  §5 空间数据的内插 ? 定义：通过已知点或多边形分区的数据，推求任意点或多边形分区数据的方法称为空间数据的内插（或插值）。? l? 点的内插： 已知点数据 ——任意点数据 l??区域的内插：已知分区数据——重新分区数据 B C A D (b)分区数据 两种空间数据 (a)离散点数据 4．1 点的内插 点的内插是用来建立具有连续变化特征现象（如高程）的数值方法。它是用未知点附近的已知点的数据，推测未知点的数据。 以建立数字高程模型为例，点的内插方法可分为几个类型： N次多项式拟合法 整体内插法 克里金法 加权平均法 移动拟合法 逐点内插法 二元样条函数内插法 双线性多项式内插法 线性内插法 分块内插 点的空间内插 点的空间内插法分类 点插值包括3步骤(数字高程模型建立为例)： 数据取样、数据内插、数据精度分析 （一）数据取样 等值线取样：数字化地形图等高线； 随机取样：地形特征点（坡度、坡向改变处）； 规则取样：理想的方法。周期性将影响结果。 取样方式 (b)随机取样 (c)规则取样 (a)等值线取样 （二）数据内插 根据已知离散点数据，选择并拟合数据模型，用模型推算任意点数据。由于地面形态复杂，用低次多项式拟合整个地形面不分切合实际，整体内插法少用。常用两种方式： 分块内插法 逐点内插法 1.分块内插法 把整体内插空间划分成若干分块(三角形、正方形)，根据每块的取样点，求各分块的曲面函数的插值方法。 取样方法： l??不规则三角网：线性内插法 l???规则格网：双线性多项式内插法、双三次多项式内插法? (1)线性内插法 线性内插法是先将所有已知数据点连接成三角网的形式，使用靠近内插点的三个已知数据点，确定三解网中的一个三角形形成的空间平面，继而求出该内插点在平面中的高程值。 x y 分块内插的线性内插法 0 P P1 P2 P3 z P 线性内插函数： z=a0+a1x+a2y 根据最靠近待定点的3个已知数据点P1（x1,y1,z1）、P2（x2,y2,z2）、P3（x3,y3,z3），代入平面方程，求系数 a0、a1、a2： z1=a0+a1x1+a2y1 z2=a0+a1x2+a2y2 z3=a0+a1x3+a2y3 矩阵形式 把系数a0、a1、a2代入方程，可求出任意给定平面点P(xP,yP)的内插高程zP zP=a0+a1xp+a2yp x y 分块内插的线性内插法 0 P P1 P2 P3 z P (2)双线性多项式内插法 双线性多项式内插法常常是在规则分布的已知数据点时，使用最靠近内插点的4个已知数据点组成一个四边形，确定一个双线性多项式来内插其中点的高程。 x y 分块内插的双线性多项式内插法 0 P2 P4 z P1 P3 P 把a0、a1、a2、a3代入方程，可求出任意待定点P(xP,yP)的zP zP=a0+a1xp+a2yp+a3xpyp 双线性多项式内插函数形式 z=a0+a1x+a2y+a3xy 设已知4个参考点P1(x1,y1,z1), P2(x2,y2,z2), P3(x3,y3,z3), P4(x4,y4,z4) (3)双三次多项式插值（二元样条函数） l正方形格网，每个分块单元拟合双三次多项式曲面； l?相邻曲面连接处的x轴方向的斜率R保持连续； l?相邻曲面连接处的y轴方向的斜率S保持连续； l?相邻曲面连接处的扭矩T连续。 1 2 3 8 4 5 7 6 C(0,1) A(0,0) B(1,0) z y x D(1,1) 样条函数的插值方法 双三次多项式插值模型，待定点高程 z=f(x,y) = a1x3y3+ a2x2y3+ a3xy3+ a4y3 + a5x3y2+ a6x2y2+ a7xy2+a8y2 + a9x3y+ a10x2y +a11xy+ a12y +a13x3 + a14x2 + a15x + a16 式中，有16个待定系数，必须列出16个线性方程，才能确定它们的数值。 按正方形格网分块。 分块边长为单位1。分块左下角为坐标原点。 样条函数的插值方法 y 15 14 13 12 7 11 8 0 4 10 1 2 3 9 A(0