ArcGIS提取流域边界以及绘制土地利用类型

1.  数据获取

获取流域边界需要用到DEM高程数据，以及GLC(globleland30)的土地利用类型数据。

本次绘图的示例数据精度都在30m，需要更高精度的土地利用类型数据可以利用SinoLC-1数据，这个精度在1m，需要更高精度识别河流边界，可以使用landsat8数据。

或者是公众号：GIS遥感数据处理应用

1.1 高程数据获取

DEM高程数据地址：

地理空间数据云 ，需要注册才能下载数据

选择高程数据

框选你要的区域，然后点击检索

逐个选择你需要的数据，然后下载，下载前要记得每个数据的相对位置（通过数据条带号和行编号记录）

1.2 土地利用类型数据获取

1.2.1 数据简介以及下载

这个数据不需要像landsat8，在识别完不同土地利用类型之后需要自己去一个个框选不同的土地利用类型。

GlobeLand30数据采用WGS-84坐标系。

数据下载地址：成果查询

这个网址只有2010年的数据，如果需要2020年的数据，可以去欢迎光临国家基础地理信息中心网站

这里看看三种方式，第一种网址是http开头的，现在好像进不去。关于数据获取的话，和DEM的步骤差不多。

南纬85°-北纬85°之间的区域，投影方式采用UTM投影， 6度分带，坐标单位为米，坐标不加带号。 南北纬85°-90°之间的区域，投影方式采用极地方位投影，投影面切于地球南北极点。按照产品定义的全球分幅，共有966个数据集；面积覆盖全球陆地表面1.49亿平方千米。

30m分辨率是面向全球尺度的，格网大小是30米。数据包括耕地、森林、草地、灌木地、湿地、水体、苔原、人造地表、裸地、冰川及永久积雪10种地表覆盖类型。这份数据不带高程，重点是看宏观、看区域性、看趋势。

1.2.2 数据分类系统

类型	内容	代码
耕地	用于种植农作物的土地，包括水田、灌溉旱地、雨养旱地、菜地、牧草种植地、大棚用地、以种植农作物为主间有果树及其他经济乔木的土地，以及茶园、咖啡园等灌木类经济作物种植地。	10
林地	乔木覆盖且树冠盖度超过30%的土地，包括落叶阔叶林、常绿阔叶林、落叶针叶林、常绿针叶林、混交林，以及树冠盖度为10-30%的疏林地。	20
草地	天然草本植被覆盖，且盖度大于10%的土地，包括草原、草甸、稀树草原、荒漠草原，以及城市人工草地等。	30
灌木地	灌木覆盖且灌丛覆盖度高于30%的土地，包括山地灌丛、落叶和常绿灌丛，以及荒漠地区覆盖度高于10%的荒漠灌丛。	40
湿地	位于陆地和水域的交界带，有浅层积水或土壤过湿的土地，多生长有沼生或湿生植物。包括内陆沼泽、湖泊沼泽、河流洪泛湿地、森林/灌木湿地、泥炭沼泽、红树林、盐沼等。	50
水体	陆地范围液态水覆盖的区域，包括江河、湖泊、水库、坑塘等。	60
苔原	寒带及高山环境下由地衣、苔藓、多年生耐寒草本和灌木植被覆盖的土地，包括灌丛苔原、禾本苔原、湿苔原、高寒苔原、裸地苔原等。	70
人造地表	由人工建造活动形成的地表，包括城镇等各类居民地、工矿、交通设施等，不包括建设用地内部连片绿地和水体。	80
裸地	植被覆盖度低于10%的自然覆盖土地，包括荒漠、沙地、砾石地、裸岩、盐碱地等。	90
冰川和永久积雪	由永久积雪、冰川和冰盖覆盖的土地，包括高山地区永久积雪、冰川，以及极地冰盖等。	100

2.  ArcGIS提取流域边界

！！！！这里注意：

a  以下每一步导出数据之前，都要提前建立好一个文件夹用于接收保存导出的数据，以免数据保存的乱七八糟，影响后续数据处理。

