第五节 DEM的分析和应用.docVIP

第五节 DEM的分析和应用 1、格网DEM应用 2、三角网DEM的应用 1．地形曲面拟合 DEM最基础的应用是求DEM范围内任意点的高程，在此基础上进行地形属性分析。由于已知有限个格网点的高程，可以利用这些格网点高程拟合一个地形曲面，推求区域内任意点的高程。曲面拟合方法可以看作是一个已知规则格网点数据进行空间插值的特例，距离倒数加权平均方法，克里金插值方法，样条函数等插值方法均可采用。 2 ．立体透视图 透视立体图能更好地反映地形的立体形态，非常直观。与采用等高线表示地形形态相比有其自身独特的优点，更接近人们的直观视觉。从一个空间三维的立体的数字高程模型到一个平面的二维透视图，其本质就是一个透视变换。将?°视点?±看作为?°摄影中心?±，可以直接应用共线方程从物点（X，Y，Z）计算?°像点?±坐标（X,Y）。透视图中的另一个问题是?°消隐?±的问题，即处理前景挡后景的问题。 3 ．通视分析 通视问题可以分为五类： 1）已知一个或一组观察点，找出某一地形的可见区域。 2）欲观察到某一区域的全部地形表面，计算最少观察点数量。 3）在观察点数量一定的前提下，计算能获得的最大观察区域。 4）以最小代价建造观察塔，要求全部区域可见。 5）在给定建造代价的前提下，求最大可见区。 两点是否可见的常用算法： a）倾角法 O(xo，yo，zo)为观察点，P(xp,yp,zp)为某一格网点，OP与格网的交点为A、B、C。 设OP的倾角为α O与各交点的倾角为βi (i＝A,B,C) 若tgα＞max(tgβi )，则OP通视 否则，不通视。 b) 剖面图 两点连线是否与剖面相交 AA A A B 4 ．流域特征地貌提取与地形自动分割 基于格网DEM自动提取流域特征地貌和进行地形自动分割技术主要包括两个方面： 1）流域地貌形态结构定义，定义能反映流域结构的特征地貌，建立格网DEM对应的微地貌特征。 2）特征地貌自动提取和地形自动分割算法。格网DEM数据是一些离散的高程点数据，每个数据本身不能反映实际地表的复杂性。为了从格网DEM数据中得到流域地貌形态结构，必须采用一个清晰的流域地貌结构模型，然后针对该结构模型设计自动提取算法。 凹陷区域识别与处理 平坦区域识别与处理 1）流域结构定义 可以使用一个具有根的树状图来描述流域结构[Shreve]，目前绝大多数算法都沿用这一描述方法。在此结构中主要包括三个部分，即结点集、界线集和汇流区集。 a.内部沟谷段 b.外部沟谷段 c.内部汇流区 d. 外部汇流区 e.沟谷结点 f.汇流源点 g. 分水线段 h.分水线源点 2）流域特征地貌自动提取和地形自动分割 特征地貌定义与提取：根据网格点高程与周围高程值的关系，将格网点分为坡地、洼地、分水线、谷地、阶地和鞍部等几类。 先计算中心点与八邻点的高程差 然后对高程差进行排序，再根据高程差序列的特性给中心点格网赋一个特征编码 最后通过一系列特征码的组合特征，用模式识别的方法，将格网点划分到已知的特征地貌类别。 山脊线和山谷线提取 ：山脊线和山谷线的自动探测实际上是凹点和凸点的自动搜索。较为简单的算子是2*2的局部算子。将算子在DEM数据中滑动，比较每个格网点与行和列上相邻格网点的高程，标出其中高程最小（探测山谷线）或高程最大（探测山脊线）的格网点。对整个DEM数据计算一遍后，剩下的未标记格网点就是山脊线或山谷线上的格网点。 流域地形自动分割：流域地形自动分割的目标是将整个流域分割成一个个子汇流区。大多数算法是利用3*3窗口计算流向和基于“溢流跟踪”算法确定汇流网络。 （1）格网点流向定义 采用3×3窗口按8方向搜索计算最大坡向为各网格点的流向。分别为8方向赋不同的代码，如右图所示。每个格网有一个从1到9的数值，代表它流向相邻象元的方向，如该象元为凹点，则其值为5 （2）凹点处理 搜索所有凹点的相邻最低点（有时可能有多个高程相等的最低点），作为凹点的溢出点，以溢出点为起点继续搜索比它的高程低或相等的邻点（已经搜索的点忽略），判断是否有比原凹点更低的格网点，如果没有则以该凹点的溢出点为起点，重复上述搜索过程；如果搜索到比原凹点低的格网点，将凹点和最低邻点的方向倒转 （3）提取汇流网络 给定一个点，搜索8邻点，记录所有流向它的格网点的位置，然后再以找到的格网点为基点继续搜索记录流向它的格网点，直到没有新的汇流点为止，所有记录的格网点构成该点的汇流区。 设定阈值，将汇流区面积大于此阈值的格网点，标识为沟谷点。 沟谷网络编码: 沟谷结点编码:从流域出口开始搜索遍历整个汇流网络，对每个沟谷段的上下游结点进行编码标识，标识值是沟谷段的编码值，并记录下这些结点的位置。 把沟谷段中的每个格网点标识为沟谷段的编码值。 根据沟谷段上游结点的类型判定