需要 3D Analyst 许可。需要 Spatial Analyst 许可。
摘要
计算或更新 LAS 数据集的统计数据并生成可选的统计数据报表。
用法
通过使用分组方法(可逐步评估文件的小区域来确定估计),计算每个 LAS 文件的平均点间距的统计数值。
统计数据对 LAS 数据集图层启用过滤选项来自动显示可用的类代码并返回在引用的 LAS 文件中找到的值。在 ArcMap 和 ArcScene 中,过滤选项可通过“图层属性”对话框指定。
以下有关 LAS 文件的信息将存储在 LAS 数据集中,并且可通过 ArcCatalog 中的“LAS 数据集属性”对话框进行查看:
- 高程值和强度值的范围
- 每个类代码和返回值中数据点的数量
- RGB 值的存在
- LAS 版本和点格式 ID
- 投影 ID 和文件源 ID
- 用来生成激光雷达文件的软件的名称
可选的统计数据报告文件概述 LAS 数据集中所有文件的 LAS 属性,或者概述可导入任意数量应用程序的文本格式的各个独立 LAS 文件的 LAS 属性。
语法
LasDatasetStatistics_management (in_las_dataset, calculation_type, {out_file}, {summary_level}, {delimiter}, {decimal_separator})
参数 | 说明 | 数据类型 |
in_las_dataset | 待处理的 LAS 数据集。 | LAS Dataset Layer |
calculation_type | 指定是否针对所有激光雷达文件或仅针对那些没有统计数据的文件计算统计数据:
| Boolean |
out_file (可选) | 包含 LAS 数据集统计数据汇总的输出文本文件。 | Text File |
summary_level (可选) | 指定包含在报表中的汇总类型。
| String |
delimiter (可选) | 用来指示文本文件表的列中输入条目间隔的分隔符。
| String |
decimal_separator (可选) | 文本文件中用于区分数字的整数部分与其小数部分的小数分隔符。
| String |
代码实例
LAS 数据集统计数据 (LasDatasetStatistics) 示例 1(Python 窗口)
下面的示例演示了如何在 Python 窗口中使用此工具。
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.LASDatasetStatistics_3d("test.lasd", "NO_FORCE", "LAS_FILE",
"DECIMAL_POINT", "SPACE", "LAS_summary.txt")
LAS 数据集统计数据 (LasDatasetStatistics) 示例 2(独立脚本)
下面的示例演示了如何在独立 Python 脚本中使用此工具。
'''*********************************************************************
Name: Modify Files in LAS Dataset& Calculate Stats for LASD
Description: Adds files & surface constraints to a LAS dataset, then
calculates statistics and generates report.
*********************************************************************'''
# Import system modules
import arcpy
try:
# Script variables
arcpy.env.workspace = 'C:/data'
lasd = 'sample.lasd'
oldLas = ['2006', '2007/file2.las']
newLas = ['2007_updates_1', '2007_updates_2']
oldSurfaceConstraints = ['boundary.shp', 'streams.shp']
newSurfaceConstraints = [['sample.gdb/boundary', '<None>',
'Soft_Clip']
['sample.gdb/streams', 'Shape.Z',
'Hard_Line']]
arcpy.management.RemoveFilesFromLasDataset(lasd, oldLas,
oldSurfaceConstraints)
arcpy.management.AddFilesToLasDataset(lasd, newLas, 'RECURSION',
newSurfaceConstraints)
arcpy.management.LasDatasetStatistics(lasd, "UPDATED_FILES",
"lasd_stats.txt",
"LAS_FILE", "DECIMAL_POINT",
"SPACE", "LAS_summary.txt")
except arcpy.ExecuteError:
print(arcpy.GetMessages())
except Exception as err:
print(err.args[0])