Я ищу метод для обработки изображения дистанционного зондирования и извлечения областей кроны отдельных деревьев из изображения.
У меня есть как визуальные снимки с использованием длин волн, так и лидарные данные из этой области. Место, о котором идет речь, является пустынным районом, поэтому древесный покров не такой плотный, как лесной массив. Разрешение аэрофотоснимков составляет 0,5 фута на 0,5 фута. Разрешение лидара составляет приблизительно 1 x 1 фут. И визуальные данные, и лидар поступают из набора данных округа Пима, штат Аризона. Образец типа аэрофотоснимков, который у меня есть, находится в конце этого поста.
Этот вопрос Обнаружение единого дерева в ArcMap? кажется, та же самая проблема, но там, кажется, нет хорошего ответа.
Я могу получить разумную классификацию типов растительности (и информацию об общем процентном покрытии) в области, используя классификацию Iso Cluster в Arcmap, но это дает мало информации об отдельных деревьях. Ближе всего к тому, что я хочу, это результаты прохождения выходных данных классификации изокластера через растр в Polygon в Arcmap. Проблема в том, что этот метод объединяет деревья рядом в один многоугольник.
Изменить: я, вероятно, должен был бы включить некоторые более подробные сведения о том, что у меня есть. У меня есть наборы необработанных данных:
- Полные данные LAS и сгенерированный из него TIFF растр.
- Визуальные образы (например, показанный образец изображения, но охватывающий гораздо более широкую область)
- Ручные прямые измерения подмножества деревьев в этом районе.
Из них я сгенерировал:
- Классификация почв / растительности.
- Растры DEM / DSM.
источник
Ответы:
По спектральным и лидарным данным имеется значительное количество литературы по обнаружению индивидуальной короны. Мудрые методы, возможно, начнем с:
Фальковский, MJ, AMS Smith, PE Gessler, AT Hudak, LA Vierling и JS Evans. (2008). Влияние покрытия полога хвойного леса на точность двух отдельных алгоритмов измерения деревьев с использованием лидарных данных. Канадский журнал дистанционного зондирования 34 (2): 338-350.
Smith AMS, EK Strand, CM Steele, DB Hann, SR Garrity, MJ Falkowski, JS Evans (2008). Создание карт пространственной структуры растительности путем анализа каждого объекта вторжения можжевельника на многовременных аэрофотоснимках. Канадский журнал дистанционного зондирования 34 (2): 268-285
Если вас интересует метод Wavelet (Smith et al., 2008), я написал его на Python, но он очень медленный. Если у вас есть опыт работы с Matlab, то здесь он реализован в производственном режиме. У нас есть две работы, в которых мы идентифицировали ~ 6 миллионов акров вторжения можжевельника в восточном Орегоне, используя вейвлет-метод с изображениями NAIP RGB-NIR, поэтому это хорошо доказано.
Baruch-Mordo, S., JS Evans, J. Severson, JD Naugle, J. Kiesecker, J. Maestas и MJ Falkowski (2013) Спасение шалфея от деревьев: упреждающее решение для снижения ключевой угрозы для кандидата Биологическая консервация видов 167: 233-241
Познанович, AJ, MJ Falkowski, AL Maclean и JS Evans (2014) Оценка точности алгоритмов обнаружения деревьев в лесах можжевельника. Фотограмметрическая инженерия и дистанционное зондирование 80 (5): 627–637
Есть несколько интересных подходов, в общем случае декомпозиции объектов, из прикладной математической литературы по пространству состояний, использующей многоразрешающие гауссовские процессы для декомпозиции характеристик объектов в масштабе. Я использую эти типы моделей для описания многомасштабного процесса в экологических моделях, но он может быть адаптирован для разложения характеристик объекта изображения. Весело, но немного эзотерично.
Gramacy, RB и HKH Lee (2008). Модели гауссовского процесса с байесовским анализом с применением к компьютерному моделированию. Журнал Американской статистической ассоциации, 103 (483): 1119–1130
Kim, HM, BK Mallick, и CC Holmes (2005). Анализ нестационарных пространственных данных с использованием кусочно-гауссовских процессов. Журнал Американской статистической ассоциации, 100 (470): 653–668
источник
Чтобы создать DHM вычесть DEM из DEM, это можно сделать в Esri Raster Calculator или GDAL_CALC . Это поместит все ваши возвышения на «ровное игровое поле».
Синтаксис (замена полных путей для DEM, DSM и DHM):
Значение DHM будет в основном 0 (или достаточно близко), что вы задаете в качестве значения нодаты. С помощью Raster Calculator или GDAL_CALC вы можете извлекать значения, превышающие произвольное значение, на основе количества шума, которое вы наблюдаете в DHM. Цель этого состоит в том, чтобы уменьшить шум и выделить только кроны растительности - в случае, когда два «дерева» соседствуют, это должно быть разделено на две отдельные капли.
Синтаксис (замените полные пути для двоичного и DHM и наблюдаемое значение для значения):
Теперь с помощью GDAL_CALC или Esri IsNull создайте двоичный растр, который можно полигонизировать с помощью GDAL_Polygonize или Esri Raster to Polygon .
Чтобы уточнить полигоны, удалите чрезмерно маленькие полигоны и затем сравните их с полосами RGB, которые ищут подписи, в Esri поможет инструмент Zonal Statistics . Затем вы можете отказаться от полигонов, которые явно не имеют правильной статистики (основываясь на экспериментах и ваших данных, я не могу дать вам значения).
Это должно дать вам около 80% точности при построении отдельных коронок.
источник
eCognition - лучшее программное обеспечение для этого, я сделал это, используя другое программное обеспечение, но eCognition лучше. Вот ссылка на литературу на эту тему:
например, http://www.mdpi.com/1424-8220/14/12/22643
Дополнительно:
например, http://www.nrcresearchpress.com/doi/abs/10.1139/x05-030#.VJmMb14gAA
источник
У меня была такая же проблема пару лет назад. У меня есть решение, которое не требует отфильтрованных данных LAS или других вспомогательных данных. Если у вас есть доступ к данным LiDAR, и вы можете генерировать DEM / DSM / DHM (далее DEM, я не обсуждаю семантику номенклатуры модели поверхности) из разных возвратов, может пригодиться следующий скрипт.
Скрипт arcpy принимает 3 матрицы высот и выплевывает лесной полигон и шейп-файлы деревьев. 3 матрицы высот должны иметь одинаковое пространственное разрешение (т.е. 1 метр) и экстенты, и представлять первые возвращения, последние возвращения и голую землю. У меня были очень специфические параметры для извлечения овощей, но параметры могут быть изменены в соответствии с другими потребностями. Я уверен, что этот процесс можно улучшить, так как это была моя первая серьезная попытка написания скриптов на python.
источник
Я отправляю это как ответ из-за ограничения длины в комментарии, никаких надежд на кредиты :). Очень широкая кисть, при условии, что у вас есть DEM.
Максимальное количество групп в процессе = количество деревьев внутри отдельного многоугольника. Дополнительные критерии, например, расстояние между «деревьями» внутри многоугольников, могут помочь ... Сглаживание матрицы высот с использованием ядра также вариант.
источник