एलएएस फाइलों से दो डिजिटल एलिवेशन मॉडल (डीईएम) की तुलना करें?

Question

मेरे पास दो LiDAR फाइलें (.las) हैं, एक मूल है जो X बिंदुओं के साथ कहता है। और दूसरी पहली .las फ़ाइल की प्रतिलिपि है, लेकिन Y बिंदुओं के साथ, जहां Y X से कम है।

अब, मैं इन दोनों .las फ़ाइलों के डिजिटल एलिवेशन मॉडल (DEMs) की तुलना करना चाहता हूं और कल्पना करता हूं कि वे कितने अलग हैं।

मैं आरएमएसई, मानक विचलन जैसी अन्य प्रकार की तुलनाओं के बीच जानकारी प्राप्त करना चाहता हूं।

मैं सराहना करूंगा, अगर कोई मुझे बता सकता है कि सॉफ्टवेयर्स, और तुलनात्मक जानकारी प्राप्त करने के तरीके क्या हैं।

Comments

शायद आप एम और एन (एक्स और वाई के बजाय) का उपयोग करने के लिए प्रश्न को संपादित कर सकते हैं। पहले पढ़ने पर मुझे लगा कि X और Y समन्वय मूल्य हैं!

— मार्क आयरलैंड

प्रासंगिक सहायता प्राप्त करने के लिए आपको वास्तव में अधिक जानकारी प्रदान करने की आवश्यकता है। आपका वर्तमान प्रश्न बहुत कम दृश्य बनाता है। आपकी लास फ़ाइल में कौन सा फ़ील्ड मानों को पकड़े हुए है। जिस तरह से ग्राउंड को लास फॉर्मेट में वर्गीकृत किया गया है वह एक क्लासिफिकेशन फील्ड है न कि अलग-अलग z (ऊंचाई) मान। किसी विक्रेता को z मानों में अंतर रखने के लिए अनसाइन किए गए फ़ील्ड का उपयोग करना होगा।

— जेफरी इवांस

Answer 1

आर में दो डिजिटल एलिवेशन मॉडल (डीईएम) की तुलना कैसे करें ।

#-------------------------------------------------------------------------
#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 का उपयोग करके बनाए गए दो TINs की तुलना करने पर @ whuber का उत्तर देखें ? इस मुद्दे के बारे में एक सैद्धांतिक अंतर्दृष्टि के लिए।

Answer 2

आप ईएसआरआई आर्कगिस भूस्थैतिक विश्लेषण एक्सटेंशन के माध्यम से ऐसा कर सकते हैं - सबसेट पर सत्यापन करने में मदद पर एक अनुभाग है ।

आप R इंटरफ़ेस के माध्यम से GRASS के माध्यम से भी ऐसा कर सकते हैं । टॉमिस्लाव हेंगल ने अपनी पुस्तक ए प्रैक्टिकल गाइड टू जियोस्टेट्रिक मैपिंग में ऐसा करने के बारे में कुछ विस्तार से वर्णन किया है । यह खुली पहुंच है, इसलिए पीडीएफ डाउनलोड करने के लिए स्वतंत्र है।

Answer 3

जहाँ तक मुझे पता है, RMSE केवल DEM के निर्माण के दौरान ही कहा जाता है, न कि आगे की वापसी के लिए एक विशेषता के रूप में, इसलिए आपको DEM के निर्माण के दौरान इसे मैन्युअल रूप से "पकड़ना" होगा (यह कहा, मैंने कभी नहीं बनाया लिडार से एक डीईएम, केवल अन्य डेटा से)।

यदि आप दो DEMS के अंदर डेटा के बीच अंतर देखना चाहते हैं, तो मैं कट / भरण का उपयोग करूँगा जो कि आर्कगिस के स्थानिक विश्लेषक विस्तार ("भूतल विश्लेषण" के तहत) में है। कट / भरण आपको एक सरल विषयगत मानचित्र में डेम के बीच के परिवर्तनों को दिखाता है।

Comments

कट / भरण बहुत क्रूड है क्योंकि यह अंतरों को निर्धारित नहीं करता है। आरएमएसई दो डेटासेट की तुलना करने का एक सामान्य तरीका है: न केवल डीईएम को जमीनी सच्चाई के आंकड़ों से तुलना करना उपयोगी है, बल्कि यह दो डीईएम के बीच अंतर को निर्धारित करने का एक तरीका है।

— व्हिबर

Answer 4

मैं अंतर का एक साधारण डीईएम करूंगा। DEM2-DEM1। यह सभी क्षेत्रों को दिखाएगा जो अलग-अलग हैं और कितने से हैं।

मेरी वेबसाइट के मुखपृष्ठ पर अंतर के एक उच्च रेस डेम पर एक छवि है। thadwester.com
रंगीन बाईं छवि पर एक नज़र डालें।