ArcGIS实验-Ex20-沟谷网络地提取及沟壑密度地计算.docVIP

实用标准文档 文案大全 第十一章 水文分析 练习3：沟谷网络的提取及沟壑密度的计算 一、背景 沟壑密度是描述地面被水道切割破碎程度的一个术语。沟壑密度越大，地面越破碎。破碎的地面必然起伏不平，多斜坡。这样一方面使地表物质稳定性降低，另一方面易形成地表径流。沟壑密度越大，地面径流和土壤冲刷越快，沟蚀发展越快。沟壑密度是地形发育阶段、降水量或地势高差、土壤渗透能力和地表抗蚀能力的重要特征值。所以，沟壑密度是气候、地形、岩性、植被等因素综合影响的反映，对于水土流失监测、水土保持规划有着重要的意义。因此，沟壑密度的提取具有在地形分析中，特别是在进行水土保持和水土流失研究中具有重要的意义。 沟壑密度也称沟谷密度或沟道密度，指单位面积内沟壑的总长度。单位一般以km/km2表示，数学表达为： DS = ΣL/A ，式中：DS指沟壑密度；ΣL指研究区域内的沟壑总长度（单位：km），A指特定研究区域的面积（单位：km2）。 二、目的 能熟悉掌握水文分析中的沟谷网络的提取原理及过程。 三、要求 1、用水文地质分析工具提取出研究区域的沟谷网络； 2、算出该研究区域的沟壑密度。 四、数据 一幅25m分辨率的黄土地貌DEM数据，数据的区域大概有59 km2。数据存放于…/ChP11 /Ex3中，请将其拷贝到E：/ChP11/Ex3。结果保存在…/ChP11/Ex3/Result文件夹中。 五、操作步骤 1、沟谷网络的提取： (1) 启动ArcToolbox，展开Analysis Tools工具箱，打开hydrology工具集。在图层管理器中加载研究区域的原始DEM数据。 (2) 原始DEM数据提取水流方向数据：双击hydrology工具集中的flow direction工具，选择原始DEM数据作为输入表面数据，将输出的水流方向数据命名为flowdir，如图1。 图1 水流方向数据的计算 (3) 洼地的计算：双击hydrology工具集中的sink工具，选择上一步计算的水流方向数据flowdire作为输入数据，将输出数据命名为sink，如图2。 图2 洼地的计算 (4) 计算得到原始DEM上有洼地，需要进行洼地填充。双击hydrology工具集中的fill工具，进行原始DEM的洼地点填充。在Input surface raster文本框中选择原始DEM数据dem，将输出数据命名为filldem，因为选择的是将所有洼地全部填充，所有在填充容限Z limit为默认值，如图3。 图3 洼地填充 (5) 基于无洼地的水流方向的计算：同步骤2）中一样，打开水流方向计算，选择的输入表面数据是无洼地DEM数据filldem，将输出的水流方向数据命名为flowdirfill，如图4。 图4 基于无洼地的水流方向的计算 (6) 汇流累积数据的计算。双击hydrology工具集中的flow accumulation工具。选择flowdirfill数据作为输入的水流方向数据；输出数据命名为flowacc；输入的权值数据不选择，利用系统默认的等权且权值为1的模式。点击OK进行运算。运算出的汇流累积数据flowacc，如图5所示。 图5 汇流累积数据的计算 (7) 栅格河网的生成。在栅格河网的生成中，需要设置一个汇流累积阈值。双击spatial analysis tools工具箱中的map algebra工具集中的multi output map工具。在对话框的文本框中输入：E:\chp11\ex3\result\streamnet = con (E:\chp11\ex3\result\flowacc 100, 1) ，如图6，计算出栅格河网数据streamnet，结果如图7所示。 图6 栅格河网的生成 （用Map Algebra总会出问题，要注意写法！） 图7 栅格河网 (8) 栅格河网矢量化。双击hydrology工具集中的stream to feature工具。选择streamnet作为栅格河网数据输入，将水流方向数据flowdirrfill作为输入的水流方向数据。将输出的数据命名为stream1，也就是矢量形式的沟谷网络。 图8 栅格河网矢量化 (9) 伪沟谷的删除。由于基于DEM的河网的提取是采用最大坡降的方法，那么在平地区域例如谷底等，在这些区域上的水流方向是随机的，那么就很容易生成平行状的河流等错误的形状，这种平行状的沟谷被称为伪沟谷，需要进行手工编辑剔除的。伪沟谷的剔除步骤为： A．加载模块，点击editor模块的下拉箭头，点击start editing菜单工具，如图9，在操作目标层（target）中选择矢量形式的河网数据stream1。 B．点击编辑工具中的选择符号，将其激活。在ArcMap地图显示窗口中选择那些沟谷网络上的平行状