'''****************************************************************************
Name: Define Data Boundary of LAS File
Description: This script demonstrates how to delineate data boundaries of
LAS files with irregularly clustered points. It is intended for
use as a script tool with one input LAS file.
****************************************************************************'''
import arcpy
inLas = arcpy.GetParameterAsText(0)
ptSpacing = arcpy.GetParameterAsText(1)
classCode = arcpy.GetParameterAsText(2)
returnValue = arcpy.GetParameterAsText(3)
outTin = arcpy.GetParameterAsText(4)
outBoundary = arcpy.GetParameterAsText(5)
try:
lasMP = arcpy.CreateUniqueName('lasMultipoint', 'in_memory')
arcpy.ddd.LASToMultipoint(inLas, LasMP, ptSpacing, class_code,
"ANY_RETURNS", "", sr, inFormat, zfactor)
arcpy.ddd.CreateTin(outTin, sr, "{0} Shape.Z masspoints"\
.format(lasMP), "Delaunay")
arcpy.ddd.CopyTin(outTin, "{0}_copy".format(outTin))
maxEdge = ptSpacing * 4
arcpy.ddd.DelineateTinDataArea(outTin, maxEdge, "PERIMETER_ONLY")
arcpy.ddd.TinDomain(outTin, outBoundary, "POLYGON")
except arcpy.ExecuteError:
print(arcpy.GetMessages())
except Exception as err:
print(err)