Hello,
I am currently attempting to create a Python script tool that takes a user-specified GPX file (GPS data), and converts this data into ArcGIS features, which need to be in a specific coordinate system. I understand that there is a 'GPX to Features' Conversion tool within ArcToolbox, but the added funbtionality of the tool will be that the user will also be able to specify if they would like the data (which is always an array of points if collected as lines or polygons) converted into line features or polygon features. I work in county government, and my intentions are to share this tool with users in different departments who regularly collect GPX data and would like the above functionality within the framework of a single script tool. I am currently having some problems with the code I have written. The script runs correctly until reaching the point in which a spatial reference object must be created with the correct projection, at which point i get the error: RuntimeError: ERROR 999999: Error executing function. Here is the code that I have so far:
import arcpy
import os
import sys
import traceback
# Set variables
GPX = arcpy.GetParameterAsText(0) # User-specified GPX file
outlocGDB = arcpy.GetParameterAsText(1) # Output workspce to store GDB
fcnm = arcpy.GetParameterAsText(2) # Name/location of the feature class that is to be created
##que = ##arcpy.GetParameterAsText(3) # "Would you like to display you input GPX file as a Line Feature Class in ArcGIS?"
##que1 = ##arcpy.GetParameterAsText(4) # "Would you like to display your input GPX file as a Polygon Feature Class in ArcGIS?"
# Create new File Geodatabase to store the output data
GDB = arcpy.CreateFileGDB_management(outlocGDB, "GPXFeatures")
outgdb = GDB.getOutput(0) # Derive output of creation of File Geodatabase above
outfile = os.path.join(outgdb, fcnm) # output feature class for GPX to Features conversion
# Convert GPX files into features
arcpy.GPXtoFeatures_conversion(GPX, outfile)
# Set variables in preparation for projection operation
outprj = os.path.join(outgdb, str(fcnm) + "_prj") # Create output name for projected version of above Feature Class
outSR = arcpy.SpatialReference("NAD_1983_StatePlane_Virginia_North_FIPS_4501_Feet")
# PROJECT DATA
arcpy.Project_management(outfile, outprj, outSR, "WGS_1984_(ITRF00)_To_NAD_1983")
Any sort of help would be greatly appreciated. Thanks so much.
-Jacob Snyder
