坡长坡度因子的生成过程

坡长坡度因子的生成过程如下：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*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_yinziPow([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.img4、坡长坡度因子LS=L*S根据rastercalculator中的公式，得到ls