LASファイルの2つのデジタル標高モデル（DEM）を比較しますか？

12

2つのLiDARファイル（.las）がありますが、1つは元のXポイントです。もう1つは、最初の.lasファイルのコピーですが、Yポイントがあります（YはXより小さい）。

次に、これら2つの.lasファイルのデジタル標高モデル（DEM）を比較し、それらの違いを視覚化します。

他の種類の比較の中で、RMSE、標準偏差などの情報を取得したい。

誰がどのソフトウェアを比較できるか、比較情報を取得する方法を教えていただければ幸いです。

— アンドレ・シルバ
ソース

1

おそらく、XとYの代わりにMとNを使用するように質問を編集できます。最初に読んだとき、私はXとYが座標値だと思っていました！

— マークアイルランド

1

関連するヘルプを受け取るためには、さらに情報を提供する必要があります。あなたの現在の質問では、シーンはほとんどありません。lasファイルのどのフィールドが値を保持していますか。分類された地面がlas形式で割り当てられる方法は、分類フィールドであり、異なるz（標高）値ではありません。ベンダーは、z値の違いを保持するために未割り当てフィールドを使用する必要があります。

— ジェフリーエヴァンス

8

Rの2つのデジタル標高モデル（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を使用して作成された2つのTINの比較に関する@whuberの回答を参照してください。この問題についての理論的洞察。

— アンドレ・シルバ
ソース

5

これは、ESRI ArcGIS Geostatistical Analysis Extensionを使用して実行できます。サブセットの検証の実行に関するヘルプのセクションがあります。

Rインターフェースを介したGRASSでも同じことができます。トミスラフ・ヘングルは、彼の著書「地球統計学的マッピングの実践ガイド」で、その方法を詳細に説明しています。オープンアクセスなので、PDFは無料でダウンロードできます。

— om_henners
ソース

2

私の知る限り、RMSEはDEMの作成中にのみ記載されており、今後の参考資料としてではなく、DEMの作成中に手動で「キャッチ」する必要があります（つまり、 LidarのDEM、他のデータのみ）。

2つのDEMS内のデータの違いを確認したい場合は、ArcGISのSpatial Analystエクステンション（「表面分析」の下）にある切り取り/塗りつぶしを使用します。切り取り/塗りつぶしは、DEM間の変更を単純な主題図で示します。

— ジョナトル
ソース

3

カット/フィルは、差異を定量化しないため、あまりにも粗雑です。RMSEは、2つのデータセットを比較する一般的な方法です。DEMを地上の真実のデータと比較するだけでなく、2つのDEMの違いを定量化する1つの方法です。

— whuber

2

単純な差分のDEMを実行します。DEM2-DEM1。これにより、異なるすべての領域とその量が表示されます。

私のウェブサイトのホームページには、高解像度の違いがあるイメージがあります。thadwester.com
カラフルな左の画像をご覧ください。

— サッド
ソース