关于arcpy获取每个面要素中的折点坐标一、获取思路二、实现方法三、总结和反思

124 阅读 0 评论 82 点赞

我是靠谱客的博主执着人生，这篇文章主要介绍关于arcpy获取每个面要素中的折点坐标一、获取思路二、实现方法三、总结和反思，现在分享给大家，希望可以做个参考。

#关于arcpy获取每个面要素中的折点坐标

一、获取思路

通过使用arcpy里面的“游标”:arcpy.da对要素类里面的各种要素进行遍历和获取部分值，其中需要注意的是，在使用arcpy.da.SearchCursor()函数时，需要对SearchCursor(infeature,[fieldsList])中的字段进行元组取值，若不取值则会导致解包Array数组十分麻烦。另外，在解包

二、实现方法

实现方法1：

复制代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
with arcpy.da.SearchCursor(inputfeature,["SHAPE@"]) as cur:
flag=1
for rows in cur:
for elems in rows:
part=elems.getPart()
numOfpart=elems.pointCount
print(part)
print(numOfpart)
print("outer boundary:")
for fea in part:
for i in range(numOfpart):
if fea[i]:
x=fea[i].X
y=fea[i].Y
else:
print("interior rings:")

实现方法2：

复制代码

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
with arcpy.da.SearchCursor(curFeature,["shape@"]) as cur:
flag=1
for rows in cur:
print(rows)
for item in rows[0]:
part=item.count
print(part)
print("outer boundary:")
for i in range(part):
if item[i]:
x=item[i].X
y=item[i].Y
print(str(x)+";"+str(y))
else:
print("interior rings:")