GRD建模与内插幻灯片.pptVIP

第五讲、 DEM表面内插方法 一、DEM地形表面重建的数学机理 一、DEM地形表面重建的数学机理 DEM表面重建通用多项式 基于点的表面建模 基于格网的表面建模 基于三角形的表面建模 回顾：格网DEM生产流程 内插和外推 空间插值的理论假设 模型内插概述（1） 模型内插概述（2） DEM内插方法分类 2、模型整体内插 模型整体内插原理 整体内插的缺点 3、分块内插 分块内插原理 分块内插原理 3.1、线性内插 3.2、双线性内插 双线性多项式曲面插值 3.3、样条函数内插 样条函数内插 3.4、多面函数内插 多面函数内插的核函数 3.5、最小二乘拟合法 3.5、最小二乘拟合法 最小二乘拟合法 自然表面的拟合要求 趋势面(trend)方法 趋势面分析是一种整体内插法，这种方法假设地形表面的一般趋势与高程值的局部变化无关，并利用曲面方程来模拟待定点附近地形表面的一般趋势。 通常使用的是一次、二次、三次趋势面，过高次的趋势面不利于反映空间趋势，并可能存在趋势面的“畸变”。其中，二次趋势面可用待定点附近的6个数据点来计算方程式系数。若取用的点数多于6个，则要用测量平差方法来拟合。 趋势面(trend)方法产生的DEM难以反映地形起伏的特性。 克里金(Kriging)插值 Kriging插值方法在空间相关范围分析的基础上，用相关范围内的采样点来估计待插点属性值。 （1）数据检验与分析，删去明显偏离实际的采样数据点。 （2）数据预处理。 （3）绘制方差图，了解空间变量的集聚范围与方向。 （4）克里金插值估计。 相对来说，克里金插值则能较好地反映各种地形变化，但克里金方法的计算量很大，因此在对大面积区域大数据量内插时，这是一个不能不考虑的因素。? 4、模型逐点移面内插 逐点移面内插与分块、剖分内插的区别 模型逐点移面内插原理 4.1、邻域及其内点的确定(1) 4.1、邻域及其内点的确定(2) 4.2、邻域点权值计算 4.3、内插函数模型(1) 4.3、内插函数模型(2) 4.4、考虑地貌特征的逐点内插 逐点内插法小结 五、模型内插的效率问题 原始采样点的索引组织 六、基于内插的GRD模型生成 6.1、不规则采样点的Grd建立 6.2、规则采样点的Grd建立 6.3、基于等高线的Grd建立 逐点移面内插：以待插点为中心，确定一个邻域范围，用该邻域内的采样点计算内插点的高程值。 剖分内插可以有不同的形状和大小，但它们的位置是固定的（如TIN）。 分块内插的分块范围在内插过程中一经确定，其形状、大小和位置都保持不变，凡落在分块上的待插点都用展铺在该分块上的唯一确定的数学面进行内插。 逐点移面法一般有相同的形状和面积，位置却随内插点位置而移动，其以待插点为中心，定义一个局部函数去拟合周围的数据点，以适当半径或边长的圆或正方形作为移动面去捕捉适当数目的数据点，并以此展铺一张数学面，内插该中心的高程。数据点的范围随待插点位置的变化而变化。 以待插点为中心，定义一个局部函数去拟合周围的数据点，以适当半径或边长的圆或正方形作为移动面去捕捉适当数目的数据点，并以此展铺一张数学面，内插该中心的高程。数据点的范围随待插点位置的变化而变化。 主要特点： 基本步骤： 定义内插点的邻域范围 确定落在邻域内的采样点 选定内插数学模型 计算内插点的高程 利用搜索圆实现对数据点的搜索 k一般取7 按方位取点图示 利用搜索矩形实现对数据点的搜索 k一般取7 若内插点(x, y)与采样点(xi, yi) 满足下式则为邻域点： 注： Pi是参考点的权，R是圆的半径，di是待插点到参考点的距离。 常用的权函数 ： 赋权依据： 地形表面具有自相关性，较接近的采样点对内插点的影响要大一些。 因此采用内插点和采样点的距离来刻画采样点对内插点的贡献程度，即距离越近贡献越大，反之就小。 是移动拟合法的特例，它是在解算待定点的高程时，使用加权平均值代替误差方程： 反距离加权平均法 Zp是待定点户的高程，Zi是第i个参考点的高程值，n为参考点的个数，只是第i个参考点的权重，权函数及参考点范围的选取与移动拟合法相同。 利用内插点邻域内的采样点，可选用各种方法来计算内插点的高程值，如多层曲面叠加法、曲面拟合法、最小二乘法、有限元法、加权平均法等。 对于每个待插的点，可选取其邻近的n个数据点(可称其为参考点)拟合一多项式曲面，拟合的曲面可选用如下的形式： Z = AX2 + BXY + CY2 + DX + EY + F 式中，X、Y、Z是各参考点的坐标值，A、B、C、D、E、F为待定的参数。多项式中的各参数可由n个选定的参考点用