第六章：数字高程模型的内插题库.pptx

第六章 数字高程模型的内插;教学目的与要求;本章重点与难点;内容提要;1.1内插方法的分类;1.1内插方法的分类;1.1内插方法的分类;1.2整体内插;1.2整体内插;1.2整体内插;1.2整体内插;2.1分块内插;2.1分块内插;2.1分块内插（线性内插）;10/18/2015;第16页;2.1分块内插（双线性内插）;10/18/2015;2.1分块内插（二元样条函数内插）;2.1分块内插（二元样条函数内插）;2.1分块内插（二元样条函数内插）;2.1分块内插（二元样条函数内插）;2.1分块内插（二元样条函数内插）;2.1分块内插（二元样条函数内插）;2.1分块内插（多面叠加内插法）;10/18/2015;第27页;10/18/2015;10/18/2015;2.1分块内插（多面叠加内插法）;10/18/2015;该函数的改进形式的函数如下： 式中(dki)min2表示数据点i与距离最近的数据点k的距离。当n＝m时，Q矩阵也不是对称矩阵，因为在每个数据点上有各自的参数(dki)min。利用该核函数可以很好地考虑地貌结构线的作用，此时只要沿地貌结构线上取一组密集数据点(或先内插出来)，就会产生很小的(dki)min值，结果在双曲面顶端产生一个大的斜率，由此保证了内插曲面上突变性的转折。 对多层叠加面的解算，可通过将m个参考点的三维坐标代入．得一误差方程组，按最小二乘法解求待定系数(m＞n)。 ;2.1分块内插（多面叠加内插法）;2.1分块内插（最小二乘配置法）;2.1分块内插（最小二乘配置法）;2.1分块内插（最小二乘配置法）;10/18/2015;2.1分块内插（最小二乘配置法）;2.1分块内插（最小二乘配置法）;n为分块扩充范围内参考点的个数。按最小二乘法相关平差方法求解，得到趋势面系数向量： 任一内插点P的信号为 用待插点在趋势面上的高程hp’加上待插点的信号sp’，即得所求待插点的高程Hp’： ; 最小二乘配置法数字高程分块内插的关键问题之一是如何建立z或s的协方差矩阵，也就是如何解决信号相关性规律的问题。 由数理统计理论得知，二维各态历经性平稳随机过程的协方差仅与不同点间的水平距离有关；最小二乘配置法内插高程时，认为信号和趋势面起算高程的协方差仅与点间的水平距离有关：距离愈近，协方差愈大，超过一定的距离．协方差趋近于零值。高斯函数正好满足函数值随距离缩短而增大的条件，所以习惯上以高斯函数作为相关函数，用来计算协方差。 ;10/18/2015;10/18/2015;10/18/2015;10/18/2015;10/18/2015;10/18/2015; 最小二乘配置法有严密的数理统计理论依据，但大量的试验结果表明，它未必能在数字地面模型内插应用中取得良好的拟合效果，原因主要是以下两点： (1)应用最小二乘配置法的前提是处理对象必须属于遍历性平稳随机过程。但实际地表起伏现象都十分复杂，各类地貌形态未必都符合各态历经性平稳随机过程的统计规律，地面点间趋势面起算高程的相关度量未必仅与距离有关。实际上，大多数地貌变化都不是各向同性的，地表起伏的相关性不仅与距离有关，也与方向有关。如果前提条件不符合，就难以保证得到良好的内插质量。 (2)确定趋势面和协方差函数的参数是一个循环迭代过程。当迭代收敛速度慢时，其计算量可能比大多数高程内插算法都大，因而此方法并不实用。 ;10/18/2015;有限元法是以离散方式处理连续量的一种数学方法，它的思路是将一定范围的连续整体分割成为有限个单元(如三角形、正方形等)的集合。相邻单元边界的端点称做结点，通过解求各个结点处的物理量来描述物理量的整体分布。 ;有限元法通常采用分片光滑的奇次样条作为各个单元的内插函数，已经用于数字高程模型内插的有双线性B样条和双三次B样条两种，其整体解是一系列基函数的线性组合，形式 为解求上述函数的全部系数．须列出与所求问题等价的二次泛函数取极小值的条件，建立并计算系数向量的线性方程组，使上式有确定解。 有限元法的计算量取决于分块内结点的个数，而不是像上述其他分块内插方法那样主要与参考点个数有关。所以单元划分越细．有限元法的计算量越大。 ;分块内插的分块范围在内插过程中一经确定，其形状、大小相位置都保持不变，凡落在分块上的待插点都用展铺在该分块上的惟一确定的数学面进行内插。逐点内插法是以待插点为中心、定义一个局部函数去拟合周围的数据点，数据点的范围随待插点位置的变化而移动，因此又称移动曲面法。 ;对于每个待插的点，可选取其邻近的n个数据点(可称其为参考点)拟合一多项式曲面 拟合的曲面可选用如下的形式： Z＝AX2十BX2十CY2十DX十EY十F 式中：X，Y，Z是各参考点的坐标值，