ArcGIS使用DEM数据划定汇水区具体步骤过程
1 前言
本文主要记录了使用ArcGIS水文分析过程划分基于地形的集水区过程,主要参照了李远祥专家有关于水文分析的博客内容。经过初步的成功实验后,特记录了整个步骤流程。
1.1 实验数据
- DEM数据:等高线数据或DEM数据,本人利用的是来自地理空间数据云的30米DEM数据。对于等高线数据,也可以参照李远祥专家在ArcGIS水文分析实战教程(3)DEM数据准备中所述方法进行DEM数据预处理。
- 研究区边界:本文所用到的研究区边界下载自阿里云平台,并使用mapshaper生成带空间属性的区域边界。
- 实验设备:操作软件为 ArcGIS 10.8.1
1.2 操作步骤简要概括
- 按研究区边界裁剪DEM
- 流向初提取
- 寻找汇
- 计算Z值限制
- 精确填洼
- 流向分析
- 流量统计
- 利用条件函数进行流量划分
- 河流分级
- 河流链接
- 栅格河网矢量化
- 绘制倾泻点
- 栅格流域提取
- 栅格流域转矢量面要素
2 具体操作步骤
2.1 按研究区边界裁剪DEM
网上有诸多教程,此步骤省略
2.2 流向初提取
输入裁剪后的DEM数据,输出流向数据,并且命名为 FlowDir_0
从流向提取结果中发现其栅格值是一组连续值,即不符合D8算法,因此可确定该DEM中存在汇。
2.3 寻找汇
输入流向数据 FlowDir_0,输出汇栅格,并命名为 Sink_FlowDir0
2.4 计算Z值限制
(1) 使用【分水岭】(即【集水区】)工具为每个汇创建汇流区域栅格
输入流向数据FlowDir_0、汇点数据Sink_FlowDir0,输出汇流区域栅格,并命名为 Watersh_Sink0
(2) 使用【分区统计】在每个汇的分水岭中创建最小高程的栅格
输入汇流区域栅格 Watersh_Sink0、区域字段为VALUE、裁剪后的DEM数据、统计类型为MINIMUM,输出每个汇分水岭的最小高程栅格,并命名为sink_min
(3) 使用【区域填充】在每个汇的分水岭中创建最大高程的栅格
输入汇流区域栅格 Watersh_Sink0、裁剪后的DEM数据、输出每个汇分水岭的最大高程栅格,并命名为sink_max
(4) 使用【栅格计算器】计算Z值
输入 “sink_max” - “sink_min”,输出Z值,并命名为 Z
计算完成后可得到一组Z值如下:
其中,最小值为0,最大值为8
2.5 精确填洼
使用【填洼】工具,输入原始裁剪后的DEM,Z值取为上一步骤计算结果中的最大值(8),并进行加1(即为9);
输出填洼后的DEM,并命名为 dem_fill
2.6 流向分析
输入填洼后的DEM数据dem_fill,输出流向数据,并且命名为 FlowDir_1
如果填洼成功,则此处的流向计算结果应该是一组关于2的0到10次幂的离散值,即:
(1) 否则即填洼失败,需要再次重复步骤2~步骤5,其中DEM数据使用上一次填洼后的DEM数据即dem_fill;
(2) 本文也是进行了1次重复寻找汇以及计算Z值,最终才输出无洼DEM (dem_fill_1)
2.7 流量统计
输入流向数据 FlowDir_1,输出流量栅格数据,并命名为 FlowAcc
2.8 利用条件函数进行流量划分
流量划分的目的是确定地标径流能汇聚成河流的最小蓄积栅格,以R表示;其中R值的确定公式为:
R = W / ( X × Y ),其中
W,为划分的集水区的汇水面积;
X、Y为DEM的最小单元格尺寸,如我这里用的说是30米分辨率的DEM,则X、Y都为30
操作为:
使用【条件函数】工具,输入流量数据 FlowAcc,表达式为VALUE>=R,条件为true时选择流量数据 FlowAcc,输出筛选栅格,并命名Con_1000
在本例中,取划分的集水区的汇水面积为0.9㎞²,(即900000m²),则取R=1000
若划分的汇水面积过小,即R,对于小研究区域来说,则提取到的河网数据少,从而使得研究意义降低,因此要选取恰当的R值
2.9 河网分级
输入流量划分后的结果 Con_1000,以及流量栅格数据 FlowDir_2,输出河网分级结果 StreamO_1000
2.10 河流链接
输入流量划分后的结果 Con_1000,以及流量栅格数据 FlowDir_2,输出河流链接结果 StreamL_1000
2.11 栅格河网矢量化
(1) 使用【栅格河网矢量化】工具对分级后的河网进行矢量化操作;
输入分级后和河流数据StreamO_1000、流量数据FlowDir_2,输出河网折线要素,并命名StreamT_Con_1000.shp
(2) 打开矢量化结果StreamT_Con_1000.shp的图层属性,选择符号系统——数量——分级色彩,其中字段值选项选择‘grid_code’,色带可自行选择,其他项为默认值
(3) 在标注选项卡中进行如下相应设置
(4) 生成河流的数字化方向
从目录中再次选择矢量化结果加入操作图层中,并对该线要素进行符号设置:
以上步骤操作结果如下:
2.12 绘制倾泻点
使用【要素折点转点】工具来捕捉其线段的终点,即倾泻点的位置;
该工具位于【数据管理】–【要素】–【要素折点转点】;
输入要素为矢量河网数据StreamT_Con_1000,点类型为END,数据输出记作FeatureVert_1000.shp
2.13 栅格流域提取
选择【集水区】工具,输入流向栅格数据FlowDir_2、河流链接数据StreamL_1000,倾泻点字段值VALUE,输出流域栅格 watersh_1000
2.14 栅格流域转矢量面要素
使用【栅格转面】工具,(【转换工具】→【由栅格转出】→【栅格转面】),将流域栅格转矢量面要素。
输入流域栅格数据 watersh_1000、字段值选择VALUE,输出流域矢量面数据,记作RasterT_1000.shp
对矢量化结果进行相关显示处理
最终结果为:
3 总结
在流量划分中,最小蓄积栅格R的确定影响着最终的集水区划分结果,本实验做了4组不同的比较,其结果如下:
(1) R值选取为100时:
(2) R值选取为500时:
(3) R值选取为1000时:
(4) 不进行流量划分,即河网分级、河流链接使用原始流量统计结果时:
由以上结果可知,是否进行流量划分,以及R值选取的大小都影响着最终结果。
ArcGIS使用DEM数据划定汇水区具体步骤过程相关推荐
- Arcgis利用dem数据生成等高线
Arcgis利用dem数据生成等高线 问题描述 如果你在网上查,如何利用dem生成等高线,那么所有教程基本都是一样的,告诉你在工作箱中找到等值线工具,填入参数即可生成等高线.贴几个网上的链接教程: a ...
- arcgis dem栅格立体感_如何使用ArcGIS从DEM数据中提取水系
1. 概述 在比较偏远的地方,往往会缺少水文信息,我们可以通过ArcGIS对高程DEM数据进行水文分析,为地表水流建立模型,进而获取到该地的水文信息,DEM数据精度越高,获取到的水文数据精度也就越高, ...
- arcmap提取dem高程_如何使用ArcGIS从DEM数据中提取水系
1. 概述 在比较偏远的地方,往往会缺少水文信息,我们可以通过ArcGIS对高程DEM数据进行水文分析,为地表水流建立模型,进而获取到该地的水文信息,DEM数据精度越高,获取到的水文数据精度也就越高, ...
- gis怎么提取水系_如何使用ArcGIS从DEM数据中提取水系
1. 概述 在比较偏远的地方,往往会缺少水文信息,我们可以通过ArcGIS对高程DEM数据进行水文分析,为地表水流建立模型,进而获取到该地的水文信息,DEM数据精度越高,获取到的水文数据精度也就越高, ...
- ArcGIS 从DEM 数据提取对应点的高程值
将DEM数据文件和一个shapefile点文件(分别命名为"DEM"和"Point")加载到ArcMap中. [Spatial Analyst 工具]-[提取分 ...
- ArcGIS:根据DEM数据提取区域水系
目录 前言 0.认识需要哪些工具 1.填洼 2.流向分析 3.流量分析 4.栅格计算器 5.河网矢量化 6.结语 前言 利用DEM数据提取区域水系是常见的水文分析应用,为了获取区域的矢量水系数据,我们 ...
- arcgis里dem填洼_ArcGIS水文分析实战教程(4)地形预处理
ArcGIS水文分析实战教程(4)地形预处理 本章导读:ArcGIS的水文分析工具是基于DEM进行地表水流动的模拟,其本身不涉及到精确数值的水流流量.在形成径流的过程中考虑的全是地形因素,D8单流向算 ...
- ArcGIS水文分析实战教程(3)DEM数据准备
ArcGIS水文分析实战教程(3)DEM数据准备 ArcGIS的水文分析工具需要DEM数据才能执行,因此在使用之前必须要准备好DEM数据并且针对DEM数据进行预处理.本章节内容介绍如何根据已有的数据制 ...
- ArcGIS中利用DEM数据生成地形图既视感的等高线;利用掩膜进行等高线注记;DEM的可视化表达总结
文章目录 前言 一.效果图展示 二.地形(DEM)可视化表达 1.一维可视化 2.二维可视化 3.三维可视化 二.ArcGIS中利用DEM数据制作等高线 1.DEM数据下载 2.镶嵌 3.提取出需要的 ...
- 【Arcgis水文分析】使用DEM数据提取流域水系
[下载数据]根据自己的需要选择分辨率精度下载DEM高程数据:地理空间数据云: [ 拼接DEM数据]打开Arcmap使用镶嵌工具对下载的DEM数据进行拼接:Data Management Tools.t ...
最新文章
- BZOJ4076 : [Wf2014]Maze Reduction
- Oracle-AWR管理包DBMS_WORKLOAD_REPOSITORY.MODIFY_SNAPSHOT_SETTINGS
- 16款新品发布,数据揭秘小米MIX荣归背后逻辑
- 环信集成 2---基于环信Demo3.0,实现单聊功能
- Delphi-IOCP学习笔记三====工作线程和Listener
- 单片机工作笔记002---认识CPU的结构_哈佛结构和冯·诺依曼结构的区别
- 计算机系统层次中 从上层,在计算机系统层次结构中下层是上层的上层是下层的一个.ppt...
- CPython 解释器 与 VC 编译器版本对应关系
- Linux线程同步读写锁 rwlock
- Android自动挂断电话
- 微信小程序图片转换成文字_涨知识!这个微信小程序,能将纸上文字转换成电子版...
- python正则匹配数字_python——正则匹配数字
- (转)PQ分区魔术师中文版分区教程
- 从零开始学python第八篇
- MyBatisCodeHelper-Pro插件
- 国科大本科生9个月设计出处理器芯片 可运行Linux系统
- 跑步用挂脖耳机好还是无线耳机、公认最好的跑步耳机推荐
- 电影暗杀了一只巨可爱的猫[转帖]
- 小米CC9系列发布会,要用文科设计打动女生的雷军真的懂女生么
- python获取路由器信息_python 3.5实现检测路由器流量并写入txt的方法实例
热门文章
- word打开html是空白,网页复制到Word之后或者出现空白或者格式变乱该怎么办
- Linux安全模块(LSM)学习——简单的LSM demo(1)
- 斯托尔兹(Stolz)定理
- 有关电脑新手入门教程,电脑入门基础知识有哪些?
- Kindle 3---原系统下的字典的使用
- QCC302X/QCC303X蓝牙对讲与蓝牙扩音器
- 交互设计起源于计算机的人机界面设计的例子,交互设计概述.ppt
- Unity AI 之 行为树 的简单介绍
- 使用jemalloc优化java_C++性能优化(十) —— JeMalloc
- 计算机桌面计划,如何用便签在电脑桌面安排日程?