2つのDEMを比較してその違いを示す方法は?


10

DEMを相互に比較したい(ASTER GDEM2およびSRTM4.1)。GDEMには標高の間違った「穴」がいくつかあることを知っています(作業エリアに3つ見つかりました)。SRTMで最終的にそれらを修正するためにそれらのより多くがある場合、私は今知っている必要があります。SAGA、GRASS、QGISを使用しています。


それらを何と比較しますか?お互いに?1つはベースラインで、もう1つは不良ですか。
ニックス

正しい!それらの違いを描きたいだけです。
MAP

回答:


5

Raster1が正しいとしましょう。Raster2は障害があると思われるものです。

また、両方とも同じ単位にあります(両方のラスターがDEMであると言うため、z値はメートル単位である必要があります-私たちは正しいと仮定し、同じものを比較します)

いずれかのパッケージ(qgis、grass、saga)のrast calcを使用して、それらを比較できます。例:

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

結果のラスタ​​の各セルには、raster2がraster1から逸脱している割合が表示されます。

その後、あなたのニーズに基づいて偏差_分類を再分類することができます:例

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

問題のある領域を特定し、介入するかどうかを選択する


8
標高が海抜に近いDEMの場合、この比率はお勧めできません。raster2== 0の場所では、多くの偽の結果と完全に未定義の結果が得られます。すべてのDEMの正しいアプローチは、他のラスターと画面から1つのラスターを差し引くことです。外れ値(多くの方法で実行できます。いくつかの方法については、stats.stackexchange.comを参照してください)。
whuber

11

最初に、2つのDEMの違いを取得します。QGISには便利なラスター計算ツールがあります。同じ領域をカバーする両方のDEMのタイルを取得し、一方のDEMの値をもう一方から差し引くだけです。これにより、2つのDEM間の標高の違いを示す素敵なラスターレイヤーが得られます。


5
正確には、一方のDEMをもう一方から差し引くだけです。次に、それをQGISで表すか、ゼロ周辺の許容誤差を使用して再分類します。If -5m <Result <5mのような場合は、DEMがその特定のセルで同じ値を持っていると考えてください。
Alexandre Neto 2012

6
(+1)それは良い出発点です-そして簡単です。DEMのエラーは少し複雑です。高度エラーと位置エラーが原因です。位置誤差がない場合は、通常、一定の許容誤差に基づく画面で問題ありませんが、あるDEMが別のDEMに対して(小さな)一定の垂直方向のシフトを持っている可能性があるため、必ずしも0である必要はありません。代わりに中央値の差の周りの許容誤差を使用してください。起伏のある地域では、傾斜が急な場所では位置誤差によりさらに変動が追加されるため、許容差は、傾斜の正接にほぼ比例して増加する必要があります。
whuber

同様の問題がありますが、ラスターが部分的にしかオーバーラップしておらず(ラスターBはラスターAの一部のみをカバーしています)、「グリッド差」ツールは機能しません。私にできることはありますか?
H.Wiener 2015年
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.