【正文】 “臨時(shí)工作空間”專門(mén)用于存放不愿保留的輸出數(shù)據(jù)。 嘗試下面代碼，展示了InsertCursor的用法，同樣的基本方法可用來(lái)讀取外部數(shù)據(jù)，比如文本文件。有些工具根據(jù)屬性選擇記錄，或復(fù)制或刪除選中記錄。建議：(“g*”)，因?yàn)槟阈枰獙ふ襝overage，在你處理列表里成員時(shí)，可以將成員名稱賦給變量f，(f+”/polygon”)查找多邊形要素集。+str()print envel ，Extent是一個(gè)Object，而clip語(yǔ)句的第二個(gè)參數(shù)Rectangle是一個(gè)Envelope類型數(shù)據(jù)，所以這里我使用了envel接收Extent的四個(gè)值。lugrid39。然后，一系列調(diào)試工具出現(xiàn)；一個(gè)黃色三角出現(xiàn)在第一行代碼前。嘗試輸入?yún)?shù)如下圖，得到右下結(jié)果：5 在PythonWin里調(diào)試地理處理腳本既然我們已經(jīng)認(rèn)真地學(xué)會(huì)了從Python中創(chuàng)建并運(yùn)行地理統(tǒng)計(jì)工具，那么現(xiàn)在需要考慮如何調(diào)試我們的程序了。216。然后，我們需要將腳本加進(jìn)工具（Making a script into a tool），那樣才能在ArcToolbox或ModelBuilder或Command line中使用。② Geoprocessor Programming Model（PDF），包含方法（左箭頭表示）、屬性（可讀寫(xiě)的表示為杠鈴形，只讀的表示為部分杠鈴形）：如何使用Geoprocessor Programming Model中的ListsLists（列表）及其屬性和方法在圖表中用紫色標(biāo)出，如下：現(xiàn)在我們?cè)囍帉?xiě)一段腳本列舉出屬性表中所有屬性值（Fields）（以hmbarea柵格土地利用為例，文件存在c:\prog\hmbarea下）import arcgisscripting, sysgp = () = c:/prog/hmbarea fieldList = (landuse, *, all)dsc=(landuse)print landuse+ ++:for field in fieldList:... print , (此時(shí)輸出結(jié)果如右圖)3使用地理處理工具——Toolboxes和Aliases總所周知，地理處理工具在腳本中的使用和ArcToolbox中相同，但是需要提供工具名稱來(lái)使用它們。當(dāng)循環(huán)一個(gè)數(shù)據(jù)集時(shí)（GIS中很常用的工作）while循環(huán)很有用?！痠f（continued）接下來(lái)，我們會(huì)探索一下另一個(gè)方便的模塊：: 開(kāi)始之前，在d:/下創(chuàng)建一個(gè)“testfolder”文件夾，然后新建一個(gè)“”文件； 嘗試以下代碼段，確保print語(yǔ)句前有縮進(jìn)。當(dāng)要引用多個(gè)模塊是，中間用逗號(hào)分隔，比如：import arcgisscripting,sys,string,os,math我們也可以自己為頻繁使用的方法創(chuàng)建Module，下面，我們開(kāi)始體驗(yàn)內(nèi)置的Modules。d:\\work\\39。.capitalize()Zhulj()即將capitalize()方法用于ss=39。我想自動(dòng)化和邏輯是關(guān)鍵，它們區(qū)別于我們多數(shù)使用計(jì)算機(jī)時(shí)的交互活動(dòng)。Python腳本使用詳解目錄寫(xiě)在前面的話 2前言 2一、PYTHON語(yǔ)言基礎(chǔ) 31數(shù)學(xué)運(yùn)算符 32字符串操作 43模塊的使用（Modules） 54使用def構(gòu)建函數(shù) 65流程控制結(jié)構(gòu)：If，While，F(xiàn)or 66簡(jiǎn)單輸入和輸出 9二、ARCGISamp。我們發(fā)Email，寫(xiě)文章或者設(shè)計(jì)地圖，都需要和計(jì)算機(jī)交互，而處理一系列數(shù)據(jù)，我們需要自動(dòng)化和利用邏輯來(lái)指導(dǎo)自動(dòng)化。zhulj39。print p2d:\work\反斜線“\”和某些字母一起有特殊用法，如\n為換行，“\”為轉(zhuǎn)義字符，如“\\”則表示“\”print 39。math和random模塊很多常用的數(shù)學(xué)計(jì)算功能都可以通過(guò)math找到，比如三角計(jì)算或?qū)?shù)計(jì)算，如果要使用復(fù)雜數(shù)字，就使用cmath模塊。import if (d:/testfolder/): print 測(cè)試文件夾存在 print txt文件存在elif (d:/testfolder): print 測(cè)試文件夾存在 print 測(cè)試文件夾存在，但txt文件不存在else: print 兩者都不存在可選探索示例接下來(lái)的例子做的事情對(duì)GIS非常重要，但是實(shí)際上不用任何地理處理代碼。在后面地理處理中我們會(huì)接觸一些例子。但是記住一個(gè)名稱可能有好幾個(gè)工具，比如，裁剪工具（Clip）在CoverageAnalysisExtract工具集里，另一個(gè)是在Data Management Tools下的Raster工具集下。如前面提到的那樣，這個(gè)腳本只能用于工具，包括輸入/輸出方法是PythonWin不能處理的，但這些是多數(shù)工具必需的。 在你的toolbox中新建一個(gè)model，將剛才創(chuàng)建的腳本工具（script tool）拖進(jìn)來(lái)。我們需要經(jīng)常在Python和添加的地理處理系統(tǒng)引用之間調(diào)試程序。按下（Step over）按鈕進(jìn)入下一行，繼續(xù)執(zhí)行；點(diǎn)擊（Watch）按鈕，出現(xiàn)，在Expression下New Item里輸入lat等可以查看當(dāng)前變量值；當(dāng)你遍歷完程序，點(diǎn)擊關(guān)閉調(diào)試；嘗試添加一個(gè)斷點(diǎn)，當(dāng)定位到某一行時(shí)，在點(diǎn)，然后運(yùn)行此腳本在“Run in the debugger”下，發(fā)現(xiàn)沒(méi)有任何影響，但是修改第一行為“l(fā)at=10”，運(yùn)行發(fā)現(xiàn)在斷點(diǎn)處停止，點(diǎn)擊按鈕繼續(xù)程序。,39。try: (c:\prog\pendata\lugrid,envel,luwood)except:print ()216。你也可以裁剪所有數(shù)據(jù)時(shí)設(shè)置通配符為“*”，或者不適用通配符。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。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)”、“寫(xiě)入”、“關(guān)閉”一個(gè)文本文件，修改上面的腳本輸出節(jié)點(diǎn)值信息。“臨時(shí)工作空間”環(huán)境的主要用途是供“模型構(gòu)建器”使用。附錄2：其他 工作空間和臨時(shí)工作空間設(shè)置支持“臨時(shí)工作空間”環(huán)境設(shè)置的工具可將指定的位置用作輸出數(shù)據(jù)集的默認(rèn)工作空間。+str() pnt=() a=a+1row=()216。首先，我們看一下一些能夠在處理表格數(shù)據(jù)（主要是包括屬性字段）時(shí)用到的工具。創(chuàng)建一個(gè)小的工作空間命名為“Sanmateo”，‘g’開(kāi)頭（“g*”）的coverage文件的polygon，存入此空間。 39。,39。選擇“Stepthrough in the debugger”。這里我們直接指定輸出文件在輸入數(shù)據(jù)文件夾內(nèi)，省去了setparameterastext()方法，當(dāng)然這個(gè)方法在PythonWin中也無(wú)法運(yùn)行。 但是它能在作為Modelbuilder工具正常運(yùn)行，下面我們將使用它的輸出參數(shù)創(chuàng)建一個(gè)hillshade。前兩個(gè)參數(shù)索引為0和1，所以接下來(lái)輸出參數(shù)的索引為2和3。比如，你想得到特定類型的文件，就找到listfeatureclasses方法：fcList=(“w*”,”polygon”)注：此方法的語(yǔ)法為：({wildCard} As String, {FeatureType} As String) As Python List{} = optional wildcard為通配符，和星號(hào)（*）組合使用，如果沒(méi)有通配符則返回工作目錄下的所有feature classes。由于while提供一種容易結(jié)束循環(huán)的方法，我們甚至用它代替if語(yǔ)句。②如果你只需做一件事情，你可以在冒號(hào)后面同一行添加簡(jiǎn)短的語(yǔ)句，比如： if x0: print ‘x 比0大’ print ‘下一行不要縮進(jìn)了。通常我們會(huì)將一行import Module名放在程序頂部，比如：import arcgisscripting當(dāng)然，這不必成為你程序的第一行，但必須在使用它里面方法之前。]你可以用split()方法解析出不同的字符串片段，并創(chuàng)建一個(gè)列表（List），我們可以使用其中不同的元素print plist[0]d:print plist[1]workp2=39。zhulj39。腳本編程的主要目的是使枯燥的處理數(shù)據(jù)工作自動(dòng)化，通過(guò)邏輯來(lái)指揮處理過(guò)程。PYTHON 101如何創(chuàng)建地理處理對(duì)象（geoprocessor object） 102獲取地理處理幫助 10：如何使用Geoprocessor Programming Model中的Lists 113使用地理處理工具——Toolboxes和Aliases 114在建模中使用腳本（Scripts in ModelBuilder） 125 在PythonWin里調(diào)試地理處理腳本 18 調(diào)試選擇和消息 19 20 216使用描述（Describe）和存在（Exists）獲取數(shù)據(jù)信息 21 22（Exists） 23 237在Python腳本中使用地圖代數(shù)（Map Algebra） 268數(shù)據(jù)管理和指針（Data Management and Cursors） 27（Data Management） 27（Cursors） 28附錄1：地理處理腳本中輸入amp。在地理處理腳本邏輯中，我們需要在允許我們做的事情中作出決定，比如，處理柵格數(shù)據(jù)不同于矢量數(shù)據(jù)，或?yàn)闆](méi)投影的數(shù)據(jù)設(shè)置投影，或處理僅在特定時(shí)間搜集的數(shù)據(jù)集。print ()print s[0]zStrings可以像一個(gè)字母列表一樣處理，第一個(gè)字母下標(biāo)為0，某個(gè)字符段可以用1:3來(lái)格式化：從第1個(gè)的開(kāi)頭到第3個(gè)的開(kāi)頭，不包括下標(biāo)為3的字母；s[1]表示倒過(guò)來(lái)第一個(gè)，相當(dāng)于s[len(s)1]s1=s[1]print s1hprint s[2:]ljprint s[2:3]uprint s[2:4]ulprint s[2:],s[:5]ulj zhuljs2=()print s2ZHULJ我們可以將字符串方法的結(jié)果賦給新的變量s3=s+s2print s3zhuljZHULJ字符串組合用“+”print s*3zhuljzhuljzhulj字符串重復(fù)用“*”，后為重復(fù)次數(shù)selstr=39。Jerry\39。和之前一樣，通過(guò)以下表格來(lái)體現(xiàn)模塊的使用：輸入結(jié)果Notesimport mathprint (100)以10為底的對(duì)數(shù)print (100)自然對(duì)數(shù)print π是一個(gè)靜態(tài)常量，所以不需要括弧pi=print pi如果不想總是輸入“”可以將其賦給一變量pi不需要print即可查看變量值print (radians)print (radians)print (radians)三角函數(shù)的計(jì)算是弧度制degrad=pi/18045*degrad度轉(zhuǎn)化為弧度sin=sin(45*degrad)sin(90*degrad)即使功能函數(shù)（像sin）都可以賦給一個(gè)變量