实验9、水文分析-DEM应用

1、实验九、水文分析-dem 应用一、实验目的水文分析使用dem 数据派生其它水文特征：提取河流网络、自动划分流域。这些是描述某一地区水文特征的重要因素。通过本实验应达到以下目的：a)理解基于 dem 数据进行水文分析的基本原理（重要）。b)掌握利用arcgis的提供的水文分析工具进行水文分析的基本方法和步骤。二、实验准备预备知识：数据基础：无洼地的dem ，被较高高程区域围绕的洼地是进行水文分析的一大障碍，因此在确定水流方向以前，必须先将洼地填充。有些洼地是在dem 生成过程中带来的数据错误，但另外一些却表示了真实的地形如采石场或岩洞等。通过填充洼地(fill sinks) 得到无洼地的dem

2、关键步骤：流向分析流向分析原理1641281632842方向约定如左图：共有八个方向，分别是2 的 n 次方。水流的流向是通过计算中心栅格与邻域栅格的最大距离权落差来确定的。距离权落差是指中心栅格与邻域栅格的高程差除以两栅格间的距离，栅格间的距离与方向有关，如果邻域栅格对中心栅格的方向值为2、8、32、128，则栅格间的距离为 sqrt(2) 1.414 ，否则距离为 1。如果高程差为正值，则为流出。为负值则为流入。基本概念dem 数据存储高程数据的栅格称为digital elevation model (dem). 每个像元只有一个高程值. 基本概念排水系统水文分析步骤：实验数据：某地区 1

3、:5 万 dem 数据软件准备： 水文分析工具 (hydrology modeling) 的使用有两种途径使用水文分析功能：(1)通过 arctoolbox ：水文分析工具位于spatial analyst toolshydrology之下如果hydrology 工具集没有出现， 可以选中某个工具箱后新建一个工具集hydrology，然后右键点新建的工具集，在出现的菜单中执行添加 工具 会出现如右图所示的对话框， 将需要的水文分析工具添加到上面新建的工具集中。(2)另一种方法是添加hydrology 工具栏到arcmap 中。在 arcmap 中执行菜单命令：工具 定制 命令点击从文件添加 按

4、钮找到 esrihydrology_v2.dll 文件注意 ：这个文件通常是在arcgis 的安装路径下，默认的情况是c:program filesarcgisdeveloperkitsamplesspatialanalysthydrologicmodelingvisual_basic hydrology modeling工具条就被加载到arcmap ，在其前面的检查框上打上勾，如下图所示。 hydrology modeling工具条就可以显示在arcmap 中注意：以下的练习基于hydrology modeling 工具三、实验内容及步骤1. 数据基础：无洼地的dem 在 arcmap中加载

5、dem数据，执行工具条hydrology modeling 中的菜单命令 hydrology fill sinks，在出现的对话框中将input surface 参数指定为“ dem ”确定后得到无洼地的dem 数据： filled sink1 2. 关键步骤：流向分析在 上 一 步 的 基 础 上 进 行 ， 执 行 工 具 条hydrology modeling 中 的 菜 单 命 令 hydrology flow direction ，在出现的对话框中将input surface参数指定为“ filled sink1”确定后得到流向栅格 flow direction1 ，了解流向栅格单元

6、的数值表示的含义是什么3. 计算流水累积量在 上 一 步 的 基 础 上 进 行 ， 执 行 工 具 条hydrology modeling 中 的 菜 单 命 令 hydrology flow accumulation ，在出现的对话框中将direction raster参数指定为“ flow direction1 ”确定后得到流水累积量栅格flow accumulation1 4. 提取河流网络(1)提取河流网络栅格：在上一步的基础上进行，打开arctoolbox ，运行工具 spatial analyst toolsmap algebra 单输出地图代数 在 地图代数表达式中输入公式：c

7、on (flow accumulation1800,1) 输出栅格 指定为： streamnet 说明：通过此操作将流水累积量栅格flow accumulation1 中栅格单元值（流水累积量）大于800 的栅格赋值为1，从而得到河流网络栅格streamnet 得到的的河流网络栅格：streamnet 关闭除 streamnet之外的其它图层(2)提取河流网络矢量数据在上一步的基础上进行，执行工具条hydrology modeling 中的菜单命令 hydrology stream network as feature ，在出现的对话框中将direction raster参 数 指 定 为“

8、flow direction1 ” ， accumulation raster 参 数 指 定 为“ flow accumulation1 ” ，minimum number of cells for a stream参数指定为1000 确定后得到河流网络矢量数据(3)平滑处理河流网络打开 编辑器 工具栏，执行工具栏中的命令编辑器 开始编辑 ，确保目标图层为河流网络图层shape1, 通过打开 shape1 属性表，并选择属性表的所有行选择图层shape1中的所有要素，也可以通过要素选择按钮选择图层中所有要素执行编辑器 工具栏中的命令 编辑器 更多的编辑工具 高级编辑 打开工具条： 高级编辑，

9、点击其上的 平滑按钮（下图中前头所指） ：在平滑 处理对话框中输入参数允许最大偏移 ：3 得到平滑后的河流网络矢量图层，执行命令: 编辑器 停止编辑 ，保存所做修改。比较平滑处理后的数据与没有进行处理过的数据5流域分析在上一步的基础上进行，执行工具条hydrology modeling 中的菜单命令 hydrology watershed ，在出现的对话框中将direction raster参数指定为“ flow direction1 ” ， accumulation raster 参数指定为 “flow accumulation1 ” ， minimum number of cells for a stream 参数指定为2000 确定后得到流域栅格watershed1 打开 空间分析 工具栏，执行命令:空间分析 转换 栅格到要素 将流域栅格转换成为矢量图层，按下图所示指定参数：得到矢量数据：watershed.shp 设置图层 watershed属性 符号 后得到类似上图的效果重复以上操作步骤并将minimum number of cells for a stream参数指定为10000，看看结果是否有所不同。6. 其它应用 -降水