ArcGIS使用DEM数据划定汇水区具体步骤过程

您所在的位置:网站首页 怎么提取流域边界 ArcGIS使用DEM数据划定汇水区具体步骤过程

ArcGIS使用DEM数据划定汇水区具体步骤过程

2023-03-31 08:05| 来源: 网络整理| 查看: 265

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到7次幂的离散值,即:在这里插入图片描述

(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值选取的大小都影响着最终结果。



【本文地址】


今日新闻


推荐新闻


CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3