坡长坡度因子的生成过程.docVIP

坡长坡度因子的生成过程如下： 要估算坡长因子，可以利用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 * 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/ 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 根据raster calculator 中的公式，得到 ls