ArcGIS土壤侵蚀模数坡度因子S

一、DEM数据下载及处理

1. 下载研究区域高程（地理空间数据云 (gscloud.cn)）

2. 栅格镶嵌

• 确定栅格文件的波段数和像素类型。（右键单击内容列表，然后依次单击属性和源选项卡。）输入必须具有相同的波段数和相同的位深度。

• 导航至ArcToolbox>数据管理工具>栅格>栅格数据集>镶嵌至新栅格工具。

二、坡度因子S计算

1.公式

S=10.8sinθ+0.03    θ<5.2°

S=16.8sinθ-0.5    5.2°≤θ≤10.2°

S=21.91sinθ-0.96   θ＞10.2°

2.在ArcMAP中加载DEM数据，选中系统工具箱——SpatialAnalystTools——表面分析——坡度工具。

3.确定z因子

        z因子是一个转换因子，当输入表面的垂直坐标（或高程）单位与水平坐标(x,y)单位不同时，可使用z因子调整垂直坐标（或高程）单位的测量单位。如果垂直单位没有被改成水平单位，则表面工具的结果将不会正确。

        在计算输出表面时，将使用z因子乘以输入表面的z值。如果x,y单位与z单位完全相同（例如都是以英尺为单位），则z因子为1。这是z因子的默认值。再例如，如果单位不同（垂直的z单位是英尺而水平的x,y单位是米），则应使用z因子0.3048将z单位从英尺转换为米（1英尺=0.3048米）。

        当输入栅格位于球面坐标系中，例如采用十进制度时，z因子的正确使用尤为重要。如果输入表面栅格不在投影坐标系中，则从山体阴影中得到的输出看起来通常会很奇怪。这是由于在水平地面单位与高程z单位之间的测量存在差异。由于经度的长度随着纬度而变化，因此需要为该纬度指定一个适当的z因子。

即：

1）如果DEM是投影坐标系，单位是米，DEM的高程单位也是米，则Z因子默认为1即可；

2）若DEM是地理坐标系，平面单位是度，高程单位是米，所以要根据所以地区纬度设置Z因子；
或者将地理坐标系转化为投影坐标系，转化投影之后单位就变成了米，此时Z因子默认为1 就可以。

arcgis计算坡度（需注意DEM是地理坐标系还是投影坐标系）_arcgis坡度单位-CSDN博客

如果x,y单位是十进制度而z单位是米，则适用于特殊纬度的正确的z因子如下：

4、系统工具箱——Spatial Analyst Tools——地图代数——栅格计算器，输入公式：

Con("dem-坡度"<5.2,10.8*Sin("dem-坡度"*3.1415926/180)+0.03,Con("dem-坡度">10.2,21.91*Sin("dem-坡度"*3.1415926/180)-0.96,16.8*Sin("dem-坡度"*3.1415926/180)-0.5))