基于python的图属互操作

team470

　　在实际的业务需要中你可能有这样的要求将空间要素的所有节点按照一定的规则存储到属性字段中，然后使用后面的业务逻辑处理这些节点的坐标信息，最后在通过这些数据字段中节点信息的变化完成图形信息的变更操作。本人使用python对该种业务方式进行了初步研究。
　　1.首先要将节点信息写入属性字段
　　import string, os, sys, locale, arcgisscripting
gp = arcgisscripting.create()
gp.overwriteoutput = 1
　　msgNotEnoughParams = "Incorrect number of input parameters."
msgUseValidDecimalPointSep = "Please use one of the valid decimal point separators."
　　try:
　　if len(sys.argv) < 2: raise Exception, msgNotEnoughParams
inputFC = sys.argv[1]
cf=sys.argv[2]
tmp=str(cf)
inDesc = gp.describe(inputFC)
inRows = gp.UpdateCursor(inputFC)
inRow = inRows.next()
outLine=''
flag=False
partflag=True
while inRow:
feat = inRow.GetValue(inDesc.ShapeFieldName)
if inDesc.ShapeType.lower() == "point":
pnt = feat.getpart()
outLine = "{"+str(pnt.x) + "," + str(pnt.y) + "," + str(pnt.z) + "," + str(pnt.m)+";"

elif inDesc.ShapeType.lower() == "multipoint":
partnum = 0
partcount = feat.partcount
outLine=outLine+"{"
while partnum < partcount:
pnt = feat.getpart(partnum)
outLine = outLine+ str(pnt.x) + "," + str(pnt.y) + "," + str(pnt.z) + "," + str(pnt.m)+";"
partnum += 1

else:
partnum = 0
partcount = feat.partcount
while partnum < partcount:
part = feat.getpart(partnum)
if partflag:
outLine ="{"+outLine
else :
outLine =outLine+"{"
pnt = part.next()
if flag:
tmp1=outLine[:-1]
outLine =tmp1+"];"
flag=False
while pnt:
outLine =outLine+ str(pnt.x) + "," + str(pnt.y) + "," + str(pnt.z) + "," + str(pnt.m)+";"
pnt = part.next()
if not pnt:
pnt = part.next()
if pnt:
if flag:
tmp3=outLine[:-1]
outLine=tmp3+"];"
flag=True
outLine=outLine+"["+str(pnt.x) + "," + str(pnt.y) + "," + str(pnt.z) + "," + str(pnt.m)+";"
pnt = part.next()
if flag:
flag=False
tmp2=outLine[:-1]
outLine =tmp2+"];"
else:
tmp5=outLine[:-1]
outLine =tmp5+"}"
if partflag:
tmp4=outLine[:-1]
outLine =tmp4+"}"
partflag=False
partnum += 1

outLine=outLine[:-1]+"}"
inRow.SetValue(tmp,outLine)
inRows.UpdateRow(inRow)
outLine=''
inRow = inRows.next()
del inRow
del inRows
except Exception, ErrorDesc:
gp.AddError(ErrorDesc[0])
gp.AddError(gp.getmessages(2))

　　2.通过变更操作修改其中的坐标点信息，然后更新图形
　　import arcgisscripting
import fileinput
import os
import string
gp = arcgisscripting.create(9.3)
infile = sys.argv[1]
tmp=sys.argv[2]
# Create update cursor for feature class
inDesc = gp.describe(infile)
rows = gp.UpdateCursor(infile)
row = rows.Next()
gp.AddMessage(inDesc.ShapeType.lower())
flag=False
while row:
lineArray = gp.CreateObject("Array")
partArray=gp.CreateObject("Array")
featureVertexArray = gp.CreateObject("Array")
if inDesc.ShapeType.lower() == "polygon":
txt=str(row.GetValue(tmp))
strlist=string.split(txt,'}')
temstr=''
for value in strlist:
arr=string.split(value,';')
for pntstr in arr:
if len(pntstr)>1:
if pntstr[0:1]=='{':
temstr=pntstr[1:]
elif pntstr[0:1]=='[':
temstr=pntstr[1:]
elif pntstr[-1]==']':
temstr=pntstr[:-1]
flag=True
else:
temstr=pntstr
a=string.split(temstr,',')
pnt = gp.createobject("Point")
if len(a)==4:
pnt.x=string.atof(a[0])
pnt.y=string.atof(a[1])
pnt.z=string.atof(a[2])
pnt.m=string.atof(a[3])
lineArray.add(pnt)
gp.AddMessage(pnt.x)
gp.AddMessage(pnt.y)
if flag==True:
flag=False
featureVertexArray.add(lineArray)
lineArray.removeAll()
gp.AddMessage("newpart")
featureVertexArray.add(lineArray)
lineArray.removeAll()
polyGeom = gp.createobject("geometry", "polygon", featureVertexArray)
row.Shape=polyGeom
rows.UpdateRow(row)
row = rows.Next()
elif inDesc.ShapeType.lower() == "point":
pnt = gp.CreateObject("Point")
txt=str(row.GetValue(tmp))
arrp=string.split(txt,",")
x=arrp[0]
m=arrp[3]
pnt.x=string.atof(x[1:])
pnt.y=string.atof(arrp[1])
pnt.z=string.atof(arrp[2])
pnt.m=string.atof(m[:-1])
row.Shape=pnt
rows.UpdateRow(row)
row = rows.Next()
elif inDesc.ShapeType.lower() == "multipoint":
txt=str(row.GetValue(tmp))
strlist=string.split(txt,'}')
temstr=''
for value in strlist:
arr=string.split(value,';')
for pntstr in arr:
if pntstr[0:1]=='{':
temstr=pntstr[1:]
else:
temstr=pntstr
a=string.split(temstr,',')
pnt = gp.createobject("Point")
if len(a)==4:
pnt.x=string.atof(a[0])
pnt.y=string.atof(a[1])
pnt.z=string.atof(a[2])
pnt.m=string.atof(a[3])
lineArray.add(pnt)
gp.AddMessage(pnt.x)
gp.AddMessage(pnt.y)
gp.AddMessage("newpart")
row.Shape=lineArray
rows.UpdateRow(row)
row = rows.Next()
elif inDesc.ShapeType.lower() == "polyline":

txt=str(row.GetValue(tmp))
strlist=string.split(txt,'}')
temstr=''
for value in strlist:
arr=string.split(value,';')
for pntstr in arr:
if pntstr[0:1]=='{':
temstr=pntstr[1:]
else:
temstr=pntstr
a=string.split(temstr,',')
pnt = gp.createobject("Point")
if len(a)==4:
pnt.x=string.atof(a[0])
pnt.y=string.atof(a[1])
pnt.z=string.atof(a[2])
pnt.m=string.atof(a[3])
lineArray.add(pnt)
gp.AddMessage(pnt.x)
gp.AddMessage(pnt.y)
gp.AddMessage("newpart")
featureVertexArray.add(lineArray)
lineArray.removeAll()
polyGeomline = gp.createobject("geometry", "polyline", featureVertexArray)
row.Shape=polyGeomline
rows.UpdateRow(row)
row = rows.Next()
# Delete cursor and row objects to remove locks on the data
del row
del rows
　　

　　测试中我采用编辑图形的方式完成节点坐标信息的修改然后将图形通过上面的代码变更回原来的要素（应为第一步我已经将要素的节点信息存储到属性表中）

　　通过工具变更后为

　　通过以上的脚本最终完成数据的变更初步测试了最负责的面包括三种情况1：只有一个部分2：包括多个部分；3包括环的情况，以上脚本还在进一步测试希望大家多多指正