坡长坡因子的生成过程.docx

坡长坡度因子的生成过程如下：
要估算坡长因子，可以利用ARCGIS中的水文模块提取流水累积量。
流水累积量的计算，要先建立无洼地的DEM。（见水文分析过程）
利用jx_srtm数据，先生产任意流向图：flowdir，再寻找洼地：sink，并计算洼地流域贡献范围：sink_watersh。再计算每个洼地的最低高程和出水口高程：zonalmin，zonalmax，并可以求出洼地深度：sinkdep=zonalmax–zonalmin。
填充洼地：fill_jx_srtm。并计算无洼地的dem的水流流向：flowdir_fill1。然后计算流水累积量：flowacc。

根据坡长因子公式（参考文献：土壤侵蚀研究中的坡长因子评价问题）：
L=(λ/22.13)m,λ=flowacc*cellsize
其中：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)根据rastercalculator中的公式，计算出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)根据rastercalculator中的公式，计算出m

(3)根据rastercalculator中的公式，计算出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因子图。

根据坡度因子公式：S=10.8*sinθ+0.03(θ<5)，
16.8*sinθ-0.50(5<=θ<10)，
21.9*sinθ-0.96(θ>=10)
利用建模工具，得到s_yinzi.img

坡长坡度因子LS=L*S

根据rastercalculator中的公式，得到ls