ArcGIS10.1 python 提取面的所有拐点坐标到属性字段中

也许你会有这样的需求,是要把面数据的所有拐点坐标提取出来存储到一个属性字段中使用,但是又不想大费周折的写一些AO、AE甚至是addin的程序,那么你可以选择是哟个python来完成你的工作。示例代码如下,在arcgis10.1中通过测试

def MySub(feat):
partnum = 0
#multipart feature
partcount = feat.partCount
pntcount = 0
str=”
# Enter while loop for each part in the feature (if a singlepart feature
# this will occur only once)
while partnum < partcount:
part = feat.getPart(partnum)
pnt = part.next()

# Enter while loop for each vertex
#
str=str+”[“
while pnt:
pntcount += 1
px=’%f’ %pnt.x
py=’%f’ %pnt.y
str=str+px+”,”+py +”;”
pnt = part.next()

# If pnt is null, either the part is finished or there is an
# interior ring
if not pnt:
str=str[:-1]
str=str+”]”
pnt = part.next()
partnum += 1
return str

转载自:https://blog.csdn.net/iteye_8029/article/details/82486208

You may also like...