arcpy利用XY创建点

参考了很多博客吧，最后把自己都改蒙了，这个博客有分析arcpy创建带高程的点图层、

# -*- coding: utf-8 -*- """ Created on Sun Apr 7 15:32:24 2019@author: """# XYTableToPoint.py # Description: Creates a point feature class from input table# import system modules import arcpy import csv # Set environment settings import os rootdir = outpath=list = os.listdir(rootdir) #列出文件夹下所有的目录与文件 for i in range(0,len(list)):path = os.path.join(rootdir,list[i])print pathif os.path.isfile(path): # Set the local variablesoutputname = path[40:57]+"_"+path[58:73]+".shp"print outputnamespatRef = arcpy.SpatialReference("WGS 1984") #"WGS 1984"in_table = path with open(in_table,'rb') as csvfile:reader = csv.reader(csvfile)column1 = [row[0] for row in reader]with open(in_table,'rb') as csvfile:reader = csv.reader(csvfile) column2 = [row[1] for row in reader]with open(in_table,'rb') as csvfile:reader = csv.reader(csvfile)end_lon = [row[2] for row in reader]with open(in_table,'rb') as csvfile:reader = csv.reader(csvfile) end_lat = [row[3] for row in reader]with open(in_table,'rb') as csvfile:reader = csv.reader(csvfile)distance = [row[4] for row in reader]with open(in_table,'rb') as csvfile:reader = csv.reader(csvfile) speed = [row[5] for row in reader]with open(in_table,'rb') as csvfile:reader = csv.reader(csvfile)angle = [row[6] for row in reader]with open(in_table,'rb') as csvfile:reader = csv.reader(csvfile) u = [row[7] for row in reader]with open(in_table,'rb') as csvfile:reader = csv.reader(csvfile) v = [row[8] for row in reader]out_feature_class = "ice vector"x_coords = column1 #"longitude"y_coords = column2 # "latitude"fc=arcpy.CreateFeatureclass_management(outpath, outputname, "POINT", "","","", spatRef)arcpy.AddField_management(outpath+"\\"+outputname, "start_lon", "DOUBLE")arcpy.AddField_management(outpath+"\\"+outputname, "start_lat", "DOUBLE")arcpy.AddField_management(outpath+"\\"+outputname, "end_lon", "DOUBLE")arcpy.AddField_management(outpath+"\\"+outputname, "end_lat", "DOUBLE")arcpy.AddField_management(outpath+"\\"+outputname, "distance", "DOUBLE")#marcpy.AddField_management(outpath+"\\"+outputname, "speed", "DOUBLE")arcpy.AddField_management(outpath+"\\"+outputname, "angle", "DOUBLE")arcpy.AddField_management(outpath+"\\"+outputname, "u", "DOUBLE")arcpy.AddField_management(outpath+"\\"+outputname, "v", "DOUBLE")cursor=arcpy.InsertCursor(fc,["SHAPE@XY"]) # Add the point geometry to the feature vertex = arcpy.CreateObject("Point") vertex.X = line[0] vertex.Y =line[1] feature.shape = vertex # Add attributes feature.leibie = "shumu" # write to shapefile cursor.insertRow(feature) del cursor del fcnum=0for x in x_coords: feature = cursor.newRow() # Add the point geometry to the feature vertex = arcpy.CreateObject("Point") vertex.X = line[0] vertex.Y =line[1] feature.shape = vertex # Add attributes feature.leibie = "shumu" # write to shapefile cursor.insertRow(feature) del cursor del fc# Add the point geometry to the featurevertex = arcpy.CreateObject("Point")vertex.X = float(x)vertex.Y = float(y_coords[num])feature.shape = vertexfeature.start_lon = float(x)feature.start_lat = float(y_coords[num])feature.end_lon = float(x)feature.end_lat = float(y_coords[num])feature.distance = float(distance[num])feature.speed = float(speed[num])feature.angle = float(angle[num])feature.u = float(u[num])feature.v = float(v[num])#print vertex.X,vertex.Y,vertex.speedcursor.insertRow(feature)num=num+1arcpy.DefineProjection_management(outpath+"\\"+outputname,spatRef)# Add attributes #z_coords = "elevation"# Make the XY event layer... # ============================================================================= # =============================================================================# newfc="newpoint.shp" # arcpy.CreateFeatureclass_management(outpath, newfc, "Point") # cursor=arcpy.da.InsertCursor(newfc, ["SHAPE@"]) # array=arcpy.Array()# # # Print the total rows # print(arcpy.GetCount_management(out_feature_class)) # ============================================================================= # =============================================================================

 

总结

以上是生活随笔为你收集整理的arcpy利用XY创建点的全部内容，希望文章能够帮你解决所遇到的问题。

如果觉得生活随笔网站内容还不错，欢迎将生活随笔推荐给好友。