8-9数字高程模型数字地形分析18.pptVIP

洼地填平算法——Moran和Vezina算法 首先用一个极大高程水面数据将原始地面的DEM数据表面淹没，然后移除DEM上多余的水，最后得到的高程数据就是填洼处理的高程数据。 Ａ、对DEM数据初始化，除了边界外，用一极大临界水面高程数据将原始地面上DEM数据表面淹没。 Ｂ、通过迭代方法将W值逐步减小，最后W会收敛于Wf ,Wf代表的高程点数据既是填洼后的数据。 平地垫高算法—— Martz和Garbrecht算法 用高程增量叠加算法处理平地。对平地范围内的单元格增加一微小增量，每个单元格的增量大小是不一样的，就可以消除平地。 数据预处理方法 D8算法： A、如果最大落差值小于0，则赋以负值以表明此格网方向未定； B、如果最大落差值大于或等于0，且最大值只有一个，则将对应此最大值的方向值作为中心格网处的方向值； C 、如果最大落差等于0，且有一个以上的0值，则以这些0值所对应的方向值相加。在极端情况下，如果8个邻域高程值都与中心格网高程值相同，则中心格网方向值赋以255； D 、如果最大落差值大于0，且有一个以上的最大值，则任选一个方向作为水流方向。 被处理栅格单元同相邻8个栅格单元之间坡降的算法为： Slope=Dz/Di 式中: Slope为两个栅格之间的坡降； Dz为两个栅格单元之间的高程差； Di为两个栅格单元中心之间的距离。 32 64 128 16 K 1 8 4 2 水文分析——水流方向 78 72 69 71 58 49 74 67 56 49 46 50 69 53 44 37 38 48 64 58 55 22 31 24 68 61 47 21 16 19 74 53 34 12 11 12 2 2 2 4 4 8 2 2 2 4 4 8 1 1 2 4 8 4 128 128 1 2 4 8 2 2 1 4 4 4 1 1 1 1 4 16 32 64 128 16 K 1 8 4 2 原始DEM 矩阵 水流方向矩阵 水流方向提取结果 区域地形每点的流水累积量，它可以用区域地形表面的流水模拟方法获得。 规则格网表示的数字地面高程模型每点处有一个单位的水量，按照自然水流从高处流到低处的自然规律，根据区域地形的水流方向数字矩阵计算每点处所流过的水量数值，便可以得到该区域汇流累积量。 水文分析——汇流累积量 2 2 2 4 4 8 2 2 2 4 4 8 1 1 2 4 8 4 128 128 1 2 4 8 2 2 1 4 4 4 1 1 1 1 4 16 0 0 0 0 0 0 0 1 1 2 2 0 0 3 7 5 4 0 0 0 0 20 0 1 0 0 0 1 24 0 0 2 4 7 35 2 32 64 128 16 K 1 8 4 2 水流方向矩阵 汇流累积量矩阵 汇流累积量提取结果 河网的提取采用的是地表径流漫流模型，通过模拟地表径流的流动来产生水系。 A、在无洼地的DEM上计算出每一个栅格的水流方向矩阵； B、根据自然水流由高处流到低处的自然规律利用水流方向矩阵计算出汇流累积量； C、当汇流累积量达到一定值得时候，就会产生地表水流，所有汇流累积量大于临界值的栅格就是潜在的水流路径，由这些水流路径构成的网络就是河网。 在实现的过程中要确定汇流能力阈值，不低于该阈值的 单元格标记为水系组成部分，该方法简单，直接产生连续的 流线段。 水文分析——提取栅格河网 单阈值20 单阈值100 提取栅格河网结果 河流的分级一般依据河流的流量、形态等因素进行，不同级别的河网所代表的汇流累积量不同，级别越高，汇流累积量越大，一般是主流，而级别低的河网则一般是支流。 两种常用的分级方法： A、Strahler 分级： B、Shreve分级： 水文分析——栅格河网分级 将所有河网中没有支流的河网定义为一级，两个一级河网汇流成二级，如此下去分别为第3级，第4级，一直到河网出水口。当且仅当同级别的两条河网汇流成一条河网时级别才会才会增加，对于低级河网汇入高级情况，高级河网不会增加。 河网分级——STRAHLER分级 将所有河网中没有支流的河网定义为一级，两个一级河网汇流成二级，对于以后更高级别的河网级别是汇入的河网级别之和。 河网分级——SHREVE分级 Shreve strahler 河网分级结果对比 （a）流域DEM三维图； （b）流域DEM平面图；（c）基于DEM的河网生成图；（d）基于DEM的子流域划分图。 谢 谢！ 叙拉古国王希伦二世发现的公式 ,但根据 Morr