इलाके की वक्रता की गणना कैसे करें?


29

मेरे पास एक ऊंचाई है और मैं ढलान, असभ्यता, पहलू या वक्रता जैसे कुछ बुनियादी इलाके चर प्राप्त करना चाहूंगा।

समस्या यह है कि मुझे लगता है कि इलाके के विश्लेषक के साथ वक्रता की गणना करने की संभावना नहीं है:

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

क्या कोई मुझे बता सकता है कि यह संकेतक सूची में क्यों नहीं है? मैं इसे अन्यथा कैसे गणना कर सकता हूं?

जवाबों:


20

वक्रता गणना करने के लिए एक जटिल भू-भाग व्युत्पन्न है, जो समीकरण आप उपयोग करते हैं वह आपके इनपुट डेटा के रिज़ॉल्यूशन पर निर्भर करता है, क्योंकि आपको यह सुनिश्चित करना है कि आपके द्वारा गणना की जाने वाली वक्रता परिणाम डेटा में शोर से अलग हो सकती है।

उच्च संकल्प LiDAR डेटा पर वक्रता पर हाल ही में बहुत सारे शोध किए गए हैं जिसमें पता चला है कि स्केलिंग ब्रेक लगभग 2 या 3 मीटर रिज़ॉल्यूशन पर मौजूद है और इस बिंदु के ऊपर और अधिक अलग-अलग एल्गोरिदम (जिनके साथ मैं परिचित नहीं हूं) का उपयोग करने की आवश्यकता है । स्थलाकृतिक वक्रता की गणना के बारे में सबसे अच्छी जानकारी संभवतः हर्स्ट एट अल 2012 और उसमें दिए गए संदर्भों से आती है।

ढलान और पहलू के साथ वक्रता गणना का मूल सिद्धांत, ऊँचाई की सतह पर एक चलती हुई खिड़की को पास करना और ऊँचाई के मानों को 6 शब्द के बहुपद समारोह में फिट करना है, जिसके गुणांक केंद्र के ढलान, पहलू और वक्रता को उत्पन्न करेंगे। चलती खिड़की का सेल।

आर्कगिस एक 3x3 खोज खिड़की का उपयोग करता है जो केवल वनस्पति से पूरी तरह से रहित क्षेत्रों में अच्छे परिणाम देगा, जो उपकरण को काफी बेकार बना देता है जब तक कि लोग इस सीमा से अवगत नहीं होते हैं, यह सुझाव दे सकता है कि यह क्यूजीआईएस में मौजूद क्यों नहीं है।

गणित मूल रूप से (मुझे लगता है) इवांस (1980) में व्युत्पन्न किया गया था और भौगोलिक सूचना प्रणाली (अमेज़ॅन लिंक) के सिद्धांतों में कुछ पन्नों में सरल किया गया था, जिसे मैं एक बुनियादी स्तर पर इस तरह के इलाके विश्लेषण के लिए एक अच्छे मार्गदर्शक के रूप में सुझा सकता हूं।

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

QGIS में ऑपरेशन करने के लिए आप GRASS प्लगइन r.slope.aspect का उपयोग कर सकते हैं जो कि 3x3 चिप विंडो द्वारा भी सीमित है।

मुझे लगता है कि यह सरल उत्तर नहीं है जिसके लिए आप निस्संदेह उम्मीद कर रहे थे, लेकिन मुझे आशा है कि आप समझते हैं कि वक्रता एक सार्थक तरीके से प्राप्त करना जटिल है। शुभकामनाएं।

Evans, I. S. (1980), An integrated system of terrain analysis and slope mapping, Z. Geomorphol., 36, 274–295.

सभी विकास के लिए धन्यवाद! हालांकि, मैं इस तथ्य से चिंतित था कि क्यूजीआईएस के रेखापुंज इलाके विश्लेषक में "वक्रता" विकल्प सूची से गायब है। सामान्य है? मैंने QGIS 1.8 के अपने संस्करण को फिर से सुनिश्चित करने के लिए पुनः स्थापित किया, लेकिन यह अभी भी नहीं है: - /
डेमियन

