LAS 파일에서 두 개의 디지털 고도 모델 (DEM)을 비교하고 있습니까?

나는 두 개의 LiDAR 파일 (.las)을 가지고 있는데, 하나는 X 포인트로 독창적입니다. 그리고 다른 하나는 첫 번째 .las 파일의 사본이지만 Y 지점은 Y이며 X보다 작습니다.

이제이 두 .las 파일의 DEM (Digital Elevation Model)을 비교하고 그 차이를 시각화하려고합니다.

다른 유형의 비교 중에서 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.

ArcGIS for Desktop을 사용하여 만든 두 개의 TIN 비교에 대한 @whuber의 답변을 참조하십시오 . 이 문제에 대한 이론적 통찰력.

ESRI ArcGIS Geostatistical Analysis Extension을 통해이 작업을 수행 할 수 있습니다 . 하위 집합에 대한 유효성 검사 수행에 대한 섹션이 있습니다.

R 인터페이스 를 통해 GRASS를 통해 동일한 작업을 수행 할 수 있습니다 . Tomislav Hengl은 자신의 저서 " 지형 통계 매핑에 대한 실용 안내서 "에서 그렇게하는 방법을 자세히 설명합니다 . 공개 액세스이므로 PDF를 무료로 다운로드 할 수 있습니다.

내가 아는 한, RMSE는 DEM 작성 중에 만 언급되며 추가 굴절의 속성으로 표시되지 않으므로 DEM 작성 중에 수동으로 "잡아"야합니다. 다른 데이터에서만 Lidar의 DEM).

두 DEMS 내부의 데이터 간 차이를 보려면 ArcGIS의 Spatial Analyst 확장 ( "Surface Analysis"아래)에있는 cut / fill을 사용합니다. 컷 / 필은 간단한 주제별 맵에서 DEM 간의 변경 사항을 보여줍니다.

나는 단순한 차이의 DEM을 할 것입니다. DEM2-DEM1. 이것은 다르고 얼마만큼 다른 모든 영역을 보여줍니다.

내 웹 사이트 홈페이지에 차이가 큰 이미지가 있습니다. thadwester.com
다채로운 왼쪽 이미지를 살펴보십시오.