एक औसत कर्नेल घनत्व मानचित्र का समर्थन करने के लिए एक त्रुटि मानचित्र कैसे बनाएं?


10

मैंने एक ही स्थानिक सीमा के भीतर ढेर किए गए बिंदुओं पर केडीई चलाकर एक औसत कर्नेल घनत्व नक्शा बनाया है। उदाहरण के लिए, मान लें कि हमारे पास एक ही आकार और आकार के तीन अलग-अलग वन अंतरालों में रोपाई का प्रतिनिधित्व करने वाले तीन बिंदु आकार आकृति हैं। मैंने प्रत्येक बिंदु आकृति के लिए एक केडीई चलाया। आर्क के रैस्टर कैलकुलेटर में औसत की गणना करने के लिए, केडीई से आउटपुट को स्थानिक सीमा के आधार पर स्टैक किया गया Float(("KDE1"+"KDE2"+"KDE3")/3)। यहाँ अंतिम उत्पाद है:

यहां छवि विवरण दर्ज करें

अब मैं औसत केडीई से जुड़ी त्रुटि को दर्शाते हुए एक मानचित्र बनाने में दिलचस्पी रखता हूं। मुझे उम्मीद है कि त्रुटि मानचित्र का उपयोग नेत्रहीन चित्रण करने के लिए हो सकता है कि हॉटस्पॉट्स के साथ कितनी त्रुटि है (उदाहरण के लिए एक अंतराल में पूरी तरह से अंक के कारण एसडब्ल्यू हॉटस्पॉट?)। मुझे औसत केडीई से संबंधित त्रुटि का नक्शा बनाने के बारे में कैसे जाना चाहिए? क्या MSE इस मामले में त्रुटि का सबसे उचित उपाय होगा ?


3
यह बहुत दिलचस्प विश्लेषण है। "मानक त्रुटि" से आपका क्या तात्पर्य है? किसी प्रकार का विचलन (अंतर) प्रत्येक घनत्व मानचित्र "मीन" परत से?
लैंडस्केप एनालिसिस

@Landscape विश्लेषण पोस्ट टिप्पणियों को स्वीकार करने के लिए संपादित किया गया। हां, मैं सोच रहा हूं कि इस मामले में एक एमएसई अनुमान सबसे उपयुक्त हो सकता है। अनिवार्य रूप से, यह दिखाते हुए कि प्रत्येक केडीई औसत केडीई से कैसे भिन्न होता है। मुझे यकीन नहीं है कि यह सब कैसे एक साथ आर्कजीआईएस और / या स्क्रिप्टिंग का उपयोग करना है।
हारून

जवाबों:


8

एक कैविएट

डेटा में कोई व्यवस्थित त्रुटि नहीं होने पर, नमूना डेटा से अनिश्चितता का अनुमान लगाने के लिए एक मानक त्रुटि एक उपयोगी तरीका है यह धारणा इस संदर्भ में संदिग्ध वैधता की है, क्योंकि (ए) केडीई मानचित्रों में स्थानीय रूप से निश्चित त्रुटियां होंगी, जो परतों के बीच व्यवस्थित रूप से बनी रह सकती हैं और (बी) कर्नेल त्रिज्या (या "बैंडविड्थ) की पसंद के कारण अनिश्चितता का एक संभावित विशाल घटक है। ") इन नक्शों के किसी भी एक संग्रह में बिल्कुल भी प्रतिबिंबित नहीं होगा।

कुछ विकल्प

फिर भी, संबंधित, ध्वस्त ("स्टैक्ड") मानचित्रों के संग्रह के बीच परिवर्तनशीलता को चित्रित करना एक महान विचार है - बशर्ते आपको वर्णित सीमाओं को याद रखें। इस सेटिंग में स्थानीय परिवर्तनशीलता के कई उपाय स्वाभाविक होंगे, जिनमें शामिल हैं:

  • मूल्यों की श्रेणी , या तो additively (अधिकतम माइनस न्यूनतम) या गुणा (अधिकतम न्यूनतम द्वारा विभाजित ) व्यक्त की गई है

  • विचरण या मानक विचलन मूल्यों की। इसका गुणात्मक संस्करण मानों के लघुगणक का विचलन या मानक विचलन होगा ।

  • फैलाव का एक मजबूत अनुमानक, जैसे कि इंटरक्वेर्टाइल रेंज (या तीसरे से पहली चतुर्थांश का अनुपात)।

