• 博客园logo
  • 会员
  • 众包
  • 新闻
  • 博问
  • 闪存
  • 赞助商
  • HarmonyOS
  • Chat2DB
    • 搜索
      所有博客
    • 搜索
      当前博客
  • 写随笔 我的博客 短消息 简洁模式
    用户头像
    我的博客 我的园子 账号设置 会员中心 简洁模式 ... 退出登录
    注册 登录

gisoracle

  • 博客园
  • 联系
  • 订阅
  • 管理

公告

View Post

ArcGIS Python任意多边形

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

 

posted on 2021-01-25 20:49  gisai  阅读(344)  评论(0)    收藏  举报

刷新页面返回顶部
 
博客园  ©  2004-2025
浙公网安备 33010602011771号 浙ICP备2021040463号-3