ArcMap

Статистика точек LAS как растра

Сводка

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

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

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

  • Формат LAS поддерживает классификацию точек на основе спецификаций, разработанных обществом American Society for Photogrammetry and Remote Sensing (ASPRS). На платформе ArcGIS применяется схема классификации, предназначенная для LAS-файлов версии 1.4:

    Значение классификации Тип классификации

    0

    Классификация не выполнялась

    1

    Не назначено

    2

    Земля

    3

    Низкая растительность

    4

    Средняя растительность

    5

    Высокая растительность

    6

    Здание

    7

    Низкий шум

    8

    Ключевые для модели / Зарезервировано

    9

    Вода

    10

    Железная дорога

    11

    Дорожное покрытие

    12

    Наложение / Зарезервировано

    13

    Проволочное ограждение

    14

    Провод

    15

    Опора ЛЭП

    16

    Изолятор

    17

    Мостовой настил

    18

    Высокий шум

    19 – 63

    Зарезервировано для Определения ASPRS (LAS от 1.1 до 1.3 поддерживают коды классов до 31)

    32 – 255

    Задается пользователем (поддерживается только в LAS 1.0 и в некоторых версиях 1.4)

  • Для определения местоположений с самым высоким отраженным сигналом, что может указывать на наличие растительности, рекомендуется использовать значение PREDOMINANT_LAST_RETURN для параметра Method.

  • Для определения местоположений с потенциальными выбросами можно использовать значения Z_RANGE для параметра Method.

Синтаксис

arcpy.management.LasPointStatsAsRaster(in_las_dataset, out_raster, {method}, {sampling_type}, {sampling_value})
ПараметрОбъяснениеТип данных
in_las_dataset

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

LAS Dataset Layer
out_raster

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

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

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

Тип собираемых статистических данных по точкам LAS в каждой из ячеек выходного растра.

  • PULSE_COUNT —Количество точек последнего отраженного сигнала.
  • POINT_COUNT —Количество точек всех отраженных сигналов.
  • PREDOMINANT_LAST_RETURN —Наиболее частое значение последнего отраженного сигнала.
  • PREDOMINANT_CLASS —Наиболее часто встречаемый код классификации.
  • INTENSITY_RANGE —Диапазон значений интенсивности.
  • Z_RANGE —Диапазон значений высот.
String
sampling_type
(Дополнительный)

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

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

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

Double

Пример кода

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

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

import arcpy
from arcpy import env

env.workspace = "C:/data"
arcpy.LasPointStatsAsRaster_3d("test.lasd", "lidar_intensity.img", 
                             "INTENSITY_RANGE", "CELLSIZE", 15)
LasPointStatsAsRaster, пример 2 (автономный скрипт)

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

'''**********************************************************************
Name: LAS Point Statistics As Raster
Description: Identifies the most frequently occurring return value for 
             each pulse in a given set of LAS files. 
             Designed for use as a script tool.
**********************************************************************'''
# Import system modules
import arcpy

# Set Local Variables
lasD = arcpy.GetParameterAsText(0)
inLas = arcpy.GetParameterAsText(1) #input las files
sr = arcpy.GetParameter(2) #spatial reference of las dataset
statsRaster = arcpy.GetParameterAsText(3)

try:
    arcpy.CheckOutExtension('3D')
    # Execute CreateLasDataset
    arcpy.management.CreateLasDataset(inLas, lasD, 'RECURSION', '', sr)
    # Execute LasPointStatsAsRaster
    arcpy.management.LasPointStatsAsRaster(lasD, statsRaster, 
                                           "PREDOMINANT_RETURNS_PER_PULSE", 
                                           "CELLSIZE", 15)
    arcpy.CheckInExtension('3D')
except arcpy.ExecuteError:
    print(arcpy.GetMessages())
except Exception as err:
    print(err.args[0])

Параметры среды

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

  • Basic: Требуется Дополнительный модуль 3D Analyst или Spatial Analyst
  • Standard: Да
  • Advanced: Да

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