DEM实验DEM实验ARCGIS等高线实验

实验四
沟谷网络的提取及沟壑密度的计算
学习并掌握利用水文分析工具提取沟谷网络的原理和方法。

1．沟谷网络的提取
<1）在ArcMap中加载原始DEM，打开Hydrology工具集；
图1
<2）运用Flow Direction工具在原始DEM上提取水流方向<flowdir）；Spatial analyst tool → hydrology→flow direction
在弹出的flow direction窗口中
Input :→C:\dem04\Ex1\dem
Output：→C:\dem04\Ex1\flowdir
图2flowdir
<3）运用Sink工具进行洼地的计算<sink）；
Spatial analyst tool → hydrology→sink
Input： C:\dem04\Ex1\flowdir
OutputC:\dem04\Ex1\Sink_flowdir1
由图像可知，且流向并不遵从D8算法，值域为1-192，证明存在洼地。

<4）计算得知原始DEM上有洼地，运用Fill工具进行洼地填充<filldem）；
Spatial analyst tool → hydrology→fill
Input： C:\dem04\Ex1\dem
Output C:\dem04\Ex1\Fill_dem1
图3fill_dem1
<5）运用Flow Direction工具在无洼地DEM上提取水流方向<flowdirfill）；
Spatial analyst tool → hydrology→flow direction
Input： C:\dem04\Ex1\fill_dem1
Output C:\dem04\Ex1\Flowdir_fill2
此时的流向就满足D8算法，分别是1、2、4、8、16、32、64、128 <6）运用Flow Accumulation工具进行汇流累积数据的计算<flowacc）；
Spatial analyst tool → hydrology→f。

ow accumulation
Input： C:\dem04\Ex1\flowdir_fill2
Output C:\dem04\Ex1\FlowDir_flow1
图4flowacc_flow1
<7）根据实际情况设置汇流累积阈值<100），生成栅格河网<streamnet）；
进入raster calculator ，输入con([flowacc_flow1] > 100 , 1>
即可得到汇流累计值超过100的河网
得到calculation
<8）运用Stream to Feature工具，进行栅格河网矢量化；
Spatial analyst tool → hydrology→stream to feature
Input：calculation
Input flow_direction raster：flowdir_fill2
Output：stream.shp
<9）伪沟谷的删除<在Editor模块下进行）。

加载上等高线和地貌晕眩进行比较,将那些伪沟谷删除。

2.沟壑密度的计算
在catalog中新建一个个人数据库
在ex1文件夹上右击→New→personal geodatabase
再在新建的personal geodatabase上右击→import→feature class
图5
Input：C:\dem04\Ex1\stream.shp
Output location：C:\dem04\Ex1\New Personal Geodatabase.mdb Output feature class：length
在Arcmap中打开length的属性表，在shape_length上右击 satistics
图6
记下这个sum=612606m 单位为m
在DEM上右击进入layer properties，选择属性source
图7
记下cell size 和 column and rows
<612606/1000）/【<25*25*696*622）/1000000】=2.26千M/平方公里
图8。