आर-पैकेज gstat के साथ स्पॉटी-टेम्पोरल ब्लॉक क्रिंगिंग?


9

spमैंने दस्तावेज़ीकरण में देखा कि फ़ंक्शन coordsका तर्क SpatialPolygonsकेवल 2 कॉलम स्वीकार करता है। क्या 3-आयामी बहुभुज बनाना संभव नहीं है? मैं gstatएक तीसरे स्थानिक आयाम के रूप में समय जोड़कर पैकेज का उपयोग करके स्पैट-टेम्पोरल क्रिपिंग कर रहा हूं SpatialPointsDataFrame। मैं विस्तार करना चाहूंगा कि ब्लॉक के रूप में बहुभुज के साथ सिंचाई करना। क्या उसे करने का कोई तरीका है?

जवाबों:


5

में sp, SpatialPoints*, SpatialPixels*और SpatialGrid*(के साथ *छोड़े गए या द्वारा प्रतिस्थापित DataFrame), 2 स्थानिक आयामों से अधिक समर्थन करते हैं ओपी किया गया है, लेकिन SpatialPolygons*और SpatialLines*नहीं है। साथ gstatआप 3 डी ब्लॉक (उपयोग करने के साथ 3 डी ब्लॉक kriging कर सकते हैं block = c(10,10,10)आप ओ पी चाहता है के रूप में गैर आयताकार ब्लॉकों के लिए ऐसा नहीं कर सकते हैं), लेकिन,। तीसरे आयाम के लिए समय बदलने के लिए यह पूरी तरह से ठीक है, लेकिन आप मीट्रिक एसटी वैरोग्राम के लिए विवश हैं।

library(gstat)
vignette("st")

आपको वैरोग्राम मॉडल के लिए अधिक विकल्प देता है, लेकिन ब्लॉक माध्य मान की भविष्यवाणी के लिए नहीं (यह FYI है, प्रश्न का उत्तर नहीं है)।

प्रश्न का एकमात्र उत्तर 3 डी सशर्त सिमुलेशन करना होगा, और आपके 3 डी (2 डी बहुभुज + समय सीमा) ब्लॉकों पर कुल बिंदु मानों को एकत्रित करना होगा। थकाऊ, लेकिन संभव; केवल 3 डी पथ के साथ, एसटी विगनेट में वर्णित पथ के साथ नहीं ( krigeSTसिमुलेशन नहीं करता है - अभी तक!)।


3

आपको वास्तव में अपनी कार्यप्रणाली पर थोड़ा और अधिक शोध करने की आवश्यकता है और एस 4 श्रेणी की वस्तुओं की संरचना और प्रासंगिक gstat कार्यों के साथ sp ऑब्जेक्ट्स की बातचीत को समझने के लिए प्रलेखन को पढ़ना चाहिए। Sp Vignette में SpatialPolygons (केवल बहुभुज टोपोलॉजी) और SpatialPolygonDataFrame (विशेषता वाले बहुभुज) वस्तुओं के बीच अंतर का विस्तृत विवरण है।

आप जो समझा रहे हैं, वह क्रिंगिंग को ब्लॉक नहीं कर रहा है और समय का एक विशेषता के रूप में उपयोग करने से स्थानिक-अस्थायी अनुमान नहीं लगता है। आपके द्वारा वर्णित वैचारिक पद्धति काफी अमान्य है। बहुभुज या बहुभुज सेंट्रोइड्स का उपयोग करना एक समान यादृच्छिक क्षेत्र, अनिसोट्रॉपी और गैरबराबरी की क्रिंगिंग मान्यताओं का उल्लंघन करता है।

यहाँ स्पेसटाइम -टेम्पोरल मॉडल पर एक अच्छा gstat vignette है जो स्पेसक्राफ्ट पैकेज के इंटरफेस का उपयोग कर रहा है। मुझे यह भी ध्यान देना चाहिए कि विवश किरिगिंग पैकेज एक गैर-मध्यम माध्य फ़ंक्शन और एक आइसोट्रोपिक कमजोर रूप से स्थिर वैरोग्राम का उपयोग करके मनमाने आकार के ब्लॉकों पर ब्लॉक क्रिगिंग का संचालन कर सकता है।

उस ने कहा, अपने प्रश्न का उत्तर देने के लिए, आप एक स्पैस्टियल पॉइंट्सटाटाफ्रेम ऑब्जेक्ट को सीधे gstat में एक वेरोग्राम / क्रिगिंग मॉडल में पास कर सकते हैं। इस प्रकार की sp ऑब्जेक्ट में, विशेषताएँ "डेटा" स्लॉट में रहती हैं और आंतरिक S4 वर्ग संरचना के माध्यम से पहले से ही निर्देशांक से जुड़ी होती हैं।

# 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

समय एक चौथा आयाम होगा, z स्तर या ऊंचाई तीसरी होगी, tins या dems या z स्तर जीआईएस में ऊंचाई को संभालते हैं, समय को संस्करण डेटाबेस द्वारा दर्शाया जा सकता है। आप बहुभुज को एक टिन के साथ मैप कर सकते हैं या एक डीम के साथ इंटरसेक्ट कर सकते हैं, या एज़ कॉलम को पॉलीगॉन शेप पॉइंट्स का प्रतिनिधित्व करने वाले कॉलम में जोड़ सकते हैं। हो सकता है कि आप कॉलम के रूप में बढ़े हुए समय के माध्यम से समय क्षेत्र के साथ-साथ मूल्यों और मूल्य परिवर्तनों की एक श्रृंखला जोड़ सकें। बेशक ESRI के पास इसके लिए एक उपकरण है।


मुझे नहीं लगता कि भू-स्थानिक विश्लेषक अंतरिक्ष-समय की सिंचाई का समर्थन करता है।
पॉल हेमस्ट्रा
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.