水文分析通过建立地表水文模型,研究与地表水流相关的各种自然现象,在城市和区域规划、农业及森林、交通道路等许多领域具有广泛的应用。我国是一个自然灾害十分频繁的国家, 洪涝灾害一直严重威胁着人民生命财产安全和社会的稳定与发展。近几十年来, 自然资源的开发利用不断扩大, 城乡经济建设飞速发展, 洪水出现的频率及其造成的损失也不断的增加。特别是今年,特大的洪水袭击了我国南方,造成了重大损失。这里,我们以水库淹没为例,从GIS的角度讲解一下如何使用ArcGIS获取洪水淹没区域和计算面积。
打开青海湖部分30m高程DEM数据
(1)填洼(Fill)
在ArcToolbox中点击“Spatial Analys工具水文分析 | 填洼”,调用填洼工具。
如果填洼不成功,打开工具框的Environments(环境)设置,找到Parallel Processing(合并处理)工具,将里面的值设置为“0”即可(下同)
在显示的填洼对话框内,输入表面栅格数据为下载的DEM数据,设置好输出表面栅格
点击确定之后可以得到填洼后的栅格
(2)流向
在ArcToolbox中点击“Spatial Analyst工具 | 水文分析 | 流向”,调用流向工具。在显示的流向对话框内,输入表面栅格数据为填洼后的栅格数据
点击确定之后可以得到流向栅格
(3)河流长度
在ArcToolbox中点击“Spatial Analyst工具 | 水文分析 | 流向长度”,调用流向长度工具。在显示的流向长度对话框内,输入流向栅格数据
点击确定之后可以得到流向长度栅格
(4)流量
在ArcToolbox中点击“Spatial Analyst工具 | 水文分析 | 流量”,调用流量工具。在显示的流量对话框内,输入流向栅格数据为上一步得到的流向数据
点击确定之后可以得到流量栅格
(5)设立阈值
在ArcToolbox中点击“Spatial Analyst工具 | 地图代数 | 栅格计算器”,调用栅格工具。在显示的对话框内,输入公式:Con(FlowAcc_Flow>3000,1)
点击确定之后可以得到计算阈值栅格
(6)栅格河网矢量化
在ArcToolbox中点击“Spatial Analyst工具 | 水文分析| 栅格河网矢量化”,调用工具。在显示的对话框内,输入河网栅格:上一步阈值处理的RC_yuzhi,输入流向数据
点击确定之后可以得到栅格数据的矢量化数据
(7)平滑河网
加载编辑器工具条,启动编辑器开始编辑。在工具条中加载高级编辑工具条。
在高级编辑工具条中单击平滑按钮,输入允许最大偏移参数,点击确定
点击确定之后可以得到平滑的矢量化数据
(8)生成河流连接
启动ArcToolbox,展开“Spatial Analyst工具 | 水文分析 | 河流连接”,弹出河流连接对话框。加载相关数据,点击确定。
查看河流连接的属性表。
(9)河网分级
河网分级是一种将级别数分配给河流网络中的连接线的方法,此级别是一种根据支流数对河流类型进行识别和分类的方法。仅需知道河流的级别,即可推断出河流的某些特征。河网分级方法分为Strahler法和Shreve法,这两种方法始终将1级分配给上游河段。
启动ArcToolbox,展开“Spatial Analyst工具 | 水文分析 | 河流分级”,弹出河流分级对话框。加载相关数据。输入流网栅格:RC_yuzhi,输入流向数据
点击确定。得到Strahler方法的河流分级图
选择Shreve方法
点击确定。得到Shreve方法的河流分级图
(10)盆域分析
流域盆地是由分水岭分割而成的汇水区域,它通过分析水流方向数据确定所有相互连接并处于同一流域盆地的栅格。
启动ArcToolbox,展开“Spatial Analyst | 水文分析 | 盆域分析”,弹出盆域分析对话框。加载相关数据,点击确定。
(11) 生成分水岭
经过上一步得到的流域盆地是一个比较大的流域盆地,在很多的水文分析中,还需要基于更小的流域单元进行分析,这就需要进行流域的分割。
启动ArcToolbox,展开“Spatial Analyst工具 | 水文分析 | 分水岭”,弹出分水岭对话框。加载相关数据,倾泻点字段设置为VALUE,点击确定
(12)分水岭栅格转面
启动ArcToolbox,展开“Conversation Tools工具 | 从栅格转出 |栅格转面”,弹出对话框。加载相关数据,倾泻点字段设置为VALUE,点击确定。输入栅格:分水岭数据
(13)盆域栅格转面
启动ArcToolbox,展开“Conversation Tools工具 | 从栅格转出 |栅格转面”,弹出对话框。加载相关数据,倾泻点字段设置为VALUE,点击确定。输入栅格:盆域数据
(14)计算DEM面积
在ArcToolbox内点击“3D Analyst工具 | 功能性表面 | 表面体积”,调用表面体积工具。在显示的表面体积对话框内,输入表面为上一步设为空函数后的DEM数据,设置输出结果的文本文件的位置,由于DEM数值都是500,所以选择ABOVE或者BELOW都一样,这里选择ABOVE
以此小块区域为例——得到的属性表即可查看所需区域的面积特征