【ArcGIS】地形特征描述

【ArcGIS】描述饮用水源保护区地形特征

问题描述

利用DEM生成不规则三角网，并基于DEM和TIN分别生成坡度并对坡度进行分级：0-2°、2-6°、6-15°、15-25°和25°以上；
分析不同分区单元的平均坡度、最大坡度、最小坡度、以及不同级别坡度面积的情况。

软件及数据

软件：

ArcGIS 10.8

案例数据：

操作步骤

1. 导入源文件：首先，将“DEM.tif”文件导入图层，点击属性表查看投影坐标系，打开工具箱，点击“数据管理工具”，找到“投影栅格”，将“DEM.tif”作为输入元素进行重投影，得到“DEM_project.tif”文件；

2. 生成TIN：打开工具箱，在“3D Analyst工具”中找到“栅格转TIN”，将“DEM_project.tif”文件作为输入元素，转换生成“DEM_to_TIN.TIN”文件；

3. DEM生成坡度文件：打开工具箱，在“3D Analyst工具”中找到“栅格表面”，点击“坡度”工具，将“DEM_project.tif”文件作为输入元素，生成坡度文件“Slope_of_DEM”；

4. TIN生成坡度文件：打开工具箱，在“3D Analyst工具”中找到“表面三角化”，点击“表面坡度”工具，将“DEM_to_TIN.TIN”文件作为输入元素，坡度单位选为“DGREE”，生成坡度文件“Slope_of_TIN.shp”文件；

5. DEM坡度分级：打开工具箱，在“3D Analyst工具”中找到“栅格重分类”，点击“重分类”工具，将“Slope_of_DEM”文件作为输入元素，系统读取默认分类后，点击右侧“分类”，依次输入“2”“6”“15”“25”，将原文件重分类为5类，得到“Reclass_dem”文件；

6. TIN坡度分级：打开工具箱，找到“转换工具”，选择“要素转栅格”，将“Slope_of_TIN”文件作为输入要素，生成“TIN_to_Rast”栅格文件，在“3D Analyst工具”中找到“栅格重分类”，点击“重分类”工具，将“Slope_of_DEM”文件作为输入元素，系统读取默认分类后，点击右侧“分类”，依次输入“2”“6”“15”“25”，将原文件重分类为5类，得到“Reclass_dem”文件。

分析不同分区单元情况：

1. 首先通过合并操作将“饮用水源保护区_project”和“饮用水源保护区_buffer”合并成同一个文件，命名为“饮用水源保护区_merge”，然后打开“饮用水源保护区_merge”文件属性表，调整不同分区颜色、图例名称等；

2. 通过掩膜：打开“Spatial Analyst工具”，选择“提出分析”中的“按掩膜提取”工具，将“slope_of_dem”作为输入栅格（即需要裁切的栅格数据），将“饮用水源保护区_merge”作为要素掩膜数据（即用以裁切的栅格或矢量数据），生成“mask_slope”文件，为后续进行不同分区坡度分析做准备；

3. 打开“Spatial Analyst工具”，选择“区域分析”中的“以表格显示分区统计”工具，将“饮用水源保护区_merge”作为要素区域数据，将“slope_of_dem”作为输入赋值栅格，选择区域字段为“FID”方便后续进行对应，生成不同分区的坡度统计表。

由表可知：

水域一级保护区（1号）的平均坡度为24.205867、最大坡度为54.317135、最小坡度为0、坡度面积为44175 m²。

水域一级保护区（2号）的平均坡度为0.09856、最大坡度为32.508808、最小坡度为0、坡度面积为469425 m²。

水域二级保护区的平均坡度为0.304726、最大坡度为56.729576、最小坡度为0、坡度面积为4619900 m²。

陆域二级保护区的平均坡度为20.397699、最大坡度为81.623764、最小坡度为0、坡度面积为29345875 m²。

区外500米的平均坡度为20.75715、最大坡度为69.545074、最小坡度为0、坡度面积为16073825 m²。

区外500-1000米的平均坡度为18.123636、最大坡度为63.854946、最小坡度为0、坡度面积为16098150 m²。