ढलान की गणना के लिए SRTM Global DEM का उपयोग करना?


17

मैंने SRTM GDEM (~ 90 किमी रिज़ॉल्यूशन) डाउनलोड किया है।

मैं आर्कगिस 10 का उपयोग कर रहा हूं।

मैंने ढलान के लिए गणना करने के लिए स्थानिक विश्लेषक का उपयोग करने की कोशिश की है।

हालांकि, मैं ढलान के लिए गणना नहीं कर सकता।

आउटपुट मानों में केवल दो रेंज 0 और 0.1-90 हैं।

मैं वास्तव में निश्चित नहीं हूं कि समस्या क्या है?


यह इस बात पर निर्भर करता है कि आप दुनिया में कहां विश्लेषण कर रहे हैं। प्रत्येक स्थान के लिए अलग-अलग अनुमान हैं। आप कहां परीक्षा दे रहे हैं?
djq

5
संकल्प वास्तव में ~ 90 मीटर है, न कि ~ 90 किमी।
अखिलेश 10

बस एक टिप्पणी, यदि आपका डेस्कटॉप के रखरखाव पर है, तो आप आर्कजीआईएस ऑनलाइन पर लॉग इन कर सकते हैं और उनकी उन्नयन सेवाओं का उपयोग कर सकते हैं (कोई एक्सटेंशन की आवश्यकता नहीं है)। ढलान परत एक संदर्भ परत के रूप में उपयोग करने के लिए स्वतंत्र है। ऑस्ट्रेलिया में, हमारे पास 1 दूसरा SRTM डेटा (~ 30 मीटर रेस) blogs.esri.com/esri/arcgis/2014/07/11/…-
साइमन

जवाबों:


29

यह विश्व स्तर पर व्यापक CEM के लिए ढलानों की गणना करने के लिए एक सरल, तेज, और यथोचित सटीक तरीके से अधिक वर्णन करने के लिए एक अच्छी जगह की तरह लगता है

सिद्धांतों

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

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

मर्केटर प्रोजेक्शन बिल को फिट करता है: मान लेना भूमध्य रेखा पर सही है, इसकी विकृति अक्षांश के सेकंड के बराबर होती है। यही है, मानचित्र पर दूरी को सेक्युलर द्वारा गुणा किया जाता है। यह वृद्धि की गणना करने के लिए किसी भी ढलान की गणना का कारण बनता है: (सेकंड (एफ) * रन) (जो एक अनुपात है), जहां एफ अक्षांश है। इसे ठीक करने के लिए, हमें गणना (ढ) द्वारा गणना की गई ढलानों को गुणा करना होगा; या, समतुल्य, उन्हें cos (f) से विभाजित करते हैं। यह हमें सरल नुस्खा देता है:

एक मर्केटर प्रोजेक्शन का उपयोग करके ढलान (वृद्धि: रन या प्रतिशत) की गणना करें, फिर परिणाम को अक्षांश के कोसाइन द्वारा विभाजित करें।

कार्यप्रवाह

दशमलव डिग्री में दी गई ग्रिड (जैसे SRTM DEM) के साथ ऐसा करने के लिए, निम्न चरणों का पालन करें:

  1. एक अक्षांश ग्रिड बनाएं। (यह सिर्फ y- समन्वित ग्रिड है।)

  2. इसकी कोज्या की गणना करें।

  3. परियोजना दोनों डीईएम और अक्षांश एक मर्केटर प्रक्षेपण जिसमें पैमाने भूमध्य रेखा पर सच है का उपयोग कर की कोज्या।

  4. यदि आवश्यक हो, तो अनुमानित निर्देशांक (आमतौर पर मीटर) की इकाइयों से सहमत होने के लिए उन्नयन इकाइयों को परिवर्तित करें।

  5. या तो शुद्ध ढलान या प्रतिशत के रूप में अनुमानित DEM की ढलान की गणना करें ( रूप में करें कोण के रूप में नहीं )।

  6. इस ढलान को अनुमानित कोसाइन (अक्षांश) ग्रिड द्वारा विभाजित करें।

  7. यदि वांछित है, तो आगे के विश्लेषण या मानचित्रण के लिए ढलान ग्रिड को किसी अन्य समन्वय प्रणाली में फिर से डालें।

