地形因子提取

一、基本地形属性

1.1 总曲率、剖面曲率、平面曲率、水平曲率

基于dem数据，采用ArcGIS中工具箱->Spatial Analyst Tools->表面分析->曲率，来获取剖面曲率、平面曲率和和总曲率，另外采用栅格计算器（raster calculator）得到水平曲率（坡度×平面曲率）。

 二、派生地形属性

2.1 收敛指数CI

（1）加载DEM数据。将需要处理的栅格文件（dem.tif）直接拖拽至软件灰色窗体。

（2）打开【Terrain Analysis】—【Compound Analyses】—【Basic Terrain Analysis】工具。设置“Grid System”和“Elevation”参数，点击“Okay”，以下所有参数都会输出。

（3）导出CI数据，将其保存为TIFF格式。打开【Import/Export】—【GDAL/OGR】—【 Export GeoTIFF】工具。设置“Grid System”和“Elevation”参数，点击“Okay”，以下所有参数都会输出。设置输出位置，选择存储的文件夹。

（4）在ArcGIS中打开SPI数据。

2.2 径流强度指数SPI

（1）打开【Terrain Analysis】—【Hydrology】—【Stream Power Index】工具。可以看到需要Slope、Catchment Area两个参数，分别求Slope、Catchment Area。

（2）求取Slope。

方法①：打开【Terrain Analysis】—【Compound Analyses】—【Basic Terrain Analysis】工具，求取Slope（求TWI时已经介绍相关步骤，这里不再赘述）。

方法②：打开【Terrain Analysis】—【Morphometry】—【Slope, Aspect, Curvature】工具，求取Slope（过程略）。

（3）求取Catchment Area。

打开【Terrain Analysis】—【Hydrology】—【Melton Ruggedness Number】工具，求取Catchment Area。

（4）打开【Terrain Analysis】—【Hydrology】—【Stream Power Index】工具，求取SPI.

（5）导出SPI数据，将其保存为TIFF格式。打开【Import/Export】—【GDAL/OGR】—【 Export GeoTIFF】工具。设置输出路径。

（6）在ArcGIS中打开SPI数据。

2.3 坡长因子LS

（1）加载DEM数据。将需要处理的栅格文件（dem.tif）直接拖拽至软件灰色窗体。

（2）打开【Terrain Analysis】—【Compound Analyses】—【Basic Terrain Analysis】工具。设置“Grid System”和“Elevation”参数，点击“Okay”，以下所有参数都会输出。

（3）导出LS数据，将其保存为TIFF格式。打开【Import/Export】—【GDAL/OGR】—【 Export GeoTIFF】工具。设置“Grid System”和“Elevation”参数，点击“Okay”，以下所有参数都会输出。设置输出位置，选择存储的文件夹。

（4）在ArcGIS中打开LS数据。

2.4 地形湿度指数TWI

（单流向算法SFD、多流向算法MFD）——SAGA 9.1.1软件

（1）加载DEM数据。将需要处理的栅格文件（dem.tif）直接拖拽至软件灰色窗体。

（2）打开【Terrain Analysis】—【Compound Analyses】—【Basic Terrain Analysis】工具。设置“Grid System”和“Elevation”参数，点击“Okay”，以下所有参数都会输出。

（3）导出TWI数据，将其保存为TIFF格式。打开【Import/Export】—【GDAL/OGR】—【 Export GeoTIFF】工具。设置“Grid System”和“Elevation”参数，点击“Okay”，以下所有参数都会输出。设置输出位置，选择存储的文件夹。

（4）在ArcGIS中打开TWI数据。

2.5 地形粗糙指数TRI

（1）加载DEM数据。将需要处理的栅格文件（dem.tif）直接拖拽至软件灰色窗体。

（2）打开【Terrain Analysis】—【Morphometry】—【 Terrain Ruggedness Index (TRI)】工具。设置“Grid System”和“Elevation”参数，点击“Okay”，以下所有参数都会输出。

（3）导出TRI数据，将其保存为TIFF格式。打开【Import/Export】—【GDAL/OGR】—【 Export GeoTIFF】工具。设置“Grid System”和“Elevation”参数，点击“Okay”，以下所有参数都会输出。设置输出位置，选择存储的文件夹。

（4）在ArcGIS中打开TRI数据。

2.6 地形位置指数TPI

（1）加载DEM数据。将需要处理的栅格文件（dem.tif）直接拖拽至软件灰色窗体。

（2）打开【Terrain Analysis】—【Morphometry】—【 Terrain Position Index (TPI)】工具。设置“Grid System”和“Elevation”参数，点击“Okay”，以下所有参数都会输出。

（3）导出TPI数据，将其保存为TIFF格式。打开【Import/Export】—【GDAL/OGR】—【 Export GeoTIFF】工具。设置“Grid System”和“Elevation”参数，点击“Okay”，以下所有参数都会输出。设置输出位置，选择存储的文件夹。

（4）在ArcGIS中打开TRI数据。