Arcpy批量通过XY坐标创建点shp

# -*- coding: utf-8 -*-

# 批量通过xy创建shp

import arcpy
import csv  
# Set environment settings
import os
rootdir ="D:\\Desktop\\diantongji"
outpath="D:\\Desktop\\diantongjishp"

list = os.listdir(rootdir) #列出文件夹下所有的目录与文件
for i in range(0,len(list)):
    path = os.path.join(rootdir,list[i])
    print path
    if os.path.isfile(path):        
        outputname = list[i].split('.')[0]+"_XY.shp"
        print outputname
        spatRef = arcpy.SpatialReference(4326) #"WGS 1984"
        in_table = path 
        with open(in_table,'r') as csvfile:
            reader = csv.reader(csvfile)
            column1 = [row[1] for row in reader] 
        with open(in_table,'r') as csvfile:
            reader = csv.reader(csvfile) 
            column2 = [row[2] for row in reader]
        with open(in_table,'r') as csvfile:
            reader = csv.reader(csvfile)
            count=[row[3] for row in reader]
        
        out_feature_class = "ice vector"
        x_coords = column1  #"longitude"
        y_coords = column2  # "latitude"

#创建shp
        fc=arcpy.CreateFeatureclass_management(outpath, outputname, "POINT", "","","", spatRef)
        arcpy.AddField_management(outpath+"\\"+outputname, "count", "DOUBLE")

        cursor=arcpy.InsertCursor(fc,["SHAPE@XY"]) 

#向shp中添加数据
        num=0
        for x in x_coords: 
            feature = cursor.newRow() 

            vertex = arcpy.CreateObject("Point")
            vertex.X = float(x)
            vertex.Y = float(y_coords[num])
            feature.shape = vertex
            
            feature.count = float(count[num])
            #print vertex.X,vertex.Y,vertex.speed

            cursor.insertRow(feature)
            num=num+1
        arcpy.DefineProjection_management(outpath+"\\"+outputname,spatRef)
        print (list[i] + "         ---- OK! ----         ")
print("   ---   All project is OK!   ---   ")
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值