ढलान की गणना में त्रुटियां 0.3% तक होंगी (क्योंकि यह प्रक्रिया एक दीर्घवृत्त के बजाय एक गोलाकार पृथ्वी मॉडल का उपयोग करती है, जिसे 0.3% द्वारा समतल किया जाता है)। यह त्रुटि अन्य त्रुटियों की तुलना में काफी कम है जो ढलान की गणना में जाती हैं और इसलिए इसे उपेक्षित किया जा सकता है।


पूरी तरह से वैश्विक गणना

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

संपूर्ण वैश्विक समाधान के लिए, स्थलीय ग्रिड को तीन भागों में तोड़ने पर विचार करें - प्रत्येक ध्रुव के चारों ओर और भूमध्य रेखा के चारों ओर -, एक उपयुक्त प्रक्षेपण का उपयोग करके प्रत्येक भाग में अलग से एक ढलान गणना का प्रदर्शन करना, और परिणामों को मोज़ाइक करना। ग्लोब को विभाजित करने के लिए एक उचित स्थान 2 * आर्कटैन (1/3) के अक्षांशों पर अक्षांशों के घेरे में है, जो लगभग 37 डिग्री है, क्योंकि इन अक्षांशों पर मर्केटर और स्टेरियोग्राफिक सुधार कारक एक दूसरे के बराबर हैं (एक सामान्य मूल्य होने पर) 5/4) और किए गए सुधार के आकारों को कम करना अच्छा होगा। गणनाओं की जांच के रूप में, ग्रिड बहुत करीबी समझौते में होना चाहिए, जहां वे ओवरलैप करते हैं (अनुमानित बिंदुओं की पुन: प्राप्ति के कारण फ्लोटिंग पॉइंट इंप्रेशन और अंतर की थोड़ी मात्रा में विसंगतियों का एकमात्र स्रोत होना चाहिए)।

संदर्भ

जॉन पी। स्नाइडर, मैप प्रोजेक्शंस - ए वर्किंग मैनुअल । यूएसजीएस प्रोफेशनल पेपर 1395, 1987।


7
मैं खुद को इस स्थिति में पाता हूं, जैसा कि मैं अक्सर करता हूं, एक समाधान का वर्णन करने के लिए एक बार फिर से व्हुबर को धन्यवाद देता हूं और साथ ही इसे बनाने वाले तर्क को भी सामने रखता हूं । मेरी टोपी तुम्हारे पास है सर।
मैट विल्की

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

आह। स्पष्टीकरण के लिए धन्यवाद। मैंने जवाब बहाल कर दिया है, लोगों को भरोसा है कि उनके सामने अब पर्याप्त जानकारी है एक सूचित विकल्प बनाने के लिए :)
मैट विल्की

2
एक फ्रेंच बैकग्राउंड से आने के बाद, मुझे इस बेहतरीन उत्तर को समझने के लिए आवश्यक शब्दावली का अनुवाद करने में थोड़ा समय
अखिलेश

यह एक महान दृष्टिकोण है और मैंने पहले से ही वैश्विक ढलान रास्टर उत्पन्न करने के लिए आपके समाधान का उपयोग किया है। व्यावहारिक अनुभव से एक संकेत: चूंकि अक्षांश मान दक्षिण के ऋणात्मक हैं, इसलिए आपको निम्नलिखित समीकरण में पूर्ण अक्षांश मान का उपयोग करना होगा: (1 + पाप (अक्षांश)) / 2
सेलिका

18

मूल उत्तर

मुझे लगता है कि आपके रेखापुंज के लिए क्षैतिज इकाइयाँ या तो डिग्री या आर्सेकंड में हैं। आपको इस रेखापुंज को एक स्थानिक प्रक्षेपण के लिए फिर से प्रस्तुत करना होगा जहां आपकी क्षैतिज और ऊर्ध्वाधर इकाइयां समान हैं (यानी, यदि ऊर्ध्वाधर इकाइयां मीटर में हैं, तो मैं यूटीएम का उपयोग करने का सुझाव देता हूं, जिसमें मीटर की क्षैतिज इकाइयां हैं)।

