ArcMap

LAS 数据集转 TIN

描述

通过 LAS 数据集导出不规则三角网 (TIN)。

插图

LAS 数据集转 TIN

使用方法

  • 您可通过在图层过滤器设置中选择分类代码、分类标记和返回值的任意组合来使 LAS 数据集图层对所显示和处理的 LAS 点进行限制。过滤器的定义可通过图层属性对话或创建 LAS 数据集图层工具实现。

  • LAS 数据集图层还可用于控制表面约束要素的强化，而该要素可由 LAS 数据集引用。当将 LAS 数据集作为三角化网格面显示或处理时，将强制执行约束。

  • 虽然 TIN 可以支持的总点数可以超过 1500 万个点，但建议您将 TIN 数据集限制为不超过 500 万个点，以确保在显示和分析数据时的响应性能。可以使用点稀疏化方法和控制输出处理范围来减少 TIN 结点计数。

语法

arcpy.ddd.LasDatasetToTin(in_las_dataset, out_tin, {thinning_type}, {thinning_method}, {thinning_value}, {max_nodes}, {z_factor}, {clip_to_extent})
参数说明数据类型
in_las_dataset

待处理的 LAS 数据集。

LAS Dataset Layer
out_tin

将要生成的 TIN 数据集。

TIN
thinning_type
（可选）

指定用于选择可导出至 TIN 的 LAS 数据点子集的技术。

  • NONE —不应用稀疏化功能。这是默认设置。
  • RANDOM —基于相应的 thinning_method 选择和 thinning_value 条目随机选择 LAS 数据点。
  • WINDOW_SIZE —LAS 数据集分为由 thinning_value 定义的方形切片，并使用 thinning_method 选择 LAS 点。
String
thinning_method
（可选）

指定用于减少 LAS 数据点的技术，该技术将对稀疏化值的解释产生影响。可用选项取决于所选的稀疏化类型。

  • PERCENT —稀疏化值将反映将在输出中保留的 LAS 点的百分比
  • NODE_COUNT —稀疏化值将反映输出中所允许的结点总数。
  • MIN —在每个窗口大小区域中选择具有最低高程的 LAS 数据点。
  • MAX —在每个自动确定的窗口大小区域中选择具有最高高程的 LAS 数据点。
  • CLOSEST_TO_MEAN —选择高程最接近自动确定的窗口大小区域中平均值的 LAS 数据点。

指定用于减少 LAS 数据点的技术，该技术将对 thinning_value 的解释产生影响。可用选项取决于所选的 thinning_type。

  • PERCENT — thinning_value 将反映 LAS 数据集中总点数的百分比。仅当选择 thinning_type="RANDOM" 时才会提供此选项。
  • NODE_COUNT —thinning_value 将反映输出中所允许的节点总数。仅当选择 thinning_type="RANDOM" 时才会提供此选项。
  • MIN —在每个窗口大小区域中选择具有最低高程的 LAS 点。仅当选择 thinning_type="WINDOW_SIZE" 时才会提供此选项。
  • MAX —在每个窗口大小区域中选择具有最高高程的 LAS 点。仅当选择 thinning_type="WINDOW_SIZE" 时才会提供此选项。
  • CLOSEST_TO_MEAN —在每个窗口大小区域中选择高程与所有 LAS 点的平均值最接近的 LAS 点。仅当选择 thinning_type="WINDOW_SIZE" 时才会提供此选项。
String
thinning_value
（可选）

如果选择 thinning_type="WINDOW_SIZE"，该值表示将要划分 LAS 数据集的采样区。

如果选择 thinning_type="RANDOM" 和 thinning_method="PERCENT"，该值表示要被导出到 TIN 的 LAS 数据集的点百分比。

如果选择 thinning_type="RANDOM" 和 thinning_method="NODE_COUNT"，该值表示可被导出到 TIN 的 LAS 点的总数。

Double
max_nodes
（可选）

输出 TIN 中允许的结点的最大数量。默认值为 5 百万。

Double
z_factor
（可选）

Z 值将乘上的系数。此值通常用于转换 Z 线性单位来匹配 XY 线性单位。默认值为 1，此时高程值保持不变。如果输入表面的空间参考具有已指定线性单位的 Z 基准，则将禁用此参数。

Double
clip_to_extent
（可选）

指定是否根据分析范围裁剪生成的 TIN。仅当分析范围是输入 LAS 数据集的子集时，该选项才有效。

  • CLIP —根据分析范围裁剪输出 TIN。这是默认设置。
  • NO_CLIP —不根据分析范围裁剪输出 TIN。
Boolean

代码示例

LasDatasetToTin 示例 1（Python 窗口）

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

import arcpy
from arcpy import env

arcpy.CheckOutExtension('3D')
env.workspace = 'C:/data'

arcpy.LasDatasetToTin_3d('se_baltimore.lasd', 'se_bmore', 'RANDOM', 15, 3.28)
LasDatasetToTin 示例 2（独立脚本）

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

'''**********************************************************************
Name: LAS Dataset to TIN Example
Description: Create a TIN using bare earth lidar measurements. This 
             script is 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
surfCons = arcpy.GetParameterAsText(2) #input surface constraints
sr = arcpy.GetParameter(3) #spatial reference of las dataset
outTin = arcpy.GetParameterAsText(4)
thinningType = arcpy.GetParameterAsText(5)
thinningMethod = arcpy.GetParameterAsText(6)
thinningValue = arcpy.GetParameter(7)
zFactor = arcpy.GetParameter(8)

try:
    arcpy.CheckOutExtension('3D')
    # Execute CreateLasDataset
    arcpy.management.CreateLasDataset(inLas, lasD, 'RECURSION', surfCons, sr)
    lasLyr = arcpy.CreateUniqueName('lasdToTin', 'in_memory')
    classCode = 2
    returnValue = 'LAST'
    # Execute MakeLasDatasetLayer
    arcpy.management.MakeLasDatasetLayer(lasD, lasLyr, classCode, returnValue)
    # Define extent of the area of interest
    env.extent(1426057, 606477, 1449836, 623246)
    # Execute LasDatasetToTin
    arcpy.ddd.LasDatasetToTin(lasLyr, outTin, thinningType, 
                              thinningMethod, thinningValue, zFactor)
    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)

环境

许可信息

  • Basic: 需要 3D Analyst
  • Standard: 需要 3D Analyst
  • Advanced: 需要 3D Analyst