Набор данных LAS в TIN (3D Analyst)
Резюме
Экспортирует нерегулярную сеть треугольников TIN на основе набора данных лазерного сканирования LAS.
Рисунок
Использование
Если набор данных LAS указан в качестве входных данных, будут обработаны все точки данных в LAS-файлах, на которые он ссылается.
Слой набора данных LAS можно использовать для фильтрации точек LAS по коду класса или возвращаемым значениям. Слой можно создать с помощью инструмента Создать слой набора данных LAS (Make LAS Dataset Layer) или загрузив набор данных LAS в ArcMap или ArcScene и указав нужные коды классов и возвращаемые значения в диалоговом окне свойств слоя.
-
Параметры выборки точек Минимум (MINIMUM), Максимум (MAXIMUM) и Среднее (AVERAGE) служат для использования автоматически рассчитанной области размера окна в целях построения набора данных TIN с числом узлов, наиболее удобным для быстродействия системы.
Синтаксис
| Параметр | Объяснение | Тип данных |
in_las_dataset |
Входной набор данных LAS. | LAS Dataset Layer |
out_tin |
Выходной набор данных TIN. | TIN |
thinning_type (дополнительно) |
Тип прореживания, используемый для уменьшения числа точек в данных LAS для получения узлов TIN.
| String |
thinning_method (дополнительно) |
Метод прореживания определяет методику уменьшения точек в данных LAS, а также влияет на трактовку Значения прореживания (Thinning Value). Доступные опции зависят от выбранного Типа прореживания (Thinning Type). Для типа Случайный (RANDOM):
Для типа Размер окна (WINDOW_SIZE):
| String |
thinning_value (дополнительно) |
Значение, связанное с выбранным Типом прореживания (Thinning Type) и Методом прореживания (Thinning Method). Для методов прореживания, доступных при типе выбора точек Случайный (RANDOM):
Для любого из методов прореживания типа Размер окна (WINDOW_SIZE) значение характеризует области, на которые делится экстент набора данных LAS для выборки точек. | Double |
max_nodes (дополнительно) |
Максимальное количество узлов, разрешенное для выходной TIN. Значение по умолчанию — 5 миллионов узлов. | Double |
z_factor (дополнительно) |
Коэффициент, на который будут умножаться значения высоты. Обычно он используется для преобразования линейных единиц измерения Z, соответствующих линейным единицам измерения XY. Значение по умолчанию равно 1, при нем значения высоты не изменяются. | Double |
Пример кода
В следующем примере показано использование этого инструмента в окне 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)
В следующем примере показано использование этого инструмента в автономном скрипте 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)