第5章空间数据处理素材.ppt

9.1 不规则三角构网 9.2 空间内插 9.3 区域内插 9.1 不规则三角构网 不规则三角构网有不同的方法（准则），如最大-最小距离准则、圆准则、最大-最小内角准则等；方法不同，构网结果也不完全相同 从理论上讲，三角网的建立应该基于最佳三角形条件：保证每个三角形是锐角三角形或三角形边长近似相等，避免出现过大的钝角和过小的锐角 书中介绍了角度判别法建立TIN、Delaunay三角构建方法 1）三角网生长算法的基本思路 1）在点集中任选一点B1,找出离该点最近的另一点B2 2）根据余弦定理计算被选第三顶点为角顶点的三角形内角的大小，其中选择最大者对应的点为该三角形的第三顶点B3,构成第一个三角形和新增两条边 3）以新增三角形新增的两条边按步骤2）的方法（角度最大）扩展三角形 4）重复进行，并进行检测（任意一条边的只能被利用两次，如果利用三次则无效；新增边不能已形成的三角形相交） 为提高邻近点的检索速度，需要对数据分块 构建的TIN不唯一 2）Delaunay三角网（简称D-TIN） 各种不同的TIN构建方法中，所形成的TIN应该满足三个基本要求： 惟一性（TIN是惟一的） 最大最小角特性（三角形最小内角尽量最大，尽量接近等边三角形） 空圆性（三角形边长和尽量最小、且三角形外接圆中不包含其它三角形顶点，也称Delaunay 法则） 由于原始数据的不同， Delaunay三角网有多种构网方法：基于Vonoroi图（泰森多边形）构建TIN；根据离散点直接构建TIN 基于离散点直接构网方法又分无约束和有约束方法 无约束法可分：逐点插入、三角网生长法、分治算法、凸包法等 3）三角网生长算法的基本思路 首先找出离散点集中相距最短的两点，连接两点成为D-TIN的初始基线 然后按D-TIN的判断法则找出包含此基线的Delaunay三角形的第3端点，该端点位于基线右测 连接新点与原来两点形成两条新边 再按按D-TIN的判断法则找出包含此两边的另外两个Delaunay三角形的第3端点 依次循环处理所有的新边，直到所有离散点均成为D-TIN的端点 4）逐点插入法的基本步骤 先定义一个包含所有数据的初始多边形 从数据集中任意选择一点（如A），插入到初始多边形中建立初始三角网 然后按以下步骤进行迭代计算，直到所有数据点被处理 插入一个离散采样点P，在初始三角网中找出包含P点的三角形T，把P与T的三个顶点相连，生成3个新的小三角性 用局部优化方法（LOP,Local Optimization Procedure）从里到外进行优化，确保所有点内插后的三角网为D-TIN 9.2 空间内插 内插(Interpolation)含义: 将离散的数据点转化为连续的数据曲面 用已知点来估算其他未知点的过程 需要插值的原因: 现有离散曲面的分辨率、像元大小、方向与要求不符 现有连续曲面的数据模型与要求不一致 现有数据不能完全覆盖所要求的区域 内插含义: 内插和外推 Sampled points Estimated points 空间插值基础：控制点 控制点分布 控制点密度 控制点的自相关程度 内插方法 内插方法根据内插时使用采样点（控制点）的范围内插可以分为：全局内插方法（整体拟合） 、局部内插方法（局部拟合） 根据具体的插值内容可以分为：点内插、区域内插 全局内插方法（整体拟合）是基于研究区域内的所有采样点建立（考虑全部控制点）拟合方程的方法，如 趋势面法、回归模型、最小二乘法 特点 不能提供内插区域的局部特性,如金矿品位富集、辐射源局部异常等 通常用于拟合大范围、长周期变化情况，如沙漠地貌、平原地貌、地下水位等 内插结果粗糙 局部内插方法（局部拟合）是指仅用邻近于位置点的少数已知采样点来估算未知点的特征值，如 样条函数法、移动平均法、克里金内插方法、线性插值法 特点 可以提供内插区域的局部特性,且不受其它区域的内插影响 通常用于金矿品位估计、高程内插等 内插结果具有精确性 1）内插方法：全局内插 z=f( x, y) Control Points 常用全局内插函数 常用全局内插函数特性分析 1st Order Trend Original surface 2nd Order Trend 3rd Order Trend 系数求解 2）内插方法：局部内插 移动拟合内插 选择的数据点数6 计算每个采样点的权，应与距离有关 根据平差理论求解参数，F值就是内插点的高程 多面函数法插值 基本思想是：任何一个圆滑的数学表面总可以用一系列有规则的数学表面的总和以任意的精度进行逼近 若有m个控制点，选择其中的n个，为核函数的中心点Pj(Xj,Yj)，令 则数据点满足 得误差方程 法方程求解得 任意一点的高程值为 克里金（Krigin