ArcMap

Make LAS Dataset Layer

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.

Syntax

arcpy.management.MakeLasDatasetLayer(in_las_dataset, out_layer, {class_code}, {return_values}, {no_flag}, {synthetic}, {keypoint}, {withheld}, {surface_constraints})
ParameterExplanationData 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.

  • 0 —Never processed by a classification method
  • 1 —Processed by a classification method but could not be determined
  • 2 —Bare earth measurements
  • 3 —Vegetation whose height is considered to be low for the area
  • 4 —Vegetation whose height is considered to be intermediate for the area
  • 5 —Vegetation whose height is considered to be high for the area
  • 6 —Structure with roof and walls
  • 7 —Erroneous or undesirable data that is closer to the ground
  • 8 —Reserved for later use, but used for model key points in LAS 1.1 - 1.3
  • 9 —Water
  • 10 —Railway tracks used by trains
  • 11 —Road surfaces
  • 12 —Reserved for later use, but used for overlap points in LAS 1.1 - 1.3
  • 13 —Shielding around electrical wires
  • 14 —Power lines
  • 15 —A lattice tower used to support an overhead power line
  • 16 —A mechanical assembly that joins an electrical circuit
  • 17 —The surface of a bridge
  • 18 —Erroneous or undesirable data that is far from the ground
  • 19 - 63 —Reserved class codes for ASPRS designation.
  • 64 - 255 —User-definable class codes.
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.

  • Last Return —Last return
  • First of Many —First of many
  • Last of Many —Last of many
  • Single Return —Single return
  • 1 —1st Return
  • 2 —2nd Return
  • 3 —3rd Return
  • 4 —4th Return
  • 5 —5th Return
  • 6 —6th Return
  • 7 —7th Return
  • 8 —8th Return
  • 9 —9th Return
  • 10 —10th Return
  • 11 —11th Return
  • 12 —12th Return
  • 13 —13th Return
  • 14 —14th Return
  • 15 —15th Return
String
no_flag
(Optional)

Specifies whether data points that do not have any classification flags assigned should be enabled for display and analysis.

  • INCLUDE_UNFLAGGED —Unflagged points will be displayed. This is the default.
  • EXCLUDE_UNFLAGGED —Unflagged points will not be displayed.
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..

  • INCLUDE_SYNTHETIC —Synthetic points will be displayed. This is the default.
  • EXCLUDE_SYNTHETIC —Synthetic points will not be displayed.
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.

  • INCLUDE_KEYPOINT —Model key points will be displayed. This is the default.
  • EXCLUDE_KEYPOINT —Model key points will not be displayed.
Boolean
withheld
(Optional)

Specifies whether data points flagged as withheld, which typically represent unwanted noise measurements, should be enabled for display and analysis.

  • EXCLUDE_WITHHELD —Withheld points will not be displayed. This is the default.
  • INCLUDE_WITHHELD —Withheld points will be displayed.
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

Code sample

MakeLasDatasetLayer example 1 (Python window)

The following sample demonstrates the use of this tool in the Python window.

import arcpy
from arcpy import env

arcpy.CheckOutExtension('3D')
env.workspace = 'C:/data'
arcpy.MakeLasDatasetLayer_management('Baltimore.lasd', 'Baltimore Layer', 
                                     2, 'LAST', 'INCLUDE_UNFLAGGED', 
                                     'EXCLUDE_SYNTHETIC', 'INCLUDE_KEYPOINT',
                                     'EXCLUDE_WITHHELD')
MakeLasDatasetLayer example 2 (stand-alone script)

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)
    if arcpy.ProductInfo == 'ArcView':
        arcpy.CheckOutExtension('3D')
    # 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)
    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

Related topics