Come confrontare due DEM e mostrare la loro differenza?

10

Voglio confrontare i miei DEM tra loro (ASTER GDEM2 e SRTM4.1). So che il GDEM ha alcuni "buchi" con elevazioni errate (ne ho trovati 3 nella mia area di lavoro). Ora devo sapere se ce ne sono altri per correggerli alla fine con SRTM. Uso SAGA, GRASS e QGIS.

— CARTA GEOGRAFICA
fonte

confrontarli con cosa? insieme? uno è la linea di base e l'altro è difettoso?

— Nickves,

destra! Voglio solo disegnare la differenza tra loro.

— MAPPA

5

Diciamo che Raster1 è quello corretto, e Raster2 è quello che sospetti abbia dei difetti.

Entrambi sono anche nelle stesse unità (dal momento che dici che entrambi i raster sono DEM, il valore z dovrebbe essere in metri - presumo che siamo corretti e confrontiamo le stesse cose)

Puoi usare rast calc da entrambi i pacchetti (qgis, grass, saga) per confrontarli, ad esempio:

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

ogni cella dal raster risultante visualizzerà la percentuale che raster2 devia dal raster1.

Successivamente è possibile riclassificare deviation_rast in base alle proprie esigenze: ad es

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

identificare le aree problematiche e scegliere se si desidera intervenire o meno

— nickves
fonte

8

Il rapporto non è una buona idea per i DEM con elevazioni vicino al livello del mare: otterrai molti risultati fasulli e risultati completamente indefiniti ovunque raster2 == 0. L'approccio giusto, per tutti i DEM, è solo quello di sottrarre un raster dall'altro e schermare per outlier (che può essere fatto in molti, molti modi: consultare stats.stackexchange.com per alcuni metodi).

— whuber

11

Potresti iniziare ottenendo la differenza tra i due DEM. QGIS ha uno strumento calcolatrice raster che dovrebbe tornare utile. Basta ottenere tessere di entrambi i DEM che coprono la stessa area e sottrarre i valori di un DEM dall'altro. Questo dovrebbe darti un bel livello raster che mostra le differenze di elevazione tra i due DEM.

— RK
fonte

5

Esatto, sottrai solo un DEM dall'altro. Lo rappresenterei quindi in QGIS o lo riclassificherei, usando una certa tolleranza intorno a Zero. Qualcosa come Se -5m <Risultato <5m, allora considera che i DEM hanno lo stesso valore in quella particolare cella.

— Alexandre Neto,

6

(+1) È un buon inizio, e anche semplice. Gli errori nei DEM sono un po 'complicati: derivano da errori di elevazione e da errori di posizione. Se non si verifica alcun errore di posizione, in genere una schermata basata su una tolleranza costante è buona, ma non dovrebbe necessariamente essere intorno allo zero, poiché un DEM potrebbe avere uno spostamento verticale (piccolo) costante rispetto a un altro. Utilizzare invece una tolleranza attorno alla differenza mediana. Nelle aree collinari, l'errore di posizione aggiunge ulteriori variazioni nei casi in cui le pendenze sono più ripide, quindi la tolleranza deve essere aumentata approssimativamente in proporzione alla tangente della pendenza.

— whuber

Ho un problema simile, ma i miei raster si sovrappongono solo parzialmente (il raster B copre solo una parte del raster A) e lo strumento "differenza griglia" non funzionerà. Qualche idea su cosa posso fare?

— H.Wiener il