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