# -*- coding: cp936 -*-
#本脚以最左边、Y值最大的点为起始点按顺时针为多边形节点编码,生成一个包含记录编码值和多边形FID字段的点要素类

#注意:
#1.本脚本作为arcgis脚本工具使用、脚本测试版本为10.0中文,未装sp5补丁
#2.暂不支持带环多边形,处理带环多边形程序直接崩溃
#3.本工具脚本用到系统工具"多部分(Multipart)至单部分(Singlepart)",因该工具无法完美处理所有的CAD多边形要素类,
#  所以本脚本工具只支持简单CAD多边形要素类。

#Creater: ych
#Create Time:2013年5月1日

import arcpy
import os

#获取单部分多边形最右边且Y值最大的点(边线相交的多边形不适合)
def getRightupperPoint(fc): 
    points = [point for point in fc.getPart(0)]
    XValues = [point.X for point in points]
    Xmax_indexs = [XValues.index(x) for x in XValues if x == max(XValues)]
    Ymax_index = [i for i in Xmax_indexs if points[i].Y == max([points[j].Y for j in Xmax_indexs ])][0]
    RightupperPoint = points[Ymax_index]
    return   RightupperPoint,Ymax_index

#判断单部分多边形的方向(True 为顺时针,False为逆时针 圆、椭圆返回-1)     
def isclockwise(fc): 
     #找出最右且Y值最大的节点
     points = [point for point in fc.getPart(0)]
     if len(points)!= 2:
        rightupperPoint,rightupperPoint_index = getRightupperPoint(fc)
    
        #如果最右且Y值最大的节点的后一个点(以在points列表中的顺序为参考顺序 )的Y值,比最右且Y值 \
        #最大的节点的前一个点的Y值大,则多边形为顺时针,否则为逆时针     
        if rightupperPoint_index == len(points)-1:
           BehindPt = points[1]
        else:
           BehindPt = points[rightupperPoint_index+1]
       
        prePt = points[rightupperPoint_index-1]
    
        if prePt.X == rightupperPoint.X:
            return True
        else:
            y = (prePt.Y - rightupperPoint.Y)*(BehindPt.X - rightupperPoint.X)/(prePt.X - rightupperPoint.X ) + rightupperPoint.Y

if y <BehindPt.Y:
             return True
        else:
             return False
     else:
        return -1
         
           
       
#将单部分多边形要素的起始点只设置为最右且Y值最大的点
def changeStartPoint(fc):
      points = [point for point in fc.getPart(0)]
      if len(points)!= 2:
         #计算多边形最右边且Y值最大的点在Points中索引
         index_RightupperPoint = getRightupperPoint(fc)[1]
         #计算更改了起始点的多边形要素的newpoints
         newpoints =[]
         newpoints.extend(points[index_RightupperPoint:-1])
         newpoints.extend(points[:index_RightupperPoint+1])
         #创建新要素
         newPolygon = arcpy.Polygon(arcpy.Array(newpoints))
         return newPolygon
      #len(points) ==2为圆、椭圆 
      else:
         return fc

#主程序
def main():
  #设置参数
  fcl_input = arcpy.GetParameterAsText(0)  #输入多边形要素类
  out_PointClass = arcpy.GetParameterAsText(1) #输出点要素类

#创建空的、包含用于存储节点编号的子段的点要素类
  name = os.path.basename(out_PointClass)
  path = out_PointClass[:-(len(name)+1)]
  desc = arcpy.Describe(fcl_input)
  pointClass = arcpy.CreateFeatureclass_management(path,name,"POINT","","","",desc.spatialReference)
 
  #添加字段
  arcpy.AddField_management(pointClass,"NodeCode","SHORT")
  arcpy.AddField_management(pointClass,"PolygonFID","LONG")

#多部分至单部分
  scrathName = arcpy.CreateScratchName ("TEMP","","",path)
  fcl_input=arcpy.MultipartToSinglepart_management (fcl_input,scrathName)

#创建cursor
  rows = arcpy.InsertCursor(pointClass)
  rows_fc = arcpy.SearchCursor(fcl_input)
 
 
  #输入多边形要素类的FID字段
  fcID = arcpy.ListFields(fcl_input,"","OID")

for row_fc in rows_fc:
        fc = row_fc.shape
        #将多边形要素的起始点只设置为最右且Y值最大的点
        fc = changeStartPoint(fc)
        #获取多边形节点列表
        points = [point for point in fc.getPart(0)]
        #如果多边形为逆时针,则翻转节点列表,使多边形变成顺时针方向
        if isclockwise(fc)==False: 
           points.reverse()   
        else:
           pass
        #将顺时针方向的多边形点列表写到多边形要素类中
        for point in  points[:-1]:    
           newrow = rows.newRow()
           newrow.shape = point
           newrow.NodeCode =points.index(point)+1
           newrow.PolygonFID = eval("row_fc."+fcID[0].name)
           rows.insertRow(newrow)
        arcpy.AddMessage("Nodes of Polygon %s has been successfully decoded."% eval("row_fc."+fcID[0].name))    
  #删除游标
  del rows,rows_fc,row_fc
  #删除临时文件
  arcpy.Delete_management (fcl_input)

if __name__ == "__main__":

main()

来自:http://bbs.esrichina-bj.cn/esri/viewthread.php?tid=127216

