使用坐标点生成shp图斑

以下代码由PHP站长网 52php.cn收集自互联网

现在PHP站长网小编把它分享给大家，仅供参考

?#coding:utf-8
import arcpy
from arcpy import env
env.workspace = 'g:/python'    ## 设置工作空间
file = open('11.txt')        ## 读取text文件
def create_polygon(coord_l):
    point = arcpy.Point()
    array = arcpy.Array()
    featureList = []
    for feature in coord_l:
        for coord in feature:
            point.X = coord[0]
            point.Y = coord[1]
            array.add(point)

        array.add(array.getObject(0))
        polygon = arcpy.Polygon(array)
        array.removeAll()
        featureList.append(polygon)
    return featureList

li = []
for line in file:
    li.append(map(float,line.split()))  ##将坐标导入列表
coord_l = []
code_l = []
for i in range(len(li)):
    if len(li[i]) == 2:
        code_l.append(str(int(li[i][0]))+str(int(li[i][1])))
        coord_l.append([])
    else:
        coord_l[len(coord_l)-1].append([li[i][1],li[i][2]])

poly = create_polygon(coord_l)
arcpy.CopyFeatures_management(poly,'g:/python/polygons.shp')  ##设置生成的shp文件名
arcpy.AddField_management('polygons.shp','code','TEXT')
cur = arcpy.UpdateCursor('polygons.shp')
i = 0
for row in cur:
    row.code = code_l[i]
    cur.updateRow(row)
    i = i + 1
del cur,row



    
    

以上内容由PHP站长网【52php.cn】收集整理供大家参考研究

如果以上内容对您有帮助，欢迎收藏、点赞、推荐、分享。

（编辑：李大同）

【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容!