İki DEM nasıl karşılaştırılır ve farkları nasıl gösterilir?

10

DEM'lerimi birbirleriyle karşılaştırmak istiyorum (ASTER GDEM2 ve SRTM4.1). Biliyorum, GDEM yanlış yükseklik ile bazı "delikleri" olduğunu biliyorum (3 tanesi çalışma alanımda buldum). Şimdi bilmeliyim ki, sonunda SRTM ile düzeltmek için daha fazlası varsa. SAGA, GRASS ve QGIS kullanıyorum.

— MAP
kaynak

onları ne ile karşılaştır? birbirleriyle? biri temel ve diğeri hatalı?

— nickves

sağ! Sadece aralarındaki farkı çizmek istiyorum.

— HARİTA

5

Diyelim ki Raster1 doğru, Ve Raster2 hatalarından şüphelendiğiniz.

İkisi de aynı birimlerde (her iki raster de DEM olduğunu söylediğiniz için z değeri metre cinsinden olmalıdır - doğru olduğumuzu ve aynı şeyleri karşılaştırdığımızı varsayıyorum)

Bunları karşılaştırmak için her iki paketten (qgis, grass, saga) rast calc kullanabilirsiniz:

deviation_rast  = abs( (raster1 - raster2) / raster1 * 100 )

sonuçta elde edilen taramadan her hücre, tarama2'nin taramadan sapan yüzdesini gösterecektir.

Bundan sonra, ihtiyaçlarınıza göre deviation_rast'ı yeniden sınıflandırabilirsiniz: örn.

red -> deviation_rast_value > 75
orange -> deviation_rast_value > 50
yellow ->  deviation_rast_value > 25
green -> deviation_rast_value  > 0

sorunlu alanları belirlemek - ve müdahale etmek isteyip istemediğinizi seçmek için

— nickves
kaynak

8

Oran, deniz seviyesine yakın yüksekliklere sahip DEM'ler için iyi bir fikir değildir : raster2 == 0 olduğunda birçok sahte sonuç ve tamamen tanımlanmamış sonuçlar elde edersiniz. Tüm DEM'ler için doğru yaklaşım, bir raster diğerinden ve ekrandan çıkarmaktır. aykırı değerler için (bu birçok yolla yapılabilir: bazı yöntemler için stats.stackexchange.com adresine bakın ).

— whuber

11

İki DEM'in farkını alarak başlayabilirsiniz. QGIS, kullanışlı olması gereken bir raster hesaplama aracına sahiptir. Sadece aynı alanı kapsayan her iki DEM'in döşemelerini alın ve bir DEM'in değerlerini diğerinden çıkarın. Bu, iki DEM arasındaki yükseklik farklarını gösteren güzel bir raster katmanı elde etmelidir.

— RK
kaynak

5

Kesinlikle, sadece bir DEM'i diğerinden çıkarın. Daha sonra onu QGIS'de temsil edebilirim, hatta Sıfır etrafında hoşgörü kullanarak yeniden sınıflandırabilirim. Eğer -5m <Sonuç <5m gibi bir şey DEM'lerin o hücrede aynı değere sahip olduğunu düşünün.

— Alexandre Neto

6

(+1) Bu iyi bir başlangıç - ve basit. DEM'lerde hatalar biraz karmaşıktır: yükseklik hatalarından ve konumsal hatalardan kaynaklanır. Konum hatası yoksa, genellikle sabit bir toleransa dayalı bir ekran iyidir - ancak mutlaka sıfır civarında olmamalıdır, çünkü bir DEM diğerine göre (küçük) sabit dikey kaymaya sahip olabilir. Bunun yerine medyan farkının etrafında bir tolerans kullanın. Tepelik alanlarda, konum hatası eğimlerin daha dik olduğu yerlerde daha fazla varyasyon ekler, bu nedenle toleransın eğimin teğetiyle orantılı olarak kabaca artırılması gerekir.

— whuber

Benzer bir sorunum var, ancak rasterlerim sadece kısmen çakışıyor (raster B raster A'nın sadece bir kısmını kaplıyor) ve "ızgara farkı" aracı çalışmaz. Ne yapabileceğime dair bir fikrin var mı?

— H.Wiener