多边形节点编码python脚本

# -*- 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