用作栅格的 LAS 点统计数据 (数据管理)

摘要

创建栅格,栅格的像元值反映的是 LAS 数据集所引用 LAS 文件的测量值的相关统计信息。

使用情况

  • LAS 数据集图层还可用于控制表面约束要素的强化,而该要素可由 LAS 数据集引用。 当将 LAS 数据集作为三角化网格面显示或处理时,将强制执行约束。

  • ArcGIS 使用美国摄影测量及遥感协会 (ASPRS) 定义的 LAS 分类方案。 了解有关激光雷达点分类的详细信息

  • 可考虑对方法参数应用 PREDOMINANT_LAST_RETURN 选项,以识别回波值较高(可表示存在植被)的位置。

  • 可考虑对方法参数应用 Z_RANGE 选项,以确定可能包含异常值的位置。

参数

标注说明数据类型
输入 LAS 数据集

待处理的 LAS 数据集。

LAS Dataset Layer
输出栅格

输出栅格的位置和名称。 将栅格数据集存储到地理数据库或文件夹(例如 Esri GRID)时,请勿向栅格数据集的名称添加文件扩展名。 在将栅格存储到文件夹中时,可提供文件扩展名以定义栅格的格式,例如 .tif(生成 GeoTIFF)或 .img(生成 ERDAS IMAGINE 格式文件)。

如果栅格存储为 TIFF 文件或存储在地理数据库中,可使用地理处理环境设置指定其栅格压缩类型和质量。

Raster Dataset
方法
(可选)

对输出栅格每个像元内的 LAS 点所采集的统计数据的类型。

  • 脉冲计数最后回波点的数量。
  • 点计数所有回波中的点的数量。
  • 最常见的最后回波最常见的最后回波值。
  • 最常见类代码最常见的类代码。
  • 强度值范围强度值范围。
  • 高程值范围高程值范围。
String
采样类型
(可选)

指定将用于解译采样值参数值以定义输出栅格分辨率的方法。

  • 观测将使用分割 LAS 数据集范围的最长边的像元数。
  • 像元大小将使用输出栅格的像元大小。 这是默认设置。
String
采样值
(可选)

采样类型结合使用以定义输出栅格分辨率的值。

Double

arcpy.management.LasPointStatsAsRaster(in_las_dataset, out_raster, {method}, {sampling_type}, {sampling_value})
名称说明数据类型
in_las_dataset

待处理的 LAS 数据集。

LAS Dataset Layer
out_raster

输出栅格的位置和名称。 将栅格数据集存储到地理数据库或文件夹(例如 Esri GRID)时,请勿向栅格数据集的名称添加文件扩展名。 在将栅格存储到文件夹中时,可提供文件扩展名以定义栅格的格式,例如 .tif(生成 GeoTIFF)或 .img(生成 ERDAS IMAGINE 格式文件)。

如果栅格存储为 TIFF 文件或存储在地理数据库中,可使用地理处理环境设置指定其栅格压缩类型和质量。

Raster Dataset
method
(可选)

对输出栅格每个像元内的 LAS 点所采集的统计数据的类型。

  • PULSE_COUNT最后回波点的数量。
  • POINT_COUNT所有回波中的点的数量。
  • PREDOMINANT_LAST_RETURN最常见的最后回波值。
  • PREDOMINANT_CLASS最常见的类代码。
  • INTENSITY_RANGE强度值范围。
  • Z_RANGE高程值范围。
String
sampling_type
(可选)

指定将用于解译采样值参数值以定义输出栅格分辨率的方法。

  • OBSERVATIONS将使用分割 LAS 数据集范围的最长边的像元数。
  • CELLSIZE将使用输出栅格的像元大小。 这是默认设置。
String
sampling_value
(可选)

采样类型结合使用以定义输出栅格分辨率的值。

Double

代码示例

LAS 点统计栅格 (LasPointStatsAsRaster) 示例 1(Python 窗口)

下面的示例演示了如何在 Python 窗口中使用此工具。

import arcpy
from arcpy import env

env.workspace = "C:/data"
arcpy.LasPointStatsAsRaster_3d("test.lasd", "lidar_intensity.img", 
                             "INTENSITY_RANGE", "CELLSIZE", 15)
LAS 点统计栅格 (LasPointStatsAsRaster) 示例 2(独立脚本)

下面的示例演示了如何在独立 Python 脚本中使用此工具。

'''**********************************************************************
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)

许可信息

  • Basic: 需要 3D Analyst 或 Spatial Analyst
  • Standard: 是
  • Advanced: 是

相关主题