摘要:地形是影響土壤侵蝕的重要因子,,在侵蝕估算模型中常用坡度和坡長(zhǎng)(LS)來衡量,,在大區(qū)域上常基于數(shù)字高程模型(DEM)提取,。SRTM作為大區(qū)域尺度上質(zhì)量高,、易獲取的高程數(shù)據(jù),在全球土壤侵蝕評(píng)價(jià)中得到廣泛應(yīng)用,。但現(xiàn)有地形因子提取算法要求高程和柵格單元的單位(通常為m)一致,,需對(duì)SRTM進(jìn)行坐標(biāo)變換才能使用。針對(duì)大區(qū)域上SRTM坐標(biāo)轉(zhuǎn)換時(shí)間開銷大的問題,,提出了一種直接基于SRTM提取地形因子的算法(LSA-SRTM),。該算法利用地理坐標(biāo)系下的經(jīng)緯度信息計(jì)算柵格單元長(zhǎng)度及單元坡長(zhǎng),結(jié)合最陡坡降策略獲得坡度及流向,,進(jìn)而提取匯水面積,,根據(jù)坡度設(shè)置坡度截?cái)帱c(diǎn),根據(jù)匯水面積閾值設(shè)置溝道截?cái)帱c(diǎn),,經(jīng)正,、反遍歷后獲得累積坡長(zhǎng),采用CSLE的分段公式計(jì)算LS因子,。以Himmelblau-Orlandini數(shù)學(xué)曲面,、5個(gè)典型樣區(qū)的1″SRTM作為數(shù)據(jù)源,將LSA-SRTM,、投影坐標(biāo)系下的LS算法(LSA-DEM)與手工測(cè)量的結(jié)果進(jìn)行對(duì)比,。LSA-SRTM方法與測(cè)量值在數(shù)學(xué)曲面和典型樣區(qū),坡長(zhǎng)的R2分別為0.8552,、0.7788,、0.7269、0.7024,、0.6909,、0.7255,LS因子的R2分別為0.8907,、0.8209,、0.8213,、0.7142、0.7145,、0.8212,。在運(yùn)行時(shí)間方面,LSA-SRTM方法具有較高的效率,。結(jié)果表明,,LSA-SRTM算法計(jì)算精度、效率更高,,可為大區(qū)域地形因子提取的研究提供支撐,。