LAS 点统计栅格 (Data Management)

摘要

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

用法

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

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

• LAS 数据集图层可用于按类代码或返回值过滤 LAS 点。创建图层的方法是使用创建 LAS 数据集图层工具、通过在 ArcMap 或 ArcScene 中加载 LAS 数据集以及通过图层属性对话框指定所需类代码和返回值。

语法

LasPointStatsAsRaster_management (in_las_dataset, out_raster, {method}, {sampling_type}, {sampling_value})

参数	说明	数据类型
in_las_dataset	输入 LAS 数据集。	LAS Dataset Layer
out_raster	输出栅格的位置和名称。在向地理数据库或文件夹（如 Esri Grid）中存储栅格数据集时，不应向栅格数据集的名称添加文件扩展名。可提供文件扩展名，以在向文件夹存储时定义栅格的格式： .bil—Esri BIL .bip—Esri BIP .bsq—Esri BSQ .dat—ENVI DAT .img—ERDAS IMAGINE .png—PNG .tif—TIFF 如果栅格存储为 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
import exceptions, sys, traceback

# 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:
    # Get the traceback object
    tb = sys.exc_info()[2]
    tbinfo = traceback.format_tb(tb)[0]
    # Concatenate error information into message string
    pymsg = 'PYTHON ERRORS:\nTraceback info:\n{0}\nError Info:\n{1}'\
          .format(tbinfo, str(sys.exc_info()[1]))
    msgs = 'ArcPy ERRORS:\n {0}\n'.format(arcpy.GetMessages(2))
    # Return python error messages for script tool or Python Window
    arcpy.AddError(pymsg)
    arcpy.AddError(msgs)

环境

当前工作空间, 临时工作空间, 范围, 输出坐标系, 地理变换, 像元大小, 金字塔, 栅格统计, 输出配置关键字, 自动提交

许可信息

ArcGIS for Desktop Basic:需要 3D Analyst 或者Spatial Analyst

ArcGIS for Desktop Standard:是

ArcGIS for Desktop Advanced:是