कई मामलों में, गुणन उपाय घनत्व के लिए अधिक उपयुक्त हो सकते हैं, क्योंकि प्रति एकड़ 100 और 101 पेड़ों के बीच का अंतर असंगत हो सकता है, जबकि प्रति एकड़ 2 और 1 पेड़ों के बीच का अंतर अपेक्षाकृत महत्वपूर्ण हो सकता है। दोनों १०१ - १०० = २ - १ - १ = १ की एक ही (योजक) सीमा को प्रदर्शित करते हैं, लेकिन उनकी १.०१ और २.२ की गुणात्मक सीमा काफी भिन्न होती है। (ध्यान दें कि एक गुणक सीमा हमेशा 1 से अधिक होती है, इसलिए कि 1 से 1.01 से 2.00 एक सौ गुना आगे है।)

गणना

इन उपायों की गणना के लिए स्थानीय आंकड़ों के कुछ रूप की आवश्यकता होती है सेल आँकड़े स्थानिक विश्लेषक में कार्यक्षमता प्रसरण, पर्वतमाला, और मानक विचलन की गणना करेगा। स्थानीय मात्राएँ रैंक के साथ मिल सकती हैं । उधम मचाते रहने के बजाय, जो रैंकों का उपयोग करने के लिए सुविधाजनक है, का उपयोग करने के लिए रैंक करता है। उन्हें खोजने के लिए, n ढेर में ग्रिड की संख्या हो। माध्यिका की एक रैंक होती है (n + 1) / 2 - जो कि पूरी संख्या नहीं हो सकती है, यह दर्शाता है कि इसे n / 2 और n / 2 + 1 रैंक के औसत से गणना की जानी चाहिए, जिसमें से कोई भी माध्यिका अनुमानित होगी। चतुर्थकों को अनुमानित करने के लिए, फिर, (n + 1) / 2 को निकटतम पूर्ण संख्या में नीचे करें, फिर 1 जोड़ें और 2 से विभाजित करें। इस संख्या को r होने दें । उपयोगआर और एन + 1 - आर क्वार्टल्स के रैंकों के लिए।

एक उदाहरण के रूप में, यदि स्टैक में n = 6 ग्रिड हैं, (n + 1) / 2 राउंड डाउन 3 है और (3 + 1) / 2 = 2 राउंडिंग की आवश्यकता नहीं है। उपयोग आर = 2 और आर = 6 + 1 - 2 = 5 रैंकों के लिए। वास्तव में, यह प्रक्रिया प्रत्येक सेल में छह मूल्यों के दूसरे सबसे कम ( आर = 2) और दूसरे उच्चतम ( आर = 5) मान लौटाएगी। आप या तो उनके अंतर या उनके अनुपात का नक्शा बना सकते हैं।


2
मैं कल्पना करूंगा कि भिन्नता का गुणांक इस संदर्भ में उपयोगी होगा।
जेफरी इवांस

1
@ जेफ्री धन्यवाद, यह भी एक विचार है। सीवी (स्थानीय) माध्य ग्रिड द्वारा स्थानीय (स्थानीय) मानक विचलन ग्रिड को विभाजित करके प्राप्त किया जाता है। मैंने इसका उल्लेख नहीं किया, लेकिन इस तरह के गुणात्मक सारांशों के लिए, कुछ क्षेत्रों को उन क्षेत्रों से बाहर निकालने के लिए ध्यान रखा जाना चाहिए जहां भाजक (मतलब या न्यूनतम, जैसा भी मामला हो) शून्य के करीब है: परिणाम वहां अविश्वसनीय हो सकते हैं और संभव है गुठली को सन्निकट करने में संख्यात्मक अभेद्यता और छोटे अशुद्धियों के अलावा कुछ भी नहीं दर्शाएगा।
whuber

क्या आप अपने पहले पैराग्राफ में (ए) पर विस्तृत जानकारी दे सकते हैं? उदाहरण के लिए, क्या आप डेटा संग्रहण तकनीकों में त्रुटियों का उल्लेख कर रहे हैं जो प्रत्येक वन गैप (और इसलिए व्यवस्थित रूप से प्रत्येक केडीई रेखापुंज में खुद को प्रकट कर सकते हैं), या फोकल फ़ंक्शन के कार्यान्वयन के लिए बंधी त्रुटियों को संदर्भित कर सकते हैं ?
मन्नीग

@ शुभंकर चारों ओर महान विचार - बहुत धन्यवाद!
आरोन
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.