ভূখণ্ড বক্রতা গণনা কিভাবে?


29

আমার একটি এলিভেশন রাস্টার রয়েছে এবং আমি কিছু বেসিক টেরিন ভেরিয়েবলগুলি পেতে চাই, যেমন ,াল, রাগাদ্বল্য, দিক বা বক্রতা।

সমস্যাটি হ'ল ভূখণ্ড বিশ্লেষকের সাথে আমার বক্রতা গণনা করার সম্ভাবনা আছে বলে মনে হয় না:

এখানে চিত্র বর্ণনা লিখুন

কেউ আমাকে বলতে পারেন কেন এই সূচকটি তালিকায় নেই? অন্যথায় আমি এটি কীভাবে গণনা করতে পারি?

উত্তর:


20

বক্রতা গণনা করার জন্য একটি জটিল অঞ্চলভিত্তিক ডেরাইভেটিভ, আপনি যে সমীকরণটি ব্যবহার করেন তা আপনার ইনপুট ডেটার রেজোলিউশনের উপর নির্ভর করে, কারণ আপনাকে নিশ্চিত করতে হবে যে আপনি যে বক্রতা ফলাফলগুলি গণনা করেছেন তার তথ্যগুলিতে গোলমাল থেকে আলাদা করা যায়।

হাই রেজোলিউশন লিজার ডেটাতে বক্রতা গণনার উপর সম্প্রতি অনেক গবেষণা করা হয়েছে যা দেখিয়েছে যে প্রায় 2 বা 3 মিটার রেজোলিউশনে একটি স্কেলিং বিরতি বিদ্যমান এবং এই বিন্দুটির উপরে আরও বিভিন্ন অ্যালগোরিদম (যা আমি তেমন পরিচিত নই) ব্যবহার করা দরকার । টপোগ্রাফিক বক্রতা গণনা সম্পর্কে সর্বোত্তম তথ্য সম্ভবত হুরস্ট এট আল 2012 এবং এর উল্লেখগুলি থেকে এসেছে ।

Opeালু এবং দিকের মতো বক্রতা গণনার মূল নীতিটি হ'ল উচ্চতা পৃষ্ঠের উপর দিয়ে চলমান উইন্ডোটি পাস করা এবং উচ্চতা মানগুলিকে একটি 6 টি বহুবচনীয় ফাংশনে ফিট করা, যার গুণফলগুলি কেন্দ্রের opeাল, দিক এবং বক্রতা অর্জন করবে চলমান উইন্ডোর সেল

আর্কজিআইএস একটি 3x3 অনুসন্ধান উইন্ডো ব্যবহার করে যা কেবলমাত্র উদ্ভিদবিহীন অঞ্চলে ভাল ফলাফল দেয়, যা লোকেরা এই সীমাবদ্ধতা সম্পর্কে অবগত না হলে এই সরঞ্জামটি মোটামুটি অকেজো করে তোলে, এটি কিউজিআইএসে উপস্থিত না হওয়ার কারণেই এটি প্রস্তাব দিতে পারে।

গণিতগুলি মূলত ইভান্স (1980) সালে উদ্ভূত হয়েছিল এবং ভৌগলিক তথ্য সিস্টেমের নীতিমালা (অ্যামাজন লিঙ্ক) এর কয়েকটি পৃষ্ঠায় সরলীকৃত হয়েছিল যা আমি এই স্তরের বিশ্লেষণের একটি মৌলিক স্তরের একটি ভাল গাইড হিসাবে সুপারিশ করতে পারি।

কোনও ডেমের বক্রতা গণনা করার একটি উপায় হ'ল ডিইএমকে একটি এসিআই রাস্টার হিসাবে রূপান্তর করা, এটি একটি অলস অ্যারেতে পড়ুন এবং তারপরে ডেটাগুলির মধ্য দিয়ে প্রবাহিত একটি চলমান উইন্ডোতে বহুবর্ষীয় ফিটিং সম্পাদন করা। এটি করা মোটামুটি সহজ, তবে এটি কার্যকর করতে খুব ধীর এবং অপরিশোধিত পরিমাণের যথেষ্ট পরিমাণের প্রয়োজন (এই ধরণের ক্রিয়াকলাপগুলি প্রায়শই তাদের গতি বাড়ানোর জন্য সি ++ এ পোর্ট করা হয়)।

কিউজিআইএস-এ অপারেশন করতে আপনি জিআরএসএস প্লাগইন 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
নোট করুন যে গ্রাসে আপনি 'r.param.scale' ব্যবহার করতে পারেন যা উইন্ডো অপারেটরের জন্য ব্যবহারকারী-সংজ্ঞায়িত আকারের সাথে বক্ররেখার মতো ভূখণ্ডের ডেরাইভেটিভগুলি গণনা করে (কেবলমাত্র ডিফল্ট 3x3 নয়)।
ব্যবহারকারী 1013346

