地理信息系统原理（高勇）09 数字地形模型与地形分析.pptVIP

格网DEM计算地形属性 单要素属性 可由高程数据直接计算得到 如坡度、坡向 复合属性 由几个单要素属性按一定关系组合成的复合指标 用于描述某种过程的空间变化 这种组合关系通常是经验关系，也可以使用简化的自然过程机理模型 坡度、坡向 坡度为水平面与局部地表之间的正切值。包含两个成分： 斜度——高度变化的最大值比率（常称为坡度） 地表单元的法向与Z轴的夹角，即切平面与水平面的夹角 坡向——变化比率最大值的方向 地表单元的法向量在水平面上的投影与X轴之间的夹角 地貌分析还可能用到二阶差分凹率和凸率 比较通用的度量方法是 斜度用百分比度量 坡向按从正北方向起算的角度测量 凸度按单位距离内斜度的度数测量 制图 坡度图 对不同的坡度设定不同的灰度级 坡向图 分为东、南、西、北、东北、西北、 东南、西南8类，再加上平地，共9类 用不同的色彩显示 格网DEM计算坡度、坡向 使用3*3窗口，窗口在DEM高程矩阵中连续移动后，完成整幅图的计算 坡度计算 坡向计算 二阶差分计算方法：提高运算速度 有限二阶差分法：点(i,j)在x方向上的斜度 式中σx是格网间距，沿对角线时应乘以sqrt(2) 计算八个方向的斜度，运算速度也快得多 但地面高程得局部误差将引起严重得坡度计算误差 数字分析方法 用数字分析方法计算东西方向得坡度公式 二阶差分法计算坡度坡向的详细说明 对于规则格网DEM，栅格尺寸为D 使用3×3规则窗口 窗口内每点的高程为ei 其中中心点高程为e0 窗口中心点的坡度S、坡向A： 二阶差分法 e5 e2 e6 e1 e0 e3 e8 e4 e7 格网DEM计算其他地形属性 剖面积 根据工程设计的线路，计算其与DEM各格网边交点Pi(Xi, Yi, Zi)，则线路剖面积为 n为交点数 Di,i+1为Pi与Pi+1之距离 同理可计算任意横断面及其面积 体积 DEM体积由四棱柱（无特征的格网）与三棱柱体积进行累加得到 四棱柱体上表面用抛物双曲面拟合 三棱柱体上表面用斜平面拟合 下表面均为水平面或参考平面 可计算工程中的挖方、填方及土壤流失量 表面积 对于含有特征的格网，将其分解成三角形 对于无特征的格网，可由4个角点的高程取平均即中心点高程，然后将格网分成4个三角形 由每一三角形的三个角点坐标（xi, yi, zi）计算出通过该三个顶点的斜面内三角形的面积 最后累加得到实地的表面积 剖面分析 意义 以线代面，研究区域的地貌形态、轮廓形状、地势变化、地质构造、斜坡特征、地表切割强度等 可以提供土地利用规划、工程选线和选址等的决策依据 分析方法 可在格网DEM或三角网DEM上进行。 已知A、B两点的坐标，则可求出两点连线与格网或三角网的交点，并内插交点上的高程，以及各交点之间的距离。然后按选定的垂直比例尺和水平比例尺，按距离和高程绘出剖面图。 剖面图不一定必须沿直线绘制，也可沿一条曲线绘制 开挖与回填 开挖土方量与回填土方量（挖方与填方） 开挖 回填 通视分析 通视分析 指以某一点为观察点，研究某一区域通视情况的地形分析 观察哨所的设定，显然观察哨的位置应该设在能监视某一感兴趣的区域，视线不能被地形挡住。 森林中火灾监测点的设定，无线发射塔的设定等 有时还可能对不可见区域进行分析，如低空侦察飞机在飞行时，要尽可能躲避敌方雷达的捕捉，飞行显然要选择雷达盲区飞行 通视问题 已知一个或一组观察点，找出某一地形的可见区域 欲观察到某一区域的全部地形表面，计算最少观察点数量 在观察点数量一定的前提下，计算能获得的最大观察区域 以最小代价建造观察塔，要求全部区域可见 在给定建造代价的前提下，求最大可见区 分类 点的通视：计算视点与待判定点之间的可见性问题 线的通视：已知视点，计算视点的视野问题 区域的通视：已知视点，计算视点能可视的地形表面区域集合的问题 基于格网模型与基于TIN模型的DEM计算通视的方法差异很大 通视示例 观察点 不通视 通视 雷达盲区飞行 可视域的三维显示 通视算法（1） 点对点通视 格网点作为计算单位。这样点对点的通视问题简化为离散空间直线与某一地形剖面线的相交问题 算法 已知视点V的坐标，以及P点的坐标。DEM为二维数组Z[M][N] 生成V到P的投影直线点集， 并得到直线点集对应的高程数据{Z[k], ( k=1,...K-1 )}，这样形成V到P的DEM剖面曲线。 以V到P的投影直线为X轴，V的投影点为原点，求出视线在X-Z坐标系的直线方程H 比较数组H[k]与数组Z[k]中对应元素的值，如果存在Z[k]H[k]，则V与P不可见，否则可见 V P 通视算法（2） 点对线通视 点对线的通视，实际上就是求点的视野 对于视野线之外的任何一个地形表面上的点都是不可见的，但在视野线内的点有可