如何比较两个DEM并显示它们的差异?


10

我想相互比较我的DEM(ASTER GDEM2和SRTM4.1)。我知道GDEM的标高有一些“洞”(我在工作区域发现了3个)。现在,我需要知道,是否还有更多的产品可以最终通过SRTM进行更正。我使用SAGA,GRASS和QGIS。


将它们与什么进行比较?彼此?一个是您的基准,另一个是错误的?
尼克2012年

对!我只想画出它们之间的区别。
MAP 2012年

Answers:


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,正确的方法是从另一个栅格中减去一个栅格并进行筛选离群值(可以通过许多方法完成:有关某些方法,请参阅stats.stackexchange.com)。
ub

11

您可以先获取两个DEM的差异。QGIS有一个应该很方便的栅格计算器工具。只需获得覆盖相同区域的两个DEM的切片,然后从另一个DEM中减去一个DEM的值即可。那应该给您一个很好的栅格图层,该栅格图层显示两个DEM之间的高程差异。


5
确实,只需从另一个中减去一个DEM。然后,我将在QGIS中表示它,甚至使用零附近的一些容忍度对其进行重新分类。如果-5m <结果<5m,则认为该特定单元格中的DEM具有相同的值。
亚历山大·内托

6
(+1)这是一个好的开始-而且也很简单。DEM中的错误有点复杂:它们是由高程错误和位置错误引起的。如果没有位置误差,则通常基于恒定公差的屏幕是好的-但它不一定必须约为零,因为一个DEM相对于另一个DEM可能有一个(较小的)恒定垂直偏移。请在中间差附近使用公差。在丘陵地区,在坡度较陡的地方,位置误差会增加更多的变化,因此公差需要与坡度的切线成比例地增加。
ub

我有一个类似的问题,但是我的栅格仅部分重叠(栅格B仅覆盖栅格A的一部分),并且“栅格差异”工具将不起作用。有什么想法我能做什么?
H.Wiener
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.