LAS 数据集转 TIN (3D Analyst)

许可等级:BasicStandardAdvanced

摘要

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

插图

LAS Dataset to TIN

用法

语法

LasDatasetToTin_3d (in_las_dataset, out_tin, {thinning_type}, {thinning_method}, {thinning_value}, {max_nodes}, {z_factor})
参数说明数据类型
in_las_dataset

输入 LAS 数据集。

LAS Dataset Layer
out_tin

输出 TIN 数据集。

TIN
thinning_type
(可选)

用于减少生成的 TIN 中被保存为节点的 LAS 数据点的细化类型。

  • NONE - 未应用细化,并且不需要细化方法或细化值。这是默认设置。
  • RANDOM - 基于相应的细化方法选择和细化值条目对 LAS 数据点进行随机选择。
  • WINDOW_SIZE - 通过对每个由细化值定义的方形区域进行评估并使用细化方法选择 LAS 点的方式来减少 LAS 数据点。
String
thinning_method
(可选)

细化方法定义了用于减少 LAS 数据点的具体方法,并会对细化值的解释方式产生影响。可用选项取决于所选的细化类型

对于 RANDOM

  • PERCENT细化值将反映 LAS 数据集全分辨率的节点的百分比。
  • NODE_COUNT细化值将反映输出中所允许的节点总数。

对于 WINDOW_SIZE

  • MINIMUM在每个自动确定的窗口大小区域中选择具有最低高程的 LAS 数据点。
  • MAXIMUM在每个自动确定的窗口大小区域中选择具有最高高程的 LAS 数据点。
  • CLOSEST_TO_MEAN选择高程最接近自动确定的窗口大小区域中平均值的 LAS 数据点。
String
thinning_value
(可选)

与所选细化类型细化方法相关联的值。

RANDOM 点选择方法可用的细化方法:

  • PERCENT - 值将表示 LAS 数据集全分辨率的数据点的百分比。
  • NODE_COUNT - 值将表示输出 TIN 中所允许的节点总数。

对于任何 WINDOW_SIZE 细化方法来说,值都表示 LAS 数据集范围所分成的用于对数据点进行采样的区域。

Double
max_nodes
(可选)

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

Double
z_factor
(可选)

高程值将乘上的系数。此系数通常用于转换 Z 线性单位,以匹配 XY 线性单位的值。默认值为 1,此时高程值保持不变。

Double

代码实例

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)

环境

许可信息

ArcGIS for Desktop Basic:需要 3D Analyst
ArcGIS for Desktop Standard:需要 3D Analyst
ArcGIS for Desktop Advanced:需要 3D Analyst
9/15/2013