（158点）数据集制作1——点转shp+地形特征提取

搞笑的秀儿

已于 2025-03-18 09:10:58 修改

阅读量1k

点赞数 30

分类专栏：数据集制作文章标签： arcgis

于 2025-03-17 16:49:16 首次发布

本文链接：https://blog.csdn.net/m0_52460096/article/details/146308568

版权

数据集制作专栏收录该内容

8 篇文章

订阅专栏

一、158灾害点转.shp

158个灾害点，度分秒——>小数，xlsx——>.csv，XY表转点

二、DEM提取

根据DEM提取高程、坡向、坡度、平面曲率、剖面曲率、地表起伏度、地形粗糙度、径流强度指数（SPI）、地形湿度指数（TWI）、地形崎岖指数（TRI）

2.1 浙江省DEM获取

数据名称	格式	来源	网址	注意
30mDEM	.tif	地理空间数据云	地理空间数据云	下载的为.img，需要拉入arcgis建立金字塔
浙江省行政区域	.shp	国家地理信息公共服务平台	天地图 https://cloudcenter.tianditu.gov.cn/administrativeDivision/	下载的为geojson格式，需要转换

批量下载解压显示文件损坏？？？19个文件，一个个下载……

下载完后，先把19块拼在一起，然后按照浙江省边界进行裁剪：【教程可见12.18-19 GIS学习——基于GIS和DEM的地形地貌特征提取与分析 + 多值提取到点_径流强度指数-CSDN博客】

2.2 地形特征提取

1、投影转换

感觉原来是地理坐标系的话，应该先转换成投影坐标系，然后再进行坡度、坡向等地形特征的提取。
WGS1984——>51N 如下图

2、等高线提取

采用【等值线】工具，结合DEM数据进行等值线分析

3、山体阴影

4、坡度

5、剖面曲率

根据坡度计算破面曲率，输入是坡度！！

曲率分为平面曲率和剖面曲率。

剖面曲率：是对一个研究区的dem进行坡度分析后的结果再一次进行坡度分析。

平面曲率：是对研究区的dem进行坡向分析后的结果再一次进行坡度分析。

6、坡向

7、平面曲率

根据坡向计算平面曲率，注意输入是坡向！！！

曲率分为平面曲率和剖面曲率。

剖面曲率：是对一个研究区的dem进行坡度分析后的结果再一次进行坡度分析。

平面曲率：是对研究区的dem进行坡向分析后的结果再一次进行坡度分析。

8、地形起伏度

采用【焦点统计】，计算DEM数据，11×11窗口内的最大值

局部最大高程计算：

局部最小高程计算：

栅格计算器 ——地形起伏度计算：max-min

9、地形粗糙度

采用【栅格计算器】，对坡度进行弧度制转换

1/Cos("坡度.tif"*3.1415926/180)

10、径流强度指数（SPI）

SPI值越高，表明径流集中可能导致土壤侵蚀。SPI公式形式较多，常见公式如下：

SPI=SCA * Tan(slope)

SPI=ln[SCA * Tan(slope)]

SPI=ln[SCA * Tan(slope)*100]

其中，SCA和slope意义与上述相同。本文基于公式2进行计算。所以要先计算SCA和slope

ArcToolbox→Spatial Analyst Tools (Spatial Analyst工具)→Hydrology (水文分析)→Fill (填洼)

要计算SCA：填洼——流向——流量——栅格计算单位面积汇流量

填洼（并行处理设计为0）

流向（并行处理设计为0）

流量（并行处理设计为0）

SCA栅格计算器
公式为：

Con("流量_丽水" ==0,1,"流量_丽水") *900 / Con("流向_丽水" ==1,30, Con("流向_丽水" ==4,30,Con("流向_丽水" ==16,30, Con("流向_丽水" ==64,30, Con("流向_丽水" ==2,30 * SquareRoot(2) , Con("流向_丽水" ==8,30 * SquareRoot(2), Con("流向_丽水" ==32,30 * SquareRoot(2), Con("流向_丽水" ==128,30 * SquareRoot(2)))))))))

其中，900为栅格大小，即栅格分辨率的平方，需要根据自己的数据像元大小进行修改（本文为30m*30m,因此为900），且当"Flow Accumulation"为0时没意义，至少应为1。对于流向为1、4、16、64 的栅格，除以栅格尺寸，这里是30；对于流向为2、8、32、128的栅格，除以栅格尺寸的√2倍，这里是√2倍的30。