【正文】 使用工具對話框時，輸出數(shù)據(jù)集名稱會按照當前工作空間環(huán)境和臨時工作空間環(huán)境設(shè)置自動生成。最后一個建議是在模型中使用你的腳本工具，這是一個非常值得和有用的工作。,+,str()+\n),pnt=(),a=a+1,row=() ()挑戰(zhàn)2：寫一個可以讀取剛才輸出的文本文件并且創(chuàng)建一個新的shp文件的腳本。 嘗試下面的代碼，：import,arcgisscriptinggp=()=c:/prog/surf_bldcur=()row=()while,row:,print,row=()216。 創(chuàng)建新腳本（），添加字段“elevm”到“”（surf_bld文件夾內(nèi)），并計算其值為“[elev]*”。在一些情況下，我們需要處理數(shù)據(jù)字段時有一系列工具可以使用，比如，“Add,Fields”、“Calculate,values,for,fields”、“delete,fields”以及通過一個公共字段鏈接表格獲取附加數(shù)據(jù)字段。),wsName=San_Mateo,newWS=workspaceLoc+/+wsName,if,not,(newWS):,(workspaceLoc,wsName),fList=(g*),for,f,in,fList:,if,(f+/polygon):,(f+/polygon,newWS+/+f)except:,print,()探索2：創(chuàng)建兩個Python列表：cov=[“geology”,”landuse”,”publands”,”flood”,”streams”,”roads”,”vegetation”]fld=[“TYPEID”,”LUCODE”,”PUBCODE”,”FLOODCODE”,”STCODE”,”ROADCODE”,”VEGCODE”]將每個coverage根據(jù)對應(yīng)的字段轉(zhuǎn)換為柵格，分辨率設(shè)置為60。在本例中，我們使用了ArcInfo,Workspace，并把格網(wǎng)柵格存在這里，當然你也可以將他們存在Geodatabase中。測試一個數(shù)據(jù)是否存在也可幫助我們避免錯誤。+str()+39。比如，在柵格運算中，我們可能想用一個數(shù)據(jù)的邊界去界定另個數(shù)據(jù)集，很像裁剪操作，或者檢查一個數(shù)據(jù)的拓撲錯誤。landuse/polygon39。 轉(zhuǎn)換腳本（conversion,script）既然我們已經(jīng)會在幫助系統(tǒng)里尋找答案，也掌握了幾種程序調(diào)試方法，現(xiàn)在讓我們建立一個腳本。這些都可以在PythonWin的幫助系統(tǒng)里找到相關(guān)示例教程。216。178。打開“Sun,Angle”和“Azimuth”發(fā)現(xiàn)它們還是默認值，說明此腳本工具還沒有運行。那么，它是干什么的呢？還記得結(jié)果是輸出兩個數(shù)字參數(shù)，那么，這些數(shù)字哪去了呢？很好的問題，這僅能說明你能創(chuàng)建一個工具，但是不能想ArcToolbox那樣運行。當然也可以使用之前創(chuàng)建的toolbox。每次你運行這個工具時，都會看到一個對話框，提示輸入?yún)?shù)，這個方法允許你在接下來的程序中使用。每一個工具都有自己的別名，我們可以通過右鍵屬性來查看：AliaseToolbox“conversion”Conversion“3d”3D,Analyst“geocoding”Geocoding“analysis”Analysis“ga”Geostatistical,Analyst“arc”Coverage“l(fā)r”Linear,Referencing“management”Data,Management“sa”Spatial,Analyst“cartography”Cartography“stats”Spatial,Statistics現(xiàn)在我們用一段腳本來解釋：import,arcgisscripting,sysgp=()=”c:/prog/hmbarea”=1,OverWriteOutput:Boolean,為1表示允許覆蓋已存在文件,將streams/arc轉(zhuǎn)換為shp文件(streams/arc,),利用轉(zhuǎn)換后的shp文件，做200米的緩沖(,200),用做過緩沖的shp裁剪(,)注：上面腳本用“”注釋的中文內(nèi)容不要出現(xiàn)在腳本文件中，否則會出現(xiàn)錯誤結(jié)果截圖：如果你一次使用一個工具集中的若干工具，可以通過環(huán)境設(shè)置省下一些文字：,=,Analysis,(,200),(,)4在建模中使用腳本（Scripts,in,ModelBuilder）首先，需要記住的很重要的一點是，ArcToolbox里相當數(shù)量的工具實際上都是腳本。比如，特征數(shù)據(jù)可能是①單個shp文件；②geodatabase（地理數(shù)據(jù)庫，我們可能指定地理數(shù)據(jù)庫為工作空間）；③多邊形、弧或點要素的coverage數(shù)據(jù)。 下面的代碼創(chuàng)建并輸出指定文件夾內(nèi)shp文件名列表（每個都以‘.shp’結(jié)尾）import,osws=c:/prog/hmbareailist=(ws)創(chuàng)建一個列表保存工作文件夾內(nèi)的文件fcs=[]創(chuàng)建一個空列表，保存結(jié)尾為‘.shp’的文件for,i,in,ilist:,if,(.shp):,(i)for,fc,in,fcs:print,fc（輸出結(jié)果如右圖所示）216。下面例子簡單體會一下布爾型表達式：x=1while,x10: print,x x=x+1表達式“x10”結(jié)果是true或false，所以這樣允許我們在計算完一種情況時運行一系列代碼。將會出現(xiàn)錯誤結(jié)果。 有三種流程控制操作：if,僅在一個特定情形下才執(zhí)行語句；while,當一種情形存在下，持續(xù)執(zhí)行語句for,遍歷一系列值這些語法和def有些相似：初始語句后加頓號、需要執(zhí)行的語句塊有縮進。y2=4152963不同賦值語句間用“；”分隔xr=x2x1yr=y2y1(xr,yr)(xr**2+yr**2)(xr*xr+yr*yr)**不同的方式，相同的結(jié)果import,random()random()方法，每次結(jié)果都不同，值域為[,)rnd=rnd()mu=50s=10print,(mu,s)4使用def構(gòu)建函數(shù)有點像Module，但更簡單，函數(shù)是一個自己定義功能，用在之后的代碼中，并且提供任何你想要使用的參數(shù)。Python安裝時自帶了大量Modules，最常用的有math,sys,random,。Jerry39。work39。print,selstrelev1000字符串可以使用單引號或雙引號，跨行時用雙引號。下面的表格顯示了整型（Integer）和浮點型（Float）各種組合運算的結(jié)果，記住一條規(guī)則，只要參與運算的有浮點型，則結(jié)果為浮點型；全為整型時，結(jié)果才為整型。在接下來的聯(lián)系中，我們會探索Python的使用以及創(chuàng)建腳本來使用ArcGIS里眾多的地理處理工具。值得一提的是ArcGIS的在線幫助文檔，一個實時更新的GIS寶庫，很多專業(yè)性知識都可以找到答案，點擊鏈接ArcGIS10中文幫助、。茫茫互聯(lián)網(wǎng)，終于找到了舊金山州立大學(xué)Jerry,Davis教授的個人主頁，對其中《Geoprocessing,Scripts,With,Python》如獲至寶，獨樂樂不如眾樂樂，現(xiàn)在將其教程翻譯并結(jié)合自己的學(xué)習(xí)情況給出總結(jié)。希望能夠給更多想學(xué)習(xí)Python的同學(xué)一個參考。 獲取更過腳本例子來學(xué)習(xí),：ESRI的地理處理模型和腳本工具庫是個不錯的選擇。所有你能在ArcToolbox或Model中使用的工具都能夠用在Python腳本中，這些腳本可以生成腳本工具，像其他地理處理工具一樣使用。輸入表達式結(jié)果Notes2+35整型結(jié)果2.+3，結(jié)果浮點型2312*36整型結(jié)果2.*3浮點型5/22整型5./25%21取模Az=270Newaz=az+180Print,newaz%36090取模的用途之一——方位角加180后逆轉(zhuǎn)方向5**22525**沒有sqrt()功能，除非添加math模塊2字符串操作注：使用Python幫助：有超過30種內(nèi)置方法來處理字符，請到Sequence,Types下的String,Methods尋找?guī)椭∽址且淮帜?，比如’San,Francisco’，字符串下標從0開始。othersel=”’elev’1000”print,othersel‘elev’1000print,()False一些方法返回值為布爾型（True或False），一些返回索引值（下標值）print,()Truep=39。,39。s,Kidsprint,39。當然還有好多Modules可以下載，比如數(shù)字處理（Numeric）——numpy。這個函數(shù)從此可像變量那樣在程序中使用，結(jié)合例子更容易理解。這三個結(jié)構(gòu)的一些重要的公共特征：①if、while、for語句均以冒號結(jié)尾，接下來是縮進的代碼塊，用于if、while、for定義的情形。但是不幸的是，你可能不知道DEM文本文件的垂直單位是英尺還是米。許多情況下我們需要使用條件代碼。 下面這個例子的循環(huán)較多次數(shù)。當我們想遍歷工作空間里的coverage文件時，應(yīng)使用ListDatasets而不是ListFeatureClasses。腳本都有一個圖標。索引0和1指第一個和第二個參數(shù)。216。比如輸出一種數(shù)據(jù)，柵格或特征數(shù)據(jù)（.shp）之類的。右鍵單擊工具，選擇Run，然后發(fā)現(xiàn)兩個輸出參數(shù)已經(jīng)改變！需要注意的是：latitude范圍是90~90，~。 探索2：如何在PythonWin里運行這個腳本？首先我們得明確幾個點：①我們將把hillshade作為腳本的一部分使用，并為其提供輸入?yún)?shù)：一個高程柵格（elevation,raster）；②GetParameterastext僅在用作工具時起作用。 打印語句（Print,statements）一開始就養(yǎng)成良好的調(diào)試方法是：將變量的當前值或腳本的處理過程打印在屏幕上。這里我們演示單步執(zhí)行代碼的過程。① 添加引用import,arcgisscripting② 創(chuàng)建“sendmsg”函數(shù) def,sendmsg(msg): print,msg (msg)③ =1④ 使用Data,Management,toolbox中的Workspace工具集中的CreateFolder工具在“c:/prog”中創(chuàng)建新文件夾“woodside”(c:\prog,woodside)⑤ 設(shè)置gp的工作文件夾為woodside=c:\prog\woodside⑥ 使用Conversion里的To,Raster工具集中的DEMtoRaster工具將c:/prog/，不要提供任何可選參數(shù)值(39。,39。39。，比如在使用數(shù)據(jù)之前添加一個判斷是否存在的語句：if,(“”) (“”,””,200)在下面的代碼中我們將使用這種方法判斷數(shù)據(jù)是否存在，其可以成為overwriteoutput設(shè)置的替代語句。代碼如下：import,arcgisscriptinggp=()(spatial)=c:/prog/hmbareaif,not,(+/): ,(c:/prog/pendata/cities/polygon,),，所以我使用pendata文件夾下的cities里的polygon轉(zhuǎn)換的shp文件取代。代碼如下：import,arcgisscriptinggp=()=c:/prog/pendatacov=[geology,landuse,publands,flood,streams,roads,vegetation]fld=[TYPEID,LUCODE,PUBCODE,FLOODCODE,STCODE,ROADCODE,VEGCODE]