Skip To Content

Classify LAS Ground

Summary

Classifies ground points in aerial lidar data.

Illustration

Classify LAS Ground

Usage

  • This tool requires the input LAS dataset to have a projected coordinate system. Data stored in a geographic coordinate system can be reprojected by using the Extract LAS tool with a projected coordinate system specified in the output coordinate system environment setting.

  • Only the last return of LAS points with class code values of 0, 1, or 2 will be considered for reclassification as ground. If your LAS files use different class code values to represent unclassified or ground measurements, consider using the Change LAS Class Codes tool to reassign them accordingly. The classification process will also ignore points that are assigned with either the overlap or noise classification flag.

  • If the ground classification is being done solely to create a ground raster surface of a specific resolution, consider using the DEM Resolution parameter with the desired raster resolution to create the requisite classification result for interpolating the surface.

  • When classifying LAS returns over a terrain with divergent slope characteristics, such as relatively flat areas alongside locations with steep slope profiles, consider running the ground classifier once with the standard method and a second time with the aggressive detection method and the reuse ground option enabled. Apply a processing extent or specify a polygon boundary to limit this operation to just the region it is needed.

  • Locations with bridges and freeway on-ramps should be reviewed as they may be misclassified as ground.

Syntax

ClassifyLasGround_3d (in_las_dataset, method, {reuse_ground}, {dem_resolution}, {compute_stats}, {extent}, boundary, {process_entire_files})
ParameterExplanationData Type
in_las_dataset

The LAS dataset to process. Only the last return of LAS points with class code values of 0, 1, and 2 will be evaluated.

LAS Dataset Layer
method

The method used for detecting ground points.

  • STANDARDThis method has a tolerance for slope variation that allows it to capture gradual undulations in the ground's topography that would typically be missed by the conservative option but not capture the type of sharp reliefs that would be captured by the aggressive option. This method is the default.
  • CONSERVATIVE When compared to other options, this method employs a tighter restriction on the variation of the ground's slope, which allows it to differentiate the ground from low-lying vegetation such as grass and shrubbery. It is best suited for topography with minimal curvature.
  • AGGRESSIVEThis method detects ground areas with sharper reliefs, such as ridges and hill tops, that may be ignored by the STANDARD method. This method is best used in a second pass of this tool with the reuse_ground parameter set to REUSE_GROUND. Avoid using this method in urban areas or flat, rural areas, as it may result in the misclassification of higher objects—such as utility towers, vegetation, and portions of buildings—as ground.
String
reuse_ground
(Optional)

Specifies if existing ground points should be reclassified or reused.

  • RECLASSIFY_GROUND Existing ground points will be reclassified. Points that are not found to be a part of the ground will be re-assigned a class code value of 1, which represents unclassified points. This is the default.
  • REUSE_GROUND Existing ground points will be accepted without scrutiny and contribute to the determination of unclassified points.
Boolean
dem_resolution
(Optional)

Specifies a distance that will result in only a subset of points being evaluated for classification as ground, thereby making the process faster. Consider using this option when a faster method for generating a DEM surface is needed. The minimum distance is 0.3 meters, but the specified distance must be at least 1.5 times the average point spacing of the lidar data for this process to take effect.

Linear Unit
compute_stats
(Optional)

Specifies whether statistics should be computed for the LAS files referenced by the LAS dataset. Computing statistics provides a spatial index for each LAS file, which improves analysis and display performance. Statistics also enhance the filtering and symbology experience by limiting the display of LAS attributes, like classification codes and return information, to values that are present in the LAS file.

  • COMPUTE_STATSStatistics will be computed.
  • NO_COMPUTE_STATSStatistics will not be computed. This is the default.
Boolean
extent
(Optional)

Specifies the extent of the data that will be evaluated by this tool.

Extent
boundary

A polygon feature that defines the area of interest to be processed by this tool.

Feature Layer
process_entire_files
(Optional)

Specifies how the processing extent is applied.

  • PROCESS_EXTENTOnly LAS points that intersect the area of interest will be processed. This is the default.
  • PROCESS_ENTIRE_FILESIf any portion of a LAS file intersects the area of interest, all the points in that LAS file, including those outside the area of interest, will be processed.
Boolean

Derived Output

NameExplanationData Type
out_las_dataset

The LAS dataset that was modified.

LAS Dataset Layer

Code sample

ClassifyLasGround example 1 (Python window)

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

arcpy.env.workspace = 'C:/data'
arcpy.ClassifyLasGround_3d('metro.lasd', 'CONSERVATIVE', 
                           boundary='study_area.shp', 
                           process_entire_files='PROCESS_ENTIRE_FILES')
ClassifyLasGround example 2 (stand-alone script)

The following sample demonstrates the use of this tool in a stand-alone Python script.

'''****************************************************************************
Name:        Classify Ground & Vegetation in Forest Environment
Description: Classify points representing vegetation with LAS class code values
             of 3, 4, and 5. The code is designed for use as a script tool.
****************************************************************************'''
# Import system modules
import arcpy

# Set Local Variables
inLas = arcpy.GetParameterAsText(0)
recursion = arcpy.GetParameterAsText(1)
lasd = arcpy.GetParameterAsText(2)

try:
    arcpy.CheckOutExtension('3D')
    # Execute CreateLasDataset
    arcpy.management.CreateLasDataset(inLas, lasd, folder_recursion=recursion)
    # Make an initial pass of ground classifier
    arcpy.ddd.ClassifyLasGround(lasd, method="Conservative")
    # Make a secondary pass to capture ridges
    arcpy.ddd.ClassifyLasGround(lasd, method="Aggressive", 
                                reuse_ground="REUSE_GROUND")
    # Classify vegetation
    arcpy.ddd.ClassifyLasByHeight(lasd, ground_source='GROUND', 
                                  height_classification=[[3, 5], 
                                                         [4, 17], 
                                                         [5, 120]], 
                                  noise='HIGH_NOISE', compute_stats="COMPUTE_STATS")
    arcpy.CheckInExtension('3D')

except arcpy.ExecuteError:
    print(arcpy.GetMessages())

Licensing information

  • Basic: Requires 3D Analyst
  • Standard: Requires 3D Analyst
  • Advanced: Requires 3D Analyst

Related topics