LAS Point Statistics As Raster (Data Management Tools)
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.
In addition to filtering points by certain attributes, you can also filter LAS points by elevation using the Elevation Range control in the Filters group of the ribbon for a LAS dataset layer. This tool will honor the elevation range set on this control.
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.
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 Count—The number of last return points will be collected.
Point Count—The number of points from all returns will be collected.
Most Frequent Last Return—The most frequent last return value will be collected.
Most Frequent Class Code—The most frequent class code will be collected.
Range of Intensity Values—The range of intensity values will be collected.
Range of Elevation Values—The 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.
Observations—The 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 Size—The 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.
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_COUNT—The number of last return points will be collected.
POINT_COUNT—The number of points from all returns will be collected.
PREDOMINANT_LAST_RETURN—The most frequent last return value will be collected.
PREDOMINANT_CLASS—The most frequent class code will be collected.
INTENSITY_RANGE—The range of intensity values will be collected.
Z_RANGE—The 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.
OBSERVATIONS—The 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.
CELLSIZE—The 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:
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)
# Run CreateLasDataset
arcpy.management.CreateLasDataset(inLas, lasD, 'RECURSION', '', sr)
# Run LasPointStatsAsRaster
arcpy.management.LasPointStatsAsRaster(lasD, statsRaster,
"PREDOMINANT_RETURNS_PER_PULSE",
"CELLSIZE", 15)