数字地面模析.pptVIP

2、 平均高程 设某区域具有等间距的n个规则格网，每个格网的高程为h，就可计算区域平均高程。 1）区域的平均高程 其中n为区域内格网个数；h（Pk）为第k个格网的（平均）高程。 每个格网的平均高程 h(pk) =(h00 + h01 + h10 + h11)/4 h00、 h01 、 h10、 h11分别是格网的四个角点的实际高程值。 1）坡度的定义 地面的坡度是表示地表面上该点倾斜的程度。 从数学说，地面的坡度是一个矢量，它既有大小，又有方向，坡度定义为地形表面和水平面之间夹角的正切值。 用矢量来描述坡度时，坡度是地面的法向量与Z轴之夹角slope 。 1、坡度和坡向的定义 二、坡度和坡向的计算 Y Z X ( i, j ) ( i, j+1) ( i+1,j ) ( i+1, j+1) n 坡度计算 2）坡向的定义（Aspect） 地面的坡向就是坡面的朝向，粗略的可分为向南、向北、向东、向西。更细的还可进一步分为南、北、东、西、东南、西南、东南和西北向。 DEM模型中坡向定义为地形表面法线在水平面上的投影与X轴之间夹角。在地学领域中，通常根据法线在水平面上的投影的位置，将其分成阳坡、阴坡、半阳坡、和半阴坡。 阴 阴 半阴 半阴 半阳 半阳 阳 阳 N S E W dir Y Z X ( i, j ) ( i, j+1) ( i+1,j ) ( i+1, j+1) n 坡向计算 坡向分类 坡度、坡向分析 根据坡度、坡向的不同，对grid进行再分类计算，得到一个新的grid 每一个单元，计算与其相邻的８个单元的之间的坡度、坡向 坡度：高程差／d（单元的大小） 坡向：每一个单元，通过与相邻８个单元的关系，计算平均方向或最大的方向。在ArcView 中是最大方向，以方位角0－360度表示坡向 坡度、坡向分析 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: a b c d e f g h i (dz/dx) = ((a + 2d + g) - (c + 2f + i)) / (8 * x_mesh_spacing) (dz/dy) = ((a + 2b + c) - (g + 2h + i)) / (8 * y_mesh_spacing) 根据相邻网格点上的坡度和坡向之间的逻辑关系，可以判断坡形的凹凸变化情况，确定沟谷线、山脊和鞍部的位置，划分流域范围。 三、地面形态的计算  1）谷点 满足下面4个条件 ( Zi,j-1 - Zi,j )* ( Zi,j+1 - Zi,j) 0 AND ( Zi-1,j - Zi,j )* ( Zi+1,j - Zi,j) 0 AND Zi,j+1 Zi,j AND Zi+1,j Zi,j 2) 脊点 满足下面4个条件 ( Zi,j-1 - Zi,j )* ( Zi,j+1 - Zi,j) 0 AND ( Zi,j-1 - Zi,j )* ( Zi+1,j - Zi,j) 0 AND Zi,j+1 Zi,j AND Zi+1,j Zi,j i,j-1 i,j i,j+1 i-1,j-1 i-1,j i-1,j+1 i+1,j-1 i+1,j i+1,j+1 1、山谷点、山脊点的计算 2、沟谷密度 沟谷密度是表征地面破碎程度的一种指标，它是沟谷总长度(∑L)与地表单元总面积(∑A)之比。提取谷点和脊点，将地表单元内所有谷点在单元区域内的延伸长度累加，便获得单元的沟谷长度。沟谷密度为: 反映某一面积单元内地势起伏变化与侵蚀程度的复杂程度，是地表面积与投影面积之比： 3、地表粗糙度 计算辐照度