ArcCatalog / ArcGIS के साथ एक रेखापुंज को पुन: प्रकट करने के लिए, इसमें देखें:

आर्कटूलबॉक्स> डेटा प्रबंधन उपकरण> अनुमान और रूपांतरण> रेखापुंज> परियोजना रेखापुंज

एक अनुमानित स्थानिक संदर्भ चुनें, जो आपकी रुचि के क्षेत्र को कवर करता है, उदाहरण के लिए, UTM ज़ोन का प्रयास करें। कई अन्य विकल्प हैं जो मैनुअल में सबसे अच्छे दस्तावेज हैं । ध्यान दें, आप पूरी पृथ्वी के लिए एक ढलान डेटासेट नहीं बना सकते हैं (यदि ऐसा है जो आप करने की कोशिश कर रहे हैं)।

बेहतर उत्तर, स्केल के साथ GDAL का उपयोग करना

अब जबकि SRTM डेटा विश्व स्तर पर उपलब्ध है मैं वास्तव में फाइलों के साथ देख और काम कर सकता हूं। gdaldemGDAL से उपयोगिता ढलान की गणना और एक का उपयोग कर सकते हैं hillshade पैमाने क्षैतिज खड़ी इकाइयों के एक अनुपात के लिए विकल्प। मैनुअल SR11M टाइल्स की तरह कुछ के लिए 111120 मीटर / ° की सिफारिश करता है। उदाहरण के लिए, एक OSGeo4W शेल से:

$ gdaldem slope -s 111120 -compute_edges N44E007.hgt N44E007_slope.tif

-compute_edgesअगर आप कुछ टाइल्स साथ पिरो करना चाहते हैं विकल्प, किनारों अधिक सहज बना देता है। या एक बड़े क्षेत्र के लिए टाइल की गणना करें। "स्केल" तकनीक के साथ नुकसान यह है कि ईडब्ल्यू और एनएस दिशाओं में दूरी भूमध्य रेखा को छोड़कर समान नहीं है, इसलिए ध्रुवों के करीब टाइलों के लिए ढलान के कुछ अजीब गलत चित्रण हो सकते हैं।


यह आपकी अंतिम टिप्पणी पर जोर देने के लायक है: यह भूमध्य रेखा के पास नहीं अंक के लिए एक खराब समाधान है। यह "विषम गलत व्याख्याओं" का एक छोटा मामला नहीं है: परिणाम सकल गलत होंगे, विशेष रूप से भूमध्य रेखा की तुलना में ध्रुवों के करीब स्थानों में। gdaldemराज्यों के लिए प्रलेखन "भूमध्य रेखा के पास नहीं स्थानों के लिए, यह gdaldem का उपयोग करने से पहले अपने ग्रिड को gwarwarp का उपयोग करने से रोकने के लिए सबसे अच्छा होगा।" दुर्भाग्य से यह ग्लोब को कवर करने वाले डेटासेट के लिए काम नहीं करेगा, जब तक कि आप उन्हें छोटे टुकड़ों (74 UTM ज़ोन, शायद?) में नहीं तोड़ते, उन प्रोजेक्ट्स को करें, जो ढलान की गणना करते हैं, और परिणामों को मोज़ेक करते हैं।
whuber

7

सीधे शब्दों में कहें, वहाँ एक नहीं है। परिभाषा के अनुसार, डिग्रियों पर आधारित एक समन्वय प्रणाली अन-अनुमानित है। आम बोलचाल में हम कहते हैं कि WGS84 एक "भौगोलिक" प्रक्षेपण है, लेकिन यह असत्य है, बस सुविधा के लिए।

मुझे लगता है कि मुझे एक सॉफ्टवेयर या प्रक्रिया के बारे में पढ़ने की याद है, जो गैर-अनुमानित भौगोलिक स्थान में ऊंचाई मॉडल के साथ सटीक रूप से काम करने के लिए है, लेकिन मैं अभी इसका पता नहीं लगा सकता। किसी भी मामले में यह प्रायोगिक रहा होगा या इसे कोड प्रक्रिया से स्वयं निर्मित करेगा।


