首页 > 编程语言 > 详细

ArcGIS Python任意多边形

时间:2021-01-25 23:46:07      阅读:33      评论:0      收藏:0      [点我收藏+]
# -*- coding: utf-8 -*-
#by gisoracle 2021.01.23
import arcpy
import math

#保存数据
def Save(polygon,rows):

    row = rows.newRow()
    #row.setValue(shapefieldname,featureList)
    row.shape=polygon
    rows.insertRow(row)
    del row
outFeature=arcpy.GetParameterAsText(0)
x0=500000
y0=0
r=500.0
y2=500
n=9
h=int(360/n)
arcpy.AddMessage("h="+str(h))
pi=3.1415926
rows = arcpy.InsertCursor(outFeature)
array = arcpy.Array()
for i in range(0,360,h):
    x1=x0+r*math.cos(i*pi/180)
    y=y0+r*math.sin(i*pi/180)
    Z=n*100+360-i
    point = arcpy.Point(x1,y,Z)
    #point.X=x1
    #point.Y=y
    #point.Z=Z
    arcpy.AddMessage("Z="+str(Z))
    array.add(point)
polygon = arcpy.Polygon(array, None, True)
Save(polygon, rows)
polygon = None

 

ArcGIS Python任意多边形

原文:https://www.cnblogs.com/gisoracle/p/14327358.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!