nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Create XY Event Layer / Save to Layer file / Copy features to shapefile

BenSciance

I am attempting this simple workflow.

1) Adding an XY event layer of lat/lon points with z values from a text file.
2) Save the event layer to a layer file.
3) Copy the saved layer file to a new shapefile or feature class in geodatabase.

Now, I know the event layer that is created when you use the "create XY features" tool is a temporary file. I am having issues at the 2nd and 3rd step. The script completes, BUT the new layer file and shapefile are nowhere to be found. I have defined the
workspace path to my local machine but the files are not there after the script completes. Any advice would be greatly
appreciated.

Here is my code :

# MakeXYLayer.py
# Description: Creates an XY Event layer and saves it to a layer file then copies
# the layer file to a shapefile

# import system modules
import arcpy
from arcpy import env

# Set environment settings
env.workspace = "C:\Users\mbs7038\Documents\ArcGIS\matDEM"

try:
        # Set the local variables
        in_Table = "LatLonElev.txt"
        x_coords = "Lon"
        y_coords = "Lat"
        z_coords = "Elev"
        out_Layer = "ElevPts_layer"
        saved_Layer = r"C:\Users\mbs7038\Documents\ArcGIS\matDEM\ElevPts.lyr"

        # Set the spatial reference
        spRef = r("GCS_WGS_1984")

        # Make the XY event layer
        arcpy.MakeXYEventLayer_management(in_Table, x_coords, y_coords, out_Layer, spRef, z_coords)

        # Save to a layer file
        arcpy.SaveToLayerFile_management(out_Layer, saved_Layer)

        #===============================================================================
        # Copy features to shapefile

        arcpy.CopyFeatures_management("C:\Users\mbs7038\Documents\ArcGIS\matDEM\ElevPts.lyr", "C:\Users\mbs7038\Document         s\ArcGIS\matDEM\ElevPts.shp")
        #===============================================================================

except:
    # If an error occurred print the message to the screen
    print arcpy.GetMessages()

print 'script complete!!!!!'

Finde mehr Beiträge mit dem Tag

Accepted answers

JakeSkinner

The problem was with the coordinate system. Replace:

spRef = r("GCS_WGS_1984")

With:

spRef = arcpy.SpatialReference("WGS 1984")

Also, a good suggestion is to add print statements throughout the code to pinpoint where the code is failing.

All comments

JoshuaAhmed

Hi all,

I have a similar script where I am going through a file and converting all .txt files into events layers, and finally into points files. The points are X and Y coords. I have tried changing all manner of things yet I still get a Runtime error: object: error in executing tool message. Any ideas where the script is going wrong. I know it's something in the MakeXYEventLayer_management but cannot figure out where... Any help will be welcomed.

-------------------------------------------------------------------

import arcpy, sys, os

from arcpy import env

from arcpy.sa import *

arcpy.CheckOutExtension("Spatial")

arcpy.env.overwriteOutput = True #Set the script to be able to overwrite files

Path_table = "J:/Annual_time_series/NDVI_images/Beni1/Centrelines/" #path

for root, dirs, files in os.walk(Path_table):

for name in files:

beni3_reclass = os.path.join(root,name) #name

if beni3_reclass.endswith('_CL.txt') == True:

tilename = name[0:-4] #changed to keep my filenames

print tilename + " being processed."

table = Path_table + name

in_x_field = "cl_x"

in_y_field = "cl_y"

out_layer = Path_table + tilename + "_layer"

#saved_layer = Path_table + tilename + ".lyr"

#Convert txt file table to an events layer

arcpy.MakeXYEventLayer_management (table, in_x_field, in_y_field, out_layer, {"J:/Annual_time_series/NDVI_images/coordinate_system/WGS 1984 UTM Zone 19N.prj"}, {})

print arcpy.GetMessages()

#save to featureclass - this makes it a permenant shape points file

Path_points = "J:/Annual_time_series/NDVI_images/Beni1/Centrelines/CL_coords"

in_features = out_layer

out_featureclass = Path_points + tilename + "_CL_pts.shp"

arcpy.CopyFeatures_management (in_features, out_feature_class, {}, {}, {}, {})

-------------------------------------------------------------------------

The error message is:

Traceback (most recent call last):

File "E:/CLtxt_to_points.py", line 22, in <module>

arcpy.MakeXYEventLayer_management (table, in_x_field, in_y_field, out_layer, {"J:/Annual_time_series/NDVI_images/coordinate_system/WGS 1984 UTM Zone 19N.prj"}, {})

File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\management.py", line 6348, in MakeXYEventLayer

raise e

RuntimeError: Object: Error in executing tool

Thanks,

Josh

BenSciance

There seems to be a problem when defining the spatial reference as an input variable

spRef = arcpy.SpatialReference(4326)

AND

spRef = arcpy.SpatialReference("WGS 1984")

Since the "Make XY event layer" script has a direct input for the spatial reference, I wrote the code for WGS 1984 coordinate
system directly in the script (4326) instead.

 arcpy.MakeXYEventLayer_management(in_Table, x_coords, y_coords, out_Layer, 4326, z_coords)

After doing this, the entire script ran perfectly. The print statements kept the processing steps more organized as well. Thanks so much for your input!

import arcpy
from arcpy import env

# Set environment settings
env.workspace = r"C:\Users\mbs7038\Documents\ArcGIS\floodtest"

