坡長(zhǎng)坡度因子的生成過(guò)程

《坡長(zhǎng)坡度因子的生成過(guò)程》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。

1、坡長(zhǎng)坡度因子的生成過(guò)程如下：1、要估算坡長(zhǎng)因子，可以利用ARCGIS中的水文模塊提取流水累積量。流水累積量的計(jì)算，要先建立無(wú)洼地的DEM。（見(jiàn)水文分析過(guò)程）利用jx_srtm數(shù)據(jù)，先生產(chǎn)任意流向圖：flowdir，再尋找洼地：sink，并計(jì)算洼地流域貢獻(xiàn)范圍：sink_watersh。再計(jì)算每個(gè)洼地的最低高程和出水口高程：zonalmin，zonalmax，并可以求出洼地深度：sinkdep=zonalmax–zonalmin。填充洼地：fill_jx_srtm。并計(jì)算無(wú)洼地的dem的水流流向：flowdir_fill1。然后計(jì)算流水累積量：flowacc。2、根據(jù)坡

2、長(zhǎng)因子公式（參考文獻(xiàn)：土壤侵蝕研究中的坡長(zhǎng)因子評(píng)價(jià)問(wèn)題）：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)根據(jù)rastercalculator中的公式，計(jì)算出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

3、)根據(jù)rastercalculator中的公式，計(jì)算出m(3)根據(jù)rastercalculator中的公式，計(jì)算出l_yinziPow([flowacc]*87.72085638/22.13,[m])對(duì)得到的l_yinzi（l_yinzi1.img）進(jìn)行去除異常值：利用ERDAS中的圖像增強(qiáng)工具，對(duì)異常值進(jìn)行平滑處理。這里，我們選擇以下的運(yùn)算模板：因此，得到圖：l_yinzi_convolution.img。再利用建模工具，對(duì)l_yinzi1.img和l_yinzi_convolution.img求值。（參考文獻(xiàn)：GIS支持下的土壤侵蝕量估算_以江西省泰和縣灌溪鄉(xiāng)為例

4、）得到圖：l_yinzi_model.img，為最后的L因子圖。1、根據(jù)坡度因子公式：S=10.8*sinθ+0.03(θ<5)，16.8*sinθ-0.50(5<=θ<10)，21.9*sinθ-0.96(θ>=10)利用建模工具，得到s_yinzi.img1、坡長(zhǎng)坡度因子LS=L*S根據(jù)rastercalculator中的公式，得到ls