У меня есть следующий процесс для автоматического создания контуров / контуров зданий из файлов Las, и мне было интересно, если кто-нибудь может помочь улучшить его или есть предложения для других процессов?
ЭКСПЕРИМЕНТ:
A. Инструменты ArcGIS - Lidardataset : хотя вы можете видеть только точки класса здания, вы не можете создавать многоугольник автоматически. Существует минимальная ограничивающая геометрия, но, поскольку она смотрит на набор данных в целом, требуется вручную выбирать точки для каждого здания.
B. Инструмент внешнего лидарного анализа от Overwatch - извлечение здания является основной функцией. http://www.featureanalyst.com/lidar_analyst.htm (Больше не активен - 28 июля 2014 г.)
Потребовалось около 8 минут, чтобы создать голую землю (обязательно), а затем извлечь здания. Это дает лучшие результаты, даже захватывая здания под деревьями (1 и 4). Это даже не смотря на учебник, поскольку это очень интуитивно понятный интерфейс.
Некоторые проблемы: он пропустил крытые коридоры (3), но правильно пропустил открытую крышу (2).
C. LASTOOLS - используйте lasmerge, чтобы объединить 4 соответствующих тайла, а затем отскок, чтобы извлечь только зону AOI.
arcpy.lasboundary("AOI.las","1","buildings","true","true","shp","AOI_buildings.shp","","#","true") 
Вывод был очень неровным, поэтому я использовал:
arcpy.SimplifyBuilding_cartography("AOI","AOI_buildings_simplified.shp","4 Meters","0 SquareMeters","NO_CHECK")
а потом побежал:
arcpy.SimplifyPolygon_cartography("AOI_buildings_simplified","AOI_s.shp","POINT_REMOVE","1 Meters","0 SquareMeters","NO_CHECK","KEEP_COLLAPSED_POINTS")
1 и 2 показывают проблемные области. Отчасти это связано с плохой лидарной классификацией.
Это дает хорошие строительные линии для длинных зданий с прямыми краями, но если они перекрываются ветвями или это сложное здание с несколькими краями, крышами на нескольких уровнях близко друг к другу и т. Д., Ручное создание необходимо для любой подробной работы на уровне участка. Для общей оценки формы крыши процесс достаточно точен. Процесс занимает менее 10 минут после завершения процесса.
Мои исследования : /programming/5037588/simplified-or-smooth-polygons-that-contain-the-original-detailed-polygon и процесс, предложенный @gromix, и что такое определение, алгоритмы и практические решения для вогнутого корпуса?