多边形节点编码python脚本的更多相关文章

  1. 关于Python脚本开头两行的:#!/usr/bin/python和# -*- coding: utf-8 -*-的作用 – 指定文件编码类型

    #!/usr/bin/python指定用什么解释器运行脚本以及解释器所在的位置 # -*- coding: utf-8 -*-用来指定文件编码为utf-8的PEP 0263 -- Defining P ...

  2. 【转】关于Python脚本开头两行的:#!/usr/bin/python和# -*- coding: utf-8 -*-的作用 – 指定文件编码类型

    原文网址:http://www.crifan.com/python_head_meaning_for_usr_bin_python_coding_utf-8/ #!/usr/bin/python 是用 ...

  3. 【转载】关于Python脚本开头两行的:#!/usr/bin/python和# -*- coding: utf-8 -*-的作用 – 指定文件编码类型

    1.#!/usr/bin/python 是用来说明脚本语言是 python 的 是要用 /usr/bin下面的程序(工具)python,这个解释器,来解释 python 脚本,来运行 python 脚 ...

  4. arcgis python脚本工具实例教程—栅格范围提取至多边形要素类

    arcgis python脚本工具实例教程-栅格范围提取至多边形要素类 商务合作,科技咨询,版权转让:向日葵,135-4855_4328,xiexiaokui#qq.com 功能:提取栅格数据的范围, ...

  5. win下python脚本以unix风格换行保存将会报错为编码问题 SyntaxError: encoding problem:gbk

    utf-8与gbk编码都报错 从别人的github拉下来一个python脚本. 直接运行,python报错如下: File ".\drag_files_do_event.py", ...

  6. Python脚本开头两行的:#!/usr/bin/python和# -*- coding: utf-8 -*-的作用

    #!/usr/bin/Python指定用什么解释器运行脚本以及解释器所在的位置 # -*- coding: utf-8 -*-用来指定文件编码为utf-8的 估计有不少人注意过一些python脚本开头 ...

  7. python脚本实现集群检测和管理

    python脚本实现集群检测和管理 场景是这样的:一个生产机房,会有很多的测试机器和生产机器(也就是30台左右吧),由于管理较为混乱导致了哪台机器有人用.哪台机器没人用都不清楚,从而产生了一个想法-- ...

  8. Python脚本控制的WebDriver 常用操作 <十七> 获取测试对象的属性及内容

    测试用例场景 获取测试对象的内容是前端自动化测试里一定会使用到的技术.比如我们要判断页面上是否显示了一个提示,那么我们就需要找到这个提示对象,然后获取其中的文字,再跟我们的预期进行比较.在webdri ...

  9. python 脚本查看微信把你删除的好友--win系统版

    PS:目测由于微信改动,该脚本目前不起作用 下面截图来自原作者0x5e 相信大家在微信上一定被上面的这段话刷过屏,群发消息应该算是微信上流传最广的找到删除好友的方法了.但群发消息不仅仅会把通讯录里面所 ...

随机推荐

  1. Oracle数据库安装后,登录Database Control时密码错误

    解决方案1(实测可行): sys 和 system用户的用户名和密码还记得不?试试看如果不行,用sqlplus 在服务器本地登录%sqlplus / as sysdbaSQL>alter use ...

  2. SOAP Services for Python

    转自:https://pypi.python.org/pypi/SOAPpy/ 1.下载SOAPpy 0.12.5 2.解压tar -xvzf SOAPpy-$VERSION$.tar.gz 3.安装 ...

  3. HDU 1561-The more, The Better(树状背包)

    题意: n个城堡,每个有一定的财富,有些城堡不能直接攻克,要攻克这些城堡必须先攻克其他某一个特定的城堡,求应攻克哪m个城堡,获得最大财富. 分析:dp[i][j],以i为根的子树,攻克j个城堡,获得的 ...

  4. HDU 1520-Anniversary party(树形dp入门)

    题意: n个人参加party,已知每人的欢乐值,给出n个人的工作关系树,一个人和他的顶头上司不能同时参加,party达到的最大欢乐值. 分析:dp[i][f],以i为根的子树,f=0,i不参加,f=1 ...

  5. Loadrunner通过sitescope监控mysql

    1.   Loadrunner通过sitescope监控mysql 1.1. Sitescope下载 第一次在网上下载的是sitescope11.30,安装之后无法破解,所以最后卸载了,重新下载了si ...

  6. 高质量、处于持续更新的R包

    本文在Creative Commons许可证下发布 自由软件的问题是开发人员没有稳定的资金来源支持,可能更新上做不到持续.经过考证和圈内朋友的帮助,现在把R包中高质量.持续更新的跟大数据事业相关的R包 ...

  7. 【DOM】学习笔记

    三. 一份文档就是一颗节点树 节点类型:元素节点——属性节点.文本节点 getElementById()返回一个对象,对应一个元素节点 getElementByTagName()返回一个对象数组,分别 ...

  8. bzoj 2716 天使玩偶(CDQ分治,BIT)

    [题目链接] http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=29234 [题意] 询问当前点与已知点的最小曼哈顿距离. [思路 ...

  9. INPUT输入框灰体提示

    INPUT输入框灰体提示 <input type="text" value='15 words limit' style="color:#999999" ...

  10. FIREDAC数据引擎

    以前使用过BDE.ADO.DBX等数据引擎,后来发现它们都没有UNIDAC好用, 所以在很长的一段时间内中间件都使用UNIDAC作为数据引擎. 偶然的机会,使用了DELPHI XE5自带的FIREDA ...