अहह, यह पाया गया: भूकंप (यूएसजीएस) से भूस्खलन घटना के परिणाम के आकलन के लिए एक वैश्विक ढलान डेटासैट का विकास । पेज 4 समस्या का अच्छी तरह से वर्णन करता है

... एक डिग्री की लंबाई इसके अक्षांशीय स्थान के आधार पर भिन्न होती है। भूमध्य रेखा पर, एक-डिग्री ब्लॉक द्वारा एक-डिग्री काफी वर्ग है जब मीटर की इकाइयों में परिवर्तित किया जाता है (x- दिशा में 111,321 मीटर x- दिशा में 110,567 मीटर से y-दिशा में ... लेकिन ध्रुवों के करीब दूरी। x- दिशा अक्षांश के कोसाइन के एक समारोह के रूप में छोटी हो जाती है, मध्याह्न के अभिसरण के कारण। अधिकांश जीआईएस पैकेज, आर्कगिस शामिल हैं, केवल वर्ग पिक्सेल पर काम करते हैं, और इसलिए x, y, या z आयामों को समायोजित करने के लिए एक कारक का उपयोग करते हैं। एक सामान्य इकाई संभव नहीं है।

पेपर विशिष्ट गणनाओं और सॉफ्टवेयर टूल्स ( , , ) का वर्णन करने के लिए जाता है जो वे इस मूल मुद्दे को हल करने के लिए उपयोग करते थे। कागज में कोड शामिल नहीं है, लेकिन अगर उनसे पूछा जाए तो वे साझा कर सकते हैं। किसी भी मामले में हालांकि मैं शायद सिर्फ यह पूछूंगा कि परिणाम कहां हैं, यूएसजीएस होने के नाते यह शायद कहीं पहले से ही ऑनलाइन है। :)


1
उस कागज का सुझाव है कि ढलानों की गणना करने के लिए अज़ीमुथल समतुल्य प्रक्षेपण का इस्तेमाल किया जा सकता है और यह गलत है। यह वास्तव में प्रक्षेपण के मूल के पास सही ढलान देगा, लेकिन वे भी, उत्तरोत्तर कम सटीक प्राप्त करेंगे क्योंकि मूल बढ़ जाती है।
व्हीबर

इशारा करने के लिए धन्यवाद। पाठकों, कृपया संतुलन के लिए gis.stackexchange.com/a/40464/108 और साथ ही पढ़ने के लिए सुनिश्चित रहें
मैट विल्की

2

ग्लोबल डीईएम पैरामीटर (जहां अधिकांश सूत्र यूक्लिडियन स्थान की धारणा पर आधारित होते हैं) को एक्वि 7 जीआरआईडी प्रणाली (बाउर-मार्सक्लिंजर एट अल। 2014) का उपयोग करके कुशलता से प्राप्त किया जा सकता है। EQUI7 GRID 7 भूमि क्षेत्रों में दुनिया को विभाजित करता है, सभी को सटीक के न्यूनतम नुकसान के साथ समान प्रक्षेपण प्रणाली में पेश किया जाता है। EQUI7 GRID में 250 मीटर रिज़ॉल्यूशन पर वैश्विक DEM का उदाहरण देखें । यहां आप कुछ सैंपल कोड पा सकते हैं जो दिखाता है कि एसएजीए जीआईएस का उपयोग करके वैश्विक डीईएम मापदंडों को कैसे प्राप्त किया जाए। एक बार जब आप EQUI7 GRID सिस्टम में DEM मापदंडों को पूरा कर लेते हैं, तो आप सभी मानचित्रों को WGS84 longlatनिर्देशांक में वापस बदल सकते हैं और फिर GDAL का उपयोग करके एक वैश्विक मोज़ेक बना सकते हैं।


