第7章+数字地面模型5.ppt

Horn 算法 The actual algorithm that is used to calculate slope is: rise_run = SQRT(SQR(dz/dx)+SQR(dz/dy)) degree_slope = ATAN(rise_run) * 57.29578 degree_aspect = ATAN((dz/dy)/(dz/dx)) * 57.29578 where the deltas are calculated using a 3x3 roving window. a through i represent the z_values in the window: (dz/dx) = ((a + 2d + g) - (c + 2f + i)) / (8 * x_mesh_spacing) (dz/dy) = ((a + 2b + c) - (g + 2h + i)) / (8 * y_mesh_spacing) 根据相邻网格点上的坡度和坡向之间的逻辑关系，可以判断坡形的凹凸变化情况，确定沟谷线、山脊和鞍部的位置，划分流域范围。 三、地面形态的计算 由于实际地形无一定数学规律可循，因此，影响DEM精度的主要因素是原始数据的获取，其中主要包括: 数据采集的密度; 数据采集位置（选点）。 注意:任何一种内插方法均不能弥补由于取样不当所造成的信息损失。 2、DEM数据源的采集 1) 数据采集的精度 2) 数据采集的方法 (1)人工网格法 将地形图蒙上格网，逐格读取中心或角点的高程值、构成数字高程模型。 (2)立体像对分析 通过遥感立体像对，根据视差模型、自动选配左右影像的同名点，可建立数字高程模型。 (3)三角网方法（TIN） (4)曲面拟合 根据有限个离散点的高程、采用多项式或样条函数求得拟合公式，再逐一计算各点的高程，可得到拟合的DEM。 对有限个离散点，每三个邻近点联接成三角形，每个三角形代表一个局部平面，再根据每个平面方程，可计算各格网点高程，生成DEM。 A（x1,y1，z1） B（x2,y2,z2） C（x3,y3,z3） D（x4,y4,z4） P（x,y） — 三角网上的任意数据获取 点数据的生成与叠加 (5) 等值线插值 根据各局部等值线上的高程点，通过插值公式计算各点的高程，得到DEM。 等值线插值法是比较常用的方法，输入等值线后，可在矢量格式的等值线数据基础上进行，插值效果较好。 （1）采集前根据建立DEM模型的精度要求确立合理的取样密度，如对单调地形可均匀采集，密度不必过大； （2）地形变化处采用密集采样，确保地形转折处的数据如山谷线，山脊线，谷缘线，崖线，山坡转折线等； （3）不应出现大的空白区，即使大片平坦地区也应有最底的采集密度。 3）数据采集的原则 采集的DEM数据必需进行后处理方可应用，后处理包括： (1) 格式转换 由于采集数据所存储格式包括数据内容、数据类型等各异，因此必须进行数据格式转换，使数据成为自己所要的格式。 (2) 坐标转换 有时采集数据要进行坐标系转换，如将像片坐标转换成大地坐标。 4) 数据采集的后处理 (3) 数据编辑处理 任何数据获取后总要编辑，编辑过程通常是一种交互过程，主要包括剔除错误的、过密的、重复的点；加密要加密区域的点。 (4) 数据分块 由于采集数据的不同，数据常具有不同排列顺序，如利用地形图采集的等高线数据通常按条采集和存储查寻到，而对等高线进行区域插值时要用待插点周周的数据为依据，为迅速查询到待插点周周的数据，需要将数据重新分块存储，分块方法通常是将整个区域分成等间隔的格网，每个格网之间有一定的重叠度，相互之间可用链指针连接。 1) 网格尺寸的确定 网格尺寸的确定是DEM数据网格化的重要问题，它关系到数据模型的精度。 网格尺寸的确定首先应分析地形的形态特征。 一般地说，采样点的密度基本上确定了网格点的密度。 在采集优选点的情况下，网格点数应大于或接近于样点数，取nN2n，N 为网格点数，n为样点数。 3、数据点的加密和插值 数据插值的核心问题是确定内插方法，任何一种内插方法都是基于原始函数光滑性，也就是说认为临近数据点之间存在很大的相关性。 2) 内插方法 DEM常用的插值算法： 1）距离加权法 距离愈远的点对待插点的影响愈小； 2）双线性插值法 假定待插点的高程在直线上呈比例变化，用双线性内插函数求出； 3）趋势面插值法 用通过采样点的空间曲面来模拟地形表面，常用二次多项式来拟合 4）样条插值法 用距离函数样条法，或分片样条法进行插值 5）克吕