LAS Point Statistics As Raster (Data Management)

Summary

Creates a raster whose cell values reflect statistical information about LAS points.

Usage

  • You can filter the points that are processed by this tool using any combination of classification codes, classification flags, and return values through the LAS dataset layer's point filters. The filters can be defined on the Layer Properties dialog box or using the Make LAS Dataset Layer tool. LAS point filters will reflect in the results obtained for all statistical outputs except the Pulse Count method, which evaluates all last return points.

  • The Point Count option of the Method parameter evaluates the point density of the LAS dataset. This information is helpful when detecting irregular hot spots of high density points, because this type of data can produce degraded performance when processing the LAS files using other tools. If high density hot spots are encountered, the Thin LAS tool can be used to produce a new point cloud with more consistent spatial distribution.

  • The Most Frequent Class Code option of the Method parameter can be used to identify the number of unique objects defined by a specific class code, particularly if those objects are spatially separated in the horizontal direction. To do this, first apply a point filter to isolate the specific class code. Then create the statistics raster using a cell size that is large enough to ensure that the connectivity of cells representing points belong to the same object and that gaps between points from different objects are preserved. For example, to determine the cell size appropriate for connecting points belonging to a street light, consider the overall spacing of the points and the distance between nearby street lights. The resulting raster can be used to assess the footprint of each object. The number of unique objects can be obtained by either converting the raster to a polygon feature using the Raster To Polygon tool or assigning a unique value to each cluster of cells using the Region Group tool.

Parameters

LabelExplanationData Type
Input LAS Dataset

The LAS dataset that will be processed.

LAS Dataset Layer
Output Raster

The location and name of the output raster. When storing a raster dataset in a geodatabase or in a folder such as an Esri Grid, do not add a file extension to the name of the raster dataset. A file extension can be provided to define the raster's format when storing it in a folder, such as .tif to generate a GeoTIFF or .img to generate an ERDAS IMAGINE format file.

If the raster is stored as a .tif file or in a geodatabase, the raster compression type and quality can be specified using geoprocessing environment settings.

Raster Dataset
Method
(Optional)

Specifies the type of statistics that will be collected about the LAS points in each cell of the output raster.

  • Pulse CountThe number of last return points will be collected.
  • Point CountThe number of points from all returns will be collected.
  • Most Frequent Last ReturnThe most frequent last return value will be collected.
  • Most Frequent Class CodeThe most frequent class code will be collected.
  • Range of Intensity ValuesThe range of intensity values will be collected.
  • Range of Elevation ValuesThe range of elevation values will be collected.
String
Sampling Type
(Optional)

Specifies how the Sampling Value parameter will be interpreted to define the output raster's cell size.

  • ObservationsThe Sampling Value will define the number of columns or rows in the output raster based on whichever is longest. The cell size will be derived by dividing the longest side of the output's extent with the input in the Sampling Value parameter. If an observation value of 3000 is used on a dataset whose longest side is 23.67 kilometers, the output raster's resolution will be 7.89 meters. This method offers a helpful way of creating an output with a predictable size that can be generated rapidly.
  • Cell SizeThe cell size will be directly defined by the Sampling Value parameter. This is the default.
String
Sampling Value
(Optional)

The value used in conjunction with the Sampling Type parameter to define the output raster's cell size.

Double

arcpy.management.LasPointStatsAsRaster(in_las_dataset, out_raster, {method}, {sampling_type}, {sampling_value})
NameExplanationData Type
in_las_dataset

The LAS dataset that will be processed.

LAS Dataset Layer
out_raster

The location and name of the output raster. When storing a raster dataset in a geodatabase or in a folder such as an Esri Grid, do not add a file extension to the name of the raster dataset. A file extension can be provided to define the raster's format when storing it in a folder, such as .tif to generate a GeoTIFF or .img to generate an ERDAS IMAGINE format file.

If the raster is stored as a .tif file or in a geodatabase, the raster compression type and quality can be specified using geoprocessing environment settings.

Raster Dataset
method
(Optional)

Specifies the type of statistics that will be collected about the LAS points in each cell of the output raster.

  • PULSE_COUNTThe number of last return points will be collected.
  • POINT_COUNTThe number of points from all returns will be collected.
  • PREDOMINANT_LAST_RETURNThe most frequent last return value will be collected.
  • PREDOMINANT_CLASSThe most frequent class code will be collected.
  • INTENSITY_RANGEThe range of intensity values will be collected.
  • Z_RANGEThe range of elevation values will be collected.
String
sampling_type
(Optional)

Specifies how the Sampling Value parameter will be interpreted to define the output raster's cell size.

  • OBSERVATIONSThe Sampling Value will define the number of columns or rows in the output raster based on whichever is longest. The cell size will be derived by dividing the longest side of the output's extent with the input in the Sampling Value parameter. If an observation value of 3000 is used on a dataset whose longest side is 23.67 kilometers, the output raster's resolution will be 7.89 meters. This method offers a helpful way of creating an output with a predictable size that can be generated rapidly.
  • CELLSIZEThe cell size will be directly defined by the Sampling Value parameter. This is the default.
String
sampling_value
(Optional)

The value used in conjunction with the Sampling Type parameter to define the output raster's cell size.

Double

Code sample

LasPointStatsAsRaster example 1 (Python window)

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

import arcpy
from arcpy import env

env.workspace = "C:/data"
arcpy.LasPointStatsAsRaster_3d("test.lasd", "lidar_intensity.img", 
                             "INTENSITY_RANGE", "CELLSIZE", 15)
LasPointStatsAsRaster example 2 (stand-alone script)

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

'''**********************************************************************
Name: LAS Point Statistics As Raster
Description: Identifies the most frequently occurring return value for
             each pulse in a given set of LAS files.
             Designed for use as a script tool.
**********************************************************************'''
# Import system modules
import arcpy

# Set Local Variables
lasD = arcpy.GetParameterAsText(0)
inLas = arcpy.GetParameterAsText(1) #input las files
sr = arcpy.GetParameter(2) #spatial reference of las dataset
statsRaster = arcpy.GetParameterAsText(3)

# Execute CreateLasDataset
arcpy.management.CreateLasDataset(inLas, lasD, 'RECURSION', '', sr)
# Execute LasPointStatsAsRaster
arcpy.management.LasPointStatsAsRaster(lasD, statsRaster,
                                       "PREDOMINANT_RETURNS_PER_PULSE",
                                       "CELLSIZE", 15)

Related topics