实验六ArcGISHydrology剖析.docx
- 文档编号:13377945
- 上传时间:2023-06-13
- 格式:DOCX
- 页数:25
- 大小:713.29KB
实验六ArcGISHydrology剖析.docx
《实验六ArcGISHydrology剖析.docx》由会员分享,可在线阅读,更多相关《实验六ArcGISHydrology剖析.docx(25页珍藏版)》请在冰点文库上搜索。
实验六ArcGISHydrology剖析
实验五ArcGISHydrology
一、实验准备
1.水文分析处理步骤:
2.处理数据:
前期实验得到的tingrid.imgDEM数据
3.处理工具:
ArcToolBox工具箱中的spatialanalysttools下的Hydrology
二、无洼地DEM生成
1.水流方向提取
1)在ArcMap中用左键单击ArcToolbox图标,启动ArcToolbox。
2)打开水文分析模块。
启动ArcToolbox,展开AnalysisTools工具箱,打开hydrology工具集。
3)双击FlowDirection工具,打开如下图所示水流方向(FlowDirection)计算对话框。
a)Inputsurfacedata文本框中选择输入数据tingrid。
b)Outputflowdirectionraster文本框中命名计算出来的水流方向文件名
为FlowDir_ting3,并选择保存路径。
c)在Forcealledgecellstoflowoutward(Optional)前的复选框前打钩,所有在DEM数据边缘的栅格的水流方向全部是流出DEM数据区域。
默认为不选择。
d)dropraster输出。
dropraster是该栅格在其水流方向上与其临近的栅格之间的高程差与距离的比值,以百分比的形式记录。
它反映了在整个区域中最大坡降的分布情况。
这一步为可选步骤。
e)单击OK,进行水流方向计算过程。
f)计算完成后,计算出的水流方向数据如下图所示:
2.洼地计算
1)洼地计算
a)双击hydrology工具集中的Sink工具,弹出洼地计算对话框,如下图
b)在Inputflowdirectionraster文本框中,选择水流方向数据flowdir_ting3。
c)在Outputraster文本框中,选择存放的路径以及重新命名输出文件为sink_ting3。
d)单击OK进行洼地计算。
结果如下图所示,深色的区域是洼地。
2)洼地深度计算
a)双击hydrology工具集中的watershed工具,弹出流域计算对话框,如下图所示,用来计算洼地的贡献区域。
b)在Inputflowdirectionraster文本框中选择水流方向数据flowdir_ting3,在Inputrasterorfeaturepourpoint文本框中输入洼地数据sink_ting3,在pourpointfield文本框中选择value。
c)在Outputraster文本框中设置输出数据的名称为Watersh_Flow1。
d)单击OK,进行洼地贡献区域的计算。
显示结果如下图所示。
e)计算每个洼地所形成的贡献区域的最低高程。
双击spatialanalysistools工具箱中zonal工具集下的zonalstatistic工具,弹出如下图所示的分区统计对话框。
在Inputrasterorfeaturezonaldata文本框中,选择洼地贡献区域数据Watersh_Flow1;
在Inputvalueraster文本框中输入希望进行统计分析的数据层,现在需要统计洼地贡献区域的最低高程,选tingrid1作为valueraster。
在Outputraster文本框中将输出数据文件命名为zonalmin,存放路径保持不变。
统计类型选择。
在统计类型选择的下拉菜单中有软件所提供的一些统计类型:
分别是在分带区域中统计的每一个分带的平均值(mean)、最大值(maximum)、最小值(minimum)、分带中的属性值的变化值(rang)、标准差(std)以及总和(sum)。
这里选择最小值作为统计类型。
以上设置完毕之后,单击OK,完成计算。
f)计算每个洼地贡献区域出口的最低高程即洼地出水口高程。
双击spatialanalysistools工具箱中zonal工具集下的zonalfill工具,弹出如图所示的分区统计对话框。
在Inputzoneraster文本框中选择Watersh_Flow1,在Inputweightraster文本框中选择tingrid1,在Outputraster文本框中将输出数据文件名改为zonalmax,然后单击OK,进行运算。
g)计算洼地深度。
加载SpatialAnalyst模块,点击SpatialAnalyst模块的下拉箭头,点击rastercalculator菜单工具,如下图,在文本框里面输入sinkdep=([zonalmax]-[zonalmin]),然后点击evaluate进行计算。
经过以上七步的运算,就可到如下图所示的所有洼地贡献区域的洼地深度。
通过对研究区地形的分析,可以确定出哪些是由数据误差而产生的洼地,哪些洼地区域又是真实的反映地表形态,从而根据洼地深度来设置合理的填充阈值,使得生成的无洼地DEM更准确的反映地表形态。
3.洼地填充
1)双击hydrology工具集中的fill工具,弹出如图所示的洼地填充对话框。
2)在Inputsurfaceraster文本框中,选择需要进行洼地填充的原始DEM数据。
3)在Outputsurfaceraster文本框中设置输出文件名为Fill-tingrid3。
4)Zlimit——填充阈值,当设置一个数之后,在洼地填充过程中,那些洼地深度大于阈值的地方将作为真实地形保留,不予填充;系统默认情况是不设阈值,也就是所有的洼地区域都将被填平。
5)单击OK,进行洼地填平计算,计算后的无洼地DEM见下图所示。
三、汇流累积量
1.基于无洼地DEM的水流方向的计算。
计算方法同上有洼地的DEM水流方向的计算一样。
将生成的水流方向文件命名为fdirfill。
2.汇流累积量的计算,如图所示。
在得到水流方向之后,可以利用水流方向数据来计算汇流累积量。
双击hydrology工具集中的fillaccumulation工具,打开汇流累积量计算对话框。
1)在Inputflowdirectionraster文本框中,选择由无洼地DEM生成的水流方向栅格数据fdirfill。
2)在Outputaccumulationraster文本框中,将数据文件名修改为flowacc。
3)在Inputweightraster文本框中输入配权数据,配权数据一般是表示降水、土壤以及植被等对径流影响的因素分布不平衡而形成的,更能详细模拟该区域的地表特征。
如果无数据,系统默认为所有的栅格配以相同的权值1,那么计算出来的汇流累积量的数值就代表着该栅格位置流入的栅格数的多少。
4)点击OK,完成计算。
结果见下图所示。
四、水流长度
1.双击hydrology工具集中的flowlength工具,弹出水流长度的计算对话框,如图所示,用来计算水流长度的大小。
2.在Inputflowdirectionraster文本框中选择基于无洼地DEM提取出的水流方向数据fdirfill。
3.在Outputraster文本框中选择并命名输出的水流长度栅格数据文件名称。
分别进行顺流计算和朔流计算,输出的数据文件命名为Flowlendown和Flowlenup。
4.计算方向提供了两种选择,分别为Downstream(顺流计算)和Upstream(朔流计算)。
5.输入计算配权栅格数据。
那么对于flowlength来说,Downstream记录着其沿着水流方向到下游流域出水口中最长距离所流经的栅格数;Upstream则记录着其沿着水流方向到上游栅格的最长的距离的栅格数。
6.当设置完成后,点击OK完成。
两种方向计算出的结果图像如下图所示。
FlowlendownFlowlenup
五、河网的提取
1.河网的生成
1)河网的生成是基于汇流累积矩阵的,首先应该计算出研究区域的汇流累积矩阵。
汇流。
这里用上面计算的汇流累计栅格数据flowacc作为基础数据。
2)设定阈值。
阈值的设定在河网的提取过程是很重要的,并且直接影响到河网的提取结果。
阈值的设定应遵循科学、合理的原则。
首先应该考虑到研究的对象,研究对象中的沟谷的最小级别,不同级别的沟谷所对应的不同的阈值;其次考虑到研究区域的状况,不同的研究区域相同级别的沟谷需要的阈值也是不同的。
所以,在设定阈值时,应充分对研究区域和研究对象进行分析,通过不断的实验和利用现有地形图等其它数据辅助检验的方法来确定能满足研究需要并且符合研究区域地形地貌条件的合适的阈值。
3)栅格形式的河网的形成。
栅格河网的生成是利用mapalgebra工具集中的multimapoutput工具中的con命令或者setnull命令计算的,他们都是基于栅格进行有条件的查询并将查询结果赋予新的栅格数据中。
计算的思想是利用所设定的阈值进行整区域的分析并生成一个新的栅格图层,在生成栅格数据中将那些汇流量大于设定阈值的栅格的属性值设定为1,而小于或等于设定阈值的栅格的属性值设定为无数据。
栅格河网的生成是利用设定一个河网生成阈值,也可以利用ArcMap中的SpatialAnalysis分析模块下的RasterCalculator来计算出所有大于设定阈值的栅格,这些栅格就是河网的潜在位置。
将计算出来的栅格河网命名为streamnet。
方法一:
mapalgebra生成阈值100的河网
方法二:
RasterCalculator生成阈值1000的河网
4)栅格河网矢量化。
在hydrology工具集中提供了将上一步生成的栅格河网进行矢量化的工具streamtofeature,通过streamtofeature就可以得到矢量形式的河网图。
栅格数据的矢量化时,选择的是最短的路线绘制成线。
双击hydrology工具集中的streamtofeature工具,如下图
5)在Inputstreamraster文本框中,选择streamnet1;在Inputflowdirectionraster文本框中,输入由无洼地计算出来的水流方向数据fdirfill;在Outputpolylinefeatures文本框中将输出的数据命名为streamfea3。
生成的矢量数据如图所示。
2.streamlink的生成
1)Streamlink基于水流方向数据和栅格河网数据计算,首先在ArcMap里将水流方向数据fdirfill和栅格河网数据streamnet3打开。
2)双击hydrology工具集中的streamlink工具,弹出如图所示的streamlink计算的对话框。
在Inputstreamraster文本框中选择streamnet3,在Inputflowdirectionraster文本框中选择fdirfill。
在Outputraster文本框中将输出数据名称设为StreamLink。
然后点击OK进行运算。
3)经过Streamlink计算之后,它将栅格河网分成不包含汇合点栅格河网片段,并将片断进行记录,在属性表中除了记录该片段的ID号之外,还记录着每个片段所包含的栅格数。
如图所示。
3.河网分级的生成
1)同Streamlink的计算一样,streamorder的计算也是基于水流方向数据和栅格河网数据的基础上,首先在ArcMap里将水流方向数据fdirfill和栅格河网数据streamras打开。
2)双击hydrology工具集中的streamorder工具,弹出streamorder计算的对话框。
在Inputstreamraster文本框中选择streamnet3,在Inputflowdirectionraster文本框中选择fdirfill。
分别用Strahler分级和Shreve分级对河网进行分级,改输出数据名称分别设为Streamostr和Streamoshr,单击OK完成。
计算出的两种河网分级分别如下图所示:
Strahler分级
Shreve分级
六、流域的分割
1.流域盆地的确定
在ArcGIS中,流域盆地的计算是利用hydrology工具集中的basin工具来进行计算的。
1)双击hydrology工具集中的basin工具,打开流域盆地计算的对话框。
如下图所示。
2)输入数据为水流方向数据fdirfill,设置输出数据文件名为basin。
3)点击OK完成。
结果如下
2.汇水区出水口的确定
经过上一步得到的流域盆地是一个比较大的流域盆地,在很多的水文分析中,还需要基于更小的流域单元进行分析,那么就需要将这些流域从大的流域中分解出来,那么就需要进行流域的分割。
而流域的分割首先是要确定小级别的流域的出水口的位置,小级别的流域出水口的位置可以用spatialanalysistools工具箱下的hydrology工具集中的snappourpoint的工具寻找。
它的思想是利用一个记录着point点栅格数据,在这个数据层中,那些属性值存在的点作为潜在的出水点,在该点位置上在指定距离内在汇流累积量的数据层上搜索那些具有较高汇流累积量栅格点的位置,这些搜索到的栅格点就是小级别的流域的出水点。
也可以利用已有的出水点的矢量数据。
如果没有出水点的栅格或矢量数据,可以用已有河网数据进一步生成的streamlink数据作为汇水区的出水口数据。
因为streamlink数据中隐含着河网中每一条河网弧段的联结信息,包括弧段的起点和终点等,相对而言,弧段的终点就是该汇水区域的出水口所在位置.
3.集水流域的生成
对于低级的集水区的生成,可以使用hydrology工具集中的watershed工具生成。
其思想如下:
先确定一个出水点,也就是该集水区的最低点,然后结合水流方向数据,分析搜索出该出水点上游所有流过该出水口的栅格,直到所有的该集水区的栅格都确定了位置,也就是搜索到流域的边界,分水岭的位置。
1)首先在ArcMap中打开水流方向数据fdirfill和流域出口点数据streamlink,因为集水流域的计算需要流域出口点的位置数据和水流方向数据。
水流方向数据使用基于无洼地DEM计算出的fdirfill数据。
对于流域出口点数据,可以使用streamlink数据,因为在它的属性中具有每条河网弧段信息,而这些河网弧段的形成则是有其附近的栅格构成的贡献区域的汇流结果,这个贡献区域也就是该河网弧段的集水区域。
2)双击hydrology工具集中的watershed工具,打开集水区域(贡献区域)计算的对话框。
分别在水流方向数据和出水口数据输入的文本框中选择fdirfill和streamlink数据,在输出数据中将输出的文件名改为watershed。
3)点击OK,进行集水区域的计算,计算结果如图所示,为了更好的表现流域的分割效果,在此窗口中还将上面计算的流域盆地和矢量河网的数据打开,从结果中可以知道,通过streamlink作为流域的出水口数据所得到的集水区域是每一条河网弧段集水区域,也就是要研究的最小沟谷的集水区域,它将一个大的流域盆地按照河网弧段将其分为一个个的小的集水盆地。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 ArcGISHydrology 剖析