12

বক্রতা গণনা করার জন্য রাস্টার বিশ্লেষণের ESRI- র সংস্করণটি QGIS এর জন্য একটি প্লাগইন বিকাশ করতে সহায়ক হতে পারে।

প্রতিটি কক্ষের জন্য, ফর্মের একটি চতুর্থ-আদেশের বহুপদী:

এখানে চিত্র বর্ণনা লিখুন

জেড = আক্সিয়া + বিক্সিয়া + স্যাক্সি² + ডিএক্স² + আইআই + এফসি + জিএক্স + হাই + আই

একটি 3x3 উইন্ডো দ্বারা গঠিত একটি পৃষ্ঠের সাথে ফিট করে। এ, বি, সি এবং আরও সহগের গুণকগুলি এই পৃষ্ঠ থেকে গণনা করা হয়।

ডায়াগ্রামে প্রদর্শিত প্রতিটি কোণের জন্য সহগ এবং উচ্চ স্তরের নয় টি মানের মধ্যে সম্পর্কগুলি নিম্নরূপ: বক্ররেখা মানগুলি ডায়াগ্রাম বক্ররেখা মানগুলি চিত্র

এ = [(জেড 1 + জেড 3 + জেড 7 + জেড 9) / 4 - (জেড 2 + জেড 4 + জেড 6 + জেড 8) / 2 + জেড 5] / এল 4

বি = [(জেড 1 + জেড 3 - জেড 7 - জেড 9) / 4 - (জেড 2 - জেড 8) / 2] / এল 3

সি = [(-Z1 + জেড 3 - জেড 7 + জেড 9) / 4 + (জেড 4 - জেড 6)] / 2] / এল 3

ডি = [(জেড 4 + জেড 6) / 2 - জেড 5] / এল 2

ই = [(জেড 2 + জেড 8) / 2 - জেড 5] / এল 2

এফ = (-জেড 1 + জেড 3 + জেড 7 - জেড 9) / 4 এল 2

জি = (-জেড 4 + জেড 6) / 2 এল

এইচ = (জেড 2 - জেড 8) / 2 এল

আই = জেড 5

বক্রতা সরঞ্জামের আউটপুটটি পৃষ্ঠের দ্বিতীয় ডেরাইভেটিভ - উদাহরণস্বরূপ, opeালের opeাল — যেমন:

বক্রতা = -2 (ডি + ই) * 100

সম্পূর্ণ তথ্য এবং উত্স:

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


2
এটি গণিতগুলির একটি দুর্দান্ত সংক্ষিপ্তসার, তবে আপনার ডেটার 3 মিটারের নিচে রেজোলিউশন থাকলে কোয়ালিটির কোনও বিবেচনামূলক ক্ষতি ছাড়াই 6 টার্মের বহুবচন ব্যবহার সহজ করা যায়।
সাগরিভ

5
কারণ ইএসআরআই "বক্রতা" কেবলমাত্র ডি + ই নির্ভর করে, বাকি সমস্ত অপ্রয়োজনীয়। কম্পিউটিং -2(D+E)*100যেমন (100/L2)*(3*Z5 - [Z2+Z4+Z6+Z8+Z5])শো মূল মান তিনবার যেমন এই মান প্রাপ্ত করার কিভাবে, Z5, বিয়োগ একটি পাড়া সমষ্টি Z2+Z4+Z6+Z8+Z5(ক ব্যাসার্ধ 1 টি চেনাশোনা ব্যবহার করে), আর সবাই rescaled 100/L2। এটি কেবল তিনটি সাধারণ গ্রিড অপারেশন।
whuber

@ শুভ: এক্স এবং ওয়াইজ কি কেউ আমাকে বলতে পারবেন? এবং জেড কি বক্রতার মান উপস্থাপন করে? তাহলে -2 (ডি + ই) * 100 কি? এটা কি সম্ভব যে আমার কাছে যদি ভ্যাক্টর ফাইল থাকে যাতে বহু বহুভুজ থাকে এবং আমি আবহাওয়া জানতে চাইলে বহুভুজের অভ্যন্তরের অঞ্চল অবতল বা উত্তল (পরিকল্পনা এবং প্রোফাইল বক্রতা) হয়?
হার্ডিক

জেড = উচ্চতার মূল্য এবং -2 (D + কে ই) * 100 = বক্রতার মূল্য en.wikipedia.org/wiki/Curvature (Osculating_Circle)
Mapperz

6

