segmentGroundSMRF

Земля сегмента из данных о лидаре с помощью Простого морфологического фильтра (SMRF) алгоритм

свернуть все на странице

Синтаксис

groundPtsIdx = segmentGroundSMRF(ptCloud)
groundPtsIdx = segmentGroundSMRF(ptCloud,gridResolution)
[groundPtsIdx,nonGroundPtCloud,groundPtCloud] = segmentGroundSMRF(___)
[___] = segmentGroundSMRF(___,Name,Value)

Описание

пример

groundPtsIdx = segmentGroundSMRF(ptCloud) сегментирует облако точки ввода, ptCloud в землю и неземлю указывает и возвращает логический матричный или векторный groundPtsIdx. Функция устанавливает наземные индексы точки на true и false для неназемных точек.

groundPtsIdx = segmentGroundSMRF(ptCloud,gridResolution) дополнительно задает размерность элемента сетки.

[groundPtsIdx,nonGroundPtCloud,groundPtCloud] = segmentGroundSMRF(___) дополнительно возвращает наземные точки и неназемные точки как отдельный pointCloud объекты. Используйте этот синтаксис с любой из комбинаций входных аргументов в предыдущих синтаксисах.

пример

[___] = segmentGroundSMRF(___,Name,Value) задает опции с помощью одних или нескольких аргументов name-value. Например, 'ElevationThreshold',0.4 устанавливает порог вертикального изменения для идентификации, что неземля указывает на 0,4.

Примеры

свернуть все

Скрипт Open Live Script

Сегментируйте землю в неорганизованном воздушном облаке точек.

Создайте lasFileReader возразите, чтобы получить доступ к данным о файле LAS.

fileName = fullfile(toolboxdir('lidar'), 'lidardata', 'las', ...
    'aerialLidarData2.las');
lasReader = lasFileReader(fileName);

Считайте данные об облаке точек из файла LAS с помощью readPointCloud функция.

ptCloud = readPointCloud(lasReader);

Наземные данные о сегменте из облака точек.

[groundPtsIdx,nonGroundPtCloud,groundPtCloud] = segmentGroundSMRF(ptCloud);

Визуализируйте наземные и неназемные точки.

figure
pcshowpair(groundPtCloud, nonGroundPtCloud)

Figure contains an axes object. The axes object contains 2 objects of type scatter.

Скрипт Open Live Script

Сегмент и удаляет землю из организованного облака точек. Облако точек было получено в магистральном сценарии.

Загрузите данные об облаке точек в рабочую область. 

ld = load('drivingLidarPoints.mat');

Отобразите облако точки ввода.

pcshow(ld.ptCloud)
xlim([-40 40])
ylim([-50 50])

Figure contains an axes object. The axes object contains an object of type scatter.

Наземные данные о сегменте из облака точек.

[~,nonGroundPtCloud,groundPtCloud] = segmentGroundSMRF(...,
  ld.ptCloud,'MaxWindowRadius',5,'ElevationThreshold',0.1,'ElevationScale',0.25);

Визуализируйте неназемные точки.

figure
pcshow(nonGroundPtCloud)
xlim([-40 40])
ylim([-50 50])

Figure contains an axes object. The axes object contains an object of type scatter.

Входные параметры

свернуть все

Данные об облаке точек в виде pointCloud объект.

Размерность каждого элемента сетки в виде положительной скалярной величины.

Типы данных: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

Аргументы name-value

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

Пример: 'ElevationThreshold',0.4 устанавливает порог вертикального изменения, чтобы идентифицировать, что неземля указывает на 0,4.

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

Примечание

Значение по умолчанию работает эффективно на воздушные данные о лидаре. Для лучшей эффективности на наземных данных, набор MaxWindowRadius к меньшему значению как 8.

Типы данных: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

Наклонный порог, чтобы идентифицировать неназемные элементы сетки в минимальном вертикальном изменении появляется карта в виде неотрицательного скаляра. Функция классифицирует элемент сетки как неземлю, если ее наклон больше SlopeThreshold. Увеличьте это значение, чтобы классифицировать крутые склоны как землю.

Типы данных: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

Порог вертикального изменения, чтобы идентифицировать неземлю указывает в виде неотрицательного скаляра. Функция классифицирует точку как неземлю, если различие в вертикальном изменении между точкой и оцененной земной поверхностью больше ElevationThreshold. Увеличьте это значение, чтобы охватить больше точек от ухабистой земли.

Примечание

Значение по умолчанию работает эффективно на воздушные данные о лидаре. Для лучших результатов на наземных данных, набор ElevationThreshold к меньшему значению как 0.1.

Типы данных: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

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

Типы данных: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

Выходные аргументы

свернуть все

Бинарная карта сегментированного облака точек, возвращенного как логическая матрица для организованных облаков точек, и как логический вектор для неорганизованных облаков точек. Функция устанавливает местоположения наземных точек в матрице к true и не оснуйте точки к false.

Облако точек неназемных точек, возвращенных как pointCloud объект.

Облако точек наземных точек, возвращенных как pointCloud объект.

Алгоритмы

Простой морфологический фильтр (SMRF) алгоритм [1] данные об облаке точек сегментов в наземные и неназемные точки. Алгоритм разделен на три этапа:

  1. Создайте минимальную поверхность вертикального изменения из данных об облаке точек.

  2. Сегментируйте поверхность на наземные и неназемные элементы сетки.

  3. Сегментируйте исходные данные об облаке точек.

Минимальное поверхностное создание

  1. Разделите данные об облаке точек на сетку по измерению xy- (вид с высоты птичьего полета). Задайте размерность элемента сетки с помощью gridResolution.

  2. Найдите самое низкое вертикальное изменение (Zmin) значение для каждого элемента сетки (пиксель).

  3. Объедините все значения Zmin в 2D матрицу (растровое изображение), чтобы создать минимальную карту поверхности вертикального изменения.

Поверхностная сегментация карты

  1. Примените морфологическую вводную операцию на минимальную поверхностную карту. Для получения дополнительной информации о морфологическом открытии, смотрите Типы Морфологических Операций.

  2. Используйте дискообразный элемент структурирования с радиусом 1 пикселя. Для получения дополнительной информации смотрите Элементы Структурирования.

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

  4. Выполните шаги 1 - 3 итеративно. Увеличьте радиус элемента структурирования на 1 пиксель в каждой итерации, пока это не достигнет максимального радиуса, заданного MaxWindowRadius.

  5. Конечный результат процесса итерации является бинарной маской, где каждый пиксель классифицируется как являющийся или землей или неземлей.

Сегментация облака точек

  1. Примените бинарную маску на исходную минимальную поверхностную карту, чтобы удалить незаземляющие сетки.

  2. Заполните незаполненные сетки с помощью методов интерполяции изображения, чтобы создать предполагаемую модель вертикального изменения.

  3. Вычислите разность в вертикальном изменении между каждой точкой в исходном облаке точек и предполагаемой моделью вертикального изменения. Если различие больше ElevationThreshold, классифицируйте пиксель как неземлю.

Ссылки

[1] Pingel, Томас Дж., Кит К. Кларк и Уильям А. Макбрайд. “Улучшенный простой морфологический фильтр для классификации ландшафтов бортовых данных о ЛИДАРЕ”. Журнал ISPRS фотограмметрии и дистанционного зондирования 77 (март 2013): 21–30.

Расширенные возможности

Смотрите также

Функции

Объекты

Введенный в R2021a

Документация Lidar Toolbox

Поддержка

Для просмотра документации необходимо авторизоваться на сайте