এলএএস ফাইল থেকে দুটি ডিজিটাল এলিভেশন মডেল (ডেম) তুলনা করছেন?

Question

আমার দুটি লিডার ফাইল (.las) রয়েছে, একটিটি এক্স পয়েন্টের সাহায্যে আসল। এবং অন্যটি হ'ল প্রথম .las ফাইলের অনুলিপি কিন্তু ওয়াই পয়েন্ট সহ, যেখানে ওয়াই এক্স এর চেয়ে কম।

এখন, আমি এই দুটি .las ফাইলের ডিজিটাল এলিভেশন মডেলগুলির (ডিইএম) তুলনা করতে এবং সেগুলি কতটা আলাদা তা কল্পনা করতে চাই।

আমি অন্যান্য ধরণের তুলনার মধ্যে আরএমএসই, স্ট্যান্ডার্ড বিচ্যুতির মতো তথ্য পেতে চাই।

আমি প্রশংসা করব, যদি কেউ আমাকে কোন সফ্টওয়্যার এবং তুলনা তথ্য পাওয়ার উপায়গুলি বলতে পারে।

Comments

সম্ভবত আপনি এম এবং এন (এক্স এবং ওয়াইয়ের পরিবর্তে) ব্যবহার করতে প্রশ্নটি সম্পাদনা করতে পারেন। প্রথম পড়ার সময় আমি ভেবেছিলাম এক্স এবং ওয়াই সমন্বিত মান!

— আয়ারল্যান্ড

প্রাসঙ্গিক সহায়তা পেতে আপনার আরও সত্যিকারের আরও তথ্য সরবরাহ করতে হবে। আপনার বর্তমান প্রশ্নটি খুব সামান্য দৃশ্য তৈরি করে। আপনার লাস ফাইলের কোন ক্ষেত্রটি মান ধারণ করে। লাস ফর্ম্যাটে যেভাবে স্থল শ্রেণিবদ্ধ করা হয়েছে তা শ্রেণিবদ্ধকরণ ক্ষেত্র এবং বিভিন্ন জেড (উচ্চতা) মান নয়। 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.

ডেস্কটপের জন্য আর্কজিআইএস ব্যবহার করে তৈরি দুটি টিআইএন তুলনা করার বিষয়ে @ হোবারের উত্তর দেখুন ? এই সমস্যা সম্পর্কে একটি তাত্ত্বিক অন্তর্দৃষ্টি জন্য।

Answer 2

আপনি এটি ইএসআরআই আর্কজিআইএস জিওস্টাটিস্টিকাল অ্যানালাইসিস এক্সটেনশনের মাধ্যমে করতে পারেন - সাবসেটগুলিতে বৈধতা প্রদর্শনের জন্য সহায়তার একটি বিভাগ রয়েছে ।

আরআর ইন্টারফেসের মাধ্যমে আপনি গ্রাসের মাধ্যমেও এটি করতে পারেন । টমিস্লাভ হেনগেল তাঁর প্রযোজনীয় গাইড টু জিওস্ট্যাটাস্টিকাল ম্যাপিং বইটিতে এটি কীভাবে করবেন তা কিছুটা বিশদে বর্ণনা করেছেন । এটি উন্মুক্ত অ্যাক্সেস, তাই পিডিএফ ডাউনলোড করতে বিনামূল্যে।

Answer 3

যতদূর আমি জানি, আরএমএসই কেবল ডিএম তৈরির সময় বলা হয়েছে, এবং আরও বেশি পরিমাণে শ্রদ্ধার কারণ হিসাবে নয়, সুতরাং আপনাকে ডিএম তৈরির সময় নিজে হাতে "ধরতে" হবে (যে বলেছিল, আমি কখনই তৈরি করি নি লিডার থেকে একটি ডিইএম, কেবলমাত্র অন্যান্য ডেটা থেকে)।

আপনি যদি দুটি ডিএমএসের মধ্যে ডেটা-র মধ্যে পার্থক্যগুলি দেখতে চান তবে আমি কাট / ফিল ব্যবহার করব যা আর্কজিআইএসের স্পেসিয়াল অ্যানালিস্ট এক্সটেনশনে রয়েছে ("সারফেস অ্যানালাইসিস" এর অধীনে)। কাটা / ভরাট আপনাকে একটি সাধারণ বিষয়ভিত্তিক মানচিত্রে ডেমের মধ্যে পরিবর্তনগুলি দেখায়।

Comments

কাটা / ফিল করা অত্যন্ত অপরিশোধিত কারণ এটি পার্থক্যগুলিকে মাপ দেয় না। আরএমএসই হ'ল দুটি ডেটাসেটের তুলনা করার একটি সাধারণ উপায়: একটি ডেমকে স্থল-সত্যের ডেটার সাথে তুলনা করা কেবল কার্যকর নয়, এটি দুটি ডিইএমের মধ্যে পার্থক্যের পরিমাণ নির্ধারণের এক উপায়।

— হোবার

Answer 4

আমি পার্থক্যটির একটি সাধারণ ডিএম করব। DEM2-DEM1। এটি সমস্ত অঞ্চল পৃথক এবং কত দ্বারা প্রদর্শিত হবে।

আমার ওয়েবসাইটের হোমপৃষ্ঠাতে একটি চিত্রকে উচ্চতর পার্থক্যে ফেলেছে। thadwester.com
রঙিন বাম ইমেজ একবার দেখুন।