বক্রতা সাগা'র মডিউল 'টেরেইন বিশ্লেষণ - মরফোমেট্রি ---> opeাল, দিক, বক্রতা' ব্যবহার করে গণনা করা যেতে পারে

এই অ্যালগরিদমের কোনওটির ভিত্তিতে গণনা করা যেতে পারে:

  • সর্বাধিক opeাল (ট্র্যাভিস এট আল। 1975)
  • সর্বাধিক ত্রিভুজ opeাল (টারবোটন 1997)
  • সর্বনিম্ন স্কোয়ারগুলি লাগানো প্লেন (হর্ন 1981, কোস্টা-ক্যাব্রাল এবং বার্গেস 1996)
  • ফিট 2. ডিগ্রি পলিনম (বাউর, রোহডেনবুর্গ, বর্ক 1985)
  • ফিট 2 ডিগ্রি পলিনম (হেরডিজেন এবং বেরান 1982)
  • ফিট 2 ডিগ্রি পলিনম (জেভেনবার্গেন এবং থর্ন 1987)
  • ফিট 3. ডিগ্রি পলিনম (হারালিক 1983)

3

ল্যান্ডস্পারফ এটি করতে পারে You আপনি উইন্ডোজের আকার নির্ধারণ করতে পারেন (3 * 3, 5 * 5, 7 * 7, 11 * 11, ...) তবে এটি অবশ্যই বিজোড় সংখ্যা হতে হবে। মাল্টিস্কেল বিশ্লেষণ। আপনি opeাল, দিক এবং বক্রতার স্কেল নির্ভরতা বিবেচনা করতে পারেন। http://www.landserf.org/ ল্যান্ডস্পার্ফ নির্দিষ্ট পূর্বনির্ধারিত উইন্ডোতে সহ চতুর্ভুজ বহুমাত্রিক সমীকরণ ফিট করবে, তবে আপনি যদি 50 * 50 এর মতো বড় উইন্ডোজের আকারটি সংজ্ঞায়িত করেন তবে এটি দীর্ঘ সময় নিতে পারে। এটি আপনার রাস্টার আকার এবং উইন্ডোজ আকারের উপর নির্ভর করে যা আপনি সংজ্ঞায়িত করেছেন। জো উড তাঁর পিএইচডি থিসিসের জন্য ল্যান্ডসার্ফ লিখেছিলেন। এটি জাভাতে লেখা আছে।

  • প্রোফাইল বক্রতা
  • পরিকল্পনা বক্রতা
  • দ্রাঘিমাংশ বক্রতা
  • ক্রস বিভাগীয় বক্রতা
  • গড় বক্রতা
  • ন্যূনতম বক্রতা
  • সর্বাধিক বক্রতা

ল্যান্ডস্পার্পে বিভিন্ন স্কেলে গণনা করা যায়


2

আপনি নিখরচায় সাগা জিআইএস (http://sourceforge.net/apps/trac/saga-gis/wiki) বা টিএএস (http://www.uoguelph.ca/~hydrogeo/TAS/index.html) ব্যবহার করে দেখতে পারেন।


এই প্যাকেজগুলি বক্রতা গণনা করতে কীভাবে পরিচালিত হয় সে সম্পর্কে আপনি আরও বিশদ সরবরাহ করতে পারেন? আপনার পোস্ট করা লিঙ্কগুলি দুটি জিআইএস প্যাকেজের সাধারণ লিঙ্ক এবং জিজ্ঞাসা করা প্রশ্নের সাথে সরাসরি প্রাসঙ্গিকতা নেই।
সাগরিভ

1
খুব সংক্ষিপ্ত উত্তরের জন্য দুঃখিত। সাগায় আপনি মডিউল টেরেইন বিশ্লেষণ - যৌগিক বিশ্লেষণ -> স্ট্যান্ডার্ড টেরেইন বিশ্লেষণগুলি পেতে পারেন। সেখানে আপনি বক্রতা, প্রোফাইল বক্রতা, পরিকল্পনা বক্রতা গণনা করতে পারেন (সফ্টওয়্যার সাহায্যের ভিতরে কোনও বিশদ বিবরণ নেই)। আপনি এসএজি, ফ্লাট, ... থেকে সহজ সাগায় ডেটা আমদানি করতে পারেন ...
রক

@ এসগ্রিভে টিএএস-এ আপনি প্রোফাইল, পরিকল্পনা এবং স্পর্শকাতর বক্রতা গণনা করতে পারেন (মেনু: টেরিন বিশ্লেষণ -> প্রাথমিক ভূখণ্ডের বৈশিষ্ট্য -> সারফেস ডেরাইভেটিভস)। আবার, এটি কীভাবে কাজ করে তা খুব সামান্য বর্ণনা রয়েছে।
রক
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.