try:
    # Set the local variables
    #INPUTS
    in_Table = "TEST.txt"
    x_coords = "Lon"
    y_coords = "Lat"
    z_coords = "wlevel"

    #OUTPUTS
    out_Layer = "TEST_Layer"
    #saved_Layer = "TEST1.lyr"
    fc=r"C:\Users\mbs7038\Documents\ArcGIS\floodtest\AWESOME.gdb\TEST1"

    #===========================================================================
    # Set the spatial reference----DID NOT WORK WHEN DEFINING spRef.
    # USE COORDINATE SYSTEM CODE FOR "WGS 84" (4326) IN THE SCRIPT FOR MAKE XY EVENT LAYER
    #print "Set spatial reference"
    #spRef = arcpy.SpatialReference(4326)
    #===========================================================================

    #===========================================================================
    # Make the XY event layer
    # MakeXYEventLayer_management (table, in_x_field, in_y_field, out_layer, {spatial_reference}, {in_z_field})
    print "Make XY event layer"
    arcpy.MakeXYEventLayer_management(in_Table, x_coords, y_coords, out_Layer, 4326, z_coords)
    #===========================================================================

    #===========================================================================
    # Save to a layer file----DO NOT NEED TO PROCESS THIS STEP.  CAN SKIP TO COPY FEATURES
    #print "Saving layer file"
    #arcpy.SaveToLayerFile_management(out_Layer, saved_Layer)
    #===========================================================================

    #===========================================================================
    # Copy features to feature class
    print "Copying features to feature class"
    arcpy.CopyFeatures_management(out_Layer,fc)
    #===========================================================================


except:
    print arcpy.GetMessages()

JakeSkinner

Yes, this worked for me. Add some print statements to your code so you can see where it is failing at. Ex:

import arcpy
from arcpy import env

# Set environment settings
env.workspace = r"C:\Users\mbs7038\Documents\ArcGIS\matDEM"

try:
    # Set the local variables
    in_Table = "LatLonElev.txt"
    x_coords = "Lon"
    y_coords = "Lat"
    z_coords = "Elev"
    out_Layer = "ElevPts_layer"
    saved_Layer = "ElevPts.lyr"

    # Set the spatial reference
    print "Set spatial reference"
    spRef = arcpy.SpatialReference("WGS 1984")

    # Make the XY event layer
    print "Make XY event layer"
    arcpy.MakeXYEventLayer_management(in_Table, x_coords, y_coords, out_Layer, spRef, z_coords)

    # Save to a layer file
    print "Saving layer file"
    arcpy.SaveToLayerFile_management(out_Layer, saved_Layer)

    #===============================================================================
    # Copy features to shapefile
    print "Copying features to shapefile"
    arcpy.CopyFeatures_management(saved_Layer, "ElevPts.shp")
    #===============================================================================


except:
    print arcpy.GetMessages()

BenSciance

Still no luck. Did it work when you ran it?

JakeSkinner

The problem was with the coordinate system. Replace:

spRef = r("GCS_WGS_1984")

With:

spRef = arcpy.SpatialReference("WGS 1984")

Also, a good suggestion is to add print statements throughout the code to pinpoint where the code is failing.

DanPatterson_Retired

try saving it to a comma-separated values file (.csv)

BenSciance

Sure.

Columns are as follows: serial date number, date, Lat, Lon, and Z value.

TEST.txt.zip

GetImages.txt.zip

JakeSkinner

Can you upload a sample of your .txt file you are working with?

BenSciance

Hi Jake, I changed the workspace to the temp folder, still no luck. I guarantee the code is missing something very simple, as is the usual case from my experience. I can give it a go in ArcMap, just will take longer than desired.

If you have any other suggestions I'll give it a shot.

Thanks

Ben

JakeSkinner

As a test, try setting your env.workspace to another directory, i.e. r"C:\temp".

BenSciance

Hi Jake Thank you for the input.

I included the full filepath for variables in the original code just to double check everything was in order, I know it makes the code messy and unorganized, just a temporary thing while working this issue out.

I made the changes you suggested, specifically, adding the "r" in front of the filepath for defining the env.workspace.

The script completed, yet no new layer files and shapefiles were added to the workspace.

Any other possible workarounds or ideas as to why this is occurring?

Thank you

Ben

JakeSkinner

Hi Ben,

When specifying absolute paths, be sure to specify an 'r' before hand. Ex:

env.workspace = r"C:\Users\mbs7038\Documents\ArcGIS\matDEM"

However, once you have the env.workspace set, you can read/write data from this workspace without specifying the entire path. Ex:

arcpy.CopyFeatures_management(saved_Layer, "ElevPts.shp")

Try the following:

import arcpy
from arcpy import env

# Set environment settings
env.workspace = r"C:\Users\mbs7038\Documents\ArcGIS\matDEM"

try:
    # Set the local variables
    in_Table = "LatLonElev.txt"
    x_coords = "Lon"
    y_coords = "Lat"
    z_coords = "Elev"
    out_Layer = "ElevPts_layer"
    saved_Layer = "ElevPts.lyr"

    # Set the spatial reference
    spRef = r("GCS_WGS_1984")

    # Make the XY event layer
    arcpy.MakeXYEventLayer_management(in_Table, x_coords, y_coords, out_Layer, spRef, z_coords)

    # Save to a layer file
    arcpy.SaveToLayerFile_management(out_Layer, saved_Layer)

    #===============================================================================
    # Copy features to shapefile

    arcpy.CopyFeatures_management(saved_Layer, "ElevPts.shp")
    #===============================================================================


except:
# If an error occurred print the message to the screen
    print arcpy.GetMessages()

Nimm an der Unterhaltung teil

Angemeldete Mitglieder können Beiträge verfassen, Updates folgen und mehr. Neu hier? Registriere ein kostenloses Konto.

Schnellzugriffe

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.