分享

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

 昵称QAb6ICvc 2014-07-18

也许你会有这样的需求,是要把面数据的所有拐点坐标提取出来存储到一个属性字段中使用,但是又不想大费周折的写一些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

    本站是提供个人知识管理的网络存储空间,所有内容均由用户发布,不代表本站观点。请注意甄别内容中的联系方式、诱导购买等信息,谨防诈骗。如发现有害或侵权内容,请点击一键举报。
    转藏 分享 献花(0

    0条评论

    发表

    请遵守用户 评论公约

    类似文章 更多