【正文】 “臨時(shí)工作空間”專(zhuān)門(mén)用于存放不愿保留的輸出數(shù)據(jù)。 嘗試下面代碼，展示了InsertCursor的用法，同樣的基本方法可用來(lái)讀取外部數(shù)據(jù)，比如文本文件。有些工具根據(jù)屬性選擇記錄，或復(fù)制或刪除選中記錄。建議：(“g*”)，因?yàn)槟阈枰獙ふ襝overage，在你處理列表里成員時(shí)，可以將成員名稱(chēng)賦給變量f，(f+”/polygon”)查找多邊形要素集。+str()print,envel,，Extent是一個(gè)Object，而clip語(yǔ)句的第二個(gè)參數(shù)Rectangle是一個(gè)Envelope類(lèi)型數(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中相同，但是需要提供工具名稱(chēng)來(lái)使用它們。當(dāng)循環(huán)一個(gè)數(shù)據(jù)集時(shí)（GIS中很常用的工作）while循環(huán)很有用。’if（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ě)在前面的話(huà) 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\反斜線(xiàn)“\”和某些字母一起有特殊用法，如\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è)名稱(chēng)可能有好幾個(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)值信息?！芭R時(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。比如，你想得到特定類(lèi)型的文件，就找到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*deg