类似于上图这样

b 导出数据时，如果报错，进行如下设置，重新设置导出，在导出之前，点击环境，打开环境设置，把并行处理这里填为0。——————这一步很重要，不然渲染填洼文件可能会失败，后续处理没设置的话，也有可能会失败。

c 我用的是中文版的ArcGIS，可能一些步骤在翻译上存在差别，在按照步骤操作的时候，自行进行中英文替换哈~

举例，watershed 有的翻译成分水岭，有的翻译成集水区，是一个意思哈，点他就行。

2.1 影像拼接

首先导入DME数据

注意！！！ ---------------在拼接之前要把所有DEM数据放到一个文件夹里面，不然很麻烦的。

导入之后就会变成这样，但是需要将图像拼接，如果只有一张图就忽略这一步。

拼接开始！！！！！！

可以参考这篇博客：ArcGIS拼接多张dem影像（含裁剪、dem负值处理）_arcgis dem拼接-CSDN博客

拼接好后，未拼接的数据可以删除了，只保留DEM_mosaic.tif这个图层。

2.2 填洼设置

拼接好后就要开始填洼啦！！！

步骤：ArcToolBox——>spatial Analyst 工具——>水文分析——>填洼 （fill）

然后点击环境，打开环境设置，把并行处理这里填为0

这一步很重要，不然渲染填洼文件可能会失败，后续处理没设置的话，也有可能会失败。

填洼完成后，如下图所示，会自动生成一个新的fill

2.3 流向设置

关于ArcGIS流向设置，涉及到的三种算法，我的ArcGIS是10.+的版本，老版本的可能只有D8算法，新版本还增加了MFD和DIFD两种算法，具体介绍可以参考这个介绍，不感兴趣的可以跳过此链接：流向算法的一些总结 - 简书

按照以下说明导出流向文件就可以了

步骤：ArcToolBox——>spatial Analyst 工具——>水文分析——>流向

导出后结果，如下所示：

2.4 流量设置

步骤：ArcToolBox——>spatial Analyst 工具——>水文分析——>流量

注意：这一步需要等比较久，根据研究区域大小，时间不定。研究区域越大，出图的时间越长，这一步需要耐心等待。

数据处理完后的效果，得到流量数据（flow.tif）：

2.5 将流量数据栅格化

步骤：ArcToolBox——>spatial Analyst 工具——>地图代数——>栅格计算器

这一步是要把前面获得的流量数据进行栅格化

以下就获得了raster5000效果图

2.6 河网链接

步骤：ArcToolBox——>spatial Analyst 工具——>水文分析——>河流链接

成功之后效果图，生成了raster_link文件

2.7 捕捉倾斜点

        前面生成的河网链接（raster_link），是非常重要的一步，是用于后面提取流域面积的关键步骤。在设置倾斜点之前，可以根据流向图层（flow_direction）放大筛选最下游的点作为倾斜点，然后如果下游有多条支流，可以设置多个倾斜点。流域边界提取的原理就是，根据自己在下游设置的倾斜点以上的部分，包括以上的干流以及支流所涵盖的流域边界都会通过算法计算在内。

第一步：先设置点图层

然后这样操作

点图层就创建好了

接下来通过流向数据图层，找河流最下游点，用于绘制倾斜

步骤和思路大致如下图：

注意：创建点的时候，注意以下步骤

编辑器——>编辑窗口——>创建要素——>构造工具——>选择点——>结束编辑

