Comparaison de deux modèles numériques d'élévation (DEM) à partir de fichiers LAS?

J'ai deux fichiers LiDAR (.las), l'un est original disons avec X points. Et l'autre est une copie du premier fichier .las mais avec des points Y, où Y est inférieur à X.

Maintenant, je veux comparer les modèles numériques d'élévation (DEM) de ces deux fichiers .las et visualiser à quel point ils sont différents.

Je veux obtenir des informations comme RMSE, l'écart-type, entre autres types de comparaison.

J'apprécierais, si quelqu'un pouvait me dire quels logiciels et comment obtenir les informations de comparaison.

Comment comparer deux modèles numériques d'élévation (DEM) dans R.

#-------------------------------------------------------------------------
#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.

Voir la réponse de @ whuber sur la comparaison de deux TIN créés à l'aide d'ArcGIS for Desktop? pour un aperçu théorique de cette question.

Vous pouvez le faire via l'extension d'analyse géostatistique d'ESRI ArcGIS - il y a une section dans l'aide pour effectuer la validation sur des sous-ensembles .

Vous pouvez faire la même chose par GRASS à travers l' interface R . Tomislav Hengl décrit en détail comment procéder dans son livre A Practical Guide to Geostatistical Mapping . C'est un accès ouvert, donc le PDF est téléchargeable gratuitement.

Pour autant que je sache, RMSE n'est indiqué que lors de la fabrication du DEM, et non comme un attribut pour une réfraction ultérieure, vous devez donc "l'attraper" manuellement pendant la fabrication du DEM (cela dit, je n'ai jamais fait un DEM de Lidar, uniquement à partir d'autres données).

Si vous voulez voir les différences entre les DONNÉES à l'intérieur des deux DEMS, j'utiliserais couper / remplir qui se trouve dans l'extension Spatial Analyst d'ArcGIS (sous "Analyse de surface"). Le cut / fill vous montre dans une simple carte thématique les changements entre les DEM.

Je ferais un simple DEM de différence. DEM2-DEM1. Cela montrera tous les domaines qui sont différents et de combien.

Theres une image à une haute résolution dem de différence sur la page d'accueil de mon site Web. thadwester.com
Jetez un œil à l'image colorée de gauche.