LAS 数据集转 TIN (3D Analyst)
摘要
通过 LAS 数据集导出不规则三角网 (TIN)。
插图
用法
将 LAS 数据集指定为输入时,将处理它引用的 LAS 文件中的所有数据点。
LAS 数据集图层可用于按类代码或返回值过滤 LAS 点。创建图层的方法是使用创建 LAS 数据集图层工具、通过在 ArcMap 或 ArcScene 中加载 LAS 数据集以及通过图层属性对话框指定所需类代码和返回值。
-
MINIMUM、MAXIMUM 和 AVERAGE 点选择选项利用自动计算的窗口大小区域,该区域旨在生成 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 数据点的细化类型。
| String |
thinning_method (可选) |
细化方法定义了用于减少 LAS 数据点的具体方法,并会对细化值的解释方式产生影响。可用选项取决于所选的细化类型。 对于 RANDOM:
对于 WINDOW_SIZE:
| String |
thinning_value (可选) |
与所选细化类型和细化方法相关联的值。 RANDOM 点选择方法可用的细化方法:
对于任何 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