ArcGIS专题操作之-水文分析

地理信息系统试验水文分析-DEM应用一、试验目的与要求1.试验目的水文分析：依据DEM提取河流网络，进展河网分级，计算流水累积量、流向、水流长度、依据指定的流域面积大小自动划分流域。通过本试验应到达以下目的：①理解基于DEM②把握利用ArcGIS2.试验要求①了解水文分析工具②DEM③流向分析④计算流水累积量⑤计算水流长度〔流程〕⑥提取河流网络⑦流域分析二、试验原理水文分析根本步骤①无凹地的DEM1地理信息系统试验DEM〔如采石场或喀斯特地貌〕DEM外表存在着一些凹陷的区域。这些区域在进展地表水流模拟时，由于低高程栅格的存在，从而使得在进展水流流向计算时得到不合理的或错误的水流方向，因此，在进行水流方向的计算之前，应当首先对原始DEM数据进展凹地填充，得到无凹地的DEM。②关键步骤：流向分析―――流向分析原理ArcGIS中通过将中心栅格的8个邻域栅格编码〔D8算法〕，来确定水流方向。方向商定如左图：共有八个方向，分别是2的n次方。水流的流向是通过计算中心栅格与邻域栅格的最大距离权落差来确定的。距离权落差是指中心栅格与邻域栅格的高程差除以两栅格间的距离，栅格间的距离与方向有关，假设邻域栅格对中心栅格的方向值为2、8、32、128，则栅格间的距离为SQRT(2)≈1.414，否则距离为1。假设高程差为正值，则为流出；负值则为流入。③汇流累积量在地表径流模拟过程中，汇流累积量是基于水流方向数据计算而来的。对每一个栅格来说，其汇流累积量的大小代表着其上游有多少个栅格的水流方向最终汇流经过该栅格，汇流累积的数值越大，该区域越易形成地表径流。图有些地方的计算不是太理解④水流长度〔流程〕水流长度通常是指在地面上一点沿水流方向到其流向起点〔终点〕间的最大地面距离在水平面上的投影长度。目前水流长度的提取方式主要有两种，一种是顺流计算〔Downstream〕，一种是朔流计算〔Upstream〕。顺流计算是计算地面上每一点沿水流方向到该点所在流域出水口最大地面距离的水平投影；朔流计算者是计算地面上每一点沿水流方向到其流向起点间的最大地面距离的水平投影。2地理信息系统试验⑤提取河流网络DEM上利用最大坡降的方法得到每一个栅格的水流方向；然后利用水流方向栅格数据计算出每一个栅格在水流方向上累积的栅格数，即汇流累积量，所得到的汇流累积量则代表在一个栅格位置上有多少个栅格的水流方向流经该栅格；假设每一个栅格处携带一份水流，那么栅格的汇流累积量则代表着该栅格的水流量。基于上述思想，当汇流量到达肯定值的时候，就会产生地表水流，那么全部那些汇流量大于那个临界数值的栅格就是潜在的水流路径，由这些水流路径构成的网络，就是河网。⑥流域分析流域又称集水区域，是指流经其中的水流和其它物质从一个公共的出水口排出从而形成一个集中的排水区域。用来描述流域还有例如：流域盆地〔basin〕长江流域、集水盆地〔catchment〕或水流区域〔contributingarea〕。Watershed数据显示了区域内每个流域汇水面积的大小。汇水面积是指从某个出水口〔或点〕流出的河流的总面积。出水口〔或点〕即流域内水流的出口，是整个流域的最低处。流域间的分界限即为分水岭，流域分水线所包围的区域面积就是流域面积。三、试验数据及软件预备试验数据：某地区1:5万DEM数据。软件预备：启动ArcToolbox，开放[SpatialAnalysisTools]工具箱，翻开[Hydrology]工具集。四、试验内容及步骤无凹地的DEM生成原始DEM数据提取水流方向执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[FlowDirectionFlowDirectionForcealledgecellstoflowoutward(Optional3地理信息系统试验选框前打钩，则全部在DEM数据边缘的栅格的水流方向全部流出DEM数据区域〔默认为不选择〕。“dropraster”是该栅格在其水流方向上与其接近的栅格之间的高程差与距离的比值，以百分比的形式记录，它反映了在整个区域中最大坡降的分布状况〔可选步骤〕。凹地的计算执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[Sink局部截图凹地填充Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[Fill对De进展填充[Fill]“Zlimi”为填充阈值，当设置一个数值之后，在凹地填充过程中，那些凹地深度大于阈值的地方将作为真实地形保存，不予填充；系统默认状况是不设阈值，也就是全部4地理信息系统试验的凹地区域都将被填平。DEM〔对填充后的dem〕执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[FlowDirection理解流向栅格单元的数值表示的含义。ValueCount计算汇流累积量〔对流向数据进展分析〕执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[FlowAccumulation在[FillInputweightraste”文本框中输入权重数据，权重数据一般是表示降水、土壤以及植被等对径流影响的因素分布不平衡而形成的，更能具体模拟该区域的地表特征。假设很多据，系统默认为全部的栅格配以一样的权值1，那么计算出来的汇流累积量的数值就代表着该栅格位置流入的栅格数的多少。计算水流长度〔流程〕执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[FlowLength计算方向分别为Downstream〔顺流计算Upstream〔朔流计算“Inputweightraster”文本框中输入权重数据，Downstream记录着其沿着水流方向到下游流域出水口中最长距离所流经的栅格数；Upstream5地理信息系统试验提取河流网络提取河流网络栅格①设定阈值。阈值的设定在河网的提取过程是很重要的，并且直接影响到河网的提取结果。阈值的设定应遵循科学、合理的原则。首先应当考虑到争论的对象，争论对象中的沟谷的最小级别，不同级别的沟谷所对应的不同的阈值；其次考虑到争论区域的状况，不同的争论区域一样级别的沟谷需要的阈值也是不同的。②提取栅格河流网络方式一：翻开Arctoolbox，运行工具[SpatialAnalystTools]>>[MapAlgebra]>>[单输出地图代数]。在[地图代数表达式]中输入公式：con(FlowAccumulation1>800,1)，[输出栅格]指定为：StreamNet。计算的思想：利用所设定的栅格阈值进展条件查询并将查询结果赐予的栅格[FlowAccumulation1]中栅格单元值〔流水累积量〕大于800的栅格赋值为1，而小于或等于设定阈值的栅格属性值赋为很多据，从而得到河流网络栅格[StreamNet]。方式二：在ArcMap中，点击空间分析模块>>计算器,在计算器中输入表达式：con(FlowAccumulation1>800,1)。关闭除[Streamnet]之外的其它图层。栅格河流网络矢量化执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[streamtofeature]伪沟谷的删除和平滑处理河流网络①伪沟谷的删除6地理信息系统试验由于基于DEM的河网的提取是承受最大坡降的方法，那么在平地区域〔例如谷底等〕水流方向是随机的，就很简洁生成平行状的河流等错误外形〔伪沟谷[Editor辑剔除的。争论区域边缘的很短的沟谷也需进展删除。②平滑处理河流网络利用[Editor器]>>[更多的编辑工具]>>[高级编辑]翻开工具条：[高级编辑]，点击其上的[平滑]按钮〔以下图中箭头所指：在[平滑]处理对话框中输入参数[允许最大偏移]：3，得到平滑后的河流网络矢量图层。或执行[Arctoolbox]>>[DataManagementTools]>>[Generalization]>>[SmoothLine]比较平滑处理后的数据与没有进展处理过的数据。streamlinkStreamlink〔河网的构造信息。StreamlinkStreamlink的计算，即得到每一个河网弧段的起始点和终止点。同样，也可以得到该汇水区域〔流域〕的出水口。执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[Streamlink]经过计算，它将栅格河网在集合点栅格处分割成河网片段，并将片断进展记录，在属性表中除了记录该片段的ID[Streamlink]结果可利用[streamtofeature]转换为矢量数据。河网分级ArcGISStrahler分级和ShreveStrahler117地理信息系统试验234级中，当且仅当同级别的两条河网弧段汇流成一条河网弧段时，该弧段级别才会增加，对于那些低级弧段汇入高级弧段的状况，高级弧段的级别不会转变，这也是比较常用的一种河网分级方法。对Shreve1Strahler河网弧段，其级别的定义是由其汇入河网弧段的级别之和，执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[StreamOrder[StreamOrder]结果可利用[streamtofeature]转换为矢量数据。流域分析〔1〕流域盆地确实定流域盆地是由分水岭分割而成的汇水区域。它通过对水流方向数据的分析确定出全部相互连接并处于同一流域盆地的栅格。流域盆地确实定首先是要确定分析窗口边缘的出水口的位置，也就是说，在进展流域盆地的划分中，全部的流域盆地的出水口均处于分析窗口的边缘。当确定了出水口的位置之后，也就是找出全部流入出水口的上游栅格的位置。执行[Arctoolbox]>>[SpatialAnalystTools]>>[Hydrology]>>[Basin为了使计算结果更简洁理解，可以将上述计算出的矢量河网数据在同一个窗口中翻开，进展关心分析。全部的流域盆地的出口都在争论区域的边界上。利用流域盆地分析，可以从很大的一个争论区域中选择感兴趣的流域并将该流域从整个争论区域分割出来进展单独的分析。翻开[空间分析]工具栏，执行命令:[空间分析]>>[转换]>>[栅格到要素]将流域栅格转换成为矢量图层，并进展符号设置。〔2〕集水流域的生成经过上一步得到的流域盆地是一个比较大的流域盆地，在很多的水文分析中，还需要基于更小8地理信息系统试验的流域单元进展分析，那么就需要进展流域的分割。而流域的分割首先是要确定小级别的流域的出水口的位置。①汇水区出水口确实定小级别的流域出水口的位置，可以用[Hydrology]工具集中的[SnapPourPoint]工具查找。思想:利用point点栅格数据查找潜在的出水点，并赋属性值。在该点位置上在指定距离内在汇流累积量的数据层上搜寻那些具有较高汇流累积量栅格点的位置，这些搜寻到的栅格点就是小级别的流域的出水点。也可以利用已有的出水点的矢量数据。假设没有出水点的栅格或矢量数据,可以用上述生成的streamlink数据作为汇水区的出水口数streamlink对而言，弧段的终点就是该汇水区域的出水口所在位置。②集水流域的生成低级的集水区的生成，可以使用[Hydrology]工具集中的[Watershed]工具生成。其思想如下：先确定一个出水点，也就是该集水区的最低点，然后结合水流方向数据，分析搜寻出该出水点上游全部流过该出水口的栅格，直到全部的该集水区的栅格都确定了位置，也就是搜寻到流域的边界——分水岭的位置。通过streamlink作为流域的出水口数据所得到的集水区域是每一条河网弧段集水区域，也就是要争论的最小沟谷的集水区域，它将一个大的流域盆地依据河网弧段将其分为一个个的小的集水盆