1
ध्यान दें कि GRASS में आप 'r.param.scale' का उपयोग कर सकते हैं जो विंडो ऑपरेटर (न केवल डिफ़ॉल्ट 3x3) के लिए उपयोगकर्ता द्वारा परिभाषित आकार के साथ, टेरामीटर जैसे टेरा डेरिवेटिव की गणना करता है।
user1013346

12

वक्रता की गणना के लिए रैस्टर एनालिसिस का ESRI का संस्करण QGIS के लिए एक प्लगइन विकसित करने में मददगार हो सकता है।

प्रत्येक कोशिका के लिए, फार्म का एक चौथा क्रम बहुपद:

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

Z = Ax Cy = + Bx²y + Cxy D + Dx² + Ey F + Fxy + Gx + Hy + I

एक 3x3 खिड़की से बना सतह पर फिट है। गुणांक ए, बी, सी, और इसी तरह, इस सतह से गणना की जाती है।

गुणांक और प्रत्येक कोशिका के लिए उन्नयन के नौ मानों के बीच के संबंध, जो चित्र पर दिखाए गए हैं, निम्नानुसार हैं: वक्रता मान आरेख वक्रता मान आरेख

A = [(Z1 + Z3 + Z7 + Z9) / 4 - (Z2 + Z4 + Z6 + Z8) / 2 + Z5] / L4

B = [(Z1 + Z3 - Z7 - Z9) / 4 - (Z2 - Z8) / 2] / 3

C = [(-Z1 + Z3 - Z7 + Z9) / 4 + (Z4 - Z6)] / 2] / L3

D = [(Z4 + Z6) / 2 - Z5] / L2

E = [(Z2 + Z8) / 2 - Z5] / L2

F = (-Z1 + Z3 + Z7 - Z9) / 4L2

जी = (-Z4 + Z6) / 2L

एच = (जेड 2 - जेड 8) / 2 एल

म = Z5

वक्रता उपकरण का आउटपुट सतह का दूसरा व्युत्पन्न है - उदाहरण के लिए, ढलान का ढलान - जैसे:

वक्रता = -2 (D + E) * 100

पूरी जानकारी और स्रोत:

http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//00q90000000t000000


2
यह गणित का एक बड़ा सारांश है, लेकिन गुणवत्ता में कोई नुकसान नहीं होने के साथ 6 शब्द बहुपद का उपयोग करने के लिए इसे सरल बनाया जा सकता है, यदि आपके डेटा का रिज़ॉल्यूशन 3 मीटर से कम है।
12 sgrieve

5
क्योंकि ESRI "वक्रता" केवल D + E पर निर्भर करता है, शेष सभी अनावश्यक है। कम्प्यूटिंग के -2(D+E)*100रूप में (100/L2)*(3*Z5 - [Z2+Z4+Z6+Z8+Z5])दिखाता है कि इस मूल्य को तीन बार मूल मान के रूप में कैसे प्राप्त किया जाए Z5, एक पड़ोस राशि Z2+Z4+Z6+Z8+Z5(एक त्रिज्या 1 सर्कल का उपयोग करके) घटाकर , सभी द्वारा rescaled 100/L2। यह सिर्फ तीन साधारण ग्रिड संचालन है।
whuber

@ वाउचर: क्या कोई मुझे बता सकता है कि x और y क्या है? और Z वक्रता के मान को दर्शाता है? फिर क्या -2 (डी + ई) * 100 है? क्या यह संभव है कि अगर मेरे पास इतने सारे बहुभुजों वाली एक वेक्टर फ़ाइल है और मैं जानना चाहता हूं कि बहुभुज के अंदर का क्षेत्र अवतल या उत्तल (योजना और प्रोफ़ाइल वक्रता) है?
हार्दिक

जेड = ऊंचाई मूल्य और -2 (डी + ई) * 100 = वक्रता मूल्य en.wikipedia.org/wiki/Curvature (Osculating_Circle)
Mapperz

