স্পেসিও-টেম্পোরাল ব্লক আর প্যাকেজ জিস্ট্যাট সহ ক্রিগিং?


9

এর ডকুমেন্টেশনে spআমি দেখেছি যে ফাংশনের coordsযুক্তিটি SpatialPolygonsকেবল 2 টি কলাম গ্রহণ করে accep ত্রি-মাত্রিক বহুভুজ তৈরি করা কি সম্ভব নয়? আমি gstatএকটি এর তৃতীয় স্থানিক মাত্রা হিসাবে সময় যুক্ত করে প্যাকেজ ব্যবহার করে স্প্যাটিও-টেম্পোরাল ক্রিগিং করছি SpatialPointsDataFrame। বহুগুণের সাথে ক্রিগিংকে ব্লক হিসাবে ব্লক করতে আমি এটি প্রসারিত করতে চাই। এটা করার কোন উপায় আছে?

উত্তর:


5

ইন sp, SpatialPoints*, SpatialPixels*এবং SpatialGrid*(সঙ্গে *বাদ দেওয়া বা দ্বারা প্রতিস্থাপিত DataFrame), 2 স্থানিক মাত্রা চেয়ে বেশি সমর্থন কি হিসাবে ওপি করেছেন, কিন্তু SpatialPolygons*এবং SpatialLines*না। সঙ্গে gstatআপনি 3-ডি অবরোধের (ব্যবহারে 3-ডি ব্লক kriging কি করতে পারেন block = c(10,10,10)আপনি ওপি চায় অ-আয়তক্ষেত্রাকার ব্লক জন্য এটা করতে পারব না), কিন্তু। তৃতীয় মাত্রার জন্য সময় বিকল্প হিসাবে নিখুঁতভাবে ঠিক আছে, তবে আপনি মেট্রিক এসটি ভেরোগ্রামে সীমাবদ্ধ।

library(gstat)
vignette("st")

আপনাকে ভ্যারোগ্রাম মডেলগুলির জন্য আরও বিকল্প দেয়, তবে ব্লক গড় মানগুলির পূর্বাভাস দেওয়ার জন্য নয় (এটি এফওয়াইআই, প্রশ্নের উত্তর নয়)।

এই প্রশ্নের একমাত্র উত্তর হ'ল 3 ডি শর্তসাপ্য সিমুলেশনগুলি করা এবং আপনার স্বেচ্ছাসেবী 3 ডি (2 ডি বহুভুজ + সময় ব্যাপ্তি) ব্লকগুলির উপরে সামগ্রিক পয়েন্ট মান। ক্লান্তিকর, তবে সম্ভব; এছাড়াও কেবল 3D পথ বরাবর, এসটি ভিগনেটে বর্ণিত পথটির krigeSTসাথে নয় ( এখনও সিমুলেশন করে না -!))


3

আপনার এসপি এস 4 শ্রেণির অবজেক্টের কাঠামো এবং প্রাসঙ্গিক জিস্ট্যাট ফাংশনগুলির সাথে এসপি অবজেক্টের মিথস্ক্রিয়া বুঝতে আপনার প্রকৃতপক্ষে আপনার আরও কিছুটা গবেষণা করা উচিত এবং ডকুমেন্টেশন পড়তে হবে। স্পে ভিগনেটে স্প্যাটিয়ালপলিগনস (কেবলমাত্র বহুভুজ টপোলজি) এবং স্পেশিয়ালপলিজোনডাটা ফ্রেম (বৈশিষ্ট্যযুক্ত বহুভুজ) বস্তুর মধ্যে পার্থক্যের বিশদ ব্যাখ্যা রয়েছে।

আপনি যা ব্যাখ্যা করছেন তা ক্রিগিংকে ব্লক করা নয় এবং সময় হিসাবে একটি গুণ হিসাবে ব্যবহার করা স্থানিক-অস্থায়ী অনুমানের ফলস্বরূপ নয়। আপনার বর্ণিত ধারণাগত পদ্ধতিটি বেশ অবৈধ। বহুভুজ বা বহুভুজ সেন্ট্রয়েড ব্যবহার করা অভিন্ন র্যান্ডম ক্ষেত্র, অ্যানিসোট্রপি এবং ননস্টেশনারিটির ক্রিগিং অনুমান লঙ্ঘন করে।

স্পেসটাইম প্যাকেজের ইন্টারফেসটি ব্যবহার করে স্থানিক-টেম্পোরাল মডেলগুলির জন্য এখানে একটি চমৎকার জিস্ট্যাট ভিগনেট is আমার এও লক্ষ্য করা উচিত যে সীমাবদ্ধ ক্রিগিং প্যাকেজটি ননস্টেশনারি গড় ফাংশন এবং আইসোট্রপিক দুর্বলভাবে স্থিতিশীল ভেরোগ্রাম ব্যবহার করে নির্বিচারে আকারের ব্লকগুলিতে ক্রিগিং ব্লক পরিচালনা করতে পারে।

এটি বলেছিল, আপনার প্রশ্নের উত্তর দেওয়ার জন্য, আপনি স্পেসিয়ালপয়েন্টস ডেটাফ্রেম অবজেক্টটি সরাসরি জিস্টেটের একটি ভেরোগ্রাম / ক্রিগিং মডেলটিতে পাস করতে পারেন। এসপি অবজেক্টের এই ধরণের বৈশিষ্ট্যগুলি "ডেটা" স্লটে থাকে এবং ইতিমধ্যে অভ্যন্তরীণ এস 4 শ্রেণি কাঠামোর মাধ্যমে স্থানাঙ্কগুলিতে সংযুক্ত থাকে।

# COERCE meuse DATAFRAME TO sp SpatialPointsDataFrame OBJECT
require(gstat)
data(meuse)
coordinates(meuse) <- ~ x + y
head(meuse@data)

# CREATE SEMIVARIOGRAM USING THE zinc ATTRIBUTE
# NOTE: THERE IS NO ARGUMENT FOR A "4th DIM"
v <- variogram(log(zinc) ~ 1, meuse)
plot(v, type = "l")

0

সময়টি একটি চতুর্থ মাত্রা হবে, জেড স্তর বা উচ্চতা তৃতীয় হবে, টিনস বা ডেমস বা জেড স্তরগুলি জিআইএসে উচ্চতা হ্যান্ডেল করবে, সময়টি সংস্করণযুক্ত ডাটাবেস দ্বারা প্রতিনিধিত্ব করা যেতে পারে। আপনি বহুভুজকে টিনে মানচিত্র তৈরি করতে পারেন বা একটি ডেমের সাথে ছেদ করতে পারেন, বা বহুভুজ আকৃতির পয়েন্টগুলি উপস্থাপনকারী এক্স কলামগুলিতে অ্যাজ কলাম যুক্ত করতে পারেন। কলাম হিসাবে বর্ধিত হতে পারে আপনি সময়ের সাথে মান এবং মান পরিবর্তনের একটি সময় ক্ষেত্র যোগ করতে পারেন। অবশ্যই ইএসআরআই এর পক্ষে সম্ভবত একটি সরঞ্জাম রয়েছে।


ভূতাত্ত্বিক বিশ্লেষক স্পেস-টাইম ক্রিগিংকে সমর্থন করে বলে আমি মনে করি না।
পল হিমস্ট্র্রা
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.