क्या आप बता सकते हैं कि यह सवाल का जवाब कैसे देता है? यदि आप ढलान की गणना के लिए समवर्ती अनुमानों का उपयोग करने का प्रस्ताव कर रहे हैं, तो कृपया ध्यान रखें कि बड़े सापेक्ष मीट्रिक विकृतियों के कारण खराब समाधान होता है जो प्रक्षेपण के केंद्र से बाहर निकलता है। यद्यपि भूमि जनता पर सात ऐसे अनुमानों को केंद्रित करना उस समस्या को कम करने में मदद करता है, फिर भी यह सबसे अच्छा विकल्प नहीं है।
whuber

कागज Bauer-Marschallinger एट अल द्वारा। (2014) बताता है कि इन अनुमानों को वैश्विक भूमि जनता का प्रतिनिधित्व करने के लिए क्यों चुना गया था (उन्हें सटीक रूप से न्यूनतम नुकसान हुआ है)। मैं मानता हूं कि किसी भी 2D प्रोजेक्शन से अंतत: विकृति पैदा होगी, लेकिन जहां तक ​​मुझे पता है, सटीक और कमोडिटी (2D बीजगणित) में नुकसान के बीच EQUI7 एक अच्छा समझौता है। यह कहने के बाद कि, हेक्सागोन्स का उपयोग फिर से वैश्विक भूमि सतहों का प्रतिनिधित्व करने के लिए किया जा रहा है (हालांकि 3 डी हेक्सागोन के साथ डीईएम विश्लेषण अभी भी एक बोझिल है)।
टॉम हेंगल

संदर्भ के लिए धन्यवाद। इसका सार बताता है कि यह एक बहुत अलग समस्या का हल करता है, "सामान्य डेटा को नियमित रैस्टर ग्रिड के प्रक्षेपण के दौरान दिखने वाले स्थानीय डेटा ओवरसम्पलिंग को कम करना।" इसका मतलब यह नहीं है कि यह अन्य उद्देश्यों के लिए अच्छा प्रदर्शन करेगा, जैसे ढलानों का अनुमान लगाना।
whuber

बेशक EQUI7 स्थानीय ढलानों का सही आकलन करने की समस्या को पूरी तरह से हल नहीं करता है, लेकिन यह संभवतः उपर्युक्त मर्केटर प्रोजेक्शन का उपयोग करने की तुलना में अधिक सुरुचिपूर्ण समाधान है। अंत में, अगर कोई सही सटीकता के साथ ढलानों का अनुमान लगाना चाहता है, तो एकमात्र विकल्प संभवत: (1) छोटे आकार की टाइलों (जैसे 100 से 100 किमी) के प्रति स्थानीय (विषुवतीय) अनुमानों का उपयोग 10-20% ओवरलैप के साथ करते हैं, जैसा कि इसमें भी उल्लेख किया गया है। वरदीन एट अल। (2007) , या (2) हेक्सागोन ग्रिड ( dggridR पैकेज ) का उपयोग करने के लिए ।
टॉम हेंगल

समस्या सटीक नहीं है - यह व्यवस्थित रूप से पक्षपाती ढलानों और पहलुओं के उत्पादन में निहित है। क्योंकि भूमध्यवर्ती अनुमान अपने केंद्रों पर उत्पन्न होने वाले भू-आकृतियों के लिए रूढ़िवादी दिशा-निर्देशों को आंशिक रूप से विकृत करते हैं, इसलिए पहलू हमेशा गलत होंगे (हालांकि उन केंद्रों के बिल्कुल समीप जहां सभी विकृति कम है) और ढलानों में त्रुटियां तेजी से दूरी के साथ बढ़ेंगी। बेशक कई स्थानीय अनुमानों का उपयोग काम करेगा, लेकिन आपके द्वारा दिए गए लालित्य के बहुत विपरीत है।
whuber

-2

ढलान उठना / दौड़ना है। संगणना उठो और संगणना भागो और तुम्हारे पास तुम्हारा जवाब है। भौगोलिक निर्देशांक के बीच की दूरी की गणना करना सरल है। यह UTM में रूपांतरण की तुलना में कम रेज़मैपलिंग त्रुटि का परिचय देगा, आदि।

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