LAS dosyalarından iki Dijital Yükseklik Modeli (DEM) karşılaştırılsın mı?

İki LiDAR dosyam (.las) var, biri orijinal X puanı ile diyelim. Diğeri ise ilk .las dosyasının kopyası, ancak Y'nin X'in altında olduğu Y noktaları ile.

Şimdi, bu iki .las dosyasının Dijital Yükseklik Modellerini (DEM) karşılaştırmak ve ne kadar farklı olduklarını görselleştirmek istiyorum.

Diğer karşılaştırma türleri arasında RMSE, standart sapma gibi bilgileri almak istiyorum.

Herkes bana hangi yazılımlar ve karşılaştırma bilgi almak için yollar söyleyebilirdi takdir ediyorum.

R'de iki Dijital Yükseklik Modeli (DEM) nasıl karşılaştırılı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.

ArcGIS for Desktop kullanılarak oluşturulan iki TIN'i karşılaştırma hakkında @ whuber'ın cevabına bakınız. bu konuyla ilgili teorik bir kavrayış için.

Bunu ESRI ArcGIS Geostatistics Analysis Extension üzerinden yapabilirsiniz - alt kümelerde doğrulama gerçekleştirme konusunda yardımda bir bölüm vardır .

Aynı şeyi R arayüzü üzerinden GRASS üzerinden de yapabilirsiniz . Tomislav Hengl bunu nasıl ayrıntılı olarak Jeoistatistiksel Haritalama için Pratik Kılavuz kitabında açıklayacağını açıklar . Açık erişime sahiptir, bu nedenle PDF'yi indirmek ücretsizdir.

Bildiğim kadarıyla, RMSE sadece DEM'in yapımı sırasında belirtilir, daha fazla referans için bir özellik olarak değil, bu yüzden DEM'in yapımı sırasında manuel olarak "yakalamanız" gerekir (yani, asla yapmadım Lidar'dan bir DEM, sadece diğer verilerden).

İki DEMS içindeki DATA arasındaki farkları görmek istiyorsanız, ArcGIS'in Mekansal Analist uzantısında bulunan ("Yüzey Analizi" altında) kesme / dolgu kullanacağım. Kes / doldur, basit bir tematik haritada DEM'ler arasındaki değişiklikleri gösterir.

Basit bir DEM farkı yapardım. DEM2-DEM1. Bu, tüm farklı alanları ve ne kadar olduğunu gösterecektir.

Web sitem ana sayfasında fark yüksek res dem bir görüntü var. thadwester.com
Sol renkli resme bir göz atın.