दो DEM की तुलना कैसे करें और अपना अंतर कैसे दिखाएं?


10

मैं अपने DEMs की एक-दूसरे (ASTER GDEM2 और SRTM4.1) से तुलना करना चाहता हूं। मुझे पता है, कि GDEM में गलत ऊँचाई के साथ कुछ "छेद" हैं (मैंने अपने कार्य क्षेत्र में उनमें से 3 को पाया)। अब मुझे यह जानने की जरूरत है, अगर उनमें से अधिक को अंततः SRTM के साथ ठीक करना है। मैं SAGA, GRASS और QGIS का उपयोग करता हूं।


उनकी तुलना किससे करें? एक दूसरे के साथ? एक आप आधारभूत हैं और दूसरा दोषपूर्ण है?
13

सही! मैं सिर्फ उनके बीच का अंतर निकालना चाहता हूं।
एमएपी

जवाबों:


5

मान लें कि Raster1 सही है, और Raster2 वह है जिसमें आपको संदेह है कि इसमें दोष हैं।

वे दोनों भी एक ही इकाई में हैं (चूंकि आप कहते हैं कि दोनों आपदाएं डीईएम हैं, जेड मान मीटर में होना चाहिए - मुझे लगता है कि हम सही हैं और हम समान चीजों की तुलना करते हैं)

आप उन्हें तुलना करने के लिए पैकेज (qgis, घास, गाथा) से रैस्ट कैल्क का उपयोग कर सकते हैं:

deviation_rast  = abs( (raster1 - raster2) / raster1 * 100 ) 

परिणामी रेखापुंज से प्रत्येक सेल प्रतिशत प्रदर्शित करेगा जो raster2 raster1 से विचलित करता है।

उसके बाद आप अपनी आवश्यकताओं के आधार पर विचलन_धार को पुनर्वर्गीकृत कर सकते हैं: उदा

red -> deviation_rast_value > 75
orange -> deviation_rast_value > 50
yellow ->  deviation_rast_value > 25
green -> deviation_rast_value  > 0

समस्याग्रस्त क्षेत्रों की पहचान करना - और यह चुनना कि आप हस्तक्षेप करना चाहते हैं या नहीं


8
डीईएम के लिए समुद्र तल के पास की ऊंचाई के साथ अनुपात एक अच्छा विचार नहीं है: आपको जहां भी raster2 == 0. बहुत सारे फर्जी परिणाम और पूरी तरह से अपरिभाषित परिणाम मिलेंगे, सभी डीईएम के लिए सही दृष्टिकोण, बस दूसरे और स्क्रीन से एक ढलाईकार को घटाना है। आउटलायर्स के लिए (जो कि कई तरीकों से किया जा सकता है: कुछ तरीकों के लिए आँकड़ो से परामर्श करें ।stackexchange.com )।
whuber

11

आप दो DEMs का अंतर प्राप्त करके शुरू कर सकते हैं। क्यूजीआईएस में एक रास्टर कैलकुलेटर टूल है जो काम में आना चाहिए। बस दोनों डेम की टाइलें प्राप्त करें जो एक ही क्षेत्र को कवर करती हैं और एक डेम के मूल्यों को दूसरे से घटाती हैं। आपको एक अच्छी रेखापुंज परत मिलनी चाहिए जो दो DEMs के बीच ऊंचाई में अंतर दिखाती है।


5
वास्तव में, बस एक DEM को दूसरे से घटाएं। मैं फिर इसे QGIS में प्रस्तुत करूंगा या शून्य के चारों ओर कुछ सहिष्णुता का उपयोग करते हुए इसे पुन: सूचीबद्ध करूंगा। अगर ऐसा कुछ है तो -5m <परिणाम <5m तो विचार करें कि DEM का उस विशेष सेल में समान मूल्य है।
अलेक्जेंड्रे नेटो

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

मुझे एक समान समस्या है, लेकिन मेरे चीर-फाड़ केवल आंशिक रूप से ओवरलैप हो रहे हैं (रेखापुंज बी केवल रेखापुंज A के हिस्से को कवर करता है) और "ग्रिड अंतर" उपकरण काम नहीं करेगा। कोई भी विचार जो मैं कर सकता हूं?
एच। वीनर
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.