6

वक्रता की गणना एसएजीए के मॉड्यूल 'टेरेन विश्लेषण - मॉर्फोमेट्री ---> ढलान, पहलू, वक्रता' के उपयोग से की जा सकती है

गणना इनमें से किसी एक एल्गोरिथम के आधार पर की जा सकती है:

  • अधिकतम ढलान (ट्रैविस एट अल। 1975)
  • अधिकतम त्रिभुज ढलान (तारबोन 1997)
  • कम से कम वर्ग के मैदान (हॉर्न 1981, कोस्टा-कैब्रल और बर्गेस 1996)
  • फिट 2. डेज्री पोलिनोम (बाउर, रोहडेनबर्ग, बोर्क 1985)
  • फिट 2. डेज्री पोलिनॉम (हीरडेजेन और बेरन 1982)
  • फ़िट 2. डेग्री पोलिनोम (ज़ेवेंबरगेन और थॉर्न 1987)
  • फिट 3. डेज्री पोलिनॉम (हारलिक 1983)

3

LandSerf ऐसा कर सकता है। आप खिड़कियों के आकार (3 * 3, 5 * 5, 7 * 7, 11 * 11, ...) को परिभाषित कर सकते हैं, लेकिन यह विषम संख्या में होना चाहिए। मल्टीस्केल विश्लेषण। आप ढलान, पहलू और वक्रता के पैमाने पर निर्भरता पर विचार कर सकते हैं। http://www.landserf.org/ LandSerf विशिष्ट पूर्वनिर्धारित खिड़कियों के लिए सह द्विघात बहुपद समीकरण में फिट होगा, लेकिन यदि आप बड़ी खिड़कियों के आकार को 50 * 50 की तरह परिभाषित करते हैं तो इसमें लंबा समय लगता है। यह आपके रेखापुंज आकार और खिड़कियों के आकार पर निर्भर करता है जिसे आप परिभाषित करते हैं। जो वुड ने अपनी पीएचडी थीसिस के लिए लैंड्सर्फ लिखा। यह जावा में लिखा है।

  • प्रोफ़ाइल वक्रता
  • योजना वक्रता
  • देशांतर वक्रता
  • पार अनुभागीय वक्रता
  • मतलब वक्रता
  • न्यूनतम वक्रता
  • अधिकतम वक्रता

LandSerf में विभिन्न पैमाने में गणना की जा सकती है


2

आप मुफ्त सागा जीआईएस (http://sourceforge.net/apps/trac/saga-gis/wiki) या TAS (http://www.uoguelph.ca/~hydrogeo/TAS/index.html) भी आज़मा सकते हैं।


क्या आप इस बात पर अधिक विस्तार प्रदान कर सकते हैं कि वक्रता की गणना करने के लिए ये पैकेज कैसे संचालित होते हैं? आपके द्वारा पोस्ट किए गए लिंक दो जीआईएस पैकेजों के सामान्य लिंक हैं, और पूछे गए प्रश्न की कोई प्रत्यक्ष प्रासंगिकता नहीं है।
sgrieve

1
बहुत कम जवाब के लिए क्षमा करें। SAGA में आप मॉड्यूल टेरेन एनालिसिस - कंपाउंड एनालिसिस -> स्टैंडर्ड टेरेन एनालिसिस पा सकते हैं। वहां आप वक्रता, प्रोफ़ाइल वक्रता, योजना वक्रता की गणना कर सकते हैं (सॉफ्टवेयर सहायता के अंदर कोई विस्तृत विवरण नहीं है)। आप asGA, flt से आसान सागा में डेटा आयात कर सकते हैं ...
Rok

@sgrieve TAS में आप प्रोफ़ाइल, योजना और स्पर्शरेखा वक्रता की गणना कर सकते हैं (मेनू: मैदानी विश्लेषण -> प्राथमिक मैथुन गुण -> भूतल परिधीय)। फिर, इसमें बहुत कम विवरण है कि यह कैसे काम करता है।
रोक्क
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.