[初学入门]ArcGIS中Python脚本学习(9)

import arcgisscripting,os,sys gp=arcgisscripting.create(9.3) wspath=\txtfile=\

fpath=wspath+\

if gp.exists(fpath):os.remove(fpath) txtfile=open(fpath,\gp.workspace=wspath

rows=gp.SearchCursor(\row=rows.Next() while row:

feat=row.shape a=0

while a

geomArray=feat.getpart(a) fid=row.GetValue(\ geomArray.Reset

pnt=geomArray.Next() while pnt:

txtfile.write(str(fid) + \ pnt=geomArray.Next() a=a+1

row=rows.next() txtfile.close()

挑战2：

写一个可以读取刚才输出的文本文件并且创建一个新的shp文件的脚本。 如果你使用的是ArcGIS10，这里有最新的帮助，如果使用的是ArcGIS9.3.1或9.3，这里有相应帮助。

import arcgisscripting,sys,fileinput gp=arcgisscripting.create(9.3)

inFile=\gp.workspace=\fcName=\gp.overwriteoutput=1

gp.CreateFeatureClass_management(gp.workspace,fcName,\try:

cur=gp.InsertCursor(fcName) for line in fileinput.input(inFile): values=line.split(\ id=int(values[0]) x=float(values[1]) y=float(values[2]) feat=cur.NewRow()

31 / 33

pnt=gp.CreateObject(\ pnt.x=x pnt.y=y

feat.shape=pnt feat.id=id

cur.InsertRow(feat) print id,x,y del cur

fileinput.close() except:

print gp.getmessages() if cur: del cur

接下来该做什么呢？

我们已经做了很多，但是大多数都是创建处理特定事情的部分，比如遍历一系列数据。但是到目前为止我们并没有花太多时间在这些循环上。

另一间我们没有花太多时间处理的事情是，让脚本工具更好地工作。输入和输出是这个处理过程中有挑战性的部分，所以在这方面下手可能会更适合你。我已经发现这是最令人困惑的部分，所以我将能发现的尽可能多的输入/输出情况做了总结，请看附录1。