Summary
Creates a LAS dataset layer that can apply filters to LAS points and control the enforcement of surface constraint features.
Usage
A LAS dataset layer can be used to filter lidar data for subsequent analysis. For example, a digital elevation model (DEM) can be easily created by filtering for class code 2, which represents ground measurements, then using the resulting layer as an input for the LAS Dataset To Raster tool.
Note:
The resulting layer can be preserved as a layer file by using the Save To Layer File tool.
Syntax
arcpy.management.MakeLasDatasetLayer(in_las_dataset, out_layer, {class_code}, {return_values}, {no_flag}, {synthetic}, {keypoint}, {withheld}, {surface_constraints}, {overlap})
Parameter | Explanation | Data Type |
in_las_dataset | The LAS dataset to process. | LAS Dataset Layer |
out_layer | The name of the resulting LAS dataset layer. Any backslash or forward slash can be used to denote a group layer. | LAS Dataset Layer |
class_code [class_code,...] (Optional) |
Allows the filtering of LAS points by classification codes. The range of valid values will depend on the class codes supported by the version of LAS files referenced by the LAS dataset. All class codes will be selected by default.
| String |
return_values [return_values,...] (Optional) | The return values to be used for filtering LAS points. When no value is specified, all returns are used.
| String |
no_flag (Optional) | Specifies whether data points that do not have any classification flags assigned should be enabled for display and analysis.
| Boolean |
synthetic (Optional) | Specifies whether data points flagged as synthetic, or points that originated from a data source other than lidar, should be enabled for display and analysis..
| Boolean |
keypoint (Optional) | Specifies whether data points flagged as model key points, or significant measurements that should not be thinned away, should be enabled for display and analysis.
| Boolean |
withheld (Optional) | Specifies whether data points flagged as withheld, which typically represent unwanted noise measurements, should be enabled for display and analysis.
| Boolean |
surface_constraints [surface_constraints,...] (Optional) | The name of the surface constraint features that will be enabled in the layer. All constraints are enabled by default. | String |
overlap (Optional) | Specifies whether data points flagged as overlap should be enabled for display and analysis.
| Boolean |
Code sample
The following sample demonstrates the use of this tool in the Python window.
arcpy.env.workspace = 'C:/data'
arcpy.MakeLasDatasetLayer_management('Baltimore.lasd', 'Baltimore Layer',
class_code=[2, 6], return_values=['LAST', 'SINGLE'])
The following sample demonstrates the use of this tool in a stand-alone Python script.
'''*********************************************************************
Name: Export Elevation Raster from Ground LAS Measurements
Description: This script demonstrates how to export
ground measurements from LAS files to a raster using a
LAS dataset. This sample is designed to be used as a script
tool.
*********************************************************************'''
# Import system modules
import arcpy
try:
# Set Local Variables
inLas = arcpy.GetParameterAsText(0)
recursion = arcpy.GetParameterAsText(1)
surfCons = arcpy.GetParameterAsText(2)
classCode = arcpy.GetParameterAsText(3)
returnValue = arcpy.GetParameterAsText(4)
spatialRef = arcpy.GetParameterAsText(5)
lasD = arcpy.GetParameterAsText(6)
outRaster = arcpy.GetParameterAsText(7)
cellSize = arcpy.GetParameter(8)
zFactor = arcpy.GetParameter(9)
# Execute CreateLasDataset
arcpy.management.CreateLasDataset(inLas, lasD, recursion, surfCons, sr)
# Execute MakeLasDatasetLayer
lasLyr = arcpy.CreateUniqueName('Baltimore')
arcpy.management.MakeLasDatasetLayer(lasD, lasLyr, classCode, returnValue)
# Execute LasDatasetToRaster
arcpy.conversion.LasDatasetToRaster(lasLyr, outRaster, 'ELEVATION',
'TRIANGULATION LINEAR WINDOW_SIZE 10', 'FLOAT',
'CELLSIZE', cellSize, zFactor)
print(arcpy.GetMessages())
except arcpy.ExecuteError:
print(arcpy.GetMessages())
except Exception as err:
print(err.args[0])
finally:
arcpy.management.Delete(lasLyr)
Environments
Licensing information
- Basic: Yes
- Standard: Yes
- Advanced: Yes