Eu tenho o seguinte processo para criar automaticamente contornos / pegadas de construção a partir dos arquivos Las e fiquei imaginando se alguém pode ajudar a melhorá-lo ou tem sugestões para outros processos usarem?
EXPERIMENTOS:
A. ArcGIS - Ferramentas Lidardataset : enquanto você pode ver apenas os pontos da classe de construção, não é possível criar polígono automaticamente. Existe a geometria mínima do limite, mas, ao examinar o conjunto de dados como um todo, é necessário escolher manualmente os pontos para cada construção.
B. Ferramenta de análise Lidar externa da Overwatch - a extração de edifícios é uma funcionalidade essencial. http://www.featureanalyst.com/lidar_analyst.htm (Não está mais ativo - 28 de julho de 2014)
Foram necessários aproximadamente 8 minutos para criar uma terra nua (essencial) e depois extrair os edifícios. Dá os melhores resultados de longe, até capturando edifícios sob as árvores (1 e 4). Isso sem sequer olhar para um tutorial, pois é uma interface muito intuitiva.
Alguns problemas são que ele perdeu os corredores cobertos (3), mas deixou de fora o teto aberto (2).
C. LASTOOLS - lasmerge usado para mesclar 4 peças relevantes e, em seguida, clipe para extrair apenas a zona da AOI.
arcpy.lasboundary("AOI.las","1","buildings","true","true","shp","AOI_buildings.shp","","#","true") 
A saída foi muito irregular, então eu usei:
arcpy.SimplifyBuilding_cartography("AOI","AOI_buildings_simplified.shp","4 Meters","0 SquareMeters","NO_CHECK")
e depois correu:
arcpy.SimplifyPolygon_cartography("AOI_buildings_simplified","AOI_s.shp","POINT_REMOVE","1 Meters","0 SquareMeters","NO_CHECK","KEEP_COLLAPSED_POINTS")
1 e 2 mostram as áreas temáticas. Parte disso é devido à classificação ruim do Lidar.
Isso fornece boas linhas de construção para edifícios com arestas longas e retas, mas se houver sobreposição de galhos ou for um edifício complexo com várias arestas, telhados em vários níveis próximos uns dos outros etc, a criação manual é necessária para qualquer trabalho detalhado no nível do local. Para uma estimativa geral da forma do telhado, o processo é suficientemente preciso. O processo leva menos de 10 minutos no total após a finalização do processo.
Minha pesquisa : /programming/5037588/simplified-or-smooth-polygons-content-the-original-detailed-polygon e o processo sugerido por @gromix e o que são definição, algoritmos e soluções práticas para casco côncavo?