基于arcpy的rusle模型ls因子計(jì)算

1、114基于基于ArcPy的RUSLE模型模型LS因子計(jì)算因子計(jì)算（FromAMLarcgising）By：半荒漠豬毛菜屬飛天小豬于2011年將Hickey計(jì)算RUSLE模型LS因子的AML代碼移植為python腳本，基于平臺為ArcGIS9.3。目前，ArcGIS最新版本已到10.3，且自10.0之后出現(xiàn)ArcPy包，內(nèi)置模塊與函數(shù)與之前版本差異較大。為方便在現(xiàn)版本中運(yùn)行，參考Hickey原文算法的基礎(chǔ)上，將飛天小豬代碼轉(zhuǎn)譯為10.X腳

2、本形式。附以下為“飛天小豬”版使用說明：附以下為“飛天小豬”版使用說明：越陡越長的坡比地面上產(chǎn)生更高的流速，但是土壤侵蝕對S（坡度）比對L（坡長）有更強(qiáng)的敏感度。不規(guī)則和分割的坡型對RUSLE的影響沒有體現(xiàn)在本程序中。LS計(jì)算的算法時(shí)基于McCool等人在1997年的研究成果，校正了坡長的水平投影；這對GIS中通過格網(wǎng)或地圖（xy）而不是字段（xyz）計(jì)算坡長很有用。程序中坡長的計(jì)算從高點(diǎn)（highpoints如平緩的脊頂ridgeto

3、ps）到低點(diǎn)（lowpoints，比如，流域傾點(diǎn)watershedpourpoint或其他出口點(diǎn)outlet）。行政上定義的流域可能并不合適，除非其也是水文上定義的集水區(qū)域。計(jì)算LS因子的理想輸入數(shù)據(jù)是已經(jīng)裁剪過或者已經(jīng)擴(kuò)大到感興趣的區(qū)域及任何相關(guān)的集水區(qū)的DEM數(shù)據(jù)。為了避免任何由于尺度導(dǎo)致的邊界影響，最大范圍應(yīng)該稍微比目標(biāo)區(qū)域大一點(diǎn)。確保DEM的高程單位和水平距離單位是一致的（默認(rèn)是米）。L和S因子計(jì)算的輸出數(shù)據(jù)需要仔細(xì)檢查以確保和

4、DEM相符以及沒有重大格式問題。如果在處理中因?yàn)楦↑c(diǎn)型格式出現(xiàn)問題，截短或取整或許可以解決。314程序如下：imptarcpyfromarcpyimptfromarcpy.saimpt#即不需要arcgising中導(dǎo)入，用arcpy包即可arcpy.env.wkspace=“E:temp“dem_input=“E:tempDEM.tif“#輸入柵格數(shù)據(jù)wshed=“E:tempBoundary.shp“#輸入流域邊界數(shù)據(jù)demunits

5、=“meters“scf_lt5=0.7scf_ge5=0.5#定義信息提示函數(shù)defsendmsg(msg):printmsgarcpy.AddMessage(msg)#定義一個(gè)函數(shù)，輸入字符型坐標(biāo)、cellsize、倍數(shù)，返回平移后的字符型坐標(biāo)值，目的為保留原始小數(shù)位數(shù)不變defStoS(scellsizemult):stri=s.split(.)inte=float(stri[0])multcellsizereturnstr(in

6、t(inte)).stri[1]#可覆蓋文件arcpy.env.overwriteOutput=1#判斷輸入DEM數(shù)據(jù)的水平和垂直方向的單位是否一致ifdemunits==Nonedemunits.strip()==““:demunits=“meters“sendmsg(“使用默認(rèn)單位：meters“)elifdemunits!=“meters“demunits!=“feet“:demunits=“meters“sendmsg(“DEM