ArcGIS10.1 python 提取面的所有拐点坐标到属性字段中
分享
也许你会有这样的需求,是要把面数据的所有拐点坐标提取出来存储到一个属性字段中使用,但是又不想大费周折的写一些AO、AE甚至是addin的程序,那么你可以选择是哟个python来完成你的工作。示例代码如下,在arcgis10.1中通过测试
文章来源:http://blog.csdn.net/esrichinacd/article/details/8772380
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
文章来源:http://blog.csdn.net/esrichinacd/article/details/8772380
0 个评论
推荐内容
相关问题
- 安装完arcgis10.3后,为何没有Envi tools这个工具箱?
- arcgis10.2有时候无法复制粘贴要素
- ArcGIS中如何用Python给字段值补充0
- 如何使用脚本或其他方法,将图层属性中属性域原值内容替换为属性域描述?
- 属性表联接excel为什么是空的?
- ArcGIS10.2线图层分割面图层
- arcgis10.4 离线部署js api 4.3 版本 无法显示
- 如何在EXCEL里把指定时间范围内值分字段求和
- 市区择房分析时,需要添加字段并将其赋值1或者-1,怎么做? 在开始编辑的时候,总是出现空间参考与数据框不匹配提示,原因是什么?会影响下面的赋值吗?
- arcgis10.2打开arcgis10.6的mxd文档
- 如何将117°35′35″这种格式的经纬度坐标导入arcmap?