27-曲面的内插和逼近探索.ppt

曲面内插方法与DEM数据管理 * 山东理工大学 建筑工程学院 测绘工程 TIN高程内插 Grid DEM内插 双线性内插 移动曲面内插方法 多面函数内插方法 有限元内插方法 主要内容 曲面的内插方法 DEM内插就是根据参考点上的高程求出其它待定点上的高程， 整体函数内插 局部函数内插 逐点内插法 用一个整体函数拟合整个区域 逐点内插方法 以每一待定点为中心，定义一个局部函数去拟合周围的数据点。逐点内插法十分灵活，精度较高，计算方法简单又不需很大的计算机内存，但计算速度可能比其它方法慢 （l）建立局部坐标 移动曲面拟合法 对DEM每一个格网点，将坐标原点移至该DEM格网点P（Xp，Yp） （3）列出误差方程式 （2）选取邻近数据点 P x y d i 误差方程式 由n个数据点列出的误差方程为 内插参数解算 （4）计算每一数据点的权 （5）法化求解 系数F是待定点内插高程值ZP 对点的选择除满足n＞6外，应保证各个象限都有数据点， 当地形起伏较大时，半径R不能取得很大。 当数据点较稀或分布不均匀时，利用二次曲面移动拟合可能产生很大的误差 移动曲面拟合法注意事项 采用了多个邻近点之加权平均水平面移动拟合法内插： 加权平均水平面移动拟合法 多面函数法DEM内插 “任何一个圆滑的数学表面总是可以用一系列有规则的数学表面的总和，以任意的精度进行逼近。”也就是一个数学表面上某点（X，Y）处高程Z的表达式为: 核函数 可任选其中n个为核函数的中心点Pj(Xj,Yj) 各数据点应满足 i=1,2…m 误差方程 法化求解得 任意一点上的高程Zk(Kn)为 其中 m=n 全部数据点取为核函数的中心 多面函数法解算 有限元法DEM内插 为了解算一个函数，把它分成为许多适当大小的“单元”，在每一单元中用一个简单的函数，例如多项式来近似地代表它。 一次样条有限元DEM内播 ?x, ?y是以格网边长为单位时点A相对于点Pij的坐标增量 A ?x ?y 若 A点是已知高程点，作为观测值，以格网高程Zi,j…作为待定的未知数 误差方程式 虚拟观测值误差方程式 断裂线的处理 HIFI内插过程中考虑计算单元中的断裂线的基本要点如下 （l）作线性内插，加密断裂线点，特别是断裂线与DEM格网线交点之平面坐标与高程 （2）将计算单元按断裂线划分成子区，确定每个子区由哪几条断裂线与边界线组成 （4）分子区进行内插计算。 （3）分子区内插的原则是：不属于该子区的数据点不参加该子区的平差计算 , 判断的方法“跌落法” 和符号判断法 。 断裂线的处理 DEM的精度 DEM精度与应用有密切关系，必须对DEM的精度进行估计 １．由地形功率谱与内插方法的传递函 数估计DEM 精度 采样间隔和地形的复杂程度 均方误差 *