ArcGIS_9_教程_第11章_水文分析
第十一章 水文分析
水文分析是DEM 数据应用的一个重要方面。利用DEM 生成的集水流域和水流网络，成为大多数地表水文分析模型的主要输入数据。表面水文分析模型应用于研究与地表水流有关的各种自然现象如洪水水位及泛滥情况，或者划定受污染源影响的地区，以及预测当某一地区的地貌改变时对整个地区将造成的影响等，应用在城市和区域规划、农业及森林、交通道路等许多领域，对地球表面形状的理解也具有十分重要的意义。这些领域需要知道水流怎样流经某一地区，以及这个地区地貌的改变会以什么样的方式影响水流的流动。 基于DEM 的地表水文分析的主要内容是利用水文分析工具提取地表水流径流模型的水流方向、汇流累积量、水流长度、河流网络（包括河流网络的分级等）以及对研究区的流域进行分割等。通过对这些基本水文因子的提取和基本水文分析，可以在DEM 表面之上再现水流的流动过程，最终完成水文分析过程。
本章主要介绍ArcGIS 水文分析模块的应用。ArcGIS 提供的水文分析模块主要用来建立地表水的运动模型，辅助分析地表水流从哪里产生以及要流向何处，再现水流的流动过程。同时，通过水文分析工具的应用，也可以有助于了解排
水系统和地表水流过程的一些基本的概念和关键的过程，以
及怎样通过ArcGIS 水文分析工具从DEM 数据上获取更多的
水文信息。
图11.1 ArcToolBox 中的 水文分析模块
ArcGIS9将水文分析中的地表水流过程集合到
ArcToolbox 里，如图11.1所示。主要包括水流的地表模拟过
程中的水流方向确定、洼地填平、水流累计矩阵的生成、沟
谷网络的生成以及流域的分割等。
本章1至5节主要是依据水文分析中的水文因子的提取
过程对ArcGIS 中的水文分析工具逐一介绍。文中所用的
DEM 数据在光盘中chp11文件夹下的tutor 文件夹里面，每
个计算过程以及每一节所产生的数据存放在tutor 文件夹的
result 文件夹里面，文件名与书中所命名相同，读者可以利用
该数据进行参照联系。第6节主要是提供了三个使用水文分
析工具以及水文分析思想的实例。 11.1 无洼地DEM 生成
DEM 被认为是比较光滑的地形表面的模拟，但是由于内插的原因以及一些真实地形
（如喀斯特地貌）的存在，使得DEM 表面存在着一些凹陷的区域。那么这些区域在进行地表水流模拟时，由于低高程栅格的存在，从而使得在进行水流流向计算时在该区域的得到不合理的或错误的水流方向，因此，在进行水流方向的计算之前，应该首先对原始DEM 数据进行洼地填充，得到无洼地的DEM 。
11.1.1 水流方向提取
水流方向是指水流离开每一个栅格单元时的指向。在ArcGIS 中通过将中心栅格的8个邻域栅格编码，水流方向便可以其中的某一值来确定，栅格
方向编码如图11.2所示。
图11.2 水流流向编码 例如：如果中心栅格的水流流向左边，则其水流方向被赋
值为16。输出的方向值以2的幂值指定是因为存在栅格水流
方向不能确定的情况，此时须将数个方向值相加，这样在后续
处理中从相加结果便可以确定相加时中心栅格的邻域栅格状
况。
水流的流向是通过计算中心栅格与邻域栅格的最大距离权落差来确定。距离权落差是指中心栅格与邻域栅格的高程差除以两栅格间的距离，栅格间的距离与方向有关，如果邻域栅格对中心栅格的方向值为2、8、32、128，则栅格间的距离为2的开平方根，否则距离为1。
ArcGIS 中的水流方向是利用D8算法，也就是最大距离权落差（最大坡降法）来计算水流方向的。具体计算步骤如下：
1. 在ArcMap 中用左键单击ArcToolbox 图标，启动ArcToolbox 。
2. 打开水文分析模块。启动ArcToolbox ，展开Analysis Tools 工具箱，打开hydrology 工
具集。
3. 双击Flow Direction 工具，打开如图
11.3所示水流方向（Flow Direction ）
计算对话框。
（1） I nput surface data 文本框中选择输
入数据dem。
（2） O utput flow direction raster 文本
框中命名计算出来的水流方向文件名
为flowdir，并选择保存路径。
（3） 在Force all edge cells to flow outward(Optional)前的复选框前打
钩，所有在DEM 数据边缘的栅格的水流方向全部是流出DEM 数据区域。默认为不选择。这一步为可选步骤。
图11.3 水流方向Flow Direction 计算对话框
（4） d rop raster 输出。drop raster 是该栅格在其水流方向上与其临近的栅格之间的高
程差与距离的比值，以百分比的形式记录。它反映了在整个区域中最大坡降的分布情况。这一步为可选步骤。
（5） 单击OK，进行水流方向计算过程。
（6） 计算完成后，计算出的水流方向数据如图11.4所示。
11.1.2 洼地计算
洼地区域是水流方向不合理的地方，可以通过水流方向来判断那些地方是洼地，然后再对洼地进行填充。有一点必须清楚的是，并不是所有的洼地区域都是由于数据的误差造成的，有很多洼地区域也是地表形态的真实反映，因此，在进行洼地填充之前，必须计算洼地深度，判断哪些地区是由于数据误差造
成的洼地而哪些地区又是真实的地表形态，
然后在进行洼地填充的过程中，设置合理的
填充阈值。
图11.4 利用Flow Direction 工具计算出来的水流方向图
图11.5 洼地计算对话框
1. 洼地计算
（1） 双击hydrology 工具集中的Sink 工
具，弹出洼地计算对话框，如图11.5
所示。
（2） 在Input flow direction raster 文
本框中，选择水流方向数据flowdir。
（3） 在Output raster 文本框中，选择存放的路径以及重新命名输出文件为sink。
（4） 单击OK 进行洼地计算。结果如图11.6所示，深色的区域是洼地。
图11.6 计算出来的洼地区域
2. 洼地深度计算
（1） 双击hydrology 工具集中的watershed 工具，弹出流域计算对话框，
如图11.7所示，用来计算洼地的贡献区域。
图11.7 洼地贡献区域计算对话窗口（watershed ）
（2） 在Input flow direction raster 文
本框中选择水流方向数据flowdir，
在Input raster or feature pour
point 文本框中输入洼地数据sink，
在pour point field 文本框中选择
value。
（3） 在Output raster 文本框中设置输出
数据的名称为watershsink。
（4） 单击OK，进行洼地贡献区域的计算。
显示结果如图11.8所示。
图11.8 计算出来的洼地贡献区域
（5） 计算每个洼地所形成的贡献区域的最低高程。双击spatial analysis tools 工具箱
中zonal 工具集下的zonal statistic 工具，弹出如图11.9所示的分区统计对话框。
1） 在Input raster or feature zonal data
文本框中，选择洼地贡献区域数据
watershsink ；
2） 在Input value raster 文本框中输入
希望进行统计分析的数据层，现在
需要统计洼地贡献区域的最低高
程，选dem 作为value raster 。
3） 在Output raster 文本框中将输出数
据文件命名为zonalmin ，存放路径
保持不变。
4） 统计类型选择。在统计类型选择的
下拉菜单中有软件所提供的一些
统计类型：分别是在分带区域中统计的每一个分带的平均值（mean ）、
最大值（maximum ）、最小值（minimum ）、分带中的属性值的变化值（rang ）、标
准差（std ）以及总和（sum ）。这里选择最小值作为统计类型。
以上设置完毕之后，单击OK ，完成计算。
图11.9 分区统计对话框
（6） 计算每个洼地贡献区域出口的最低高程即洼地出水口高程。双击spatial analysis
tools 工具箱中zonal 工具集下的zonal fill 工具，弹出如图11.10所示的分区统计对话框。在Input zone raster 文本框中选择watershsink，在Input weight raster
文本框中选择dem，在Output raster 文本框中将输出数据文件名改为zonalmax，然后单击OK，进行运算。
（7） 计算洼地深度。加载Spatial Analyst 模块，点击Spatial Analyst 模块的下拉箭
头，点击raster calculator 菜单工具，如图11.11，在文本框里面输入sinkdep = ( [zonalmax] - [zonalmin])，然后点击evaluate 进行计算。
对于以上（4）、
（5）步的计算，可以利用spatial analysis tools 工具箱中的map algebra 工具集的multi map output 工具进行计算。双击spatial analysis tools 工具箱中的map algebra 工具集下的multi map output 工具，弹出图11.12的对话框。
对于第（4）步，在文本框中输入：E:\chp11\tutor\result\zonalmin = zonalmin
(E:\chp11\tutor\result\watershsink, E:\chp11\tutor \dem)；
对于第（5）步，在文本框中输入：E:\chp11\tutor\result\zonalmax = zonalfill
(E:\chp11\tutor\result \watershsink, E:\chp11\tutor \dem)；
经过以上七步的运算，就可到如图11.13所示的所有洼地贡献区域的洼地深度。通过对研究区地形的分析，可以确定出哪些是由数据误差而产生的洼地，哪些洼地区域又是真实的反映地表形态，从而根据洼地深度来设置合理的填充阈值，使得生成的无洼地DEM 更准确的反映地表形态。
图11.11 洼深度计算
图11.10 洼地贡献区域边缘最低高程计算对话框
图11.12 map algebra 计算窗口
11.1.3 洼地填充
洼地填充是无洼地DEM
生成的最后一个步骤。在通过洼地计算之后，知道了原始的DEM 上是否存在着洼地，如果没有存在着洼地，那么原始DEM 数据就直接可以用来进行以后的河网的生成、流域的分割等。而洼地深度的计算又为在填充洼地时设置填充阈值提供了很好的参考。
图11.13 计算出的洼地深度图
1. 双击hydrology 工具集中的fill 工具，弹出如图11.14所示的洼地填充对话框。
2. 在Input surface raster 文本框中，
选择需要进行洼地填充的原始DEM 数据。
在Ou 3. tput surface raster 文本框中设置
4. 设置一个数
5. 洼地填平计算，计算后的无洼地DEM 见图11.15所示。
图11.14 洼地填充对话框
输出文件名为filldem 。
Z limit ——填充阈值，当值之后，在洼地填充过程中，那些洼
地深度大于阈值的地方将作为真实地
形保留，不予填充；系统默认情况是
不设阈值，也就是所有的洼地区域都将被填平。
单击OK ，进行
图11.15 经过洼地填充生成的无洼地DEM
洼地填充的过程是一个反复的过程。当一个洼地区域被填平之后，这个区域与附近区域再进行洼地计算，可能有会形成新的洼地，所以呢，洼地填充是一个不断反复的过程，直到最后所有的洼地都被填平，新的洼地不再产生为止。因此，当数据量很大时，这个过程会持续一段时间。
11.2 汇流累积量
在地表径流模拟过程中，汇流累积量是基于水流方向数据计算而来的。对每一个栅格来说，其汇流累积量的大小代表着其上游有多少个栅格的水流方向最终汇流经过该栅格，汇流累积的数值越大，该区域越易形成地表径流。由水流方向数据到汇流累积量计算的过程如图11.16所示。
水流方向数据汇流累积数据
图11.16 流水累积量的计算
1. 基于无洼地DEM 的水流方向的计算。计算过程同上一节水流方向的计算一样，使用
的DEM 数据是无洼地DEM 。将生成的水流方向文件命名为fdirfill 。
2. 汇流累积量的计算，如图11.17所示。
在得到水流方向之后，可以利用水流方
向数据来计算汇流累积量。双击
hydrology 工具集中的的fill
accumulation 工具，打开汇流累积量计
算对话框。
图11.17 汇流累积量计算窗口
（1） 在Input flow direction raster 文
本框中，选择由无洼地DEM 生成的水
流方向栅格数据fdirfill。 （2） 在Output accumulation raster 文本
框中，将数据文件名修改为flowacc。
（3） 在Input weight raster 文本框中输入配权数据，配权数据一般是表示降水、土壤
以及植被等对径流影响的因素分布不平衡而形成的，更能详细模拟该区域的地表特征。如果无数据，系统默认为所有的栅格配以相同的权值1，那么计算出来的汇流累积量的数值就代表着该栅格位置流入的栅格数的多少。
（4） 点击OK，完成计算。结果见图11.18所示。
图11.18 通过计算生成的汇流累积量数据
11.3水流长度
水流长度通常是指在地面上一点沿水流方向到其流向起点（终点）间的最大地面距离在水平面上的投影长度。水流长度是水土保持上的重要因子之一，当其他条件相同时，水力侵蚀的强度依据坡的长度来决定，坡面越长，汇聚的流量越大，其侵蚀力就越强，水流长度直接影响地面径流的速度，从而影响对地面土壤的侵蚀力。因此，对于水流长度的提取和分析，在水土保持工作中有很重要的作用。目前水流长度的提取方式主要有两种，一种是顺流计算，一种是朔流计算。顺流计算是计算地面上每一点沿水流方向到该点所在流域出水口最大地面距离的水平投影；朔流计算者是计算地面上每一点沿水流方向到其流向起点间的最大地面距离的水平投影。
对于水流长度的提取，在ArcGIS中，水流长度的提取提供了以上两种方向上的水流长度的提取方式。
1.双击hydrology工具集中的的flow
length工具，弹出水流长度的计算对话
框，如图11.19所示，用来计算水流长
度的大小。
2.在Input flow direction raster文本框中
选择基于无洼地DEM提取出的水流方
向数据fdirfill。
3.在Output raster文本框中选择并命名输
出的水流长度栅格数据文件名称。分别进行顺流计算和朔流计算，输出的数据文件命名为Flowlendown和Flowlenup。图11.19 flow length的计算窗口
4.计算方向提供了两种选择，分别为Downstream（顺流计算）和Upstream（朔流计算）。
5.输入计算配权栅格数据。那么对于flow length来说，Downstream记录着其沿着水流
方向到下游流域出水口中最长距离所流经的栅格数；Upstream则记录着其沿着水流方向到上游栅格的最长的距离的栅格数。
6.当设置完成后，点击OK完成。两种方向计算出的结果如图11.20和图11.21所示。
图11.21逆流方向上的水流长度图11.20 顺流方向上的水流长度
11.4河网的提取
基于DEM的水文分析，其中一个目的就是以DEM数据为基础，经过一些计算和处理而得到地表的水流网络。目前常用的河网提取方法是采用地表径流漫流模型计算：首先是在无洼地DEM上利用最大坡降的方法得到每一个栅格的水流方向；然后利用水流方向栅格数据计算出每一个栅格在水流方向上累积的栅格数，即汇流累积量，所得到的汇流累积量则代表在一个栅格位置上有多少个栅格的水流方向流经该栅格；假设每一个栅格处携带一份水流，那么栅格的汇流累积量则代表着该栅格的水流量。基于上述思想，当汇流量达到一定值的时候，就会产生地表水流，那么所有那些汇流量大于那个临界数值的栅格就是潜在的水流路径，由这些水流路径构成的网络，就是河网。
11.4.1 河网的生成
1.河网的生成是基于汇流累积矩阵的，首先应该计算出研究区域的汇流累积矩阵。汇流
累积矩阵的计算见11.2节，这里用11.2节计算的汇流累计栅格数据flowacc作为基础数据。
2.设定阈值。阈值的设定在河网的提取过程是很重要的，并且直接影响到河网的提取结
果。阈值的设定应遵循科学、合理的原则。首先应该考虑到研究的对象，研究对象中的沟谷的最小级别，不同级别的沟谷所对应的不同的阈值；其次考虑到研究区域的状况，不同的研究区域相同级别的沟谷需要的阈值也是不同的。所以，在设定阈值时，应充分对研究区域和研究对象进行分析，通过不断的实验和利用现有地形图等其它数据辅助检验的方法来确定能满足研究需要并且符合研究区域地形地貌条件的合适的阈值。
3.栅格形式的河网的形成。栅格河网的生成是利用map algebra工具集中的multi map
output工具中的con命令或者setnull命令计算的，他们都是基于栅格进行有条件的查询并将查询结果赋予新的栅格数据中。计算的思想是利用所设定的阈值进行整区域的分析并生成一个新的栅格图层，在生成
栅格数据中将那些汇流量大于设定阈
值的栅格的属性值设定为1，而小于或
等于设定阈值的栅格的属性值设定为
无数据。栅格河网的生成是利用设定一
个河网生成阈值，也可以利用ArcMap
中的Spatial Analysis分析模块下的
Raster Calculator来计算出所有大于设
定阈值的栅格，这些栅格就是河网的潜
在位置。将计算出来的栅格河网命名为图11.22 栅格河网转换成矢量结构对话框
streamnet。
4.栅格河网矢量化。在hydrology工具集中提供了将上一步生成的栅格河网进行矢量化
的工具stream to feature，通过stream to feature就可以得到矢量形式的河网图。栅格数据的矢量化时，选择的是最短的路线绘制成线。双击hydrology工具集中的stream to feature工具，如图11.22。
在Input stream raster文本框中，选择streamnet；在Input flow direction raster文本框中，输入由无洼地计算出来的水流方向数据fdirfill；在Output polyline features 文本框中将输出的数据命名为streamfea。生成的矢量数据如图11.23所示。
图11.23 栅格河网转换成的矢量河网框
11.4.2 stream link的生成
Stream link是记录着河网中的一些节点之间
的连接信息，主要记录着河网的结构信息。如图
11.24所示，Stream link的每条弧段连接着两个作
为出水点或汇合点的结点，或者连接着作为出水
点的结点和河网起始点。
因此通过Stream link的计算，即得到每一个
河网弧段的起始点和终止点。同样，也可以得到
该汇水区域的出水点。这些出水点具有很重要的
弧段结点
水文作用，对于水量、水土流失等研究具有重要
意义。而且，这些出水口点的确定，也为进一步图11.24 Stream link示意图
的流域分割准备了数据。
1. Stream link 基于水流方向数据和栅格河网数据计算，首先在ArcMap 里将水流方向数
据fdirfill 和栅格河网数据streamnet 打开。
2. 双击hydrology 工具集中的stream link 工具，弹出如图11.25所示的stream link 计算的
对话框。在Input stream raster 文本框中选择streamnet ，在Input flow direction raster 文本框中选择fdirfill 。在Output raster 文本框中将输出数据名称设为StreamLink 。然后点击OK 进行运算。
经过Stream link 计算之后，它将栅格河网分成不包含汇合点栅格河网片段，并将片断进行记录，在属性表中除了记录该片段的ID 号之外，还记录着每个片段所包含的栅格数。如图11.26所示。
11.4.3 河网分级的生成
河网分级是对一个线性的河流网络进行分级别的数字标识。在地貌学中，对河流的分级是河网分级方法：
Strahler 分级和Shreve 分级。如图根据河流的流量、形态等因素进行河流的分级。而基于DEM 提取的河网的分支具有一定的水文意义。利用地表径流模拟的思想，不同的级别的河网首先是它们所代表的汇流累积量也不同，级别越高的河网，其汇流累积量也越大，那么在水文研究中，这些河网往往是主流，而那些级别较低的河网则是支流。河网的分级同样可以研究水流的运动、汇流模式，对于水土保持等也具有重要的意义。
在ArcGIS 的水文分析中，提供两种常用的11.27所示，
对于Strahler 分级来说，它是将所有河网弧段中没有支流河网弧段分为第1级，两个1级河网弧段汇流成的河网弧段为第2级，如此下去分别为第3级，第4级，一直到河网出水口。在这种分级中，当且仅当同级别的两条河网弧段汇流成一条河网弧段时，该弧段级别才会增加，对于那些低级弧段汇入高级弧段的情况，高级弧段的级别不会改变，这也是比较常用的一种河网分级方法。对于Shreve 分级而言，其第1级河网的定义与Strahler 分级是相同的，所不同的是以后的分级，两条1级河网弧段汇流而成的河网弧段为2级河网弧段，那么对于以后更高级别的河网弧段，其级别的定义是由其汇入河网弧图11.26 StreamLink 的属性框
图11.25 StreamLink 计算对话框
段的级别之和，如图所示，当一条3级河网弧段和一条4级河网弧段汇流而成的新的河网弧段的级别就是7，那么这种河网分级到最后出水口的位置时，其河网的级别数刚好是该河网中所有的1级河网弧段的个数。
Strahler 分级 Shreve 分级
图11.27 Strahler 分级和Shreve 分级示意图
在ArcGIS 中对河网分级的步骤如下：
1. r 的计算也是基于水流方向数据和栅格河网数
2. k 和stream order 计算出的栅格数据同样可以利用hydrology 工具集中的strea 图11.28 河网的Strahler 分级结果 同Stream link 的计算一样，stream orde 据的基础上，首先在ArcMap 里将水流方向数据fdirfill 和栅格河网数据streamras 打开。
双击hydrology 工具集中的stream order 工具，弹出stream order 计算的对话框。在Input
stream raster 文本框中选择streamnet ，在Input flow direction raster 文本框中选择fdirfill 。分别用Strahler 分级和Shreve 分级对河网进行分级，改输出数据名称分别设为Streamostr 和Streamoshr ，单击OK 完成。计算出的两种河网分级分别如图11.28和图11.29所示。
对于stream lin m to feature 工具将其转化成矢量数据进行进一步的研究和分析。
图11.29 河网的Shreve 分级结果
11.5流域的分割
流域又称集水区域，是指流经其中的水流和其它物质从一个公共的出水口排出从而形成一个集中的排水区域，如图11.30所示。用来描述流域还有例如：流域盆地（basin ）、集水盆地（catchment ）或水流区域（contributing area ）。Watershed 数据显示了区域内每个流域汇水面积的大小。汇水面积是指从某个出水口（或点）流出的河流的总面积。出水口（或点）即流域内水流的出口，是整个流域的最低处。流域间的分界线即为分水岭。
分水线包围的区域称为一条河流或水系的流域，流域分水线所包围的区域面积就是流域面积。即流域是指一条河流或水系的集水区域，河流从这个集水区域获得水量的补给。任何一个天然的河网，都是由大小不等的、各种各样的水道所联合组成的，而每一个水道都有自己的特征，自己的汇水范围，即自己的流域面积，较大的流域往往是由若干较小的流域所联合组成的。
11.5.1 流域盆地的确定
流域盆地是由分水岭分割而成的汇水区域。它通过对水流方向数据的分析确定出所有相互连接并处于同一流域盆地的栅格。流域盆地的确定首先是要确定分析窗口边缘的出水口的位置，也就是说，在进行流域盆地的划分中，所有的流域盆地的出水口均处于分析窗口的边缘。当确定了出水口的位置之后，其流域盆地集水区的确定类同于洼地贡献区域的确定，也就是找出所有流入出水口的上游栅格的位置。
在ArcGIS 中，流域盆地的计算是利用hydrology 工具集中的basin 工具来进行计算的。
1. 双击hydrology 工具集中的basin 工具，打开流域盆地计算的对话框。如图11.31所示。
2. 输入数据为水流方向数据fdirfill ，设置输出数据文件名为basin 。
3. 点击OK 完成。
为了使计算结果更容易理解，可以将上
一节计算出的矢量河网数据在同一个窗口
中打开，进行辅助分析。如图11.32所示。
所有的流域盆地的出口都在研究区域的边
界上，利用流域盆地分析，可以从很大的一
个研究区域中选择感兴趣的流域并将该流
域从整个研究区域分割出来进行单独的分
析。图11.31 流域盆地计算的对话框
图11.32 计算出的流域盆地(线状图形为矢量河网数据)
11.5.2 汇水区出水口的确定
经过上一步得到的流域盆地是一个比较大的流域盆地，在很多的水文分析中，还需要基于更小的流域单元进行分析，那么就需要将这些流域从大的流域中分解出来，那么就需要进行流域的分割。而流域的分割首先是要确定小级别的流域的出水口的位置，小级别的流域出水口的位置可以用spatial analysis tools工具箱下的hydrology工具集中的snap pour point的工具寻找。它的思想是利用一个记录着point点栅格数据，在这个数据层中，那些属性值存在的点作为潜在的出水点，在该点位置上在指定距离内在汇流累积量的数据层上搜索那些具有较高汇流累积量栅格点的位置，这些搜索到的栅格点就是小级别的流域的出水点。也可以利用已有的出水点的矢量数据。
如果没有出水点的栅格或矢量数据,可以用已有河网数据进一步生成的stream link数
据作为汇水区的出水口数据。因为stream link数据中隐含着河网中每一条河网弧段的联结信息，包括弧段的起点和终点等，相对而言，弧段的终点就是该汇水区域的出水口所在位置。
11.5.3 集水流域的生成
对于低级的集水区的生成，可以使用hydrology工具集中的watershed工具生成。其思想如下：先确定一个出水点，也就是该集水区的最低点，然后结合水流方向数据，分析搜索出该出水点上游所有流过该出水口的栅格，直到所有的该集水区的栅格都确定了位置，也就是搜索到流域的边界，分水岭的位置。
1.首先在ArcMap中打开水流方向数据fdirfill和流域出口点数据streamlink，因为集水
流域的计算需要流域出口点的位置数据和水流方向数据。水流方向数据使用11.2节基于无洼地DEM计算出的fdirfill数据。对于流域出口点数据，可以使用11.4节计算出的streamlink数据，因为在它的属性中具有每条河网弧段信息，而这些河网弧段的形成则是有其附近的栅格构成的贡献区域的汇流结果，这个贡献区域也就是该河网弧段的集水区域。
2.双击hydrology工具集中的watershed工具，打开集水区域（贡献区域）计算的对话框。
分别在水流方向数据和出水口数据输入的文本框中选择fdirfill和streamlink数据，在输出数据中将输出的文件名改为watershed。
3.点击OK，进行集水区域的计算，计算结果如图11.33所示，为了更好的表现流域的
分割效果，在此窗口中还将上面计算的流域盆地和矢量河网的数据打开，从结果中可以知道，通过streamlink作为流域的出水口数据所得到的集水区域是每一条河网弧段
11.33
图集水区域的计算结果
集水区域，也就是要研究的最小沟谷的集水区域，它将一个大的流域盆地按照河网弧段将其分为一个个的小的集水盆地。
(完整word版)Arcgis操作第九章水文分析
第九章 水文分析 水文分析是DEM 数据应用的一个重要方面。利用DEM 生成的集水流域和水流网络，成为大多数地表水文分析模型的主要输入数据。表面水文分析模型研究与地表水流有关的各种自然现象例如洪水水位及泛滥情况，划定受污染源影响的地区，预测当某一地区的地貌改变时对整个地区将造成的影响等。 基于DEM 地表水文分析的主要内容是利用水文分析工具提取地表水流径流模型的水流方向、汇流累积量、水流长度、河流网络（包括河流网络的分级等）以及对研究区的流域进行分割等。通过对这些基本水文因子的提取和分析，可再现水 流的流动过程，最终完成水文分析过程。 本章主要介绍ArcGIS 水文分析模块的应用。ArcGIS 提供 的水文分析模块主要用来建立地表水的运动模型，辅助分析地 表水流从哪里产生以及要流向何处，再现水流的流动过程。同 时，通过水文分析工具的应用，有助于了解排水系统和地表水 流过程的一些基本概念和关键过程。 ArcGIS 将水文分析中的地表水流过程集合到ArcToolbox 里，如图11.1所示。主要包括水流的地表模拟过程中的水流 方向确定、洼地填平、水流累计矩阵的生成、沟谷网络的生成 以及流域的分割等。 本章1至5节主要是依据水文分析中的水文因子的提取过 程对ArcGIS 中的水文分析工具逐一介绍。文中所用的DEM 数据在光盘中chp11文件夹下的tutor 文件夹里面，每个计算 过程以及每一节所产生的数据存放在tutor 文件夹的result 文件 夹里面，文件名与书中所命名相同，读者可以利用该数据进行 参照练习。本章最后一节还提供了三个水文分析应用的实例。 9.1 无洼地DEM 生成 DEM 一般被认为是比较光滑的地形表面的模拟，但是由于内插的原因以及一些真实地形（如喀斯特地貌）的存在，使得DEM 表面存在着一些凹陷的区域。这些区域在进行地表水流模拟时，由于低高程栅格的存在，使得在进行水流流向计算时在该区域得到不合理的或错误的水流方向。因此，在进行水流方向的计算之前，应该首先对原始DEM 数据进行洼地填充，得到无洼地的DEM 。 洼地填充的基本过程是先利用水流方向数据计算出DEM 数据中的洼地区域，然后计算出这些的洼地区域的洼地深度，最后以这些洼地深度为参考而设定填充阈值进行洼地填充。 9.1.1 水流方向提取 水流方向是指水流离开每一个栅格单元时的指向。在ArcGIS 中通过 将中心栅格的8个邻域栅格编码，水流方向便可由其中的某一值来确定， 图11.2 水流流向编码 图11.1 ArcToolBox 中的 水文分析模块
空间分析之水文分析
空间分析之水文分析 一、目的与要求： 1.学习目的 水文分析：根据DEM提取河流网络，进行河网分级，计算流水累积量、流向、水流长度、根据指定的流域面积大小自动划分流域。 通过本次学习应达到以下目的： ①理解基于DEM数据进行水文分析的基本原理。 ②掌握利用ArcGIS提供的水文分析工具进行水文分析的基本方法和步骤。 2.学习要求 ①了解水文分析工具 ② DEM的预处理：填洼 ③流向分析 ④计算流水累积量 ⑤计算水流长度 ⑥提取河流网络 ⑦流域分析（盆域、分水岭）
二、水文分析基本操作步骤 1.填充洼地 对原始DEM数据进行洼地填充，得到无洼地的DEM： 在【ArcToolbox】中，双击【SpatialAnalyst工具】→【水文分析】→【填洼】，弹出“填洼”对话框，如下图： 点击确定，得无洼地的DEM【fill_dem】，结果图如下：
2.流向分析 在上一步的基础上进行，在【ArcToolbox】中，双击【SpatialAnalyst 工具】→【水文分析】→【流向】，按下图所示指定各参数： 点击确定，得到无洼地DEM生成的水流方向栅格【Flowdir_dem】，
注意：在ArcGIS中通过将中心栅格的8个邻域栅格编码（D8算法），来确定水流方向。 3.计算汇流累积量 在上一步的基础上进行，在【ArcToolbox】中，双击【SpatialAnalyst 工具】→【水文分析】→【流量】，按下图所示指定各参数：
确定后执行完成得到汇流累积量栅格【flow_acc】,如图： 4.提取河流网络 在上一步的基础上进行，打开【Arctoolbox】，双击【Spatial Analyst 工具】→【地图代数】→【栅格计算器】，在【地图代数表达式】中输入公式：Con(“flow_acc”>800,1)，【输出栅格】指定为:【StreamNet】如图：
ArcGIS之水文分析
ArcGIS教程之DEM水文分析详细图文教程，本教程和之前的两个教程有关联的，数据上是使用上一个教程的结果，步骤相互联系！最后会提供给大家数据和教程的链接！水文分析需要： 1.理解基于DEM数据进行水文分析的基本原理。 2.利用ArcGIS的提供的水文分析工具进行水文分析的基本方法和步骤。 下面开始教程： 工具/原料 ?软件准备：ArcGIS Desktop 10.0---ArcMap(spatial Analyst模块) ?数据准备：DEM（使用由本人前面的教程【ArcGIS地形分析--TIN及DEM 的生成，TIN的显示】中使用的原始数据。 方法/步骤 1.数据基础：无洼地的DEM 在ArcMap中加载 DEM数据，右击DEM图层，点击缩放至图层，显示全部。 2.在【ArcToolbox】中，（要打开扩展模块）执行命令[SpatialAnalyst工 具]——>[水文分析]——> [填洼]，按下图所示指定各参数，其中Z限制——填充阈值，当设置一个值后，在洼地填充过程中，那些洼地深度大于阈值的地方将作为真实地形保留，不予填充；系统默认情况是不设阈值，也就是所有的洼地区域都将被填平。之后点击确定即可。 3.确定后执行结果得到无洼地的DEM数据[Fill_dem1]
4.关键步骤：流向分析 在上一步的基础上进行，在【ArcToolbox】中，执行命令[SpatialAnaly st工具]——>[水文分析]——>[流向]，按下图所示指定各参数： 5.确定后执行完成后得到流向栅格[Flowdir_fill1]，理解代表什么含义！ 6.计算流水累积量 在上一步的基础上进行，在【ArcToolbox】中，执行命令[SpatialAnaly st工具]——>[水文分析]——>[流量]，按下图所示指定各参数： 1.7 确定后执行完成得到流水累积量栅格[flowacc_flow1] 如图： 7.提取河流网络 首先，提取河流网络栅格。 在上一步的基础上进行，打开【Arctoolbox】，运行工具[Spatial Anal yst 工具]——>[地图代数]——>[栅格计算器]，在[地图代数表达式]中输入公式：Con(Flow Accumulation1>800,1)，（这里的Flow Accumulat ion1要以上一步得到的文件名为准,注意是Con,不是con,大写第一个字母，不然出错）如图： [输出栅格]指定为:StreamNet保存路径和文件名任意）
GIS水文分析
实验（五）水文分析 1.实验目的 通过Arctoolbox：水文分析工具位于[Spatial Analyst Tools]>>[Hydrology]之下。 2. 实验数据 DEM 2.实验步骤 ◆打开ArcMap，并打开dem数据。如图所示 1. 无洼地DEM生成 （1）流向分析 方法：在ArcMap中用左键单击ArcToolbox图标，启动ArcToolbox。 打开水文分析模块。启动ArcToolbox，展开Spatial Analyst Tools工具箱，打开hydrology工具集。 双击Flow Direction工具，打开水流方向（Flow Direction）计算对话框。 Input surface data文本框中选择输入数据dem。 Output flow direction raster文本框中命名计算出来的水流方向文件名为flowdir，并选择保存路径。 在Force all edge cells to flow outward(Optional)前的复选框前打钩，所有在DEM数据边缘的栅格的水流方向全部是流出DEM数据区域。默认为不选择。这一步为可选步骤。
2 .洼地计算 ①洼地计算 双击hydrology工具集中的Sink工具，弹出洼地计算对话框. 在Input surface raster文本框中，选择水流方向数据flowdir。 在Output surface raster文本框中，选择存放的路径以及重新命名输出文件为sink。
单击OK计算出洼地提取结果。 ②洼地深度计算 双击hydrology工具集中的watershed工具，弹出流域计算对话框，用来计算洼地的贡献区域。 在Input flow direction raster文本框中选择水流方向数据flowdir，在Input raster or feature pour point文本框中输入洼地数据sink，在pour point field文本框中选择value。 在Output raster文本框中设置输出数据的名称为watershsink。 单击OK，进行洼地贡献区域的计算。
河海大学水文分析与计算课程设计报告定稿版
河海大学水文分析与计算课程设计报告 HUA system office room 【HUA16H-TTMS2A-HUAS8Q8-HUAH1688】
水文分析计算课程设计报告书 学院：水文水资源 专业：水文与水资源工程 学号： 姓名： 指导老师：梁忠民、李国芳 2015年06月12日 南京 目录 1、设计任务 (1) 2、流域概况 (1) 3、资料情况及计算方案拟定 (1) 4、计算步骤及主要成果 (2) 4．1 设计暴雨X p(t)计算 (2) 4．1．1 区域降雨资料检验 (2) 4．1．2 频率分析与设计雨量计算 (3) 4．2计算各种历时同频率雨量X t,P (9) 4．3 选典型放大推求X P (t) (9) 4．4 产汇流计算 (9) 4．4．1 径流划分及稳渗μ值率定 (12) 4．4．2 地表汇流 (17) 4．5 由设计暴雨X P(t)推求Q P(t) (18) 4．5．1 产流计算 (18) 4．5．2 地面汇流 (18)
4．5．3地下汇流计算 (19) 4．5．4 设计洪水过程线 (20) 5、心得体会 (22)
1、设计任务 推求江西良田站设计洪水过程线,本次要求做P 校,即推求Q 0.01%(t)。 2、流域基本概况 良田是赣江的支流站。良田站以上控制的流域面积仅为44.5km 2，属于小流域，如右图所示。年降水均值在1500～1600mm 之内，变差系数Cv 为0.2，即该地区降雨充沛，年际变化小，地处湿润地区。暴雨集中。暴雨多为气旋雨、台风雨，季节为3～8月，暴雨历时为2～3日。 3、资料情况及计算方案拟定 3.1资料情况 设计站（良田）流量资料缺乏，邻近站雨量资料相对充分，具体如表3-1： 表3-1 良田站及邻近地区的实测暴雨系列、历时洪水、特大暴雨资料 3.2 方案拟定 本次课设采用间接法推求设计洪水，即是由推求的设计暴雨， 经过产汇流计算得到设计洪水。示意图如下： 4、设计暴雨XP(t)的计算 4.1 设计暴雨X p (t)计算 4.1.1区域降雨资料检验 站名 实测暴雨流量系列 特大暴雨、历史洪水 良田 75～78 （4年） Q=216m 3 /s ，N=80（转化成X 1日，移置峡江站） 峡江 53～80 （28年） 吉安 36～80 （45年） 桑庄 57～80 （24年） X 1日 寨头 57～80 （24年） 沙港 特大暴雨 X 1日 （移置到寨头站）
水文分析
水文分析：根据DEM提取河流网络，计算流水累积量、流向、根据指定的流域面积大小自动划分流域 水文分分析工具 有两种途径使用水文分析功能： (1)通过Arctoolbox：水文分析工具位于[Spatial Analyst Tools]>>[Hydrology]之下 如果Hydrology 工具 集没有出现，可以选中 某个工具箱后新建一 个工具集[Hydrology]， 然后右键点新建的工 具集，在出现的菜单中 执行[添加]>>[工具] 会出现如右图所示的 对话框，将需要的水文 分析工具添加到上面 新建的工具集中。 (2)另一种方法是添加[Hydrology]工具栏到ArcMap中。 在ArcMap中执行菜单命令：[工具]>>[定制] 命令
点击[从文件添加]按钮 找到esrihydrology_v2.dll 文件 注意：这个文件通常是在ArcGIS的安装路径下，默认的情况是[C:\Program Files\ArcGIS\DeveloperKit\samples\SpatialAnalyst\HydrologicModeling\Visu al_Basic] [Hydrology Modeling]工具条就被加载到ArcMap，在其前面的检查框上打上勾，如下图所示。[Hydrology Modeling]工具条就可以显示在ArcMap中
注意：以下的练习基于Hydrology Modeling 工具 1. 数据基础：无洼地的DEM 在ArcMap中加载DEM数据，执行工具条[Hydrology Modeling]中的菜单命令[ Hydrology ]>>[ Fill Sinks]，在出现的对话框中将[Input Surface]参数指定为“DEM” 确定后得到无洼地的DEM数据：[ Filled Sink1 ]
ARCGIS水文分析
ARCGIS水文分析 水文分析是DEM数据应用的一个币要方式。利用DEM生成的集水流域和水流网络，成为大多数地表水文分析模型的卞要输入数据。表ICI水文分析模型应用十研究与地表水流有关的各种自然现象如洪水水位及泛滥情况，或者一划定受污染源影响的地区，以及预测当某一地区的地貌改变时一对整个地区将造成的影响等，应用在城市和区域规划、农业及森林、交通道路等许多领域，对地球表ICI形状的理解也具有}一分要的b,义。这些领域需要知道水流怎样流经某一地区，以及这个地区地貌的改变会以什么样的方式影响水流的流动。 基十DEM的地表水文分析的卞要内容是利用水文分析土具提取地表水流径流模型的水流方向、汇流祟积量、水流长度、河流网络(包括河流网络的分级等)以及对研究区的流域进行分割等。通过对这些基木水文因子的提取和基木水文分析，可以在DEM表ICI之 上再现水流的流动过程，最终完成水文分析过程。 主要介绍ArcGIS水文分析模块的应用。ArcGIS提供的水文分析模块卞要用来建立地表水的运动模型，辅助分析地表水流从哪里产生以及要流向何处，再现水流的流动过程。同时，通过水文分析土具的应用，也可以有助了解排水系统和地表水流过程的一些基木的概念和关键的过程，以及怎样通过ArcGIS水文分析土具从DEM数据上获取更多的水文信息。 ArcGIS9将水文分析中的地表水流过程集合到ArcToolbox里，卞要包括水流的地表模拟过程中的水流方向确定、汁地填平、水流祟不}一矩阵的生成、沟谷网络 的生成以及流域的分割等。 1.无洼地DEM生成
DEM被认为是比较光滑的地形表n的模拟，但是由十内插的原因以及一些真实地形(如喀斯特地貌)的存在，使得DEM表ICI存在着一些}u}陷的区域。那么这些区域在进行地表水流模拟时一，由十低高程栅格的存在，从而使得在进行水流流向不}一算时一在该区域的得到不合理的或错误的水流方向，因此，在进行水流方向的不}一算之前，应该首先对原始DEM数据进行汁地填充，得到无洼地的DEMO 水流方向是指水流离开何一个栅格单儿时一的指向。在ArcGIS个邻域栅格编码，水流方向便可以其中的某一值来确定，栅格方向编码例如:如果中心栅格的水流流向I，边，则其水流方向被赋中通过将中心栅格的8值为160输出的方向值以2的幂值指定是因为存在栅格水流 方向不能确定的情况，此时一须将数个方向值相加，这样在后续处理中从相加结果便可以确定相加时一中心栅格的邻域栅格状己。 1.2水流流向编码 水流的流向是通过不}一算中心栅格与邻域栅格的最大距离权落差来确定。距离权落差是指中心栅格与邻域栅格的高程差除以两栅格间的距离，栅格间的距离与方向有关，如果邻域栅格对中心栅格的方向值为2, 8, 32, 128，则栅格间的距离为2的开平方根，否则距离为1。 1.1.2洼地计算 注地区域是水流方向不合理的地方，可以通过水流方向来判断那些地方是注地，然后再对注地进行填充。有一点必须清楚的是，并不是所有的注地区域都是由十数 据的误差造成的，有很多洼地区域也是地表形态的真实反映，因此，在进行洼地填充之前，必须计算 注地深度，判断哪些地区是由十数据误差造成的注地而哪些地区又是真实的地表形态，然后在进行注地填充的过程中，设置合理的
水文分析
遥感与地理信息系统上机报告 班级：地化21202 序号：15 姓名：成绩： 一、实验题目：水文分析 二、实验目的 1. 了解ArcGIs,Maplnfo的基本功能； 2. 利用水文分析工具提取地表水流径流模型的水流方向、汇流累积量、水 流长度、河流网络以及对研究区的流域进行分割等。 三、实验方法与结果： 转成矢量河流数据（Stream To Feature）Archydro模型上机实习内容； 1.数据准备：数字高程模型（DEM）， 首先导入已知DEM的矢量数据，利用ArcGIS的ArcToolBox中Data Management Tools/Raster/CLIP工具，载出所需流域DEM数据。 水文流域DEM数据图1 2.流向计算 利用原始DEM数据，采用Flow Direction工具计算，得到每个表格的流
向。 3.洼地计算 基于原始DEM的坡向，计算洼地。
经过求洼后的洼地分布图 4.填洼计算 根据Sink扫描找出洼地，用Fill工具将洼地点的高程值设为与相邻点的最小高程值，这样迭代直到填平所有的洼地。填洼后形成了新的经过修正无凹陷DEM。 填洼后无凹陷DEM图 5.重新计算流向数据利用FILL后的新的DEM重新计算Flow Direction。
填洼后的流向计算 6.汇流累积量计算 6.1 利用无凹陷的dem求坡向，得到坡向分布数据。 6.2 利用flow accumulation命令计算出每个格网上淤累积汇流数，越是 上游的格网累积量越小；越处于下游累积数越大。 汇流累积量计算
7.水流长度 水流长度指地面上一点沿水流方向到流向起点（或终点）间的最大地面距离在水平面上的投影长度。 （1）在arctoobox中选择[spatial analyst 工具]/[水文分析]/[水流长度]，打开水流长度工具； （2）[输入栅格流向数据]为fdirfill；在[输出栅格]中指定保存路径及名称；（3）[侧向方向]：选择downstream或upstream； （4）[输入权重栅格数据]； （5）单击[确定]，完成操作。
ArcGIS之水文分析
ArcGIS之水文分析
ArcGIS教程之DEM水文分析详细图文教程，本教程和之前的两个教程有关联的，数据上是使用上一个教程的结果，步骤相互联系！最后会提供给大家数据和教程的链接！水文分析需要： 1.理解基于DEM数据进行水文分析的基本原理。 2.利用ArcGIS的提供的水文分析工具进行水文分析的基本方法和步骤。 下面开始教程： 工具/原料 ?软件准备：ArcGIS Desktop 10.0---ArcMap(spatial Analyst模块) ?数据准备：DEM（使用由本人前面的教程【ArcGIS地形分析--TIN及DEM 的生成，TIN的显示】中使用的原始数据。 方法/步骤 1.数据基础：无洼地的DEM 在ArcMap中加载 DEM数据，右击DEM图层，点击缩放至图层，显示全部。
2.在【ArcToolbox】中，（要打开扩展模块）执行命令[SpatialAnalyst工 具]——>[水文分析]——> [填洼]，按下图所示指定各参数，其中Z限制——填充阈值，当设置一个值后，在洼地填充过程中，那些洼地深度大于阈值的地方将作为真实地形保留，不予填充；系统默认情况是不设阈值，也就是所有的洼地区域都将被填平。之后点击确定即可。 3.确定后执行结果得到无洼地的DEM数据[Fill_dem1]
4.关键步骤：流向分析 在上一步的基础上进行，在【ArcToolbox】中，执行命令[SpatialAnaly st工具]——>[水文分析]——>[流向]，按下图所示指定各参数： 5.确定后执行完成后得到流向栅格[Flowdir_fill1]，理解代表什么含义！
水文分析与计算
水文分析与计算 不同工程要求估算的水文设计特征值不尽相同。桥梁工程要求估算所在河段可能出现的设计最高水位和最大流量,以便合理决定桥梁的高程和跨度;防洪工程为权衡下游和自身的安全、经济和风险，要求估算工程未来运行时期可能遇到的各种稀遇的洪水；灌溉、发电、供水、航运等工程需要知道所在河流可能提供的水量和水能蕴藏量，以确定灌溉面积、发电量、城市或工矿企业供水量和航运发展规模。工程的运行时期可长达几十至几百年，不可能象水文预报那样给出该时期内某一水文特征值出现的具体时间和大小，而是用水文统计的方法，估算在该时期中可能出现的某一设计标准的水文特征值。 一般说，运用水文统计方法所依据的样本很少，抽样误差较大，往往不能满足生产需要。因此，不能单纯根据工程所在地点的水文资料进行计算，还必须对计算过程和计算结果进行充分的合理性分析，才能较可靠地求得工程所在地的设计水文数据。因此，也常称水文计算为水文分析与计算。 一、设计年径流计算 即估算符合设计标准的通过河流某一指定断面的全年和各时段的径流量及其月旬分配，为水资源开发利用的水利规划和工程设计提供科学依据。计算主要内容包括：①各种
设计标准的年最大设计洪峰流量和不同时段设计洪量;②符合设计要求的洪水过程线;③当梯级水库或单一水库下游有防洪要求时，拟定一种或几种满足设计要求的设计洪水的地区组成；④年内不同时期（如某些月份、或丰水期、枯水期和施工期等）的设计洪水。 二、设计洪水计算 即计算符合某一地点指定的防洪设计标准的洪水数值，为防洪规划或防洪工程设计提供可靠的水文数据。 计算的主要内容有：①各种历时的设计地点的雨量或流域平均面雨量；②它们的时程分配和地区分布;③大型工程和重要的中小型工程,还要求估算指定流域的可能最大暴雨，供推算可能最大洪水之用。 三、设计暴雨计算 并根据设计暴雨计算结果，推求相应的设计洪水和涝水。算主要内容有：确定某一设计标准的各年输沙量及其年内分配，以估计水库库容减少情况和工程寿命；估算水库和它的上下游河道冲淤变化，为水工建筑物布设和水库运用方式的确定提供依据。例如，通过合理布设排沙底孔和规定水库运用方式，有助于利用异重流排沙（见河流泥沙、水库淤积）。
如何使用ArcGIS进行水文分析(完整版)
如何使用ArcGIS 进行水文分析 对于做水利的朋友来说有时候需要进行水文的分析，今天给大家分享一下如何通过ArcGIS 进行水文分析，材料可以通过水经注万能地图下载器进行下载。工具/ 原料 水经注万能地图下载器ArcGIS 方法/ 步骤 1. 打开水经注万能地图下载器，框选上需要进行水文分析的地方并下载（图1） 图1 2.下载完成后会自动导出成tif 格式的高程DEM数据，将其加载到ArcGIS 内（图2）。【说明】：此处下载生成的tif 格式的图片即为大家常说的DEM数据，直接加载到ArcGIS 内即可使用。
图2 3. 点击“自定义”→“扩展模块”（图3），在弹出的对话框中将“空间分析” Spatial Analyst ）工具勾选上（图4）。 图3
图4 4. 在ArcToolbox 中点击“ Spatial Analyst 工具”→“水文分析”→“填洼” （图5），在弹出的“填洼”对话框中按图 6 进行设置。其中Z限制——填充阈值，当设置一个值后，在洼地填充过程中，那些洼地深度大于阈值的地方将作为真实地形保留，不予填充；系统默认情况是不设阈值，也就是所有的洼地区域都将被填平。【特别说明】：为了保证最终分析成功，在最终的结果之前，所有输出的数据都默认保存名称和路径，这就需要我们记清楚哪个名称是对应的哪个成果，后面会有用。
图5 图6 5. 填洼完成后得到名称为 “ Fill_tif3 的填洼成 果， 在ArcToolbox 工具中点击Spatial Analyst 工具”→“水文分析”→“流向”图7 ），在弹出的“流 向” 对话框中进行如图8 所示的设置，将上一步得到 的 Fill_tif3 ”填洼数据作为
ArcGIS空间分析报告——找出某药材地生长区域
课程：ArcGIS空间分析 实验目的：利用GIS空间分析方法，结合等高线及温度和降水数据，在充分分析某药材的生长习性的情况下，找到其生长区域，从而能够更好的保护该药材的生长环境。 数据来源：本实验所采用的数据均来自ArcGIS地理信息系统空间分析实习教程，数据有：山区等高线数据contour.shp 和山区观测点采集的年平均温度和年总降水数据climate.txt. 实验要求：根据所给条件，确定某区域适合种植这种药材的范围，求出适合种植的面积。 (1)这种药材一般生长在沟谷两侧较近的区域（不超过 500m） (2)这种药材喜阳 (3)生长气候环境为年平均温度10度-12度 (4)年总降水量为550-680mm 实验流程：利用该山区等高线数据生成DEM，基于DEM进行水文分析，提取沟谷网络；基于DEM提取坡向数据，重分类划分阴阳坡。 利用观测点采集的年平均温度和年总降水数据分别进行表面内插，生成年平均温度栅格数据和年总降水栅格数据。提取年平均温度10度-12度的区域和年总降水为
550mm-680mm的区域。 综合叠加分析满足上述4个条件的区域，得到适合该药材生长的区域，并制作专题图，计算该适合区域的面积。实验步骤： 1.利用等高线，构建DEM。首先打开ArcMap，加载等高线数据，在ArcToolbox中，选择【3D Analyst】|【Tin 管理】|【创建Tin】工具，打开工具对话框，生成tin。空间参考依然导入contour相同的坐标系统。 2.将Tin转换成格网DEM，以便于进行表面分析和与其他数据的叠加分析。选择【3D Analyst工具】|【转换】|【由Tin转出】|【Tin转栅格】工具，打开工具对话框。
ArcGIS水文分析
实验四、水文分析-DEM应用 专业年级：地信071姓名：王媛媛学号：06407024 一、实验目的与要求 1.实验目的 水文分析：根据DEM提取河流网络，进行河网分级，计算流水累积量、流向、水流长度、根据指定的流域面积大小自动划分流域。 通过本实验应达到以下目的： 1理解基于DEM数据进行水文分析的基本原理。 ②掌握利用ArcGIS提供的水文分析工具进行水文分析的基本方法和步骤。 2.实验要求 ①了解水文分析工具 2DEM的预处理：填洼与削峰 3流向分析 4计算流水累积量 5计算水流长度（流程） 6提取河流网络 7流域分析 二、实验原理 水文分析基本步骤
①无洼地的DEM DEM被认为是比较光滑的地形表面的模拟，但是由于内插的原因以及一些真实地形（如采石场或喀斯特地貌）的存在，使得DEM表面存在着一些凹陷的区域。这些区域在进行地表水流模拟时，由于低高程栅格的存在，从而使得在进行水流流向计算时得到不合理的或错误的水流方向，因此，在进行水流方向的计算之前，应该首先对原始DEM数据进行洼地填充，得到无洼地的DEM。 ②关键步骤：流向分析―――流向分析原理 水流方向是指水流离开每一个栅格单元时的指向。在ArcGIS中通过将中心栅格的8个邻域栅格编码（D8算法），来确定水流方向。 方向约定如左图：共有八个方向，分别是2的n次方。 水流的流向是通过计算中心栅格与邻域栅格的最大距离权落差来确定的。距离权落差是指中心栅格与邻域栅格的高程差除以两栅格间的距离，栅格间的距离与方向有关，如果邻域栅格对中心栅格的方向值为2、8、32、128，则栅格间的距离为SQRT(2)≈1.414，否则距离为1。如果高程差为正值，则为流出；负值则为流入。 ③汇流累积量 在地表径流模拟过程中，汇流累积量是基于水流方向数据计算而来的。对每一个栅格来说，其汇流累积量的大小代表着其上游有多少个栅格的水流方向最终汇流经过该栅格，汇流累积的数值越大，该区域越易形成地表径流。图有些地方的计算不是太理解 ④水流长度（流程） 水流长度通常是指在地面上一点沿水流方向到其流向起点（终点）间的最大地面距离在水平面上的投影长度。目前水流长度的提取方式主要有两种，一种是顺流计算（Downstream），一种是朔流计算（Upstream）。顺流计算是计算地面上每一点沿水流方向到该点所在流域出水口最大地面距离的水平投影；朔流计算者是计算地面上每一点沿水流方向到其流向起点间的最大地面距离的水平
水文分析知识
ArcGIS 水文分析介绍（详细内容见附件） 和客户交流ArcGIS的水文分析功能，搜集了些资料，顺便整理了这篇文档。从水文分析的基本原理、算法入手，分别介绍了Hydrology和 ArcHydro，并以汤国安编写ArcGIS水文分析作为实例，演示如何在ArcGIS中基于DEM执行水文分析。 1基本原理 DEM是数字高程模型的英文简称(Digital Elevation Mode)，是流域地形、地物识别的重要原始资料。自20世纪60年代以来,在利用数字高程模型DEM提取流域水文特征,模拟地表水文过程方面,国内外都开展了大量的研究。 1.1基于DEM进行流域分析的原理 从DEM提取流域特征，一个良好的流域结构模式是确定算法的前提和关键。1967年ShreveL¨描述的流域结构模式一直被后来的水文学者所引用．并设计了一些成熟的算法。Shreve使用一个具有一个根的树状图来描述流域结构(如图1所示)。在这个结构中，主要包括两个部分，一部分是结点集，一部分是界线集。沟谷结合点和沟谷源点共同组成一个沟谷结点集。所有的沟谷段组成沟谷段集，形成一个沟谷网络；所有的分水线段组成分水线段集，形成一个分水线网络；沟谷段集和分水线段集共同组成界线集。 沟谷网络中的每一段沟谷都有一个汇流区域，这些区域由流域分水线集来控制。外部沟谷段有一个外部汇流区．而内部沟谷段有两个内部汇水区，分布在内部沟谷段的两侧。整个流域被分割成一个个子流域．每个子流域好象是树状图上的一片“叶子”。 Shreve的树状图流域结构模型是简单明确的．虽然沟谷网络的结点模型和线模型与在栅格DEM中用于表示沟谷结点和沟谷线的栅格点和栅格链之间存在着拓扑不一致性。但它给出了沟谷网络、分水线网络和子汇流区的定义，明确表达了它们之间的相关关系，成为设计流域特征提取技术的基础。 1.1常用算法 流向判定建立在3×3 的DEM 栅格网的基础上，其方法有单流向法和多流向法之分，但单
ARCgiss水文提取2
利用ArcGIS水文分析工具提取河网的操作 DEM包含有多种信息，ArcToolBox提供了利用DEM提取河网的方法，但是操作比较烦琐（帮助可参看Hydrologic analysis sample applications），今天结合我自己的使用将心得写出来与大家分享。提取河网首先要有栅格DEM，可以利用等高线数据转换获得。在此基础上，要经过洼地填平、水流方向计算、水流积聚计算和河网矢量转化这几个不步骤。 1．洼地填平 DEM洼地（水流积聚地）有真是洼地和数据精度不够高所造成的洼地。洼地填平的主要作用是避免DEM的精度不够高所产生的（假的）水流积聚地。洼地填平使用ArctoolBox- >Spatial Analysis Tools->Hydrology->Fill工具。 2．水流方向计算 水流方向计算就可以使用上一步所生成的DEM为源数据了（如果使用未经洼地填平处理的数据，可能会造成精度下降）。这里主要使用ArctoolBox->Spatial Analysis Tools->Flow Direction 工具。输入的DEM采用第一步的Fill1_exam1 3．水流积聚计算 这里主要使用ArctoolBox->Spatial Analysis Tools->Flow Accumulation工具流向。栅格数据就是第二步所获得的数据（FlowDir_fill1）。可以看到，生成的水流积聚栅格已经可以看到所产生的河网了。现在所需要做的就是把这些河网栅格提取出来。可以把产生的河网的支流的象素值作为阀值来提取河网栅格。 4．提取河网栅格 使用spatial analyst中的栅格计算器，将所有大于河网栅格阀值的象素全部提取出来。至于这个阀值是多少因具体情况而定。通常是要大于积聚计算后得到栅格的最低河流象素值。这里采用的是500这个值。最后生成只有0、1值的栅格数据。其中1表示是河网，0是非河网。 5．生成河网矢量 这里主要使用ArctoolBox->Spatial Analysis Tools->Stream to Feature工具.Input Stream raster 为第四步只有0、1值的河网栅格。流向栅格使用第二步所生成的栅格数据。 6．矢量河网处理 由于Stream to Feature工具.将所有栅格象素均转为矢量线段。所以要进行处理，方法是利用属性查询的方法把所有GRID_CODE为1的全部选择出来。导出就得到了由dem所生成的河网矢量。 7．处理结果 最后,得到的河网如下图所示。但是由于是栅格转换而来。生成的河网并不是连续的矢量。可以根据需要做简单的处理。
基于ArcGIS的水文分析原理
1基本原理 DEM是数字高程模型的英文简称(Digital Elevation Mode)，是流域地形、地物识别的重要原始资料。自20世纪60年代以来,在利用数字高程模型DEM提取流域水文特征,模拟地表水文过程方面,国内外都开展了大量的研究。 1.1基于DEM进行流域分析的原理 从DEM提取流域特征，一个良好的流域结构模式是确定算法的前提和关键。1967年ShreveL¨描述的流域结构模式一直被后来的水文学者所引用．并设计了一些成熟的算法。 Shreve使用一个具有一个根的树状图来描述流域结构(如图1所示)。在这个结构中，主要包括两个部分，一部分是结点集，一部分是界线集。沟谷结合点和沟谷源点共同组成一个沟谷结点集。所有的沟谷段组成沟谷段集，形成一个沟谷网络；所有的分水线段组成分水线段集，形成一个分水线网络；沟谷段集和分水线段集共同组成界线集。 图1 流域结构模式图 (a) (b) (c) (f) (d) (e) (g) (h) 沟谷网络中的每一段沟谷都有一个汇流区域，这些区域由流域分水线集来控制。外部沟谷段有一个外部汇流区．而内部沟谷段有两个内部汇水区，分布在内部沟谷段的两侧。整个流域被分割成一个个子流域．每个子流域好象是树状图上的一片“叶子”。 Shreve的树状图流域结构模型是简单明确的．虽然沟谷网络的结点模型和线模型与在栅格DEM中用于表示沟谷结点和沟谷线的栅格点和栅格链之间存在着拓扑不一致性。但它给出了沟谷网络、分水线网络和子汇流区的定义，明确表达了它们之间的相关关系，成为设计流域特征提取技术的基础。
1.2常用算法 流向判定建立在3×3 的DEM 栅格网的基础上，其方法有单流向法和多流向法之分，但单流向法因其确定简单、应用方便而应用广泛。 1.2.1单流向法 单流向法假定一个栅格中的水流只从一个方向流出栅格，然后根据栅格高程判断水流方向。目前应用的单流向法是D8法。此外，还有Rho8 方法、DEMON 法、Lea 法和D∞法等。最常用的是D8 法：假设单个栅格中的水流只能流入与之相邻的8 个栅格中。它用最陡坡度法来确定水流的方向，即在3×3 的DEM 栅格上，计算中心栅格与各相邻栅格间的距离权落差（即栅格中心点落差除以栅格中心点之间的距离），取距离权落差最大的栅格为中心栅格的流出栅格。 所谓最陡坡度法的原理是假设地表不透水，降雨均匀．那么流域单元上的水流总是流向最低的地方“窗口滑动指以计算单元为中心，组合其相邻的若干个单元形成一个窗口”，以“窗口”为计算基本元素，推及整个DEM，求取最终结果。 目前应用最广泛的是基于流向分析和汇流分析的流域特征提取技术。Jenson and Domingue (1988)设计了应用该技术的典型算法，该算法包括3个过程：流向分析，汇流分析和流域特征提取。 1)流向分析：以数值表示每个单元的流向。数字变化范围是1～255。其中1：东；2：东南；4南；8：西南；16：西；32：西北；64：北；128：东北。除上述数值之外的其它值代表流向不确定，这是由DEM中洼地”和“平地”现象所造成的。所谓“洼地”即某个单元的高程值小于任何其所有相邻单元的高程。这种现象是由于当河谷的宽度小于单元的宽度时，由于单元的高程值是其所覆盖地区的平均高程，较低的河谷高度拉低了该单元的高程。这种现象往往出现在流域的上游。“平地指相邻的8个单元具有相同的高程，与测量精度、DEM单元尺寸或该地区地形有关。这两种现象在DEM 中相当普遍，Jenson and Domingue 在流向分析之前，将DEM进行填充；将“洼地”变成“平地”，再通过一套复杂的迭代算法确定“平地”流向。流向分析过程如图所示。
ArcGIS_9_教程_第11章_水文分析
第十一章 水文分析 水文分析是DEM 数据应用的一个重要方面。利用DEM 生成的集水流域和水流网络，成为大多数地表水文分析模型的主要输入数据。表面水文分析模型应用于研究与地表水流有关的各种自然现象如洪水水位及泛滥情况，或者划定受污染源影响的地区，以及预测当某一地区的地貌改变时对整个地区将造成的影响等，应用在城市和区域规划、农业及森林、交通道路等许多领域，对地球表面形状的理解也具有十分重要的意义。这些领域需要知道水流怎样流经某一地区，以及这个地区地貌的改变会以什么样的方式影响水流的流动。 基于DEM 的地表水文分析的主要内容是利用水文分析工具提取地表水流径流模型的水流方向、汇流累积量、水流长度、河流网络（包括河流网络的分级等）以及对研究区的流域进行分割等。通过对这些基本水文因子的提取和基本水文分析，可以在DEM 表面之上再现水流的流动过程，最终完成水文分析过程。 本章主要介绍ArcGIS 水文分析模块的应用。ArcGIS 提供的水文分析模块主要用来建立地表水的运动模型，辅助分析地表水流从哪里产生以及要流向何处，再现水流的流动过程。同时，通过水文分析工具的应用，也可以有助于了解排 水系统和地表水流过程的一些基本的概念和关键的过程，以 及怎样通过ArcGIS 水文分析工具从DEM 数据上获取更多的 水文信息。 图11.1 ArcToolBox 中的 水文分析模块 ArcGIS9将水文分析中的地表水流过程集合到 ArcToolbox 里，如图11.1所示。主要包括水流的地表模拟过 程中的水流方向确定、洼地填平、水流累计矩阵的生成、沟 谷网络的生成以及流域的分割等。 本章1至5节主要是依据水文分析中的水文因子的提取 过程对ArcGIS 中的水文分析工具逐一介绍。文中所用的 DEM 数据在光盘中chp11文件夹下的tutor 文件夹里面，每 个计算过程以及每一节所产生的数据存放在tutor 文件夹的 result 文件夹里面，文件名与书中所命名相同，读者可以利用 该数据进行参照联系。第6节主要是提供了三个使用水文分 析工具以及水文分析思想的实例。 11.1 无洼地DEM 生成 DEM 被认为是比较光滑的地形表面的模拟，但是由于内插的原因以及一些真实地形
基于DEM的ArcGIS水文分析—河网和流域的提取
基于DEM的ArcGIS水文分析 —河网和流域的提取 一、实验背景 水文分析是DEM 数据应用的一个重要方面。而利用DEM生成的集水流域和水流网络，成为大多数地表水文分析模型的主要输入数据。表面水文分析模型研究与地表水流有关的各种自然现象例如洪水水位及泛滥情况，划定受污染源影响的地区，预测当某一地区的地貌改变时对整个地区将造成的影响等。 二、实验目的 通过本实验，使读者理解基于DEM数据进行水文分析的基本原理，掌握利用ArcGIS 提供的水文分析工具进行水文分析的基本方法和步骤，并利用DEM数据提取出河网及流域。 三、实验数据 某地区栅格数据DEM，数据来源于随书光盘（…\Chp9\Ex2）。 四、实验要求 根据DEM利用水文分析工具提取地表水流径流模型的水流方向、汇流累积量、水流长度、河流网络（包括河流网络的分级等）以及对研究区的流域进行分割等。
五、实验流程图 六、实验内容及步骤 1．无洼地DEM生成 DEM 是比较光滑的地形表面模型，但由于DEM 误差以及一些真实地形或特殊地形的影响，使得DEM 表面存在一些凹陷的区域。 在进行水流方向计算时，由于这些区域的存在，往往得到不合理的甚至错误的水流方向。因此，在进行水流方向的计算之前，应该首先对原始DEM 数据进行洼地填充，得到无洼地的DEM。
洼地填充的基本过程是先利用水流方向数据计算出DEM 数据中的洼地区域，并计算洼地深度，然后，依据这些洼地深度设定填充阈值进行洼地填充。 1．1 水流方向的提取 水流的流向是通过计算中心格网与邻域格网的最大距离权落差来确定。对于每一格网的水流方向指水流离开此网格的指向。在ARCGIS 中，通过对中心栅格的1、2、4、8、16、32、64、128 等8个邻域栅格编码，中心栅格的水流方向便可有其中的某一值来确定。例如，若中心栅格的水流流向左边，则水流方向赋值16。 流向的生成是个自动的过程，可能要等一段自时间，运算的时间跟电脑性能和DEM图的精度与大小有关.。 方法是利用ArcToolbox\Spatial Analysis Tools\ Hydrology \Flow Direction，生成方向水流流向图：若从DEM中作出来的流向分析的最大数值为128则不需要填洼，否则需要填挖。
水文分析与计算
水文水利计算是工程水文的重要组成部分，分为水文计算和水利计算。 根本任务 水文计算：分析水文要素变化规律，为水利工程的建设提供未来水文情势预估。 水利计算：拟定并选择经济合理和安区可靠的工程设计方案‘规划设计参数和调度允许方式。 一、水文计算的主要研究方法 ?设计标准 ?概率预估（PMP/PMF） ?研究进展 ?基于风险理论的防洪标准研究 ?气候变化和人类活动对设计成果的影响 ?不确定性新理论 二、水利计算的主要研究方法 ?水量调节 ?洪水调节 ?枯水调节 ?水能调节 一、水文过程的随机特性 水文现象同时存在“确定性过程”和“随机性过程”。 确定性因素和随机因素共同作用下的模型，统称为“随机模型”。 二、纯随机模型对水文过程的适用性 采用随机方法解决水文计算问题时，依据的是概率统计理论中的纯随机模型，即假设所研究的水文变量是独立随机地抽自同一客观总体，而这个总体是通过概率分布函数（或概率密度函数）来描述的。 水文频率分析计算的任务，就是根据水文变量的样本对总体进行统计（如参数估计、推求制定标准的设计值等）和推断（如假设检验、推求臵信限等）。 一、洪水资料的选样
指导思想：保证纯随机模型的适用性，独立同分布。 洪水三要素：洪峰、洪量、洪水过程。 选样方法： （1）年最大值法； （2）年多次法； （3）超定量法； （4）超大值法。 二、洪水资料的审查和分析 1.可靠性审查 2.一致性审查 3.代表新审查 三、洪水资料的插补延长 1. 根据上下游测站的洪水特征相关关系进行插补延长 ?点绘相关图； ?设计站洪水由上游几个干支流测站的洪水组成，应错时叠加； ?因受洪水展开和区间来水影响，考虑能反映上述影响因素的参数； 三、洪水资料的插补延长 1. 根据上下游测站的洪水特征相关关系进行插补延长 若设计断面资料短，甚至无资料，则无法直接建立相关关系，需要修正，其做法如下： （1）两者集水面积之差小于3%，中间无天然或认为分滞洪，可直接移用； （2）面积之差大于3%，但不大于10%～20%，且暴雨分布均匀，用面积进行修正； （3）若在上下游均有参证站满足要求，则可进行内插。 三、洪水资料的插补延长 2. 利用本站峰量关系进行插补延长 通常根据调查到的历史洪峰或由相关关系法求得缺测年份洪峰流量，利用峰量关系可以推求洪水总量。 3. 利用暴雨径流关系进行插补延长 通过扣损汇流计算，推求相应于一次暴雨过程的洪水过程线，进而计算其洪峰和洪量。 4.根据相邻河流测站的洪水特征值进行插补延长 一、洪水调查的意义 河流所发生的特大洪水包括实测特大洪水与历史洪水两种。所谓历史洪水是指过去已经发生，但未通过水文站测到的特大洪水。