坡长坡度因子的生成过程

坡长坡度因子的生成过程如下：

1、要估算坡长因子，可以利用ARCGIS中的水文模块提取流水累积量。

流水累积量的计算，要先建立无洼地的DEM。（见水文分析过程）

利用jx_srtm 数据，先生产任意流向图：flowdir，再寻找洼地：sink，并计算洼地流域贡献范围：sink_watersh 。再计算每个洼地的最低高程和出水口高程：zonalmin ，zonalmax，并可以求出洼地深度：sinkdep = zonalmax – zonalmin。

填充洼地：fill_jx_srtm。并计算无洼地的dem的水流流向：flowdir_fill1。然后计算流水累积量：flowacc。

2、根据坡长因子公式（参考文献：土壤侵蚀研究中的坡长因子评价问题）：

L=(λ/22.13) m, λ=flowacc * cell size

其中：m=n/(1+n)，n=(sinθ/0.0896)/(3.0*sin0.8θ+0.56)

这里的θ用“jx_slope.img *3.1415926/180”替代

依次得到图：n、m、l_yinzi

(1) 根据raster calculator 中的公式，计算出n

(Sin([jx_slope.img] * 3.1415926 / 180) / 0.0896) / (3.0 * Pow(Sin([jx_slope.img] *

3.145926 / 180), 0.8) + 0.56)

(2) 根据raster calculator 中的公式，计算出m

(3) 根据raster calculator 中的公式，计算出l_yinzi

Pow([flowacc] * 87.72085638 / 22.13, [m])

对得到的l_yinzi（l_yinzi1.img）进行去除异常值：

利用ERDAS中的图像增强工具，对异常值进行平滑处理。

这里，我们选择以下的运算模板：

因此，得到图：l_yinzi_convolution.img。

再利用建模工具，对l_yinzi1.img和l_yinzi_convolution.img求值。（参考文献：GIS 支持下的土壤侵蚀量估算_以江西省泰和县灌溪乡为例）

得到图：l_yinzi_model.img，为最后的L因子图。

3、根据坡度因子公式：S= 10.8*sinθ+ 0.03 (θ<5)，

16.8*sinθ- 0.50 (5<=θ<10)，

21.9*sinθ- 0.96 (θ>=10)

利用建模工具，得到s_yinzi.img

4、坡长坡度因子LS=L*S

根据raster calculator 中的公式，得到ls