LAS 数据集统计数据 (Data Management)

摘要

计算或更新 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	指定是否针对所有激光雷达文件或仅针对那些没有统计数据的文件计算统计数据： SKIP_EXISTING_STATS —将跳过包含最新统计数据的 LAS 文件，仅为新添加的 LAS 文件或自初始计算后更新的 LAS 文件计算统计数据。这是默认设置。 OVERWRITE_EXISTING_STATS —为所有 LAS 文件（包括具有最新统计数据的 LAS 文件）计算统计数据。这适合于在 ArcGIS 未检测到的外部应用程序中修改了 LAS 文件的情况。	Boolean
out_file (可选)	包含 LAS 数据集统计数据汇总的输出文本文件。	Text File
summary_level (可选)	指定包含在报表中的汇总类型。 DATASET —报表将汇总整个 LAS 数据集的统计数据。这是默认设置。 LAS_FILES —报表将汇总 LAS 数据集引用的 LAS 文件的统计数据。	String
delimiter (可选)	文本文件中使用的字段定界符。 空格 —空格将用于分隔字段值。这是默认设置。 逗号 —逗号将用于分隔字段值。如果小数分隔符也是逗号，则此选项不适用。	String
decimal_separator (可选)	文本文件中用于区分数字的整数部分与其小数部分的小数分隔符。 DECIMAL_POINT —点用作小数字符。这是默认设置。 DECIMAL_COMMA —逗号用作小数字符。	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
from arcpy import env
import exceptions, sys, traceback

try:
    # Script variables
    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:
    # 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:是