आपको वास्तव में अपनी कार्यप्रणाली पर थोड़ा और अधिक शोध करने की आवश्यकता है और एस 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")