arcpy利用XY创建点

偏执的太偏执、 2022-02-23 09:56 372阅读 0赞

参考了很多博客吧，最后把自己都改蒙了，这个博客有分析[arcpy创建带高程的点图层、][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 path
        if os.path.isfile(path):        
    # Set the local variables
            outputname = path[40:57]+"_"+path[58:73]+".shp"
            print outputname
            spatRef = 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")#m
            arcpy.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 fc
            num=0
            for 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 feature
                vertex = arcpy.CreateObject("Point")
                vertex.X = float(x)
                vertex.Y = float(y_coords[num])
                feature.shape = vertex
                feature.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.speed
                cursor.insertRow(feature)
                num=num+1
            arcpy.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]: https://blog.csdn.net/A873054267/article/details/88207732