《城市地理信息系统 水文分析.doc》由会员分享，可在线阅读，更多相关《城市地理信息系统 水文分析.doc（7页珍藏版）》请在金锄头文库上搜索。
1、城市地理信息系统实验实验四、水文分析-DEM 应用 专业年级： 姓名： 学号：一、实验目的与要求1.实验目的水文分析：根据DEM提取河流网络，进行河网分级，计算流水累积量、流向、水流长度、根据指定的流域面积大小自动划分流域。通过本实验应达到以下目的： 理解基于DEM数据进行水文分析的基本原理。 掌握利用ArcGIS提供的水文分析工具进行水文分析的基本方法和步骤。2.实验要求 了解水文分析工具 DEM的预处理：填洼与削峰 流向分析 计算流水累积量 计算水流长度（流程） 提取河流网络 流域分析二、实验原理 水文分析基本步骤 无洼地的DEMDEM被认为是比较光滑的地形表面的模拟，但是由于内插的原因以
2、及一些真实地形（如采石场或喀斯特地貌）的存在，使得DEM表面存在着一些凹陷的区域。这些区域在进行地表水流模拟时，由于低高程栅格的存在，从而使得在进行水流流向计算时得到不合理的或错误的水流方向，因此，在进行水流方向的计算之前，应该首先对原始DEM数据进行洼地填充，得到无洼地的DEM。 关键步骤：流向分析流向分析原理水流方向是指水流离开每一个栅格单元时的指向。在ArcGIS中通过将中心栅格的8个邻域栅格编码（D8算法），来确定水流方向。 方向约定如左图：共有八个方向，分别是2 的n 次方。水流的流向是通过计算中心栅格与邻域栅格的最大距离权落差来确定的。距离权落差是指中心栅格与邻域栅格的高程差除以两
3、栅格间的距离，栅格间的距离与方向有关，如果邻域栅格对中心栅格的方向值为2、8、32、128，则栅格间的距离为 SQRT(2)1.414 ，否则距离为1。如果高程差为正值，则为流出；负值则为流入。 汇流累积量 在地表径流模拟过程中，汇流累积量是基于水流方向数据计算而来的。对每一个栅格来说，其汇流累积量的大小代表着其上游有多少个栅格的水流方向最终汇流经过该栅格，汇流累积的数值越大，该区域越易形成地表径流。 水流长度（流程）水流长度通常是指在地面上一点沿水流方向到其流向起点（终点）间的最大地面距离在水平面上的投影长度。目前水流长度的提取方式主要有两种，一种是顺流计算（Downstream），一种是朔
4、流计算（Upstream）。顺流计算是计算地面上每一点沿水流方向到该点所在流域出水口最大地面距离的水平投影；朔流计算者是计算地面上每一点沿水流方向到其流向起点间的最大地面距离的水平投影。 提取河流网络目前常用的河网提取方法是采用地表径流漫流模型计算：首先是在无洼地DEM上利用最大坡降的方法得到每一个栅格的水流方向；然后利用水流方向栅格数据计算出每一个栅格在水流方向上累积的栅格数，即汇流累积量，所得到的汇流累积量则代表在一个栅格位置上有多少个栅格的水流方向流经该栅格；假设每一个栅格处携带一份水流，那么栅格的汇流累积量则代表着该栅格的水流量。基于上述思想，当汇流量达到一定值的时候，就会产生地表水流
5、，那么所有那些汇流量大于那个临界数值的栅格就是潜在的水流路径，由这些水流路径构成的网络，就是河网。 流域分析流域又称集水区域，是指流经其中的水流和其它物质从一个公共的出水口排出从而形成一个集中的排水区域。用来描述流域还有例如：流域盆地（basin）、集水盆地（catchment）或水流区域（contributing area）。Watershed 数据显示了区域内每个流域汇水面积的大小。汇水面积是指从某个出水口（或点）流出的河流的总面积。出水口（或点）即流域内水流的出口，是整个流域的最低处。流域间的分界线即为分水岭，流域分水线所包围的区域面积就是流域面积。三、实验数据及软件准备1.实验数据：某
6、地区1:5 万DEM数据。2.软件准备：启动ArcToolbox，展开Spatial Analysis Tools工具箱，打开Hydrology工具集。四、实验内容及步骤1. 无洼地的DEM生成（1）原始DEM数据提取水流方向执行 Arctoolbox Spatial Analyst ToolsHydrology Flow Direction 。在 Flow Direction 对话框中，“Force all edge cells to flow outward(Optional)”的复选框前打钩，则所有在DEM数据边缘的栅格的水流方向全部流出DEM数据区域（默认为不选择）。 “drop ra
7、ster”是该栅格在其水流方向上与其临近的栅格之间的高程差与距离的比值，以百分比的形式记录，它反映了在整个区域中最大坡降的分布情况（可选步骤）。 （2）洼地的计算执行 Arctoolbox Spatial Analyst ToolsHydrology Sink 。（3）洼地填充执行 Arctoolbox Spatial Analyst ToolsHydrology Fill。在Fill对话框中，“Z limit”为填充阈值，当设置一个数值之后，在洼地填充过程中，那些洼地深度大于阈值的地方将作为真实地形保留，不予填充；系统默认情况是不设阈值，也就是所有的洼地区域都将被填平。 2. 关键步骤：无洼
8、地DEM流向分析执行 Arctoolbox Spatial Analyst ToolsHydrology Flow Direction 。理解流向栅格单元的数值表示的含义。3. 计算汇流累积量执行 Arctoolbox Spatial Analyst ToolsHydrology Flow Accumulation 。在 Fill 对话框中，“Input weight raster”文本框中输入权重数据，权重数据一般是表示降水、土壤以及植被等对径流影响的因素分布不平衡而形成的，更能详细模拟该区域的地表特征。如果无数据，系统默认为所有的栅格配以相同的权值1，那么计算出来的汇流累积量的数值就代表着
9、该栅格位置流入的栅格数的多少。 4. 计算水流长度（流程）执行 Arctoolbox Spatial Analyst ToolsHydrology Flow Length 。计算方向分别为Downstream（顺流计算）和Upstream（朔流计算）。 在“Input weight raster”文本框中输入权重数据，Downstream记录着其沿着水流方向到下游流域出水口中最长距离所流经的栅格数；Upstream则记录着其沿着水流方向到上游栅格的最长的距离的栅格数。 5. 提取河流网络(1) 提取河流网络栅格 设定阈值。阈值的设定在河网的提取过程是很重要的，并且直接影响到河网的提取结果。阈值
10、的设定应遵循科学、合理的原则。首先应该考虑到研究的对象，研究对象中的沟谷的最小级别，不同级别的沟谷所对应的不同的阈值；其次考虑到研究区域的状况，不同的研究区域相同级别的沟谷需要的阈值也是不同的。 提取栅格河流网络方式一：打开Arctoolbox，运行工具Spatial Analyst ToolsMap Algebra单输出地图代数。在 地图代数表达式中输入公式：con (Flow Accumulation1800,1)，输出栅格指定为：StreamNet。计算的思想：利用所设定的栅格阈值进行条件查询并将查询结果赋予新的栅格数据中。通过此操作将流水累积量栅格Flow Accumulation1中
11、栅格单元值（流水累积量）大于800 的栅格赋值为1，而小于或等于设定阈值的栅格属性值赋为无数据，从而得到河流网络栅格StreamNet。方式二：在Arc Map中，点击空间分析模块计算器,在计算器中输入表达式：con (Flow Accumulation1800,1)。 关闭除Streamnet之外的其它图层。(2) 栅格河流网络矢量化执行 Arctoolbox Spatial Analyst ToolsHydrology stream to feature (3) 伪沟谷的删除和平滑处理河流网络 伪沟谷的删除由于基于DEM的河网的提取是采用最大坡降的方法，那么在平地区域（例如谷底等）水流方向
12、是随机的，就很容易生成平行状的河流等错误形状（伪沟谷），需要利用Editor工具栏进行手工编辑剔除的。研究区域边缘的很短的沟谷也需进行删除。 平滑处理河流网络利用Editor工具栏进行编辑，选择图层中所有要素，执行编辑器工具栏中的命令编辑器更多的编辑工具高级编辑打开工具条：高级编辑，点击其上的平滑按钮（下图中箭头所指）：在平滑处理对话框中输入参数允许最大偏移：3，得到平滑后的河流网络矢量图层。或执行 Arctoolbox Data Management ToolsGeneralization Smooth Line 比较平滑处理后的数据与没有进行处理过的数据。（4）stream link的生成
13、 Stream link记录着河网中的一些节点之间的连接信息（河网的结构信息）。Stream link的每条弧段连接着两个作为出水点或汇合点的结点，或者连接着作为出水点的结点和河网起始点。 通过Stream link的计算，即得到每一个河网弧段的起始点和终止点。同样，也可以得到该汇水区域（流域）的出水口。执行 Arctoolbox Spatial Analyst ToolsHydrology Stream link 。经过计算，它将栅格河网在汇合点栅格处分割成河网片段，并将片断进行记录，在属性表中除了记录该片段的ID号之外，还记录着每个片段所包含的栅格数。 Stream link 结果可利用
14、stream to feature 转换为矢量数据。（5） 河网分级在ArcGIS的水文分析中，提供两种常用的河网分级方法：Strahler分级和Shreve分级。对于Strahler分级来说，它将所有河网弧段中没有支流的河网弧段分为第1级，两个1级河网弧段汇流成的河网弧段为第2级，如此下去分别为第3级，第4级，一直到河网出水口。在这种分级中，当且仅当同级别的两条河网弧段汇流成一条河网弧段时，该弧段级别才会增加，对于那些低级弧段汇入高级弧段的情况，高级弧段的级别不会改变，这也是比较常用的一种河网分级方法。对于Shreve分级而言，其第1级河网的定义与Strahler分级是相同的，所不同的是以后
15、更高级别的河网弧段，其级别的定义是由其汇入河网弧段的级别之和，执行 Arctoolbox Spatial Analyst ToolsHydrology Stream Order 。 Stream Order 结果可利用 stream to feature 转换为矢量数据。6流域分析（1）流域盆地的确定 流域盆地是由分水岭分割而成的汇水区域。它通过对水流方向数据的分析确定出所有相互连接并处于同一流域盆地的栅格。流域盆地的确定首先是要确定分析窗口边缘的出水口的位置，也就是说，在进行流域盆地的划分中，所有的流域盆地的出水口均处于分析窗口的边缘。当确定了出水口的位置之后，也就是找出所有流入出水口的上游栅格的位置。 执行 Ar
