第五章数字高程模型内插.pptVIP

* 第五章 数字高程模型内插 5.1 内插方法的分类 5.2 整体内插 5.3 分块内插 5.4 逐点内插 5.5 关于内插方法的探讨 5.1 内插方法的分类 5.2 整体内插 由研究区域内所有采样点的观测值建立的。主要通过多项式函数来实现，又 称整体函数法内插。常被用于模拟大范围的宏观变化趋势。 取研究范围内不同平面位置的n个参考点三维坐标，带入方程内，使n阶线性 方程组有唯一解。将待插点的坐标代入上式，可得到待定点的高程值。 整体内插法的优缺点。 5.3 分块内插 分块内插是把参考空间分成若干块，对各分块使用不同的函数。要考虑各 相邻分块函数间的连续问题。 5.3.1 线性内插（重点介绍内容） 首先使用最靠近插值点的三个已知数据点确定一个平面，继而求出内插点的 高程值的方法。 （1）概念 根据已知三个参考点A,B,C双线性内插p点高程值： 注意：最邻近插值一般不连续．具有连续性的最简单的插值是分片线性插值． 最邻近插值 x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? y (x1, y1) (x1, y2) (x2, y1) (x2, y2) O 二维或高维情形的最邻近插值，与被插值点最邻近的 节点的函数值即为所求． （2）Matlab算法实现 要求x0,y0单调；x，y可取为矩阵，或x取行向量，y取为列向量，x,y的值分别不能超出x0,y0的范围． z=interp2(x0,y0,z0,x,y,’method’) 被插值点 插值方法 用MATLAB作网格节点数据的插值 插值节点 被插值点的函数值 ‘nearest’ 最邻近插值； ‘linear’ 双线性插值； ‘cubic’ 双三次插值； 缺省时 双线性插值. 5.3.2 双线性多项式内插（重点介绍内容） 使用最靠近插值点的四个已知数据点组成一个四边形，进而确定一个双线性 多项式来内插待插点的高程。 正方形格网分布的双线性内插公式： （2）Matlab算法实现 双线性插值是一片一片的空间二次曲面构成． 双线性插值函数的形式如下： 其中有四个待定系数，利用该函数在矩形的四个顶点（插值节点）的函数值，得到四个代数方程，正好确定四个系数． 双线性插值 x ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? y (x1, y1) (x1, y2) (x2, y1) (x2, y2) O 要求x0,y0单调；x，y可取为矩阵，或x取行向量，y取为列向量，x,y的值分别不能超出x0,y0的范围． z=interp2(x0,y0,z0,x,y,’method’) 被插值点 插值方法 用MATLAB作网格节点数据的插值 插值节点 被插值点的函数值 ‘nearest’ 最邻近插值； ‘linear’ 双线性插值； ‘cubic’ 双三次插值； 缺省时 双线性插值. 5.3.3 反距离权重插值算法 一种局部插值方法，假设未知值的点受较近控制点的影像比较远控制点的影响更大。 Z0为点0的估计值；zi为控制点i的z值；di为控制点i与点0间的距离；s为在估算中 用到的控制点的数目；K为指定的幂。 课下练习：编写C语言程序实现该插值算法。 5.3.4 二元样条函数内插 为保证各分块曲面间的光滑性，按照弹性力学条件使所确定的n次多项式曲面 与其相邻分块的边界上所有n-1次导数都连续，这n次多项式就称为样条函数。 方程的16个待定系数，需要建立16个方程。四角点高程ZA,ZB,ZC,ZD，以及它们 的导数值RA,RB,RC,RD,SA,SB,SC,SD和TA,TB,TC,TD。 优点：保留了微地物特征；内插速度快；保证了分块间连接处为平滑联系的曲面。 缺点：采样了弹性力学条件，而地表分块不是狭义的弹性壳体。 5.3.5 多面叠加内插法（多面函数法） 基本思想是任何一个规则的或不规则的连续曲面均可以由若干个简单面（或称 单值数学面）来叠加逼近。具体做法是在每个数据点上建立一个曲面，然后在 Z方向上将各个旋转曲面按一定比例叠加成一张整体的连续曲面，使之严格的 通过各个数据点。 Q为简单数学面，又称多面函数的核函数；n为简单数学面的张数，其值与分块 扩充范围内参考点的个数相等；Ki为待定参数，代表了第i个核函数对多层叠加面 的贡献。 5.3.6 最小二乘配置法（自学） 某一个测量值包含三部分： 与某些参数有关的值。测量值是这些参数的函数，而这个函数在空间是一个 曲面，故称为趋势面。 不能简单用某个函数表达的值，称为系统的信号部分。 随机噪声。 去掉趋势面后，如果观测值包含信号和噪声两部分（且信号和噪声期望均为0， 两者的协方差也为零），则可获得信