Tengo el siguiente proceso para crear automáticamente contornos / huellas de edificios a partir de archivos Las y me preguntaba si alguien puede ayudar a mejorarlo o tiene sugerencias para otros procesos.
Experimentos
A. ArcGIS: herramientas Lidardataset : aunque puede ver solo los puntos de clase de construcción, no puede crear polígonos automáticamente. Existe la geometría de límite mínima, pero como se ve en el conjunto de datos en su conjunto, requiere que elija manualmente los puntos para cada edificio.
B. Herramienta de análisis externo Lidar de Overwatch : la extracción de edificios es una funcionalidad central. http://www.featureanalyst.com/lidar_analyst.htm (Ya no está activo - 28 de julio de 2014)
Le tomó aproximadamente 8 minutos crear una tierra desnuda (esencial) y luego extraer los edificios. Con mucho, ofrece los mejores resultados, incluso capturando edificios debajo de los árboles (1 y 4). Esto es sin siquiera mirar un tutorial ya que es una interfaz muy intuitiva.
Algunos problemas son que se ha perdido los pasillos cubiertos (3) pero que se ha salido correctamente del techo abierto (2).
C. LASTOOLS - usó lasmerge para fusionar 4 fichas relevantes y luego lasclip para extraer la zona de AOI solamente.
arcpy.lasboundary("AOI.las","1","buildings","true","true","shp","AOI_buildings.shp","","#","true") 
La salida fue muy irregular, así que usé:
arcpy.SimplifyBuilding_cartography("AOI","AOI_buildings_simplified.shp","4 Meters","0 SquareMeters","NO_CHECK")
y luego corrió:
arcpy.SimplifyPolygon_cartography("AOI_buildings_simplified","AOI_s.shp","POINT_REMOVE","1 Meters","0 SquareMeters","NO_CHECK","KEEP_COLLAPSED_POINTS")
1 y 2 muestran las áreas problemáticas. Parte de esto se debe a una mala clasificación de LIDAR.
Esto proporciona buenas líneas de construcción para edificios largos con bordes rectos, pero si hay una superposición de ramas o es un edificio complejo con múltiples bordes, techos en múltiples niveles cerca uno del otro, etc., la creación manual es necesaria para cualquier trabajo detallado a nivel del sitio. Para la estimación general de la forma del techo, el proceso es lo suficientemente preciso. El proceso tarda menos de 10 minutos en total una vez que el proceso ha finalizado.
Mi investigación : /programming/5037588/simplified-or-smooth-polygons-that-contain-the-original-detailed-polygon y el proceso sugerido por @gromix y What are Definition, Algorithms and Practical Solutions para el casco cóncavo?