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

1、实验四实验名称沟谷网络的提取及沟壑密度的计算实验目的学习并掌握利用水文分析工具提取沟谷网络的原理和方法。实验步骤1沟谷网络的提取 <1）在ArcMap中加载原始DEM，打开Hydrology工具集;+惡+苓+智+ &* S+令-II电Spat i ：=_ Analyst ToolCondi ti Dens i tyDi stanceallGr oundw at er Hydrology严 FillFlow Accumulat沪 |FL艸Flow Length产 SiBk沪 Fwur PcidStreaiTi LinkS tr e am Or derStre to FeatR Wf

2、tttrsh«dlInterpol ati on<2）运用Flow Direction工具在原始 DEM上提取水流方向vflowdir ）；Spatial an alyst tool hydrology flow direct ion 在弹出的flow direction窗口中In put : C:demO4Ex1dem图 2flowdir<3）运用Sink工具进行洼地的计算vsink ）；Spatial analyst tool f hydrology f sinkInput: C:dem04Ex1flowdirOutputC:dem04Ex1Si nk_flowdir

3、1由图像可知，且流向并不遵从 D8算法，值域为1-192，证明存在洼地。<4）计算得知原始DEM上有洼地，运用Fill工具进行洼地填充vfilldem ）；Spatial analyst tool f hydrology f fillOutput C:dem04Ex1Fill_dem1图 3fill_dem1<5)运用Flow Direction工具在无洼地 DEM上提取水流方向vflowdirfill );Spatial an alyst tool f hydrology f flow direct ionIn put: C:dem04Ex1fill_dem1Output C:d

4、em04Ex1Flowdir_fill2此时的流向就满足D8算法，分别是1、2、4、8、16、32、64、128<6)运用Flow Accumulation工具进行汇流累积数据的计算 vflowacc );Spatial analyst tool f hydrology f f。 ow accumulationOutput C:dem04Ex1FlowDir_flow1图 4flowacc_flow1<7）根据实际情况设置汇流累积阈值 <100）,生成栅格河网vstreamne）; 进入 raster calculator，输入 con（flowacc_flow1 >

5、100 , 1>即可得到汇流累计值超过100的河网得至U calculation<8）运用Stream to FeatureX具，进行栅格河网矢量化；Spatial an alyst tool hydrology stream to featureIn put: calculati onIn put flow_direct ion raster : flowdir_fill2Output: stream.shp<9）伪沟谷的删除 <在Editor模块下进行）。 加载上等高线和地貌晕眩进行比较，将那些伪沟谷删除。2.沟壑密度的计算在catalog中新建一个个人数据库在 e

6、x1 文件夹上右击 New personal geodatabase再在新建的 personal geodatabas上右击importfeature classInput Featuresl5： VdemO4Exl straam. shpOutput LocationC: dem041 IIew PerEonal Geodatabase. mdbputpiiit FEWtuireName|length|Expression (optional)IaFi皂Id Info (optici门剖)Fi el dJSameMewFi el ilWame| Vi .Spll tPo.AKirAKCIDT

7、RUENONEGRID CODEGRID CODETBUEBONEFRawjiorEFROW.NODETRUEHONETO NODETO MODETRUENONEOKCancelSitviroiunexits. Show Help In put: C:demO4Ex1stream.shpOutput location： C:demO4Ex1New Personal Geodatabase.mdbOutput feature class length在 Arcmap中打开length的属性表，在 shape_length上右击 satistics图6记下这个sum=612606m单位为m在DEM上右击进入layer prope