Набор данных LAS в растр (Конвертация)

Сводка

Создает растр, используя значения высоты, интенсивности или цвета RGB, на основе лидарных точек, на которые ссылается набор данных LAS.

Иллюстрация

Набор данных LAS в растр

Использование

  • Для слоя набора данных LAS может быть установлено ограничение количества отображаемых и обрабатываемых точек LAS, путем выбора любой комбинации кодов классификации, флагов классификации и возвращаемых значений в настройках фильтра слоя. Фильтры можно ввести в диалоговом окне Свойства слоя или с помощью инструмента Создать слой набора данных LAS.

  • Слой набора данных LAS также может быть использован для управления объектами поверхностных ограничений, которые могут находиться в наборе данных LAS. Ограничения применяются при отображении или обработке набора данных LAS в виде триангуляционной поверхности.

  • При экспорте большого растра следует указать Выходной тип данных как целое число, чтобы сэкономить место на диске, если точность ваших данных позволяет использовать целочисленные z-значения.

  • Настоятельно рекомендуется в определении набора данных LAS входные данные ограничивать изучаемой областью. Одной из причин является предотвращение интерполяции за пределами реального экстента изучаемых данных. Другой причиной является возможность значительного снижения производительности при использовании параметров естественной окрестности, если область данных определена неверно.

    Примечание:

    При использовании BINNING задействуются только ограничения на вырезание, стирание и замену. Линии перегибов и точки якоря не используются. Опция триангуляции поддерживает все типы ограничений, однако её выполнение занимает более длительное время.

Синтаксис

LasDatasetToRaster(in_las_dataset, out_raster, {value_field}, {interpolation_type}, {data_type}, {sampling_type}, {sampling_value}, {z_factor})
ParameterОбъяснениеТип данных
in_las_dataset

Обрабатываемый набор данных LAS.

LAS Dataset Layer
out_raster

Местоположение и имя выходного растра. При хранении набора растровых данных в базе геоданных или в папке, например в Esri Grid, расширение файла к имени набора растровых данных добавлять не нужно. Расширение файла, которое можно ввести, чтобы задать формат растра при его сохранении в папке, например, .tif для создания GeoTIFF или .img для создания файла формата ERDAS IMAGINE.

Если растр хранится в базе геоданных или в виде файла TIFF, тип и качество сжатия этого растра можно задать с помощью настроек среды геообработки.

Raster Dataset
value_field
(Дополнительный)

Данные лидара, используемые для создания выходного растра.

  • ELEVATIONДля создания растра будут использоваться высоты из лидарных файлов. Используется по умолчанию.
  • INTENSITYДля создания растра будет использоваться значения интенсивности из лидарных файлов.
  • RGBЗначения RGB точек лидара будут использованы для создания трехканального изображения.
String
interpolation_type
"BINNING {cell_assignment_type} {void_fill_method}" or "TRIANGULATION {interpolation_method} {point_thinning_type} {point_selection_method} {resolution}"
(Дополнительный)

Метод интерполяции, который будет использоваться для определения значений ячеек выходного растра.

Бининг включает: Метод присвоения значений ячейкам для определения каждой выходной ячейки по точкам, попадающим в ячейку, и Метод заполнения пустот для определения значений ячеек, не содержащих точек LAS.

Методы присвоения значений ячейкам

  • AVERAGE— назначение среднего значения всех точек в ячейке. Используется по умолчанию.
  • MINIMUM— назначение минимального значения, найденного в точках ячейки.
  • MAXIMUM – назначение максимального значения, найденного в точках ячейки.
  • IDW— использование интерполяции по методу взвешенных расстояний для определения значения ячейки.
  • NEAREST — использование метода назначения ближайшей окрестности для определения значения ячейки.

Методы заполнения пустот

  • NONE — ячейке присваивается значение NoData.
  • SIMPLE – служит для вычисления средних значений в ячейках данных непосредственно вокруг ячейки NoData во избежание небольших пропусков.
  • LINEAR — триангуляция незанятых областей и использование линейной интерполяции триангулированного значения для определения значения ячейки. Используется по умолчанию.
  • NATURAL_NEIGHBOR — использование интерполяции по методу естественной окрестности для определения значения ячейки.

Методы интерполяции Триангуляция вычисляют значения ячеек с использованием подхода, схожего с построением TIN, предоставляя также возможность ускорить обработку путем прореживания наборов данных LAS методом Размера окна.

Методы триангуляции

  • Линейная – для получения значений ячеек используется линейная интерполяция.
  • Естественная окрестность – для определения значения ячейки используется интерполяция по методу естественной окрестности.

Методы выбора размера окна

  • Максимум – сохраняется точка с самым высоким значением в каждом размере окна. Используется по умолчанию.
  • Минимум – сохраняется точка с самым низким значением в каждом размере окна.
  • Ближайший к среднему – сохраняется точка, находящаяся ближе всего к среднему значению всех точек в размере окна.
Interpolate
data_type
(Дополнительный)

Указывает тип числовых значений, хранящихся в выходном растре.

  • FLOATВыходной растр будет использовать 32-битовое число с плавающей точкой со значениями от -3,402823466e+38 до 3,402823466e+38. Используется по умолчанию.
  • INTВыходной растр будет использовать соответствующую целочисленную разрядность. При использовании этой опции z-значения будут округлены до целого числа, которое будет использовано в качестве значения ячейки растра.
String
sampling_type
(Дополнительный)

Указывает метод, используемый для интерпретации Значения выборки для определения разрешения выходного растра.

  • OBSERVATIONSОпределяет число ячеек, разделяющих самую длинную сторону экстента набора данных LAS.
  • CELLSIZEОпределяет размер ячейки выходного набора растровых данных. Используется по умолчанию.
String
sampling_value
(Дополнительный)

Задает значение, используемое вместе с параметром Значение выборки для определения разрешения выходного растра.

Double
z_factor
(Дополнительный)

Коэффициент, на который будут умножаться z-значения. Этот способ обычно используется для конвертации линейных Z-единиц для соответствия линейным единицам XY. Значение по умолчанию равно 1, при нем значения высоты не изменяются. Этот параметр отключен, если пространственная привязка входной поверхности имеет Z-датум с заданными линейными единицами измерения.

Double

Пример кода

LasDatasetToRaster, пример 1 (окно Python)

В следующем примере показано использование этого инструмента в окне Python.

arcpy.env.workspace = 'C:/data'
arcpy.LasDatasetToRaster_3d('baltimore.lasd', 'baltimore.tif', 'INTENSITY',
                          'TRIANGULATION LINEAR WINDOW_SIZE 10', 'FLOAT',
                          'CELLSIZE', 10, 3.28)
LasDatasetToRaster, пример 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)

    # 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)
    print(arcpy.GetMessages())

except arcpy.ExecuteError:
    print(arcpy.GetMessages())

except Exception as err:
    print(err.args[0])

finally:
    arcpy.management.Delete(lasLyr)

Информация о лицензиях

  • Basic: Требуется Spatial Analyst или 3D Analyst
  • Standard: Да
  • Advanced: Да

Связанные разделы