1 前言

本文主要记录了使用ArcGIS水文分析过程划分基于地形的集水区过程,主要参照了李远祥专家有关于水文分析的博客内容。经过初步的成功实验后,特记录了整个步骤流程。

1.1 实验数据

  1. DEM数据:等高线数据或DEM数据,本人利用的是来自地理空间数据云的30米DEM数据。对于等高线数据,也可以参照李远祥专家在ArcGIS水文分析实战教程(3)DEM数据准备中所述方法进行DEM数据预处理。
  2. 研究区边界:本文所用到的研究区边界下载自阿里云平台,并使用mapshaper生成带空间属性的区域边界。
  3. 实验设备:操作软件为 ArcGIS 10.8.1

1.2 操作步骤简要概括

  1. 按研究区边界裁剪DEM
  2. 流向初提取
  3. 寻找汇
  4. 计算Z值限制
  5. 精确填洼
  6. 流向分析
  7. 流量统计
  8. 利用条件函数进行流量划分
  9. 河流分级
  10. 河流链接
  11. 栅格河网矢量化
  12. 绘制倾泻点
  13. 栅格流域提取
  14. 栅格流域转矢量面要素

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数据划定汇水区具体步骤过程相关推荐

  1. Arcgis利用dem数据生成等高线

    Arcgis利用dem数据生成等高线 问题描述 如果你在网上查,如何利用dem生成等高线,那么所有教程基本都是一样的,告诉你在工作箱中找到等值线工具,填入参数即可生成等高线.贴几个网上的链接教程: a ...

  2. arcgis dem栅格立体感_如何使用ArcGIS从DEM数据中提取水系

    1. 概述 在比较偏远的地方,往往会缺少水文信息,我们可以通过ArcGIS对高程DEM数据进行水文分析,为地表水流建立模型,进而获取到该地的水文信息,DEM数据精度越高,获取到的水文数据精度也就越高, ...

  3. arcmap提取dem高程_如何使用ArcGIS从DEM数据中提取水系

    1. 概述 在比较偏远的地方,往往会缺少水文信息,我们可以通过ArcGIS对高程DEM数据进行水文分析,为地表水流建立模型,进而获取到该地的水文信息,DEM数据精度越高,获取到的水文数据精度也就越高, ...

  4. gis怎么提取水系_如何使用ArcGIS从DEM数据中提取水系

    1. 概述 在比较偏远的地方,往往会缺少水文信息,我们可以通过ArcGIS对高程DEM数据进行水文分析,为地表水流建立模型,进而获取到该地的水文信息,DEM数据精度越高,获取到的水文数据精度也就越高, ...

  5. ArcGIS 从DEM 数据提取对应点的高程值

    将DEM数据文件和一个shapefile点文件(分别命名为"DEM"和"Point")加载到ArcMap中. [Spatial Analyst 工具]-[提取分 ...

  6. ArcGIS:根据DEM数据提取区域水系

    目录 前言 0.认识需要哪些工具 1.填洼 2.流向分析 3.流量分析 4.栅格计算器 5.河网矢量化 6.结语 前言 利用DEM数据提取区域水系是常见的水文分析应用,为了获取区域的矢量水系数据,我们 ...

  7. arcgis里dem填洼_ArcGIS水文分析实战教程(4)地形预处理

    ArcGIS水文分析实战教程(4)地形预处理 本章导读:ArcGIS的水文分析工具是基于DEM进行地表水流动的模拟,其本身不涉及到精确数值的水流流量.在形成径流的过程中考虑的全是地形因素,D8单流向算 ...

  8. ArcGIS水文分析实战教程(3)DEM数据准备

    ArcGIS水文分析实战教程(3)DEM数据准备 ArcGIS的水文分析工具需要DEM数据才能执行,因此在使用之前必须要准备好DEM数据并且针对DEM数据进行预处理.本章节内容介绍如何根据已有的数据制 ...

  9. ArcGIS中利用DEM数据生成地形图既视感的等高线;利用掩膜进行等高线注记;DEM的可视化表达总结

    文章目录 前言 一.效果图展示 二.地形(DEM)可视化表达 1.一维可视化 2.二维可视化 3.三维可视化 二.ArcGIS中利用DEM数据制作等高线 1.DEM数据下载 2.镶嵌 3.提取出需要的 ...

  10. 【Arcgis水文分析】使用DEM数据提取流域水系

    [下载数据]根据自己的需要选择分辨率精度下载DEM高程数据:地理空间数据云: [ 拼接DEM数据]打开Arcmap使用镶嵌工具对下载的DEM数据进行拼接:Data Management Tools.t ...

最新文章

  1. BZOJ4076 : [Wf2014]Maze Reduction
  2. Oracle-AWR管理包DBMS_WORKLOAD_REPOSITORY.MODIFY_SNAPSHOT_SETTINGS
  3. 16款新品发布,数据揭秘小米MIX荣归背后逻辑
  4. 环信集成 2---基于环信Demo3.0,实现单聊功能
  5. Delphi-IOCP学习笔记三====工作线程和Listener
  6. 单片机工作笔记002---认识CPU的结构_哈佛结构和冯·诺依曼结构的区别
  7. 计算机系统层次中 从上层,在计算机系统层次结构中下层是上层的上层是下层的一个.ppt...
  8. CPython 解释器 与 VC 编译器版本对应关系
  9. Linux线程同步读写锁 rwlock
  10. Android自动挂断电话
  11. 微信小程序图片转换成文字_涨知识!这个微信小程序,能将纸上文字转换成电子版...
  12. python正则匹配数字_python——正则匹配数字
  13. (转)PQ分区魔术师中文版分区教程
  14. 从零开始学python第八篇
  15. MyBatisCodeHelper-Pro插件
  16. 国科大本科生9个月设计出处理器芯片 可运行Linux系统
  17. 跑步用挂脖耳机好还是无线耳机、公认最好的跑步耳机推荐
  18. 电影暗杀了一只巨可爱的猫[转帖]
  19. 小米CC9系列发布会,要用文科设计打动女生的雷军真的懂女生么
  20. python获取路由器信息_python 3.5实现检测路由器流量并写入txt的方法实例

热门文章

  1. word打开html是空白,网页复制到Word之后或者出现空白或者格式变乱该怎么办
  2. Linux安全模块(LSM)学习——简单的LSM demo(1)
  3. 斯托尔兹(Stolz)定理
  4. 有关电脑新手入门教程,电脑入门基础知识有哪些?
  5. Kindle 3---原系统下的字典的使用
  6. QCC302X/QCC303X蓝牙对讲与蓝牙扩音器
  7. 交互设计起源于计算机的人机界面设计的例子,交互设计概述.ppt
  8. Unity AI 之 行为树 的简单介绍
  9. 使用jemalloc优化java_C++性能优化(十) —— JeMalloc
  10. 计算机桌面计划,如何用便签在电脑桌面安排日程?