LAS 高度度量 (3D Analyst 工具)

汇总

计算有关 LAS 数据中捕获的植被点高程测量值分布的统计数据。

用法

该工具会生成地面上方点的统计衍生数据，以更好地了解森林冠层。将分析分配有类别代码值 0 或 1（代表未分类点）以及 3、4 和 5（代表植被点）的 LAS 点。
在执行该工具之前，必须对正在处理的 LAS 数据的地面点进行分类。可以使用 分类 LAS 地面点 工具满足此要求。
该工具会为指定的每个统计选项和高度百分比数输出一个栅格。每个栅格的名称以指定的基本名称和高度指标或空间统计数据开头。对于平均值栅格，文件名将追加 "mean" 一词。对于标准差，文件名将追加 "stdev"。对于中位数绝对差，文件名将追加 "mad"。
使用最小高度参数排除非常低的点，尤其是当地面点未被彻底分类时，因为它们可能会导致该工具的结果失真。
使用最小点数参数排除包含的离地点数过少而无法生成有意义结果的区域。

标注	说明	数据类型
输入 LAS 数据集	将要处理的 LAS 数据集。	LAS Dataset Layer
输出位置	输出栅格数据集所在的文件夹或地理数据库。如果输出位置为文件夹，则生成的栅格数据集将采用 TIFF 格式。	Workspace
输出基本名称（可选）	输出栅格数据集的基本名称。	String
统计数据选项（可选）	指定为输出栅格中每个像元区域内高于地面的未分类点和植被点计算的统计数据。 Mean—The average height of the LAS points. Kurtosis—The sharpness of the change in the height of the LAS points. Skewness—The direction of deviation from the nominal height of the LAS points, which indicates the level and direction of asymmetry. Standard Deviation—The variation of the height of the points. Median Absolute Deviation—The median value of the deviation from the median height.	String
高度百分比（可选）	像元中指定百分比的点所低于的高度。例如，此值为 95 意味着生成的像元值表明 95% 高于地面的点出现的高度。	Long
最小高度	将评估的点的地面上最小高度。	Linear Unit
最小点数	计算高度度量时给定像元中必须存在的最小点数。点数小于指定最小值的像元在输出中将没有数据。	Long
像元大小	输出栅格数据集的像元大小。	Linear Unit
栅格格式（可选）	指定当输出位置为文件夹时要创建的栅格格式。 GeoTiff—Output will be created in the GeoTIFF format. This is the default. Erdas Imagine (IMG)—Output will be created in the ERDAS IMAGINE format. Esri Grid—Output will be created in the Esri Grid format.	String

派生输出

标注	说明	数据类型
输出位置	输出栅格数据集的位置。	Workspace
输出栅格	输出栅格数据集。	Raster Dataset

arcpy.3d.LasHeightMetrics(in_las_dataset, out_location, {base_name}, {statistics}, {height_percentiles}, min_height, min_points, cell_size, {raster_format})

名称	说明	数据类型
in_las_dataset	将要处理的 LAS 数据集。	LAS Dataset Layer
out_location	输出栅格数据集所在的文件夹或地理数据库。如果输出位置为文件夹，则生成的栅格数据集将采用 TIFF 格式。	Workspace
base_name （可选）	输出栅格数据集的基本名称。	String
statistics [statistics,...] （可选）	指定为输出栅格中每个像元区域内高于地面的未分类点和植被点计算的统计数据。 `MEAN`—The average height of the LAS points. `KURTOSIS`—The sharpness of the change in the height of the LAS points. `SKEWNESS`—The direction of deviation from the nominal height of the LAS points, which indicates the level and direction of asymmetry. `STANDARD_DEVIATION`—The variation of the height of the points. `MEDIAN_ABSOLUTE_DEVIATION`—The median value of the deviation from the median height.	String
height_percentiles [height_percentiles,...] （可选）	像元中指定百分比的点所低于的高度。例如，此值为 95 意味着生成的像元值表明 95% 高于地面的点出现的高度。	Long
min_height	将评估的点的地面上最小高度。	Linear Unit
min_points	计算高度度量时给定像元中必须存在的最小点数。点数小于指定最小值的像元在输出中将没有数据。	Long
cell_size	输出栅格数据集的像元大小。	Linear Unit
raster_format （可选）	指定当输出位置为文件夹时要创建的栅格格式。 `TIFF`—Output will be created in the GeoTIFF format. This is the default. `IMG`—Output will be created in the ERDAS IMAGINE format. `ESRI_GRID`—Output will be created in the Esri Grid format.	String

派生输出

名称	说明	数据类型
derived_out_location	输出栅格数据集的位置。	Workspace
output_rasters	输出栅格数据集。	Raster Dataset

代码示例

LasHeightMetrics 示例 1（Python 窗口）

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

arcpy.env.workspace = 'C:/data'

arcpy.ddd.ClassifyLasBuilding('Highland.lasd', minHeight='9 feet',
                              minArea='30 Square Feet', compute_stats=True)

LasHeightMetrics 示例 2（独立脚本）

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

'''****************************************************************************
       Name: Tile & Classify LAS Files
Description: Creates & classifies tiled LAS files.
****************************************************************************'''
# Import system modules
import arcpy
import tempfile
import math

in_las = arcpy.GetParameterAsText(1) # The LAS files that need to be tiled
out_folder = arcpy.GetParameterAsText(2) # folder for LAS files
basename = arcpy.GetParameterAsText(3) # basename for output files
out_lasd = arcpy.GetParameterAsText(4) # output LAS dataset

try:
    # Create temp LAS dataset to reference LAS files that will be tiled
    temp_lasd = arcpy.CreateUniqueName('temp.lasd', tempfile.gettempdir())
    arcpy.management.CreateLasDataset(in_las, temp_lasd)
    arcpy.ddd.TileLas(temp_lasd, out_folder, basename, out_lasd, las_version=1.4,
                      point_format=7, file_size=300)
    arcpy.management.Delete(temp_lasd)
    arcpy.ddd.ClassifyLasGround(out_lasd, method='AGGRESSIVE')
    arcpy.ddd.ClassifyLasBuilding(out_lasd, min_height='3 Meters', min_area='4 Meters')
    arcpy.ddd.ClassifyLasByHeight(out_lasd, height_classification=[(3, 6), (4,20), (5,70)],
                                  noise='All Noise', compute_stats='COMPUTE_STATS')

except arcpy.ExecuteError:
    print(arcpy.GetMessages())

环境

当前工作空间, 范围, 输出坐标系, 地理变换

许可信息

基本: 需要 3D Analyst
标准: 需要 3D Analyst
高级: 需要 3D Analyst

LAS 高度度量 (3D Analyst 工具)

汇总

用法

参数

派生输出

环境

许可信息

相关主题