可以看到，我基本上捕捉了下游所有的河网，并把它命名为point_qxd(点图层）

接下来就正式开始捕捉倾斜点啦！！！！！！！！

步骤：ArcToolBox——>spatial Analyst 工具——>水文分析——>捕捉倾泻点

2.8 分水岭（流域边界提取）

这一步实际上流域边界的提取了，有些软件中文名叫集水区，和分水岭是一样的

步骤：ArcToolBox——>spatial Analyst 工具——>水文分析——>分水岭

界面长这样，按照提示输入数据就可以了。

这个就是提取好的流域边界，但是这是栅格数据，还需要把他转成shp数据作为图层使用。

转换步骤：ArcToolBox——>转换工具——>由栅格转出——>栅格转面

流域边界效果如上，是shp格式，方便其他软件也可以使用

2.9 计算流域边界面积

步骤：打开文件属性表，然后添加字段

注意：shp图层要改为投影坐标系，而不是地理坐标系，才能计算出面积。

PCS是投影坐标系，才可以计算出面积。

3. ArcGIS提取土地利用类型

接下来我们就可以用上次的来进行土地利用类型处理

3.1 导入GLC地图

用相同的方法把你提取好的流域边界也导入进来

有需求的话，可以把几张GLC的图片进行影像拼接！！

3.2 按掩膜提取（流域边界提取）

如果有多张的GLC的话，建议先把多张影响拼接起来，这个在2.1有提到。

步骤：ArcToolBox——>空间分析工具——>提取分析——>按掩膜提取

看效果图，后续还要通过颜色来调整

颜色可以按照这个来调整，这是引用知乎一个老哥的图片，以下是原文链接

【ArcGIS教程】（101）30米精度全球地表覆盖分类数据免费分享（2000 年 2010年 2020年） - 知乎

emmmmm，这个配色感觉不太好看，可以自己看文献调配色。

3.3 注意事项和一些可能会遇到的问题

3.3.1 提取图层之后，导入到新的arcgis可能会只显示一个颜色

这个时候，在内容列表一栏，右键单击图层（上图蓝色选中部分）

步骤：属性----->符号系统----->值字段----->选择value----->点击确定就可以了

3.3.2 关联属性表

这是从网上找到的对于映射关系，建立一张csv表，保存下来，用于arcgis映射使用

一样是在内容列表里面右键单击图层

步骤：连接和关联---->连接

步骤：属性----->符号系统----->值字段----->选择连接好类型----->点击确定就可以了

通过连接表格，就可以换备注

后续可以靠count的数值，计算各部分土地利用类型的占比

计算各部分土地利用类型的占比

4、站位图绘制

在ArcGIS中变换图层坐标系是一个常见的操作，具体步骤如下：
1. 查看当前坐标系统：首先需要确定图层当前使用的坐标系统。这可以通过打开ArcCatalog，找到数据集，然后右键点击数据集，选择“图层属性”，在“XY坐标系”选项卡中查看 。
2. 使用ArcToolbox进行坐标转换：打开ArcToolbox工具箱，找到“数据管理工具”下的“投影和变换”。如果处理的是影像数据，选择“栅格”目录下的“投影”；如果是矢量数据，则选择“要素”目录下的“投影” 。
3. 选择数据和目标坐标系统：在“投影”界面中，选择需要转换坐标的数据集，然后选择转换后的坐标系统。点击“选择”来确定新的坐标系统，并点击“确定”以启动转换过程 。
4. 处理多个图层的坐标系统问题：如果加载了多个图层，发现某些图层显示不正确，可能是因为坐标系统不一致。可以通过右键点击图层，选择“属性”，然后在“坐标系”选项中更改为指定的坐标系 。
5. 自定义地理坐标转换：如果需要将地理坐标系转换为投影坐标系，例如从WGS 84转换为CGCS 2000，可以在ArcToolbox中使用“自定义地理坐标转换”工具。选择输入地理坐标系和输出投影坐标系，并根据需要选择适当的转换方法 。
6. 导出转换后的图层：转换完成后，新生成的图层将显示在内容列表中。如果需要使用转换后的数据，可能需要打开一个新的ArcMap窗口，并导入转换后的图层，以便正确读取数据 。
请注意，坐标转换可能涉及到复杂的地理信息系统（GIS）概念，如地理坐标系和投影坐标系的区别，以及可能需要的坐标转换参数。务必确保在进行坐标转换前，已经充分理解了数据的当前坐标系统和目标坐标系统 。

经纬度数据需要是csv格式，或者97-03格式的xls格式（注意不是xlsx）

导入x，y轴数据，然后再把数据保存为shp格式