Home 新闻中心 技术博客

在ArcGIS中勾绘全流域方法

发布时间:2021-11-06 23:57:31   浏览量:35   作者:GIS前沿

所需数据:emidalat,一个高程栅格;emidastream.shp,一个流域shapefile。

用高程栅格数据集勾绘整个流域的操作流程。数据源是高程栅格。emidastream.shp作为参照图层。

除非特别说明,所有工具都驻留在空间分析工具/水文工具集。

  1. 打开 emidalat和emidastream.shp。在这里插入图片描述
  2. 检查中emidalat是否存在小洼地。双击流向工具。选择emidalat作为输入地面栅格网,键入temp_flowd作为输出的水流方向栅格,点击确定。

在这里插入图片描述
3. 在temp_flowd建立后,双击汇工具,选择temp_flowd作为输入的水流方向栅格,将sinks设置为输出栅格,点击确定。在这里插入图片描述
在这里插入图片描述
4. 填平emidalat中的小洼地。双击填洼工具,选择emidalat作为输出地面栅格,设定emidafill为输出地面栅格,点击确定。在这里插入图片描述
5. 用emidafill继续上述操作。双击流向工具,选择emidatfill作为输入地面栅格,设定flowdirection为输出的水流方向栅格,运行命令。在这里插入图片描述
6. 接着创建一个流量累积栅格。双击流量工具,选择flowdirection为输入的水流方向栅格,键入flowaccumu为输出累积栅格,点击确定。在这里插入图片描述

  1. 选择(或设定临界值)flowaccumu中流入值超过500个像元的单元,双击空间分析工具/条件分析工具集的条件函数工具,选择flowaccumu作为输入的条件栅格,键入表达式Value > 500,键入1作为常数值,设定net 为输出栅格,并点击确定运行命令。在接下来的分析中将用到net作为输入河流栅格。

在这里插入图片描述
8.为net中的交汇点(或交叉点)间的每个河段指定一个唯一值,返回到Hydrology工具集,双击河流链接工具,选择net作为输入的河流栅格,flowdirection为输入的流向栅格,source为输出栅格。点击确定。

在这里插入图片描述

9.至此,准备了勾绘流域所需的所有数据。双击分水岭工具,选择flowdirection作为输入的流向栅格,source为输入的栅格,watershed指定为输出格网,点击OK。将watershed的符号改变成唯一值,以便可以看出每一个流域。在这里插入图片描述
ArcMap中的Python脚本完成以上操作,用这个选项,首先在ArcMap的标准工具条中打开Python窗口。假设工作空间是d:/chap(“/”用来指定路径),其中包含emidalat,你需要在“>>>”后输入一下预计,以完成:

inport arcpy

from arcpy import env

from arcpy.sa import *

env.workspace = “ d:/chap14”

arcpy.CheckExtension(“ Spatial ”)

outflowdirection = FlowDirection (“ emidalat”)

outsink = Sink (“ outflowdirection”)

outfill = Fill (“ emidalat”)

outfd = FloeDirection (“ outfill”)

outflowac = FlowAccumulation (“ outfd”)

outnet = Con (“ outflowac”, 1, 0,”VALUE>500”)

outstreamlink = StreamLink (“ outnet”,” outfd”)

outwatershed = Watershed (“ outfd”,” outstreamlink”)

outwatershed.Save (“ outwatershed”)

脚本的前5行输入了arcpy和Spatial Analyst 工具,定义了工作空间。接下来是使用FlowDirection、Sink、Fill、Flowdirection(已填洼DEM)、FlowAccumulation、Con、StreamLink和Watershed这些工具的语句。每输入一条语句,你会在ArcMap中看到输出结果。

来源:https://blog.csdn.net/qq_43173805/article/details/118555587

成都途远GIS是一家专业致力于无人机航空摄影测绘、航空数据处理、GIS地理信息系统研发、数字孪生城市制作、数字沙盘模型等业务的创新型科技公司,为您提供一站式地理信息服务。