沟谷网络的提取及沟壑密度的计算.docxVIP

沟谷网络的提取及沟壑密度的计算1.背景:沟壑密度是描述地面被沟壑切割破碎程度的一个指标,沟壑密度是气候,地形,岩性,植被等因素综合影响的反映.沟壑密度越大,地面越破碎,平均坡度增大,地表物质稳定性降低,且易形成地表径流,土壤侵蚀加剧.因此,沟壑密度的测定对于了解区域地形发育特征,水土流失监测,水土保持规划有着重要意义.沟壑密度也成为沟谷密度或沟道密度,指单位面积内沟壑的总长度,以公里/公里2为单位,数学表达式 ,其中DS为沟壑密度,为样区内的沟壑总长度,A指特定样区的面积.2.目的:掌握利用水文分析工具提取沟壑网络的原理和方法3.要求:(1)利用水文分析工具提取样区的沟壑网络(2)计算出该样区的沟壑密度4.数据:25米分辨率的DEM数据,区域面积约为59平方公里5.操作步骤:(1)沟谷网络的提取1)在ArcMap中加载原始DEM数据2)打开Hydrology工具3)原始DEM数据提取水流方向,双击Hydrology工具集中地flowDirection工具,打开Flow Direction计算对话框.在InputSurfaceData文本框中选择原始DEM数据,指定输出水流方向数据名为flowdir,单击OK,进行水流方向数据的计算.4)洼地的计算:双击Hydrology工具集中的sink工具.如图:5)计算得知原始DEM数据上有洼地,需要进行洼地填充,双击FIll工具,进行原始DEM的洼地填充6)基于无洼地的水流方向计算:同步骤3)中一样,打开FlowDirection计算对话框,选择的输入表面数据时无洼地的DEM数据filldem,将输出的水流方向数据命名为flowdirfill7)汇流累积数据的计算.双击Hydrology工具集中的FlowAccumulation工具,如图8)栅格河网的生成:如图所示设置汇流累积阈值为100键入公式：E：\streamnet=con(E:\flowacc100,1),得到栅格河网数据streamnet9）栅格河网数据矢量化：双击Hydrology工具集中的Stream To Feature工具，选择streamnet作为河网输入数据，将flowdirfill作为水流方向输入数据，指定输出的数据名称为stream1，也就是矢量形式的沟谷网络。10）伪沟谷的删除。由于基于DEM的河网的提取时采用最大坡降的方法，在平地区域（例如谷底等地）的水流方向是随机的，容易生成平行状河流，这种平行状的沟谷称为伪沟谷，需以手工编辑剔除。伪沟谷的剔除步骤为：A.在ArcMap中加载Editor模块，单击Editor模块的下拉箭头，选择Start Editing，在操作目标层（target）中选择stream1.B单击编辑工具中的选择符号，在ArcMap地图显示窗口中选择那些平行状的沟谷，单击右键，在弹出菜单中选择Delete。样区边缘的很短的沟谷也需要进行删除。C.所有的伪沟谷删除完毕之后，单击Editor下拉菜单中的SaveEdits，保存编辑结果，并单击StopEdits停止编辑。至此，完成了伪沟谷的删除。（2）沟壑密度的计算（偶在9.3版本中没找到Xtools啊，需要的童鞋可以在这里下载安装/Soft/gxrj/2006-07-11/10.html）1）加载Xtools扩展模块。单击Xtools下拉箭头，选择Calculate area，Perimeter，Lenghth，Acres and Hectares菜单工具，弹出如图对话框，在文本框中选择需要计算的矢量数据，单击OK，进行计算，计算结果存储于属性表中的字段中。2）打开计算长度后的矢量沟壑数据的属性表，在属性字段Length上单击右键，弹出如图所示的对话框，选择Statistics命令，进行统计。3）通过统计计算得到的沟壑总长度4）计算样区的面积，在图层管理器中单击右键DEM数据层，单击属性，进入source选项卡，查询栅格数据的行、列以及栅格分辨率，即可计算出样区的面积5）得到研究区的沟壑总长度以及研究区的面积后，计算样区的沟壑密度。