【正文】 rint,id,x,y,del,cur,()except:,print,(),if,cur:, del,cur接下來該做什么呢？我們已經(jīng)做了很多，但是大多數(shù)都是創(chuàng)建處理特定事情的部分，比如遍歷一系列數(shù)據(jù)。附錄2：其他 工作空間和臨時(shí)工作空間設(shè)置支持“臨時(shí)工作空間”環(huán)境設(shè)置的工具可將指定的位置用作輸出數(shù)據(jù)集的默認(rèn)工作空間?！芭R時(shí)工作空間”環(huán)境的主要用途是供“模型構(gòu)建器”使用。另一間我們沒有花太多時(shí)間處理的事情是，讓腳本工具更好地工作。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)閉”一個(gè)文本文件，修改上面的腳本輸出節(jié)點(diǎn)值信息。代碼如下：import,arcgisscriptinggp=()=1=c:/prog/marbles ()（Cursors）指針給了一種訪問數(shù)據(jù)值的通道，允許遍歷數(shù)據(jù)表格中的所有記錄。ToolboxTool做什么？OutputAnalysisfrequency頻率統(tǒng)計(jì)Tablestatistics總結(jié)統(tǒng)計(jì)Tableselect利用where子句選擇要素新要素集table_select用SQL語句選擇并提取選中的屬性TableData,ManagementAlias:managementaddField添加字段表格中字段calculateField通過表達(dá)式為字段賦值已存在字段的值deleteField刪除字段addXY為點(diǎn)要素添加Xamp。下面的例子說明了這種用法，雖然只有一點(diǎn)語句，但包括了在處理多工具中的兩個(gè)棘手的問題。你也可以裁剪所有數(shù)據(jù)時(shí)設(shè)置通配符為“*”，或者不適用通配符。這會(huì)讓你感到腳本在處理多種數(shù)據(jù)方面的強(qiáng)大作用。try:,(c:\prog\pendata\lugrid,envel,luwood)except:print,()216。 利用描述系統(tǒng)用一個(gè)數(shù)據(jù)邊界裁剪另一個(gè)柵格數(shù)據(jù)，將裁剪后的土地利用圖存在woodside文件夾下。,39。)⑦ 在第一行前添加try:，最后一行后添加except:，最后加上，sendmsg(()).216。按下（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)沒有任何影響，但是修改第一行為“l(fā)at=10”，運(yùn)行發(fā)現(xiàn)在斷點(diǎn)處停止，點(diǎn)擊按鈕繼續(xù)程序。 獲取工具消息和（try:…except:）上面的方法很有用，但當(dāng)我們運(yùn)行我們并不了解很多信息的地理處理工具時(shí)就顯得無能為力了。我們需要經(jīng)常在Python和添加的地理處理系統(tǒng)引用之間調(diào)試程序。運(yùn)行，然后右鍵單擊輸出文件，選“Add,to,Display”在ArcMap里查看結(jié)果。 在你的toolbox中新建一個(gè)model，將剛才創(chuàng)建的腳本工具（script,tool）拖進(jìn)來。“下一步”后是參數(shù)配置頁面，如下圖設(shè)置各參數(shù)如表格所示：Display,NameData,typeDirectionDefaultLatitudeDoubleinput0DeclinationDoubleinput0Sun,AngleDoubleoutput35AzimuthDoubleoutput300216。如前面提到的那樣，這個(gè)腳本只能用于工具，包括輸入/輸出方法是PythonWin不能處理的，但這些是多數(shù)工具必需的。仍然以太陽角計(jì)算代碼為例，我們給其加上兩句引用，四句輸入輸出語句，就可以用作Modelbuilder中的一個(gè)步驟了。但是記住一個(gè)名稱可能有好幾個(gè)工具，比如，裁剪工具（Clip）在CoverageAnalysisExtract工具集里，另一個(gè)是在Data,Management,Tools下的Raster工具集下。下面的例子顯示二者區(qū)別：arcgisscripting,；通過win32創(chuàng)建的geoprocessor可以在不同的Python版本上運(yùn)行。在后面地理處理中我們會(huì)接觸一些例子。 運(yùn)行下面的代碼，說明了一種while循環(huán)：x=1while,x10:, print,x, x=x+1屏幕依次輸出1~9216。import,if,(d:/testfolder/):,print,測(cè)試文件夾存在,print,txt文件存在elif,(d:/testfolder):,print,測(cè)試文件夾存在,print,測(cè)試文件夾存在，但txt文件不存在else:,print,兩者都不存在可選探索示例接下來的例子做的事情對(duì)GIS非常重要，但是實(shí)際上不用任何地理處理代碼。山影（hillshade）工具需要有太陽高度角和方位角作為輸入?yún)?shù)。math和random模塊很多常用的數(shù)學(xué)計(jì)算功能都可以通過math找到，比如三角計(jì)算或?qū)?shù)計(jì)算，如果要使用復(fù)雜數(shù)字，就使用cmath模塊。sKidsp3=r39。print,p2d:\work\反斜線“\”和某些字母一起有特殊用法，如\n為換行，“\”為轉(zhuǎn)義字符，如“\\”則表示“\”print,39。)print,plist[39。zhulj39。警告：不要嘗試更新隨ArcGIS安裝的Python到一個(gè)新的版本！下面介紹Python的一些簡(jiǎn)單語法和規(guī)則。我們發(fā)Email，寫文章或者設(shè)計(jì)地圖，都需要和計(jì)算機(jī)交互，而處理一系列數(shù)據(jù)，我們需要自動(dòng)化和利用邏輯來指導(dǎo)自動(dòng)化。其中第一部分參考了《Python精要參考（第二版）》、《Python編程金典（讀書筆記）》等書籍文獻(xiàn)。Python腳本使用詳解目錄寫在前面的話 2前言 2一、PYTHON語言基礎(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。Python。我想自動(dòng)化和邏輯是關(guān)鍵，它們區(qū)別于我們多數(shù)使用計(jì)算機(jī)時(shí)的交互活動(dòng)。ArcGIS10中引入了全新的Python,Window來增強(qiáng)內(nèi)嵌的Python體驗(yàn)。.capitalize()Zhulj()即將capitalize()方法用于ss=39。/39。d:\\work\\39。Jerry39。當(dāng)要引用多個(gè)模塊是，中間用逗號(hào)分隔，比如：import,arcgisscripting,sys,string,os,math我們也可以自己為頻繁使用的方法創(chuàng)建Module，下面，我們開始體驗(yàn)內(nèi)置的Modules。你想要?jiǎng)?chuàng)建一系列山影柵格來代表夏天、冬天和春秋分?！痠f（continued）接下來，我們會(huì)探索一下另一個(gè)方便的模塊：: 開始之前，在d:/下創(chuàng)建一個(gè)“testfolder”文件夾，然后新建一個(gè)“”文件； 嘗試以下代碼段，確保print語句前有縮進(jìn)。,DEM?)if,unitchar==1:unit=feetif,unitchar==2:unit=metersprint,\nElevation,in+,+unit()輸出結(jié)果：Elevation,in,feetwhile（continued）216。當(dāng)循環(huán)一個(gè)數(shù)據(jù)集時(shí)（GIS中很常用的工作）while循環(huán)很有用。Python1如何創(chuàng)建地理處理對(duì)象（geoprocessor,object）所有g(shù)eoprocessing的Python腳本都可以通過import,arcgisscripting或者win32去穿件geoprocessor,object。② Geoprocessor,Programming,Model（PDF），包含方法（左箭頭表示）、屬性（可讀寫的表示為杠鈴形，只讀的表示為部分杠鈴形）：如何使用Geoprocessor,Programming,Model中的ListsLists（列表）及其屬性和方法在圖表中用紫色標(biāo)出，如下：現(xiàn)在我們?cè)囍帉懸欢文_本列舉出屬性表中所有屬性值（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中相同，但是需要提供工具名稱來使用它們。為了在ModelBuilder中使用腳本或?qū)⒛_本當(dāng)做ArcToolbox中工具使用，我們需要考慮如何給它輸入值以及讓其設(shè)置輸出值。然后，我們需要將腳本加進(jìn)工具（Making,a,script,into,a,tool），那樣才能在ArcToolbox或ModelBuilder或Command,line中使用。所以，記得移動(dòng)時(shí)將腳本工具文件和腳本本身一起拷貝。216。首先，添加hillshade工具，雙擊指定一個(gè)elevation柵格數(shù)據(jù)（這里我選擇了marbles文件夾下的elevation），用下拉條指定azimuth和altitude值為azimuth和sun,angle。嘗試輸入?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)試我們的程序了。我喜歡的做法是定義一個(gè)‘sendmsg’函數(shù)來輸出消息：def,sendmsg(msg): print,msg (msg)….sendmsg(Finished,Cliping...)216。然后，一系列調(diào)試工具出現(xiàn)；一個(gè)黃色三角出現(xiàn)在第一行代碼前。woodelev39。lugrid39。216。+str()print,envel,，Extent是一個(gè)Object，而clip語句的第二個(gè)參數(shù)Rectangle是一個(gè)Envelope類型數(shù)據(jù)，所以這里我使用了envel接收Extent的四個(gè)值。上面的方法無法測(cè)試一個(gè)數(shù)據(jù)的字段是否存在，但線面的listFields枚舉是個(gè)技巧：(listFields返回一個(gè)Python屬性對(duì)象列表)if,not,(“”,”stcode”):我們將在后面使用這句判斷字段是否存在。建議：(“g*”)，因?yàn)槟阈枰獙ふ襝overage，在你處理列表里成員時(shí)，可以將成員名稱賦給變量f，(f+”/polygon”)查找多邊形要素集。唯一的遺憾是這工具名字太長(zhǎng)，所以我已經(jīng)把它嵌入一個(gè)自定義函數(shù)中，命名為“ma”，當(dāng)然函數(shù)中包括開始和結(jié)束“Spatial”擴(kuò)展功能。有些工具根據(jù)屬性選擇記錄，或復(fù)制或刪除選中記錄。y值。 嘗試下面代碼，展示了InsertCursor的用法，同樣的基本方法可用來讀取外部數(shù)據(jù)，比如文本文件。但是到目前為止我們并沒有花太多時(shí)間在這些循環(huán)上。“臨時(shí)工作空間”專門用于存放不愿保留的輸出數(shù)