ArcGIS水文分析提取河网及流域
在进行某些研究或者一些论文插图显示的时候,有时我们会碰到在部分资料中找不到一些小的河流或者流域的数据的情况,这里讲述通过DEM数据生成河网及流域。
一、数据来源
四川省高程数据来源于中国科学院资源环境科学与数据中心(中国科学院资源环境科学与数据中心 (resdc.cn))。
可以看到该影像覆盖区域海拔高度在97m~6511m之间。
二、数据处理
1. 填洼
利用填洼工具填充四川省数字高程数据中的凹陷点来移除数据中的小缺陷。
工具位置: Spatia Analyst 工具 -> 水文分析 -> 填洼
由下图亦可看出最低点由97变为了109,说明一些凹陷点得到填充。
2. 流向分析
流向分析即确定水流方向,在ArcGIS中是采用的D8算法进行,D8算法就是假设有无限的降水,雨水不断的在地表形成径流,假设单个栅格中的水流只能流入与之相邻的8 个栅格中。它用最陡坡度法来确定水流的方向,即在3×3 的DEM 栅格上,计算中心栅格与各相邻栅格间的距离权落差(即栅格中心点落差除以栅格中心点之间的距离),取距离权落差最大的栅格为中心栅格的流出栅格。
工具位置:Spatial Analyst 工具 -> 水文分析 -> 流向
在得到的图层中像元值的值被分为8个固定的值,很多人并不了解其所代表的含义,其实是以数值表示每个单元的流向。一般数字变化范围是1~255。其中1:东;2:东南;4南;8:西南;16:西;32:西北;64:北;128:东北。除上述数值之外的其它值代表流向不确定,这是由DEM中“洼地”和“平地”现象所造成的。
所谓“洼地”即某个单元的高程值小于任何其所有相邻单元的高程。这种现象是由于当河谷的宽度小于单元的宽度时,由于单元的高程值是其所覆盖地区的平均高程,较低的河谷高度拉低了该单元的高程。这种现象往往出现在流域的上游。“平地”指相邻的8个单元具有相同的高程,与测量精度、DEM单元尺寸或该地区地形有关。这两种现象在DEM 中相当普遍。
3. 流量统计
通过流量工具统计累计每个流向栅格的总数。在ArcGIS软件中,这个流量栅格并不是水文学意义上的河流流量,与水文站实测的数据是两个完全不同意义的名词。这里的流量指的是有多少个格子的水流汇流到某一点上,是一个空间范围概念。如果要转换为水文学上的术语,以水文站为例,那就是该水文站对应的上游集水区面积。
ArcGIS的流量统计是需要流向栅格进行参考的,其统计原理如下图
工具位置:Spatial Analyst 工具 -> 水文分析 -> 流量
4. 确定最小径流
这部分是根据自己的需要来确定想要输出的最小径流量的河流,这里采用栅格计算器工具计算,本例中设定的最小流量为5000。
工具位置:Spatial Analyst 工具 -> 地图代数 -> 栅格计算器
5. 河流链接
为了将像元点串起来并做了不同线段的划分,反应河流之间的关系,这里进行河流链接操作。
工具位置:Spatial Analyst 工具 -> 水文分析 -> 河流链接
从分析结果可以看到,链接之后离散的像元会被串联起来,形成独立的河流。如下图所示,河网已经被串联成10713条河流,但仍然使用栅格的形式表示。
6. 河网分级
利用河网分级工具对支流干流进行划分。在水文学上,河流链接不是为了提取矢量河流,而是为河流分级做准备。水文学更多的是关注河流之间的关系、河流与流域的关系,而非河流本身的图形;即便关注的是单个河流的几何图形,但一般也会跟高程带上关系,例如河流比降。所以,河流链接之后,一般就是对河流进行分级了。
工具位置:Spatial Analyst 工具 -> 水文分析 -> 河网分级
7. 河网矢量化
利用栅格河网矢量化工具将河网栅格转化为栅格矢量,方便我们后续进行矢量化。
工具位置:Spatial Analyst 工具 -> 水文分析 -> 栅格河网矢量化
当然也可以使用栅格转折线工具。
工具位置:转换工具 -> 由栅格转出 -> 栅格转折线
相对来说栅格河网矢量化工具更专业一点,因为除了转化为矢量河流之外,还会自动的创建了河流链接。下图(左边)可以通过属性表看到生成的10713条河流,跟前面所说的河流链接工具的结果(右边)是一样的。
8. 流域提取
首先打开属性表,利用按属性选择工具选取3、4、5、6级的河流,再利用融合工具提取这部分河流,并融合同级别河流,取消勾选创建多部件要素,不将同级别但不接触的河流融合在一起。
工具位置:数据管理工具 -> 制图综合 -> 融合
利用要素折点转点工具从矢量河网数据中提取分水岭工具所需要的要素倾泻点(河流出水口),将点类型设置为END,即河流终点为出水口。
工具位置:数据管理工具 -> 要素 -> 要素折点转点
利用捕捉倾泻点工具进行校正。其目的是要找到这个点附近的区域流量栅格中的最大值,确保这个流量值是由上游区域贡献的流量。否则,实际输出的倾泻点数据有可能与之叠加的流量栅格不是该区域的最大值,导致出现非常细小的流域面。
工具位置:Spatial Analyst 工具 -> 水文分析 -> 捕捉倾泻点
由于上面已经得到河流链接的数据,这里通过分水岭工具进行计算获取流域范围。
工具位置:Spatial Analyst 工具 -> 水文分析 -> 分水岭
三、符号化制图
1. 河网符号化显示
右键河网矢量图层,选择属性 -> 符号系统 -> 类别 -> 唯一值,值字段设置为grid_code,选择蓝色渐变,设置河网级别由大到小的线段宽度分别为3、2.5、2、1.5、1、0.5,在标题标注栏填入河流分级。
2. 裁剪
右键数据框,选择属性 -> 数据框,在裁剪选项中选择裁剪至形状。
点击指定形状,选择要素的轮廓为四川省行政边界。
3. 更改布局
在布局工具栏中选择更改布局 -> ISO(A)Page Sizes -> ISO AO Landscape.mxd,当然也可以挑选个人喜欢的布局。
4. 添加地图必要元素
在地图中插入标题、图例、指北针、比例尺等地图必要元素准备出图。
添加比例尺时点击属性可以设置单位,这里我们选择为千米。
5. 添加经纬网
右键数据框选择属性,在格网选项卡中点击新建格网。
在弹出的对话框第二页选择仅标注,第三页每个长轴主刻度数设置为9,其它页均为默认。
设置完成后在格网选项卡中选择属性,在弹出的对话框中,轴页长轴主刻度和分刻度的刻度大小分别设置为20磅和8磅;标注页将文字大小设置为48,垂直标注勾选左和右。
四、成果展示
来源:GIS前沿
ArcGIS水文分析提取河网及流域相关推荐
- ArcGIS水文分析提取河网、流域、湖泊
本文只保留了软件操作并不涉及知识原理,为GIS小白提供的实际操作和避坑指南.本意为自己留存复盘,如有错误请各位大佬指点.(使用的软件版本为ArcGIS10.6) 一.提取河网 1.裁剪 数据管理工具 ...
- 【教程】ArcGIS水文分析提取河流及流域
在进行某些研究或者一些论文插图显示的时候,有时我们会碰到在部分资料中找不到一些小的河流或者流域的数据的情况,这里讲述通过DEM数据生成河流及流域. 一.数据准备 本文中使用的为SRTM DEM数据,主 ...
- ArcGIS水文分析实战教程(13)流域提取流程
ArcGIS水文分析实战教程(13)流域提取流程 本章导读:这里的流域是一个自然概念,而不是行政概念,或者称之为集水区更加合理一些.在流域提取中,其过程都是先定义河流然后才能通过河流定义流域.其中流域 ...
- 水文分析提取河网_基于图的河网段地理信息分析排序算法
水文分析提取河网 The topic of this article is the application of information technologies in environmental s ...
- ArcGIS水文分析实战教程(7)细说流域提取
ArcGIS水文分析实战教程(7)细说流域提取 本章导读:流域是水文领域上是频繁出现的词汇,贯穿整个水文专业.由于流域涉及到的空间范围非常的广,一般不能通过测量的方式去精确划定.而且流域这一词在水利行 ...
- ArcGIS水文分析实战教程(6)河流提取与河网分级
ArcGIS水文分析实战教程(6)河流提取与河网分级 本章导读:ArcGIS水文分析模块提供了3个工具用作河网提取.GIS和测绘专业人员无法理解,用作水文分析连DEM数据都能拿到了(精度较高的DEM数 ...
- ArcGIS水文分析实战教程(7)细说流域提取 1
ArcGIS水文分析实战教程(7)细说流域提取 本章导读:流域是水文领域上是频繁出现的词汇,贯穿整个水文专业.由于流域涉及到的空间范围非常的广,一般不能通过测量的方式去精确划定.而且流域这一词在水利行 ...
- 【Arcgis水文分析】使用DEM数据提取流域水系
[下载数据]根据自己的需要选择分辨率精度下载DEM高程数据:地理空间数据云: [ 拼接DEM数据]打开Arcmap使用镶嵌工具对下载的DEM数据进行拼接:Data Management Tools.t ...
- ArcGIS水文分析实战教程(14)提取水库湖泊等面状水系
ArcGIS水文分析实战教程(14)提取水库湖泊等面状水系 本章导读:如果说前面章节一直都是针对河流.河网.流域等进行论述,围绕的都是从基本理论到实践操作,那么湖泊.水库等面状水系的提取,就是对水文分 ...
最新文章
- openssl 从cer文件中提取公钥
- 利用序列化与反序列化对Java对象进行深度复制
- python最全知识点_史上最全的Python面向对象知识点疏理
- java.lang.OutOfMemoryError: Java heap space 错误及解决办法
- 别瞎忙活:创业公司的6条时间管理策略
- JDK1.8 中 ConcurrentHashMap源码分析(二)元素添加是线程安全的
- 使用vlc播放器做rtsp流媒体服务器
- 外贸用ERP仓储系统有那些好处?
- JS与Asp.Net的传值
- Jmeter安装步骤
- Permute 3.6.5 小巧便捷的多媒体文件格式转换器
- 【论文笔记】基于聚类特征深度LSTM的语音情感识别
- 支持向量机之思维导图
- The field file exceeds its maximum permitted size of 1048576 bytes.
- Android Studio 制作.9图片时无法拖拉,导致没法办制作.9图片
- Sentinel控制台实时监控空白
- 初次使用snaker集成项目中遇到的问题
- Luogu P2751 [USACO4.2]工序安排Job Processing
- led亮度鉴别测试软件,测试LED颜色与亮度——FeasaLED测试仪
- 【数据库】浅谈WITH CHECK OPTION的作用
热门文章
- .text‘ will not fit in region `ROM‘, region `ROM‘ overflowed by 3056 bytes问题记录
- Java后端面试(五)
- Python参数化框架parameterized+单元测试框架unittest HTMLReport nose
- 电子巡更系统服务器,一种基于iBeacon技术的电子巡更系统技术方案
- 一个大牛写的JDBC工具类
- unity开发罗技仿真驾驶环境
- 北大英语考研分数计算机,08年计算机系录取分数线及名单
- 进服务器显示系统缓冲区空间不足,由于系统缓冲区空间不足或队列已满,不能执行套接字上的操作--解决方法...
- ASEMI快恢复二极管型号大全,快恢复特性关系与选型大全
- 【网络安全】MS17-010“永恒之蓝”漏洞的利用