【正文】 。 Python和PythonWin的版本不要急于升級你的Python或PythonWin版本。盡管它主要服務(wù)于“模型構(gòu)建器”，但有時可能也需要為各工具對話框設(shè)置臨時工作空間。“臨時工作空間”環(huán)境的主要用途是供“模型構(gòu)建器”使用。附錄2：其他 工作空間和臨時工作空間設(shè)置支持“臨時工作空間”環(huán)境設(shè)置的工具可將指定的位置用作輸出數(shù)據(jù)集的默認(rèn)工作空間。附錄1：地理處理腳本中輸入amp。我已經(jīng)發(fā)現(xiàn)這是最令人困惑的部分，所以我將能發(fā)現(xiàn)的盡可能多的輸入/輸出情況做了總結(jié)，請看附錄1。另一間我們沒有花太多時間處理的事情是，讓腳本工具更好地工作。),id=int(values[0]),x=float(values[1]),y=float(values[2]),feat=(),pnt=(Point),=x,=y,=pnt,=id,(feat),print,id,x,y,del,cur,()except:,print,(),if,cur:, del,cur接下來該做什么呢？我們已經(jīng)做了很多，但是大多數(shù)都是創(chuàng)建處理特定事情的部分，比如遍歷一系列數(shù)據(jù)。如果你使用的是ArcGIS10，這里有最新的幫助，這里有相應(yīng)幫助。,+,str(),+,。import,sys,math,stringgp,=,()import,arcgisscriptinggp=(),也可以這樣引用,=,c:/prog/Marblestry:,cur,=,(),feat,=,(),=,12,=,Sky,High,Lake,Camp,pnt,=,(Point),=,485339,=,4600001,=,pnt,(feat),del,curexcept:,print,()if,cur:,del,cur挑戰(zhàn)1：下面的代碼可以“打開”、“寫入”、“關(guān)閉”一個文本文件，修改上面的腳本輸出節(jié)點值信息。+str(),pnt=(),a=a+1row=()216。 嘗試一下代碼，介紹了讀取和顯示“”內(nèi)所有頂點，不要忘記while循環(huán)內(nèi)Next()語句，否則會出現(xiàn)無限循環(huán)！也展示了幾何要素的用法！import,arcgisscriptinggp=()rows=(c:/prog/surf_bld/)row=()while,row:,feat=,a=0,while,a:,stArray=(a),pnt=(),while,pnt:,print,str()+。指針有三種類型：l SearchCursor：讀取一行中的值l InsertCursor：插入新的行l(wèi) UpdateCursor：改變行中值以及刪除行216。代碼如下：import,arcgisscriptinggp=()=1=c:/prog/marbles ()（Cursors）指針給了一種訪問數(shù)據(jù)值的通道，允許遍歷數(shù)據(jù)表格中的所有記錄。 創(chuàng)建新腳本（）使用management里的addxy工具為Marbes文件夾下的“”添加xamp。代碼如下：import,arcgisscriptinggp=()=c:/prog/surf_bld=management“management”。Y值selectLayer,By,Attribute對圖層或表格通過屬性查詢選擇、更新或移除選擇copyFeatures復(fù)制選中的要素新要素集deleteFeatures刪除選中的要素addJoin根據(jù)一個共同的字段將一表格連接到圖層（或表格）連接關(guān)系removeJoin移除已存在連接CoverageAlias:Arcselect(reselect)根據(jù)邏輯表達(dá)式提取Coverage地圖要素輸出Coverage（包含屬性表）additem為信息表添加屬性字段dropitem從信息表中移除屬性字段joinitem根據(jù)相關(guān)屬性融合數(shù)據(jù)表Tableaddxy為point、label、node表格添加x、y字段字段添加到已存在表格216。ToolboxTool做什么？OutputAnalysisfrequency頻率統(tǒng)計Tablestatistics總結(jié)統(tǒng)計Tableselect利用where子句選擇要素新要素集table_select用SQL語句選擇并提取選中的屬性TableData,ManagementAlias:managementaddField添加字段表格中字段calculateField通過表達(dá)式為字段賦值已存在字段的值deleteField刪除字段addXY為點要素添加Xamp。首先，我們看一下一些能夠在處理表格數(shù)據(jù)（主要是包括屬性字段）時用到的工具。其中有些工具會利用不同的統(tǒng)計方法將輸入數(shù)據(jù)輸出為新的表格。的那個你需要做一連串?dāng)?shù)據(jù)管理和分析時，腳本是一個很好的選擇。下面的例子說明了這種用法，雖然只有一點語句，但包括了在處理多工具中的兩個棘手的問題。Spatial,Analyst工具下的MultiOutputMapAlgebra工具允許使用完整的地圖代數(shù)賦值語句，包括嵌套結(jié)構(gòu)（?。?。代碼如下：import,arcgisscriptinggp=()=c:/prog/pendatacov=[geology,landuse,publands,flood,streams,roads,vegetation]fld=[TYPEID,LUCODE,PUBCODE,FLOODCODE,STCODE,ROADCODE,VEGCODE]=1if,not,(multiFeatRas):新建一個文件夾保存處理后數(shù)據(jù),(,multiFeatRas)for,i,in,range(len(cov)):遍歷coverage數(shù)據(jù),if,(cov[i]+/polygon):,covTop=cov[i]+/polygon,elif,(cov[i]+/arc):,covTop=cov[i]+/arc,(covTop,fld[i],multiFeatRas/+cov[i]+g,60)7在Python腳本中使用地圖代數(shù)（Map,Algebra）將AML程序移植到Python中唯一令我猶豫的是它不能像在AML中那樣編寫地圖代數(shù)語句來處理格網(wǎng)數(shù)據(jù)。CITYCODE=2639。你也可以裁剪所有數(shù)據(jù)時設(shè)置通配符為“*”，或者不適用通配符。創(chuàng)建一個小的工作空間命名為“Sanmateo”，‘g’開頭（“g*”）的coverage文件的polygon，存入此空間。代碼如下：import,arcgisscriptinggp=()(spatial)=c:/prog/hmbareaif,not,(+/): ,(c:/prog/pendata/cities/polygon,),，所以我使用pendata文件夾下的cities里的polygon轉(zhuǎn)換的shp文件取代。在這個處理過程中，我們也會在創(chuàng)建輸出數(shù)據(jù)之前用“存在工具”中一個很方面的方法檢查其是否已經(jīng)存在。這會讓你感到腳本在處理多種數(shù)據(jù)方面的強大作用。 判斷一個字段屬性是否存在。，比如在使用數(shù)據(jù)之前添加一個判斷是否存在的語句：if,(“”) (“”,””,200)在下面的代碼中我們將使用這種方法判斷數(shù)據(jù)是否存在，其可以成為overwriteoutput設(shè)置的替代語句。在柵格數(shù)據(jù)集屬性中尋找這個方法：import,arcgisscripting,sysgp,=,(),=,c:/prog/hmbareadta,=,dsc,=,(dta)print,dta,+,+,+,:,+,str(),+,bands.運行結(jié)果：“,RasterDataset:,3,bands.”（Exists）關(guān)于一個數(shù)據(jù)集的一個更基本的信息是其是否存在。try:,(c:\prog\pendata\lugrid,envel,luwood)except:print,()216。,39。39。39。 利用描述系統(tǒng)用一個數(shù)據(jù)邊界裁剪另一個柵格數(shù)據(jù)，將裁剪后的土地利用圖存在woodside文件夾下。FeatureClass屬性包括所有表個屬性和數(shù)據(jù)集屬性；柵格數(shù)據(jù)屬性包括柵格波段屬性和數(shù)據(jù)集屬性，所以數(shù)據(jù)集屬性在一定程度上是共有的（Feature,Class、Rasters、Coverages）。)except:sendmsg(())6使用描述（Describe）和存在（Exists）獲取數(shù)據(jù)信息有很多情況下我們需要使用某個GIS數(shù)據(jù)的特征去處理其他數(shù)據(jù)。,39。,39。,39。try:,import,arcgisscripting,def,sendmsg(msg):,print,msg,(msg),=1,=c:\prog\pendata,(39。)⑦ 在第一行前添加try:，最后一行后添加except:，最后加上，sendmsg(()).216。,39。① 添加引用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。216。按下（Step,over）按鈕進(jìn)入下一行，繼續(xù)執(zhí)行；點擊（Watch）按鈕，出現(xiàn)，在Expression下New,Item里輸入lat等可以查看當(dāng)前變量值；當(dāng)你遍歷完程序，點擊關(guān)閉調(diào)試；嘗試添加一個斷點，當(dāng)定位到某一行時，在點，然后運行此腳本在“Run,in,the,debugger”下，發(fā)現(xiàn)沒有任何影響，但是修改第一行為“l(fā)at=10”，運行發(fā)現(xiàn)在斷點處停止，點擊按鈕繼續(xù)程序。選擇“Stepthrough,in,the,debugger”。這里我們演示單步執(zhí)行代碼的過程。① 現(xiàn)在我們添加GetMessages查看錯誤信息：還是上面的代碼，把“”修改為“”，查看錯誤信息： =Analysis (,200) (,)運行之，查看錯誤信息，然后修改代碼如下：try:,=Analysis,(,200),(,)GetMessages():所以信息GetMessages(0):只顯示消息GetMessages(1):只顯示警告GetMessages(2):只顯示錯誤except:, ,print,()運行可以看到以下錯誤提示：順便說一下GetMessages的集中形式和代表的含義② 運行下面簡單的代碼，并改變x的值為非零，查看結(jié)果：x=0.y=15.try: ,print,y/x,except:,print,錯誤！0不能是被除數(shù)！PythonWin提供了一些工具：①單步執(zhí)行代碼②插入斷點③觀測變量，或其他事情。 獲取工具消息和（try:…except:）上面的方法很有用，但當(dāng)我們運行我們并不了解很多信息的地理處理工具時就顯得無能為力了。那么，如果想無論在工具中或Pythonwin中都可以顯示消息，就可以這兩句都寫上。 打印語句（Print,statements）一開始就養(yǎng)成良好的調(diào)試方法是：將變