Mam następujący proces automatycznego tworzenia konturów budynków / śladów na podstawie plików Las i zastanawiałem się, czy ktoś może pomóc w ulepszeniu go lub ma sugestie dotyczące innych procesów do wykorzystania?
DOŚWIADCZENIA:
A. ArcGIS - narzędzia Lidardataset : chociaż możesz zobaczyć tylko punkty klasy budynku, nie możesz automatycznie utworzyć wielokąta. Istnieje minimalna geometria graniczna, ale patrząc na zestaw danych jako całość, musisz ręcznie wybrać punkty dla każdego budynku.
B. Narzędzie zewnętrznej analizy Lidar firmy Overwatch - ekstrakcja budynku jest podstawową funkcją. http://www.featureanalyst.com/lidar_analyst.htm (Nie jest już aktywny - 28 lipca 2014 r.)
Utworzenie gołej ziemi (niezbędne) zajęło około 8 minut, a następnie wydobycie budynków. Jak dotąd daje najlepsze wyniki, nawet uchwycenie budynku pod drzewami (1 i 4). Nie trzeba nawet patrzeć na samouczek, ponieważ jest to bardzo intuicyjny interfejs.
Niektóre problemy polegają na tym, że ominął zadaszone korytarze (3), ale prawidłowo pominął otwarty dach (2).
C. LASTOOLS - używał lasmerge do łączenia 4 odpowiednich płytek, a następnie lasclip do wydobywania strefy AOI tylko.
arcpy.lasboundary("AOI.las","1","buildings","true","true","shp","AOI_buildings.shp","","#","true") 
Wyjście było bardzo postrzępione, więc użyłem:
arcpy.SimplifyBuilding_cartography("AOI","AOI_buildings_simplified.shp","4 Meters","0 SquareMeters","NO_CHECK")
a następnie pobiegł:
arcpy.SimplifyPolygon_cartography("AOI_buildings_simplified","AOI_s.shp","POINT_REMOVE","1 Meters","0 SquareMeters","NO_CHECK","KEEP_COLLAPSED_POINTS")
1 i 2 pokazują obszary problemów. Częściowo wynika to ze złej klasyfikacji lidar.
Daje to dobre linie budynków dla długich budynków o prostych krawędziach, ale jeśli gałęzie zachodzą na siebie lub jest to budynek złożony z wieloma krawędziami, dachami na wielu poziomach blisko siebie itp., Ręczne tworzenie jest niezbędne dla każdej szczegółowej pracy na poziomie witryny. W przypadku ogólnego oszacowania kształtu dachu proces jest wystarczająco dokładny. Proces trwa łącznie mniej niż 10 minut po zakończeniu procesu.
Moje badania : /programming/5037588/simplified-or-smooth-polygons-that-contain-the-original-detailed-polygon oraz proces sugerowany przez @gromix i Czym są definicja, algorytmy i praktyczne rozwiązania dla wklęsłego kadłuba?