Ho il seguente processo per creare automaticamente contorni / impronte di edifici dai file Las e mi chiedevo se qualcuno potesse aiutarlo a migliorarlo o avere suggerimenti per altri processi da utilizzare?
ESPERIMENTI:
A. ArcGIS - Strumenti Lidardataset : mentre puoi vedere solo i punti della classe di costruzione, non puoi creare automaticamente un poligono. Esiste la geometria delimitazione minima, ma poiché osserva l'insieme di dati nel suo insieme, è necessario scegliere manualmente i punti per ciascun edificio.
B. Strumento di analisi del lidar esterno di Overwatch : l'estrazione degli edifici è una funzionalità fondamentale. http://www.featureanalyst.com/lidar_analyst.htm (Non più attivo - 28 luglio 2014)
Ci sono voluti circa 8 minuti per creare una terra nuda (essenziale) e quindi estrarre gli edifici. Fornisce di gran lunga i migliori risultati, persino catturando l'edificio sotto gli alberi (1 e 4). Questo senza nemmeno guardare un tutorial in quanto è un'interfaccia molto intuitiva.
Alcuni problemi sono che ha perso i corridoi coperti (3) ma ha correttamente lasciato fuori il tetto aperto (2).
C. LASTOOLS: utilizzato il lasmerge per unire 4 tessere pertinenti e quindi scorrere per estrarre solo la zona di AOI.
arcpy.lasboundary("AOI.las","1","buildings","true","true","shp","AOI_buildings.shp","","#","true")
L'output era molto frastagliato, quindi ho usato:
arcpy.SimplifyBuilding_cartography("AOI","AOI_buildings_simplified.shp","4 Meters","0 SquareMeters","NO_CHECK")
e poi corse:
arcpy.SimplifyPolygon_cartography("AOI_buildings_simplified","AOI_s.shp","POINT_REMOVE","1 Meters","0 SquareMeters","NO_CHECK","KEEP_COLLAPSED_POINTS")
1 e 2 mostrano le aree problematiche. Parte di ciò è dovuta a una cattiva classificazione del lidar.
Ciò offre buone linee di costruzione per edifici con bordi dritti lunghi ma se vi è una sovrapposizione di rami o è un edificio complesso con più bordi, tetti su più livelli vicini tra loro ecc., La creazione manuale è necessaria per qualsiasi lavoro dettagliato a livello di sito. Per la stima generale della forma del tetto il processo è abbastanza accurato. Il processo impiega meno di 10 minuti in totale una volta finalizzato.
La mia ricerca : /programming/5037588/simplified-or-smooth-polygons-that-contain-the-original-detailed-polygon e il processo suggerito da @gromix e What are Definition, Algorithms and Practical Solutions per lo scafo concavo?