Bagaimana cara membandingkan dua DEM dan menunjukkan perbedaannya?

10

Saya ingin membandingkan DEM saya satu sama lain (ASTER GDEM2 dan SRTM4.1). Saya tahu, bahwa GDEM memiliki beberapa "lubang" dengan ketinggian yang salah (saya menemukan 3 lubang di area kerja saya). Sekarang saya perlu tahu, jika ada lebih banyak dari mereka untuk akhirnya memperbaikinya dengan SRTM. Saya menggunakan SAGA, GRASS dan QGIS.

PETA
sumber
membandingkannya dengan apa? satu sama lain? satu apakah Anda baseline dan yang lainnya salah?
Nickves
Baik! Saya hanya ingin menggambar perbedaan di antara mereka.
PETA

Jawaban:

5

Katakanlah Raster1 adalah yang benar, Dan Raster2 adalah yang Anda curigai memiliki kesalahan.

Mereka juga keduanya berada di unit yang sama (karena Anda mengatakan kedua raster adalah DEM nilai z harus dalam meter - saya berasumsi kita benar dan kita membandingkan hal yang sama)

Anda dapat menggunakan rast calc dari salah satu paket (qgis, rumput, saga) untuk membandingkannya misalnya:

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

setiap sel dari raster yang dihasilkan akan menampilkan persentase raster2 yang menyimpang dari raster1.

Setelah itu Anda dapat mengklasifikasikan ulang deviation_rast berdasarkan kebutuhan Anda: mis

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

untuk mengidentifikasi area yang bermasalah - dan untuk memilih apakah Anda ingin campur tangan atau tidak

Nickes
sumber
8
Rasio ini bukan ide yang baik untuk DEM dengan ketinggian di dekat permukaan laut: Anda akan mendapatkan banyak hasil palsu dan hasil yang sepenuhnya tidak terdefinisi dimanapun raster2 == 0. Pendekatan yang tepat, untuk semua DEM, hanya dengan mengurangi satu raster dari yang lain dan layar untuk pencilan (yang dapat dilakukan dalam banyak, banyak cara: lihat stats.stackexchange.com untuk beberapa metode).
whuber
11

Anda bisa mulai dengan mendapatkan perbedaan dari dua DEM. QGIS memiliki alat kalkulator raster yang berguna. Dapatkan ubin dari kedua DEM yang mencakup area yang sama dan kurangi nilai satu DEM dari yang lain. Itu seharusnya memberi Anda layer raster yang bagus yang menunjukkan perbedaan ketinggian antara dua DEM.

RK
sumber
5
Tepatnya, kurangi satu DEM dari yang lain. Saya kemudian akan mewakilinya di QGIS atau bahkan mengklasifikasi ulang, menggunakan toleransi di sekitar Zero. Sesuatu seperti Jika -5m <Hasil <5m maka pertimbangkan bahwa DEM memiliki nilai yang sama dalam sel tertentu.
Alexandre Neto
6
(+1) Itu awal yang baik - dan sederhana juga. Kesalahan dalam DEM sedikit rumit: mereka dihasilkan dari kesalahan elevasi dan dari kesalahan posisi. Jika tidak ada kesalahan posisi, maka biasanya layar yang didasarkan pada toleransi konstan adalah baik - tetapi seharusnya tidak harus sekitar nol, karena satu DEM mungkin memiliki (kecil) pergeseran vertikal konstan relatif terhadap yang lain. Gunakan toleransi di sekitar perbedaan median sebagai gantinya. Di daerah perbukitan, kesalahan posisi menambah lebih banyak variasi di mana lereng lebih curam, sehingga toleransi perlu ditingkatkan secara kasar sebanding dengan garis singgung lereng.
whuber
Saya memiliki masalah yang sama, tetapi raster saya tumpang tindih hanya sebagian (raster B hanya mencakup sebagian dari raster A) dan alat "perbedaan grid" tidak akan berfungsi. Ada ide yang bisa saya lakukan?
H.Wiener