【正文】 (isnull(streamsg),0,streamsg)) ma(histreams=( elev 1000) and stream2 )另外一個(gè)方法是直接將方法賦值給變量，這樣通過(guò)下面簡(jiǎn)單的一句話就可以實(shí)現(xiàn)前面所講的自定義函數(shù)了：ma=但也要確保在初始化ma變量之前檢查“Spatial”擴(kuò)展功能，安全的編程習(xí)慣是你在Except:(“spatial”)，這是為了避免在你嘗試CheckOutExtension時(shí)，其已經(jīng)被檢查過(guò)了，如下所示：(spatial)ma=ma(slopeg=slope(elev,))ma(stream2=con(isnull(streamsg),0,streamsg))ma(histreams=( elev 1000) and stream2 )(spatial)8數(shù)據(jù)管理和指針（Data Management and Cursors）（Data Management）處理表格數(shù)據(jù)、創(chuàng)建字段來(lái)儲(chǔ)存這些數(shù)據(jù)或者其他數(shù)據(jù)操作時(shí)GIS中非常重要的工作。的那個(gè)你需要做一連串?dāng)?shù)據(jù)管理和分析時(shí)，腳本是一個(gè)很好的選擇。在一些情況下，我們需要處理數(shù)據(jù)字段時(shí)有一系列工具可以使用，比如，“Add Fields”、“Calculate values for fields”、“delete fields”以及通過(guò)一個(gè)公共字段鏈接表格獲取附加數(shù)據(jù)字段。其中有些工具會(huì)利用不同的統(tǒng)計(jì)方法將輸入數(shù)據(jù)輸出為新的表格。其他情況下，我們需要處理行數(shù)據(jù)，比如柵格數(shù)據(jù)中單個(gè)要素或柵格數(shù)據(jù)的屬性值，接下來(lái)我們學(xué)習(xí)使用指針（Cursor）逐個(gè)處理行數(shù)據(jù)。首先，我們看一下一些能夠在處理表格數(shù)據(jù)（主要是包括屬性字段）時(shí)用到的工具。有些工具根據(jù)屬性選擇記錄，或復(fù)制或刪除選中記錄。ToolboxTool做什么？OutputAnalysisfrequency頻率統(tǒng)計(jì)Tablestatistics總結(jié)統(tǒng)計(jì)Tableselect利用where子句選擇要素新要素集table_select用SQL語(yǔ)句選擇并提取選中的屬性TableData ManagementAlias:managementaddField添加字段表格中字段calculateField通過(guò)表達(dá)式為字段賦值已存在字段的值deleteField刪除字段addXY為點(diǎn)要素添加X(jué)amp。Y值Xamp。Y值selectLayer By Attribute對(duì)圖層或表格通過(guò)屬性查詢選擇、更新或移除選擇copyFeatures復(fù)制選中的要素新要素集deleteFeatures刪除選中的要素addJoin根據(jù)一個(gè)共同的字段將一表格連接到圖層（或表格）連接關(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。 創(chuàng)建新腳本（），添加字段“elevm”到“”（surf_bld文件夾內(nèi)），并計(jì)算其值為“[elev]*”。代碼如下：import arcgisscriptinggp=()=c:/prog/surf_bld=management“management”。try: if not (,elevm):判斷此字段是否已經(jīng)存在 (,elevm,float) (,elevm,[elev]*)except:print (2)216。 創(chuàng)建新腳本（）使用management里的addxy工具為Marbes文件夾下的“”添加xamp。y值。代碼如下：import arcgisscriptinggp=()=1=c:/prog/marbles ()（Cursors）指針給了一種訪問(wèn)數(shù)據(jù)值的通道，允許遍歷數(shù)據(jù)表格中的所有記錄。由于記錄（行數(shù)據(jù)）可能是一個(gè)矢量要素或一柵格值，所以，指針在處理數(shù)據(jù)是十分強(qiáng)大。指針有三種類型：l SearchCursor：讀取一行中的值l InsertCursor：插入新的行l(wèi) UpdateCursor：改變行中值以及刪除行216。 嘗試下面的代碼，：import arcgisscriptinggp=()=c:/prog/surf_bldcur=()row=()while row: print ,row=()216。 嘗試一下代碼，介紹了讀取和顯示“”內(nèi)所有頂點(diǎn)，不要忘記while循環(huán)內(nèi)Next()語(yǔ)句，否則會(huì)出現(xiàn)無(wú)限循環(huán)！也展示了幾何要素的用法！import arcgisscriptinggp=()rows=(c:/prog/surf_bld/)row=()while row: feat= a=0 while a: stArray=(a) pnt=() while pnt: print str()+。+str()+。+str() pnt=() a=a+1row=()216。 嘗試下面代碼，展示了InsertCursor的用法，同樣的基本方法可用來(lái)讀取外部數(shù)據(jù)，比如文本文件。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：下面的代碼可以“打開(kāi)”、“寫入”、“關(guān)閉”一個(gè)文本文件，修改上面的腳本輸出節(jié)點(diǎn)值信息。(對(duì)“”的操作)import arcgisscripting,os,sysgp=()wspath=c:/prog/surf_bldtxtfile=fpath=wspath+/+txtfileif (fpath):(fpath)txtfile=open(fpath,w)=wspathrows=()row=()while row: feat= a=0 while a: geomArray=(a) fid=(ID) pnt=() while pnt: (str(fid) + 。 + str() + 。 + str()+\n) pnt=() a=a+1 row=() ()挑戰(zhàn)2：寫一個(gè)可以讀取剛才輸出的文本文件并且創(chuàng)建一個(gè)新的shp文件的腳本。如果你使用的是ArcGIS10，這里有最新的幫助，這里有相應(yīng)幫助。import arcgisscripting,sys,fileinputgp=()inFile=c:/prog/surf_bld/=c:/prog/surf_bldfcName==1(,fcName,point)try: cur=(fcName) for line in (inFile): values=(。) 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接下來(lái)該做什么呢？我們已經(jīng)做了很多，但是大多數(shù)都是創(chuàng)建處理特定事情的部分，比如遍歷一系列數(shù)據(jù)。但是到目前為止我們并沒(méi)有花太多時(shí)間在這些循環(huán)上。另一間我們沒(méi)有花太多時(shí)間處理的事情是，讓腳本工具更好地工作。輸入和輸出是這個(gè)處理過(guò)程中有挑戰(zhàn)性的部分，所以在這方面下手可能會(huì)更適合你。我已經(jīng)發(fā)現(xiàn)這是最令人困惑的部分，所以我將能發(fā)現(xiàn)的盡可能多的輸入/輸出情況做了總結(jié)，請(qǐng)看附錄1。最后一個(gè)建議是在模型中使用你的腳本工具，這是一個(gè)非常值得和有用的工作。附錄1：地理處理腳本中輸入amp。輸出方法指南 輸入數(shù)據(jù)集（Datasets）l 腳本：l 工具參數(shù)：（要素類、柵格數(shù)據(jù)或圖層等）必需的輸入 輸出數(shù)據(jù)集（Datasets）（如，Spatial Statistics/測(cè)量地理分布里的平均中心）l 腳本：l 工具參數(shù)：（要素類、柵格數(shù)據(jù)等）必需的輸入 輸出值和字段（如，Spatial Statistics/Utilities下的計(jì)算面積）l 腳本：輸入和輸出的數(shù)據(jù)集、復(fù)制輸入到輸出l 工具參數(shù)：輸入要素類、輸出要素類 輸出單個(gè)值（如，Spatial Statistics/Analyzing Patterns下的空間自相關(guān)）l 腳本：SetParameterAsTextl 工具參數(shù)：（Long、Double等）獲取的輸出 輸入和輸出表格或要素類（如，Management/Fields下的Transpose Time fields）l 腳本：GetParameterAsTextl 工具參數(shù)：輸入—表格視圖，必需的輸入；輸出—表格，必需的輸出 修改輸入到輸出，沒(méi)有找到簡(jiǎn)單的例子，除了Mosaic（Samples/Conversion/Raster）l 腳本：輸出用GetParameterAsText ；類似的SetParameterAsTextl 工具參數(shù)：輸出是必須的輸入，然后作為獲取的輸出，設(shè)置從原始處獲取 批處理（如，Management/Projections和Transformations/Feature下的Batch Project）l 腳本：輸入和輸出用GetParameterAsText，獲取的輸出不需要l 工具參數(shù)：輸入地理數(shù)據(jù)集，必需的輸入，允許多值；輸出工作空間，工作空間或要素?cái)?shù)據(jù)集，必需的輸入；輸出空間的單獨(dú)的獲取的輸出也是一樣 如果成功完成則返回True（如，Management/Raster下的Batch Calculate Statistics）l 腳本：nothingl 工具參數(shù)：布爾型參數(shù)，獲取的輸出默認(rèn)為True什么時(shí)候用Derived（獲取的，衍生的）？如果工具更新輸入到工具里，或者如果你需要一個(gè)輸出并用在一個(gè)模型里，但是你不想讓輸出參數(shù)在腳本對(duì)話框中出現(xiàn)。附錄2：其他 工作空間和臨時(shí)工作空間設(shè)置支持“臨時(shí)工作空間”環(huán)境設(shè)置的工具可將指定的位置用作輸出數(shù)據(jù)集的默認(rèn)工作空間?！芭R時(shí)工作空間”專門用于存放不愿保留的輸出數(shù)據(jù)?！芭R時(shí)工作空間”環(huán)境的主要用途是供“模型構(gòu)建器”使用?！澳Ｐ蜆?gòu)建器”需要使用一個(gè)工作空間來(lái)寫入中間數(shù)據(jù)集（模型運(yùn)行后便不再使用）。盡管它主要服務(wù)于“模型構(gòu)建器”，但有時(shí)可能也需要為各工具對(duì)話框設(shè)置臨時(shí)工作空間。使用工具對(duì)話框時(shí)，輸出數(shù)據(jù)集名稱會(huì)按照當(dāng)前工作空間環(huán)境和臨時(shí)工作空間環(huán)境設(shè)置自動(dòng)生成。 Python和PythonWin的版本不要急于升級(jí)你的Python或PythonWin版本。最好是使用隨ArcGIS安裝的版本，而且PythonWin不是自動(dòng)安裝的，在前文中有介紹。 68 / 6