摘要
创建引用一个或多个 LAS 文件和可选表面约束要素的 LAS 数据集。
插图
用法
-
LAS 数据集旨在配合 .las 或 .zlas 格式的激光雷达数据使用。支持的 LAS 文件版本为 1.0 - 1.4。Esri 的 EzLAS Optimizer 是一个独立的激光雷达应用程序,可用于生成 .zlas 文件或将它们转换回 .las 格式。
-
对各个 LAS 文件进行检查,以确定其内部结构是否与 LAS 规范相一致。如果 LAS 文件没能加载到 LAS 数据集中,表示可能出现了文件损坏或文件标头信息错误的情况。考虑使用 Esri 的 CheckLAS 工具来确定数据中是否存在一些常见的问题。
-
表面约束要素可用于强化 LAS 数据集中要素衍生的表示表面特征的高程值。
-
每个 LAS 文件通常在其标题(由 LAS 数据集读取)中包含空间参考信息。如果此信息丢失或未正确定义,则 LAS 文件将不会放置在正确的位置。如果已知正确的空间参考,则可考虑使用为 LAS 创建 PRJ 参数创建能够对激光雷达进行正确地理配准的 .prj 文件。.prj 文件与 LAS 文件的名称相同,位于同一文件夹,且包含 LAS 文件坐标系的熟知文本表示(类似于与 shapefile 相关联的 .prj 文件)。
-
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 的特定版本支持)
在“工具”对话框中,可将文件夹指定为输入,具体方法如下:在 Windows 资源管理器中选择所需文件夹,然后将其拖动到参数的输入框中。
语法
CreateLasDataset_management (input, out_las_dataset, {folder_recursion}, {in_surface_constraints}, {spatial_reference}, {compute_stats}, {relative_paths}, create_las_prj)
参数 | 说明 | 数据类型 |
input [input,...] | LAS 文件和包含将被 LAS 数据集引用的 LAS 文件的文件夹。该信息可通过以下形式提供:包含所有输入数据的字符串,或包含特定数据元素的字符串列表(如:"lidar1.las; lidar2.las; folder1; folder2" 或 ["lidar1.las", "lidar2.las", "folder1", "folder2"])。 有关指定输入列表的详细信息,请参阅了解工具语法。 | File; Folder |
out_las_dataset | 将创建的 LAS 数据集。 | LAS Dataset |
folder_recursion (可选) | 指定位于输入文件夹子目录中的激光雷达数据是否将添加到 LAS 数据集。
| Boolean |
in_surface_constraints [[in_feature_class, height_field, SF_type],...] (可选) | 生成三角化网格面时,LAS 数据集将引用的要素。每个要素都需要定义以下属性: in_feature_class - LAS 数据集将引用的要素。 height_field - 可以使用要素属性表中的任何数值字段来定义高度源。如果要素几何包含 Z 值,则可以通过指定 Shape.Z 进行选择。如果没有所需高度,则指定关键字 <无> 来创建 Z-less 要素,其高程由表面内插得到。 SF_type - 定义如何将要素几何加入到表面的三角网中的表面要素类型。具有硬或软标识的选项表示要素边是否表示坡度的明显中断或平缓变化。
| Value Table |
spatial_reference (可选) | LAS 数据集的空间参考。如果没有明确指定空间参考,则 LAS 数据集将使用第一个输入 LAS 文件的坐标系。如果输入文件不包含任何空间参考信息,并且没有设置输入坐标系,则 LAS 数据集的坐标系将列为未知。 | Coordinate System |
compute_stats (可选) | 指定是否计算 LAS 文件的统计数据以及是否为 LAS 数据集生成空间索引。统计数据的存在允许 LAS 数据集图层使用过滤和符号系统选项,以便仅显示 LAS 文件中存在的 LAS 属性值。将为每个 LAS 文件创建 .lasx 辅助文件。
| Boolean |
relative_paths (可选) | 指定 LAS 数据集是通过相对路径还是绝对路径来引用激光雷达文件和表面约束要素。在文件系统中使用同一相对路径将 LAS 数据集及其关联的数据重新定位到其他位置时,使用相对路径会比较方便。
| Boolean |
create_las_prj | 指定是否会为 LAS 数据集引用的 LAS 文件创建 .prj 文件。
| String |
代码示例
CreateLasDataset 示例 1(Python 窗口)
下面的示例演示了如何在 Python 窗口中使用此工具。
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.CreateLasDataset_management("folder_a; folder_b/5S4E.las",
"test.lasd", "RECURSION",
"LA/boundary.shp <None> Softclip;"\
"LA/ridges.shp Elevation hardline", "",
"COMPUTE_STATS", "RELATIVE_PATHS")
CreateLasDataset 示例 2(独立脚本)
下面的示例演示了如何在独立 Python 脚本中使用此工具。
'''*********************************************************************
Name: Export Elevation Raster from Ground LAS Measurements
Description: This script demonstrates how to export
ground measurements from LAS files to a raster using a
LAS dataset. This sample is designed to be used as a script
tool.
*********************************************************************'''
# Import system modules
import arcpy
try:
# Set Local Variables
inLas = arcpy.GetParameterAsText(0)
recursion = arcpy.GetParameterAsText(1)
surfCons = arcpy.GetParameterAsText(2)
classCode = arcpy.GetParameterAsText(3)
returnValue = arcpy.GetParameterAsText(4)
spatialRef = arcpy.GetParameterAsText(5)
lasD = arcpy.GetParameterAsText(6)
outRaster = arcpy.GetParameterAsText(7)
cellSize = arcpy.GetParameter(8)
zFactor = arcpy.GetParameter(9)
if arcpy.ProductInfo == 'ArcView':
arcpy.CheckOutExtension('3D')
# Execute CreateLasDataset
arcpy.management.CreateLasDataset(inLas, lasD, recursion, surfCons, sr)
# Execute MakeLasDatasetLayer
lasLyr = arcpy.CreateUniqueName('Baltimore')
arcpy.management.MakeLasDatasetLayer(lasD, lasLyr, classCode, returnValue)
# Execute LasDatasetToRaster
arcpy.conversion.LasDatasetToRaster(lasLyr, outRaster, 'ELEVATION',
'TRIANGULATION LINEAR WINDOW_SIZE 10', 'FLOAT',
'CELLSIZE', cellSize, zFactor)
arcpy.GetMessages()
except arcpy.ExecuteError:
print(arcpy.GetMessages())
except Exception as err:
print(err.args[0])
finally:
arcpy.management.Delete(lasLyr)
环境
许可信息
- ArcGIS Desktop Basic: 需要 3D Analyst 或 Spatial Analyst
- ArcGIS Desktop Standard: 是
- ArcGIS Desktop Advanced: 是