Kentsel LIDAR'dan temizleme blokları deseni (sokaklardan yükseltilmiş bloklar)

Biz var 1 metrelik LİDAR DEM bir şehirden.

Bu bağlantıdan küçük bir alt küme indirilebilir :

Bu ekran görüntüsü gri paletli ham DEM'i gösterir (daha koyu kayışlar sokaklardır ve grimsi ve beyazımsı dikdörtgenler bloklardır):

Bu, Santo Domingo şehrinde, bu Google ekran görüntüsünde görülebilen bir yere karşılık gelir:

Ortalama olarak, bloklar "yükseltilmiş" ca. Sokaklardan 2 metre , bu doğru değil. Akış ağı ve topografik ıslaklık indeksi (TWI) oluşturmak için temiz bir DEM'e sahip olmak istiyoruz. DEM tedarik edildiğinde (lazer tarayıcıdan orijinal bantlarımız yok), hidrolik ağ dikdörtgen bir düzen izliyor gibi görünüyordu ve TWI blok deseniyle sonuçlandı . Bu resimler sonuçları gösterir:

Bu oluşturulacak akış ağ sonucudur, r.watershediçinde Grass GIS:

Ve bu, aşağıdakilerle oluşturulan TWI sonucudur SAGA:

Bu yanlışlığı başarılı bir şekilde çözmek için bazı prosedürleri denedik :

1) Denoising aracı . r.denoiseAracı uyguladık Grass GIS, ancak modül kurulumuyla ilgili bazı sorunlar yaşadık . Windows'ta bir kabukla tekrar çalıştırdık ve yetersiz bir bellek mesajı aldık.

2) Filtreler . Biz filtreler farklı pencereler boyutları ile (ortalama düşük-pass, medyan, vs.) farklı koştu ve sokakların yönünde ağırlığını koymak için çalışıyor ( Grass GIS, SAGA, QGIS).

3) Jeoistatistik . Nokta bulutlarını kesinlikle sokaklarda ürettik (1000 ve 2000 puanları denedik), bir variogram modeli oluşturduk ve daha sonra blokları doldurmak için sıradan bir krigging yaptık. Variogram modelleme ve sıradan krigging, Rfarklı paketler kullanılarak yapıldı. Doğrusal bir variogramımız var, bu yüzden kriging sonuçlarına güvenmeyiz.

4) Diğer araçlar . Yüklü ALDPATaraç, ancak program DEM'yi okuyamadığı için çalışamadı.

Her durumda, drenaj ağı açısından sonuçlar iyi değildi , çünkü dikdörtgen akış ağından kaçınamadık; ayrıca, TWI hala bir blok paterni ile sonuçlandı .

Özellikle, Tamam enterpolasyonlu sonuç ile, ağ sonucunu etkileyen nokta benzeri bir model DEM elde ettik. Bununla birlikte, blok paterninin etkisi azalmıştır.

Ayrıca, bu soru ve cevaplara bir göz attık ...

Çıplak toprak yükselmesi için kanopileri ve binaları DSM'den filtreleme

... bizi yönlendirdi Whitebox Geospatial Analysis Tools, ancak DEM'imizi dönüştüremedi LAS format. Ayrıca, bizim için etkinliğinden emin değildik Bare-Earth DEM tool, çünkü yarı saydam nesneleri kaldırmak için tasarlandı, yanlış "yükseltilmiş" blokları değil, bizim durumumuz.

Hidrografik analizimizi yapmak için hala yüksek kaliteli bir DEM oluşturmak istiyoruz, ancak başka ne deneyebileceğimizi bilmiyoruz.

Sorununuzu çözmek için alternatif yazılım kullanmaya açıksanız, platformlar arası açık kaynaklı GIS Whitebox Jeo - Uzamsal Analiz Araçlarının Arazi Nesnelerini Kaldır aracını önerebilirim . Sorunuzda verilerinizi LAS formatına dönüştüremediğinizi söylediğinizin farkındayım, ancak araç giriş olarak LAS dosyası değil, bir raster alıyor. Bu aracı belki de bir LAS girişi alan Bare-Earth DEM aracıyla karıştırıyorsunuz . Raster DEM'inizi Whitebox'a bir GeoTIFF dosyası veya bir ArcGIS kayan noktalı ikili raster (.flt) veya herhangi bir sayıda yaygın raster formatı olarak aktarabilirsiniz.

Bir raster DEM'den binaları kaldırma yeteneğinin başka bir örneği:

Önemli olarak, algoritma bir filtre değildir ve bu nedenle ızgara hücrelerinin 'arazi dışı nesneler' (binalar) dışındaki yüksekliklerini etkilemeyecektir. Islaklık DEM'inizi ıslaklık endeksini hesaplamak için kullanmak istiyorsanız bu çok önemlidir. (Bununla birlikte, kentsel veya yoğun bir banliyö lokasyonunun ıslaklık endeksini tahmin etmenin yararlılığını sorgulayacağım. Ayrıca, görüntünüzde gösterdiğiniz yoğun kentsel manzara boyunca akış akışı yok ... akarsu menfezleri oluşturuldu.)

DÜZENLE

Aslında, yüklediğiniz örnek veri kümesine bakarak, DEM'inizin Arazi Nesnelerini Kaldır aracıyla kullanılmaya uygun olup olmadığını bilmiyorum. Yüklediğiniz görüntünün sadece zayıf sembolojiden muzdarip olduğunu düşünmüştüm, ama şimdi aslında teraslandığını görüyorum (yani DEM içinde geniş düz alanlar adımları var). DEM'inizden aşağıdaki profile bakın:

Araç doğru eğim ölçümlerine dayanır (ki bu da DEM'inizle mümkün olduğunu düşünmüyorum ... TWI'yı da hesaplamak istiyorsanız bu işaretlenmelidir) ve zemin yüzeyi ile zemin arasındaki keskin eğim sınırlarının korunması binalar. Ancak sizin durumunuzda DEM ciddi şekilde yumuşatılmıştır ve bu keskin kenarlar o kadar belirgin değildir (evleriniz de yaklaşık 2-3 m yüksekliğinde, biraz garip görünüyor). Ellerinizi orijinal işlenmemiş LiDAR DEM'e veya daha iyisi nokta bulutu verilerine alabilir misiniz? Başvurunuz için, DEM'i yeniden enterpolasyonu ciddiye almayı düşünürüm.

Orijinal LAS nokta bulutu olmadan, raster manipülasyonu yoluyla verilere yalnızca daha fazla yanlışlık ekleyeceğinizi söylemeliyim. Sağlanan DEM, yoğun bir şehirleşmiş 1m çözünürlüklü DEM için nispeten temiz görünüyor. "Yükseltilmiş kareler", binaların son zemin sınıfına dahil edilmediği veri boşlukları arasındaki üçgenlemelerin bir sonucudur. Ayrıca suyun, gerçek modelinizde olduğu gibi gerçek dünyada çıplak bir toprak DEM'den akmayacağını da düşünebilirsiniz. Binalar gerçek hidrografik süreçlere sınırlar sunacaktır.