比较LAS文件中的两个数字高程模型（DEM）？

我有两个LiDAR文件（.las），其中一个是原始的X点。另一个是第一个.las文件的副本，但是具有Y点，其中Y小于X。

现在，我想比较这两个.las文件的数字高程模型（DEM），并可视化它们之间的差异。

我想获取诸如RMSE，标准偏差之类的信息，以及其他类型的比较。

如果有人能告诉我什么软件以及获取比较信息的方式，我将不胜感激。

如何在R中比较两个数字高程模型（DEM）。

#-------------------------------------------------------------------------
#Creating a reproducible example

library(raster)

  #simulating raster_1

  f = system.file("external/test.grd", package="raster")
  DEM_1 = raster(f)

  #simulating raster_2

  DEM_2 = DEM_1
  # replacing values from raster_1 to create a new raster sample (raster_2)
    DEM_2[(DEM_2>500 & DEM_2<900)] = 550
    DEM_2[(DEM_2>200 & DEM_2<300)] = 500

#-------------------------------------------------------------------------
# Comparison 1 (DEM_3 resulted from subtracting DEM_2 from DEM_1)

  DEM_3 = DEM_1 - DEM_2

    par(mfrow=c(1,3))

    plot(DEM_1, main = "DEM_1")
    plot(DEM_2, main = "DEM_2")
    plot(DEM_3, main = "DEM_3 = DEM_1 - DEM_2")

      dev.off()

#-------------------------------------------------------------------------
#Comparison 2 (histogram)

  hist(DEM_1, prob=T, main="DEM_1", xlab="")
  hist(DEM_2, prob=T, main="DEM_2", xlab="")
  hist(DEM_3, prob=T, main="DEM_3 = DEM_1 - DEM_2", xlab="")

    par(mfrow=c(1,1))

  standard_deviation = sd(c(as.matrix(DEM_3)),na.rm=T)

    dev.off()

#-------------------------------------------------------------------------
#comparison 3 (RMSE)

  library(hydroGOF)

  DEM_1_matrix = c(as.matrix(DEM_1))
  DEM_2_matrix = c(as.matrix(DEM_2))

  rmse = rmse(DEM_1_matrix,DEM_2_matrix)
  rmse
  [1] 135.3675 # this is the root mean squared error (RMSE) result.

请参阅@whuber关于比较使用ArcGIS for Desktop创建的两个TIN的答案？对于这个问题的理论见解。

您可以通过ESRI ArcGIS Geostatistical Analysis Extension来执行此操作-帮助中有一节对子集执行验证。

您可以通过R接口通过GRASS进行相同的操作。Tomislav Hengl在他的《地统计制图实用指南》一书中详细介绍了如何进行。它是开放式访问，因此可以免费下载PDF。

据我所知，RMSE仅在制作DEM时说明，而不是进一步提炼的属性，因此您必须在制作DEM时手动“抓住”它（也就是说，我从未做过来自激光雷达的DEM，仅来自其他数据）。

如果要查看两个DEMS内部数据之间的差异，可以使用ArcGIS的Spatial Analyst扩展（在“ Surface Analysis”下）中的cut / fill。切割/填充在简单的专题图中向您显示DEM之间的变化。

我将做一个简单的差异DEM。DEM2-DEM1。这将显示所有不同的区域以及数量不同。

我的网站首页上的图片差异很大。thadwester.com
看一下彩色的左图。