标注 | 说明 | 数据类型 |
输入 LAS 数据集 | 待处理的 LAS 数据集。 | LAS Dataset Layer |
输出栅格 | 输出栅格的位置和名称。 将栅格数据集存储到地理数据库或文件夹(例如 Esri GRID)时,请勿向栅格数据集的名称添加文件扩展名。 在将栅格存储到文件夹中时,可提供文件扩展名以定义栅格的格式,例如 .tif(生成 GeoTIFF)或 .img(生成 ERDAS IMAGINE 格式文件)。 如果栅格存储为 TIFF 文件或存储在地理数据库中,可使用地理处理环境设置指定其栅格压缩类型和质量。 | Raster Dataset |
方法
(可选) | 对输出栅格每个像元内的 LAS 点所采集的统计数据的类型。
| String |
采样类型 (可选) | 指定将用于解译采样值参数值以定义输出栅格分辨率的方法。
| String |
采样值 (可选) | 与采样类型结合使用以定义输出栅格分辨率的值。 | Double |
摘要
创建栅格,栅格的像元值反映的是 LAS 数据集所引用 LAS 文件的测量值的相关统计信息。
使用情况
-
LAS 数据集图层还可用于控制表面约束要素的强化,而该要素可由 LAS 数据集引用。 当将 LAS 数据集作为三角化网格面显示或处理时,将强制执行约束。
-
ArcGIS 使用美国摄影测量及遥感协会 (ASPRS) 定义的 LAS 分类方案。 了解有关激光雷达点分类的详细信息
可考虑对方法参数应用 PREDOMINANT_LAST_RETURN 选项,以识别回波值较高(可表示存在植被)的位置。
可考虑对方法参数应用 Z_RANGE 选项,以确定可能包含异常值的位置。
参数
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 点所采集的统计数据的类型。
| String |
sampling_type (可选) | 指定将用于解译采样值参数值以定义输出栅格分辨率的方法。
| String |
sampling_value (可选) | 与采样类型结合使用以定义输出栅格分辨率的值。 | Double |
代码示例
下面的示例演示了如何在 Python 窗口中使用此工具。
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.LasPointStatsAsRaster_3d("test.lasd", "lidar_intensity.img",
"INTENSITY_RANGE", "CELLSIZE", 15)
下面的示例演示了如何在独立 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: 是