【文章內(nèi)容簡介】 t和setparameterastext只能用在腳本工具（ArcToolbox或Modelbuilder）中。import arcgisscriptinggp=()lat=float((0))decl=float((1))sunangle=90lat+declazimuth=180if sunangle90: sunangle=180sunangle azimuth=0(2,str(sunangle)) (3,str(azimuth))這段代碼中的“新面孔”：getparameterastext：是Python傳遞給geoprocessor（我們稱之為gp）的一個方法，允許工具提供輸入?yún)?shù)。每次你運行這個工具時，都會看到一個對話框，提示輸入?yún)?shù)，這個方法允許你在接下來的程序中使用。索引0和1指第一個和第二個參數(shù)。setparameterastext：和getparameterastext相反，它傳遞第二個條目（比如str(sunangle)的值）給指定的輸出參數(shù)。前兩個參數(shù)索引為0和1，所以接下來輸出參數(shù)的索引為2和3。然后，我們需要將腳本加進工具（Making a script into a tool），那樣才能在ArcToolbox或ModelBuilder或Command line中使用。如前面提到的那樣，這個腳本只能用于工具，包括輸入/輸出方法是PythonWin不能處理的，但這些是多數(shù)工具必需的。216。 在ArcCatalog里，定位到Python腳本保存的文件夾，最好和數(shù)據(jù)在一個盤里，右鍵新建toolbox。當(dāng)然也可以使用之前創(chuàng)建的toolbox。216。 在ArcToolbox里，右鍵toolbox，選添加scripts，填寫如下圖：l 注意：腳本文件是一個腳本工具的參考！非常重要的一點，你使用腳本創(chuàng)建了一個腳本工具，但是腳本本身并沒有和工具一起保存，腳本工具作為toolbox的一部分保存在“*.tbx”文件中。你也可以右鍵腳本工具點擊“編輯”，出現(xiàn)PythonWin或其他任何IDE窗口，這看起來好像是腳本存在于工具中。所以，記得移動時將腳本工具文件和腳本本身一起拷貝?！跋乱徊健焙笫菂?shù)配置頁面，如下圖設(shè)置各參數(shù)如表格所示：Display NameData typeDirectionDefaultLatitudeDoubleinput0DeclinationDoubleinput0Sun AngleDoubleoutput35AzimuthDoubleoutput300216。 現(xiàn)在可以運行此toolbox了，對話框提示你輸入latitude和declination。工具是否正確運行呢？如果是的話，你應(yīng)該可以看到“Completed”，你可能會看到有一黑色窗體一閃而過，不用擔(dān)心。那么，它是干什么的呢？還記得結(jié)果是輸出兩個數(shù)字參數(shù)，那么，這些數(shù)字哪去了呢？很好的問題，這僅能說明你能創(chuàng)建一個工具，但是不能想ArcToolbox那樣運行。比如輸出一種數(shù)據(jù)，柵格或特征數(shù)據(jù)（.shp）之類的。216。 但是它能在作為Modelbuilder工具正常運行，下面我們將使用它的輸出參數(shù)創(chuàng)建一個hillshade。216。 在你的toolbox中新建一個model，將剛才創(chuàng)建的腳本工具（script tool）拖進來。216。 雙擊工具，輸入?yún)?shù)，初始化之。打開“Sun Angle”和“Azimuth”發(fā)現(xiàn)它們還是默認值，說明此腳本工具還沒有運行。右鍵單擊工具，選擇Run，然后發(fā)現(xiàn)兩個輸出參數(shù)已經(jīng)改變！需要注意的是：latitude范圍是90~90，~。嘗試輸入latitude 70，declination ，你會得到什么？為什么？216。 確保你已經(jīng)得到值域內(nèi)的太陽角和方位角，它們將是構(gòu)建hillshade的輸入?yún)?shù)。首先，添加hillshade工具，雙擊指定一個elevation柵格數(shù)據(jù)（這里我選擇了marbles文件夾下的elevation），用下拉條指定azimuth和altitude值為azimuth和sun angle。運行，然后右鍵單擊輸出文件，選“Add to Display”在ArcMap里查看結(jié)果。178。 探索1：如果你想通過日期獲取太陽傾角，該怎么做呢？嘗試下面代碼，現(xiàn)在有五個參數(shù)，如下表進行正確設(shè)置：Display NameData typeTypeDirectionDefaultMonthLongRequiredinput1DateLongRequiredinput1LatitudeDoubleRequiredinput0Sun AngleDoubleDerivedoutput35AzimuthDoubleDerivedoutput300import ,sys,mathgp=()構(gòu)建兩個函數(shù)，首先判斷輸入月/日的合法性，然后獲取是一年中的第幾天def jdate(im, id): 通過年月日起返回一年中的第幾天 lastD = [31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31] jd = [0, 31, 59, 90, 120, 151, 181, 212, 243, 273, 304, 334] if ((im 0) and (im 13)): if ((id 0) and (id = lastD[im 1])): return jd[im 1] + id else: print Date not in month return 0 else: print Month should be between 1 and 12 inclusive return 0計算傾角的函數(shù)def declin(day): rad = / 180 xlong = + * day xanom = + * day g = xanom * rad elong = xlong + * (g) + * (2 * g) elong = elong * rad oblecl = * rad return ((oblecl) * (elong)) / rad主程序month = int([1])date = int([2])lat = float([3])decl = declin(jdate(month, date))sunangle = 90 lat + declazimuth = 180if sunangle 90: sunangle = 180sunangle azimuth = 0(3,str(sunangle))(4,str(azimuth))我們在思考遠一點，如何根據(jù)一天內(nèi)的任何時間而不是正午時間，獲取太陽角和方位角呢？感興趣的可以從本文開頭的鏈接中下砸相應(yīng)代碼學(xué)習(xí)。178。 探索2：如何在PythonWin里運行這個腳本？首先我們得明確幾個點：①我們將把hillshade作為腳本的一部分使用，并為其提供輸入?yún)?shù)：一個高程柵格（elevation raster）；②GetParameterastext僅在用作工具時起作用。下表為腳本工具的參數(shù)設(shè)置：Display nameData typeTypeDirectionDefaultMonthLongRequiredInputDateLongRequiredInputLatitudeDoubleRequiredInputWorkspaceWorkspaceRequiredInputElevationRaster DatasetRequiredInputHillshadeStringRequiredInput代碼如下：import , sys, mathgp = ()函數(shù)定義，同上，下面僅給出函數(shù)名稱：def jdate(im, id):def declin(day): 主程序，[]代替getparameterastext()month = int([1])date = int([2])lat = float([3]) = [4]輸入時注意，路徑應(yīng)為反斜杠“\”elev = [5]hillsh = [6]給輸出hillshade指定文件名decl = declin(jdate(month, date))sunangle = 90 lat + declazimuth = 180if sunangle 90: sunangle = 180sunangle azimuth = 0(about to try)try: = 1 (after overwriteoutputs setting, + + / + hillsh) (spatial) ( + / + hillsh) (elev, + / + str(hillsh), azimuth, sunangle) (done with hillshade) (spatial)except: print () (()) (spatial)閱讀代碼發(fā)現(xiàn)，第一個輸入?yún)?shù)不適用的getparameterastext(0)，[1]，這是因為getparameterastext()方法只在工具中起作用，[]有同樣的效果，索引從1而不是0開始，當(dāng)然，這需要首先引用sys模塊。這里我們直接指定輸出文件在輸入數(shù)據(jù)文件夾內(nèi)，省去了setparameterastext()方法，當(dāng)然這個方法在PythonWin中也無法運行。嘗試輸入?yún)?shù)如下圖，得到右下結(jié)果：5 在PythonWin里調(diào)試地理處理腳本既然我們已經(jīng)認真地學(xué)會了從Python中創(chuàng)建并運行地理統(tǒng)計工具，那么現(xiàn)在需要考慮如何調(diào)試我們的程序了。我們需要經(jīng)常在Python和添加的地理處理系統(tǒng)引用之間調(diào)試程序。當(dāng)一個地理處理工具運行失敗后，我們需要從Pythonwin中得到一個豐富的消息，而不是“未知錯誤”。 調(diào)試選擇和消息Python優(yōu)于AML的優(yōu)點之一是它有更好的調(diào)試方法，調(diào)試程序有很多選擇，這里不能一一列舉。216。 打印語句（Print statements）一開始就養(yǎng)成良好的調(diào)試方法是：將變量的當(dāng)前值或腳本的處理過程打印在屏幕上。比如，對之前的一段腳本加以修改：import arcgisscripting,sysgp=()=c:/prog/hmbarea=1(streams/arc, )=Analysis(, , 200)print Finished Cliping...(, , )可以看出成功運行腳本！然而當(dāng)在工具中運行時，print語句不會產(chǎn)生錯誤，但也不會輸出任何東西，因此，(Finished Cliping...)代替print語句。那么，如果想無論在工具中或Pythonwin中都可以顯示消息，就可以這兩句都寫上。我喜歡的做法是定義一個‘sendmsg’函數(shù)來輸出消息：def sendmsg(msg): print msg (msg)….sendmsg(Finished Cliping...)216。 獲取工具消息和（try:…except:）上面的方法很有用，但當(dāng)我們運行我們并不了解很多信息的地理處理工具時就顯得無能為力了。如果在Pythonwin中運行Buffer時出現(xiàn)錯誤，比如輸入文件不存在等，只能看到“未知錯誤”的提示，這就引出了GetMessage上的方便！現(xiàn)在我們要體驗兩種調(diào)試技巧：①GetMessages地理處理方法；②Python語言的錯誤處理程序：try…except。① 現(xiàn)在我們添加GetMessages查看錯誤信息：還是上面的代碼