基于数字高程模型（DEM）栅格数据，采用三阶反距离平方权差分法计算坡度与坡向。支持坡度百分比、度数与弧度三种输出格式，适用于滑坡评估、土壤侵蚀建模与工程勘察。
输入 3×3 DEM 高程矩阵与栅格分辨率，自动计算中心像元坡度与坡向
三阶反距离平方权差分，ESRI 与 GDAL 默认算法，精度最高
三阶拉格朗日插值，适用于平滑地形，减少噪声影响
四邻域平均差分，计算快速，适用于粗略估算
按行优先编号：z₁~z₃ 为第一行，z₄~z₆ 为第二行，z₇~z₉ 为第三行
像元东西方向边长，通常等于 DEM 空间分辨率
像元南北方向边长，通常与 X 方向相同
深入理解数字地形分析的核心指标及其在地理学、生态学与工程学中的应用
坡度是地表单元陡缓程度的度量，定义为坡面与水平面之间的夹角。在 DEM 栅格分析中，坡度通过中心像元与周围 8 邻域像元的高程差计算得到，反映了地表在该点的最大倾斜程度。
坡度有三种常用表达形式：度数（°）直接表示倾斜角度；百分比（%）等于 tan(坡度) × 100，常用于道路工程与土地利用规划；弧度（rad）用于数学计算与物理模型。
坡度是土壤侵蚀模型（如 USLE/RUSLE）、滑坡敏感性评价、农业适宜性分析与交通选线中的关键输入参数。中国黄土高原地区研究表明，坡度 > 15° 的坡耕地土壤侵蚀模数可达平坦耕地的 5–10 倍。
坡向是坡面法线在水平面上投影所指的方向，以正北为 0°，顺时针旋转至 360°。坡向决定了地表接收太阳辐射的时空分布，进而影响土壤温度、植被类型、积雪消融速率与微气候特征。
在北半球中纬度地区，南坡（135°–225°）接收的太阳辐射最多，通常植被覆盖较好、土壤温度较高；北坡（315°–45°）则相对阴冷湿润，常发育针叶林或草甸。在高山地区，坡向差异可导致明显的植被垂直带谱错位。
坡向也是雪崩风险评估、林火蔓延模拟与太阳能电站选址的重要指标。例如，中国青藏高原东南缘的冰川末端朝向以东南坡为主，与夏季季风带来的降水及辐射条件密切相关。
目前主流的 DEM 坡度坡向提取算法包括 Horn 三阶差分法、Zevenbergen-Thorne 法、Sharpnack-Akin 法与 Evans-Young 法。不同算法对噪声的敏感度和地形适应性存在差异。
| 算法 | 权重模式 | 适用场景 | 精度 |
|---|---|---|---|
| Horn 法 | 反距离平方权 | 通用，ESRI 默认 | 高 |
| Zevenbergen-Thorne | 拉格朗日插值 | 平滑地形 | 中高 |
| 简单差分 | 四邻域平均 | 快速估算 | 中 |
| Evans-Young | 二次曲面拟合 | 复杂地形 | 高 |
Horn 法（1981）通过对角线方向赋予 √2 倍的权重，有效减少了 DEM 高程误差对结果的影响，被 ArcGIS、GDAL 与 GRASS GIS 作为默认算法广泛采用。
DEM 分辨率、高程精度与地形复杂度直接影响坡度坡向计算结果的可靠性。研究表明，当 DEM 分辨率从 30m 降低至 90m 时，平均坡度可低估 15%–30%，尤其在沟谷深切区域误差更为显著。
主要误差来源包括：高程量化误差（DEM 垂直精度通常为 1–15m）；地形综合误差（栅格分辨率不足以表达微地形）；算法选择误差（不同算法在凸/凹地形处的系统性偏差）。
建议在实际应用中：优先使用分辨率 ≤ 30m 的 DEM（如 ASTER GDEM、SRTM1、ALOS World 3D）；对结果进行 3×3 或 5×5 中值滤波以抑制噪声；在陡坡区域（>30°）采用高精度 LiDAR DEM 替代卫星测高产品。
坡度是耕地适宜性评价的核心指标。中国《耕地质量等级》国家标准（GB/T 33469-2016）规定，>25° 坡耕地应逐步退耕还林还草。坡向则影响作物类型选择——阳坡适宜种植喜温作物（玉米、棉花），阴坡适宜耐寒作物（马铃薯、燕麦）。
坡度与坡向是滑坡敏感性区划的关键因子。统计表明，中国西南山区 70% 以上的大型滑坡发生在 20°–45° 坡度范围内。坡向通过控制岩体风化速率与地下水渗流方向，间接影响滑坡发育——顺向坡（岩层倾向与坡向一致）的失稳概率显著高于逆向坡。
在 USLE/RUSLE 土壤侵蚀模型中，坡度因子（LS）与坡长、坡度直接相关。坡度每增加 1°，径流流速约增加 3%–5%，土壤侵蚀量呈指数增长。坡向通过影响植被覆盖度与土壤含水量，间接调节侵蚀强度。
坡向决定了坡面接收的太阳辐射总量。在纬度 35°N 地区，南向 30° 坡面年太阳辐射量可比水平面高 12%–18%，是分布式光伏电站选址的优势区域。同时，坡向差异也是山区"阴阳坡"气候分异与植被分带的主因。