摘要
创建栅格，栅格的像元值反映的是 LAS 数据集所引用 LAS 文件的测量值的相关统计信息。
用法
-
LAS 数据集图层还可用于控制表面约束要素的强化，而该要素可由 LAS 数据集引用。
-
LAS 格式支持根据美国摄影测量与遥感协会 (ASPRS) 制定的规范对每个点进行分类。ArcGIS 平台应用为 LAS 文件版本 1.4 指定的分类方案：
类别值 分类类型 0
从不分类
1
未分配
2
地面
3
低植被
4
中等植被
5
高植被
6
建筑物
7
低噪音
8
模型关键点/保留
9
水域
10
铁路
11
路面
12
重叠/保留
13
钢丝护网
14
导线
15
输电塔
16
线 - 连接器
17
桥板
18
高噪音
19 – 63
为 ASPRS 定义保留（LAS 1.1 至 1.3 最高支持类代码 31）
32 – 255
可由用户定义（仅 LAS 1.0 和 1.4 的特定版本支持）
可考虑对方法参数应用 PREDOMINANT_LAST_RETURN 选项，以识别回波值较高（可表示存在植被）的位置。
可考虑对方法参数应用 Z_RANGE 选项，以确定可能包含异常值的位置。
语法
LasPointStatsAsRaster_management (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 |
代码实例
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)
try:
arcpy.CheckOutExtension('3D')
# Execute CreateLasDataset
arcpy.management.CreateLasDataset(inLas, lasD, 'RECURSION', '', sr)
# Execute LasPointStatsAsRaster
arcpy.management.LasPointStatsAsRaster(lasD, statsRaster,
"PREDOMINANT_RETURNS_PER_PULSE",
"CELLSIZE", 15)
arcpy.CheckInExtension('3D')
except arcpy.ExecuteError:
print(arcpy.GetMessages())
except Exception as err:
print(err.args[0])