#Start Space import sys import os import arcpy from arcpy import env from arcpy.sa import * arcpy.env.workspace = "C:/gis/indian" LASdata = "C:/gis/indian" #Check out extensions arcpy.CheckOutExtension("3D") arcpy.CheckOutExtension("Spatial") arcpy.CheckOutExtension("GeoStats") #Begin processing files = os.listdir(LASdata) for file in files: if file.endswith(".las"): baseNom = "C:/gis/indian/" + file[:-4] #LAS to Multipoint outname = baseNom+".shp" try: arcpy.LASToMultipoint_3d(file, outname, 1, "5", "ANY_RETURNS") print "I've finished Shapefile with " + file except: print "I had LAS to Multipoint issues with file: " + file try: coor = arcpy.SpatialReference('WGS 1984 UTM Zone 16N', 'North American Vertical Datum 1988') arcpy.DefineProjection_management(outname, coor) except: print "I had issues projecting the data" print arcpy.GetMessages()
Solved! Go to Solution.
As a first thought, why not dump the shapefile output and instead create the output of the LAS conversion in a goedatabase inside a feature dataset where the vertical datum has been defined.
Cheers,
N