আরকজিআইএস ডেস্কটপে টপোগ্রাফিক রাগডনেস সূচক গণনা করছেন?


15

আর্কজিআইএস ডেস্কটপে আর্কিআইএনফো ওয়ার্কস্টেশন অ্যাক্সেস না করে কী কী টপোগ্রাফিক রাগডনেস সূচক গণনা করতে হয় তা কি কেউ জানেন ?

"টোগোগ্রাফিক রাগডনেস ইনডেক্স (টিআরআই) হ'ল একটি ডিজিটাল এলিভেশন গ্রিডের সংলগ্ন কোষগুলির মধ্যে উচ্চতা পার্থক্যের পরিমাণ প্রকাশ করার জন্য, এট আল। (1999) দ্বারা তৈরি একটি পরিমাপ The প্রক্রিয়াটি মূলত একটি কেন্দ্রের ঘর থেকে উচ্চতার মানগুলির পার্থক্যের গণনা করে" এবং তত্ক্ষণাত এই আটটি কোষ এটি চারপাশে আটটি উচ্চতা পার্থক্যের মানগুলির প্রত্যেককে বর্গাকার করে এগুলি সমস্ত ধনাত্মক করে তোলে এবং বর্গকে গড় করে তোলে top টপোগ্রাফিক অভদ্রতা সূচকটি এই গড়ের বর্গমূল গ্রহণ করে উত্পন্ন হয় এবং গড় উচ্চতা পরিবর্তনের সাথে মিলিত হয় কোনও গ্রিডে এবং এর আশেপাশের অঞ্চলের যে কোনও বিন্দুর মধ্যে। - জেফরি ইভান্সের একটি এমএল আর্স্ক্রিপ্ট থেকে


আরকিজিআইএস আরসস্ক্রিপ্টস.এসরি.com/details.asp?dbid=12646 এর পূর্ববর্তী ফোরামে.এসরি.com/Thread.asp?c=93&f=982&t=145448 চেক না করা কিছু আলোচনা থেকে সংস্করণ নির্ভর করে তবে অনুসন্ধানে jennessent.com

উত্তর:


18

আমি আর্কজিআইএসের বাইরে সন্ধান করার পরামর্শ দিচ্ছি) ফ্রি জিডিএল সফ্টওয়্যারটি ব্যবহার করে খুব সহজ: http://www.gdal.org/gdaldem.html

gdaldem TRI input_dem output_TRI_map

বা যদি আপনি এটিকে সাগা গিসে পছন্দ করেন: http://www.saga-gis.org/saga_modules_doc/ta_morphometry/ta_morphometry_16.html


11
+1 আমি সর্বদা আর্কজিআইএস সমস্যার নন-আর্কজিআইএস সমাধানগুলি দেখার জন্য প্রশংসা করি :-)। এটি মূলত আরকজিআইএসের প্রতি বৈরিতা নয়, নীতিগত বিষয়। একক একটি সফ্টওয়্যার সমাধানে লক হওয়া এড়ানো উচিত: এটি কেবল পেশাগতভাবেই ঝুঁকিপূর্ণ নয়, এটি বৌদ্ধিকভাবে দমবন্ধ করে।
হোয়বার

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

18

আসুন একটি ছোট (কেবলমাত্র) বীজগণিত করা যাক।

X এর মধ্য বর্গক্ষেত্রে মান হওয়া যাক; x_i, i = 1, .., 8 প্রতিবেশী স্কোয়ারগুলিতে মানগুলি সূচী করা যাক; এবং টপোগ্রাফিক অভদ্রতা সূচক হতে দিন। এই রেসিপিটি বলে যে r ^ 2 সমান (x_i - x) ^ 2 এর সমান। দুটি জিনিস যা আমরা সহজেই গণনা করতে পারি তা হ'ল (i) প্রতিবেশীর মানগুলির সমষ্টি, s = Sum {x_i; + x এর সমান; এবং (ii) টি = যোগফল {x_i ^ 2} + x ^ 2 এর সমান মানগুলির বর্গের যোগফল। (এটি মূল গ্রিড এবং এর স্কোয়ারের কেন্দ্রিয় পরিসংখ্যান)

স্কোয়ারগুলি প্রসারিত করে

r ^ 2 = যোগফল x (x_i - x) ^ 2}

= সমষ্টি {x_i ^ 2 + x ^ 2 - 2 * x * x_i

= সমষ্টি {x_i ^ 2} + 8 * x ^ 2 - 2 * x * সমষ্টি {x_i

= [সমষ্টি {x_i ^ 2} + x ^ 2] + 7 * x ^ 2 - 2 * এক্স * [সমষ্টি {x_i} + x - x]

= টি + 7 * এক্স ^ 2 - 2 * এক্স * [সমষ্টি _ x_i} + x] + 2 * x ^ 2

= t + 9 * x ^ 2 - 2 * x * গুলি

উদাহরণস্বরূপ, একটি পাড়া বিবেচনা করুন

1 2 3
4 5 6
7 8 9

এখানে, x = 5, s = 1 + 2 + ... + 9 = 45, এবং টি = 1 + 4 + 9 + ... + 81 = 285. তারপরে

(1-5) ^ 2 + (2-5) ^ 2 + ... + + (9-5) ^ 2 = 16 + 9 + 4 + 1 + 1 + 4 + 9 + 16 = 60 = আর ^ 2

এবং বীজগণিত সমতুলতা বলে

60 = আর ^ 2 = 285 + 9 * 5 ^ 2 -2 * 5 * 45 = 285 + 225 - 450 = 60, যা চেক করে।

কর্মপ্রবাহ তাই হয়:

একটি ডিইএম দেওয়া হয়েছে।

  • [DEM] এর গণনা s = ফোকাল রাশি (3 x 3 বর্গ পাড়ার বেশি)।

  • গণনা DEM2 = [DEM] * [ডিইএম]।

  • [DEM2] এর গণনা টি = ফোকাল রাশি (3 x 3 বর্গ পাড়ার বেশি)।

  • গণনা r2 = [টি] + 9 * [ডেম 2] - 2 * [ডেম] * [গুলি]।

প্রত্যাবর্তন R = sqrt ([R2])।

এটি টোটোতে 9 টি গ্রিড অপারেশন নিয়ে গঠিত দ্রুত are এগুলি সহজেই রাস্টার ক্যালকুলেটর (আর্কজিআইএস 9.3 এবং তার আগের), কমান্ড লাইন (সমস্ত সংস্করণ) এবং মডেল বিল্ডার (সমস্ত সংস্করণ) এ সহজেই সঞ্চালিত হয়।

বিটিডাব্লু, এটি একটি "গড় উচ্চতর পরিবর্তন" নয় (কারণ উন্নয়নের পরিবর্তনগুলি ইতিবাচক এবং নেতিবাচক হতে পারে): এটি মূলের বর্গক্ষেত্রের উচ্চতা পরিবর্তন। এটি http://arcscriptts.esri.com/details.asp?dbid=14156 এ বর্ণিত "টপোগ্রাফিক অবস্থান সূচক" এর সমান নয় , যা ডকুমেন্টেশন অনুসারে x - (s - x) / 8 এর সমান। উপরের উদাহরণে, টিপিআই সমান 5 - (45-5) / 8 = 0 যেখানে টিআরআই, যেমন আমরা দেখেছি, স্কয়ার্ট (60)।


1
ধন্যবাদ বিল। কোনও সরঞ্জাম বা অপারেশন কীভাবে কাজ করে তার সুনির্দিষ্ট বিবরণ দেখে আমি প্রশংসা করি। এ থেকে উপযুক্ত সময় এবং বৌদ্ধিক শক্তির উপযুক্ত বিনিয়োগের সাথে যে কেউ হাতের কাছে থাকা সরঞ্জাম ব্যবহার করে এই কাজটি চালিয়ে নিতে একটি নতুন যন্ত্রপাতি তৈরি করতে পারে। এটি এর মতো তথ্য যা জিআইএস.এসকে দীর্ঘ দুরত্বের কাজে লাগবে।
ম্যাট উইলকি

1
+1 দুর্দান্ত ব্যাখ্যা। আমি অনুমান করি এটির অর্থ একটি খাড়া তবে মসৃণ পৃষ্ঠের একটি সমতল তবে গোঁফযুক্ত পৃষ্ঠের চেয়ে টিআরআই বেশি থাকতে পারে।
কर्क কুইকেনডাল

1
পছন্দ করুন আপনি যদি চান তবে "আপেক্ষিক" অভদ্রতার সূচক পেতে স্থানীয় opeালের প্রভাব সরিয়ে দেওয়ার বিভিন্ন উপায় রয়েছে। যদিও আমি বিস্তারিত কাজ করেন নি, আমি বিশ্বাস করি যে কিছু সার্বজনীন একাধিক বিয়োগ (গ * ক) ^ 2 R2 থেকে - যেখানে cellsize এবং একটি ঢাল (হয় বৃদ্ধি / চালানোর যেমন না একটি কোণের হিসাবে অথবা শতাংশ) - কৌশলটি করা উচিত।
হোবার

@ সর্বদা আপনার উত্তরগুলিতে বিস্ময়কর পরিমাণ জ্ঞান থাকে !! কেবল একটি প্রশ্ন দয়া করে: এর অর্থ কি রাস্টারের একেবারে প্রান্তে অবস্থিত কোষগুলির টিআরআই গণনা করা সম্ভব নয়? তাদের চারপাশে প্রতিবেশী কক্ষ দ্বারা ঘিরে না থাকার কারণে?
মার্কো

1
@ মার্কো টিআরআই এমনকি সীমানা কোষেও অনুমান করা যায়। প্রশ্নটিতে ইঙ্গিত হিসাবে, আমি এখানে 9 দ্বারা যে মানগুলি দিয়েছি তা ভাগ করে একটি গড় হিসাবে বরং গড় হিসাবে প্রকাশ করা উচিত সীমানা কোষগুলিতে সূত্র এবং ডিনোমিনেটরের মান "9" দ্বারা প্রতিস্থাপন করা দরকার তাদের 3X3 আশেপাশের নন-নাল মানগুলির সংখ্যা: প্রান্ত কক্ষগুলির জন্য 6, কোণার কক্ষগুলির জন্য 4। মূল মানগুলির সূচক গ্রিডের কেন্দ্রবিন্দু থেকে এই জাতীয় মানগুলির একটি গ্রিড পাওয়া যায় (এটি অন্য কোনও নোডাটা কোষে 1 এর এবং অন্য কোথাও 0 রয়েছে)। সূত্রগুলিতে ধ্রুব "9" এর জায়গায় সেই গ্রিডটি ব্যবহার করুন।
whuber

3

রিলে এট আল।, (1999) টিআরআই সংক্ষিপ্ত স্কোয়ার বিচ্যুতির স্কোয়ার রুট। এটি আনসেলড বৈকল্পিকের খুব কাছাকাছি। আপনি যদি রিলির টিআরআই বাস্তবায়ন করতে চান তবে দয়া করে @ হুবারের দ্বারা বর্ণিত পদ্ধতিটি অনুসরণ করুন (@ ব্যবহারকারী 3338736 দ্বারা প্রদত্ত পদ্ধতিটি উইন্ডোতে সর্বাধিক মেট্রিককে সাধারণীকরণ করেছে এবং কোষের প্রকরণের মাধ্যমে সেলটি উপস্থাপন করে না)।

আমাদের জিওমরফমিট্রি এবং গ্রেডিয়েন্ট মেট্রিক্স আরকজিআইএস টুলবক্সে আমার টিআরআইয়ের একটি প্রকরণ রয়েছে যা একটি নির্দিষ্ট উইন্ডোর বৈকল্পিক। আমি এটি আরও নমনীয় এবং ন্যায়সঙ্গত বলে মনে করি। রাগোসিটি এবং বিচ্ছিন্নকরণ সহ আরও কিছু পৃষ্ঠতল কনফিগারেশন মেট্রিক রয়েছে।


ধন্যবাদ জেফরি কিছু কারণে ফায়ারফক্সে শিরোনাম ব্যতীত পৃষ্ঠাটি ফাঁকা, ভেবেছিল এটি ক্রোমে ঠিক আছে; আমার এক্সটেনশনগুলির মধ্যে একটি হতে পারে। আমি কমপক্ষে প্রতিবেদন করতে পেরে খুশি যে স্ক্রিপ্টগুলি 10.2.2 এ পরিবর্তিত হয়ে কাজ করে (আমি যেভাবেই পরীক্ষা করেছি)।
ম্যাট উইলকি

1

-সম্পাদনা: নীচের তথ্যটি ভুল। সঠিক প্রক্রিয়াটি ব্যাখ্যা করে পোস্টটি দেখুন whuber .....

টিআরআই (রাইল 1999) এবং টিপিআই (জেনেস 2002) একই রকম, তবে আলাদা।

আরকিজিআইএস ১০.০ ব্যবহার করে টিআরআই এবং টিপিআই গণনা করার জন্য ...

পদক্ষেপ 1: একটি ডিইএম থেকে 2 টি নতুন রাস্টার ডেটাসেট তৈরি করতে ফোকাল স্ট্যাটিস্টিক্স সরঞ্জামটি ব্যবহার করুন।

রাস্টার 1 "MAX") প্রতিবেশী: আয়তক্ষেত্র, উচ্চতা: 3, প্রস্থ: 3, ইউনিট: সেল, পরিসংখ্যানের ধরণ: সর্বাধিক

রাস্টার 2 "এমআইএন") প্রতিবেশী: আয়তক্ষেত্র, উচ্চতা: 3, প্রস্থ: 3, ইউনিট: সেল, পরিসংখ্যানের ধরণ: সর্বনিম্ন

পদক্ষেপ 2: আপনার সবেমাত্র তৈরি করা 2 জন রাস্টার ডেটাসেটে নিম্নলিখিত ফাংশনগুলি সম্পাদন করতে রাস্টার ক্যালকুলেটরটি ব্যবহার করুন।

টিআরআইয়ের জন্য: স্কয়াররুট (অ্যাবস ((স্কোয়ার ("% ম্যাক্স%" ") - স্কয়ার ("% এমআইএন% "")))

টিপিআইয়ের জন্য: ("% ইনপুট ডেম%" - "% মিন%") / ("% ম্যাক্স%" - "% এমআইএন%")

আমি টিআরআইয়ের জন্য নির্মিত একটি মডেল থেকে পাইথন কোডটি রফতানি করা এখানে ...

# -*- coding: utf-8 -*-
# ---------------------------------------------------------------------------
# script.py
# Created on: 2014-03-06 08:56:13.00000
#   (generated by ArcGIS/ModelBuilder)
# Usage: script <Input_raster> <TRI_Raster> 
# Description: 
# ---------------------------------------------------------------------------

# Import arcpy module
import arcpy

# Check out any necessary licenses
arcpy.CheckOutExtension("spatial")

# Script arguments
Input_raster = arcpy.GetParameterAsText(0)

TRI_Raster = arcpy.GetParameterAsText(1)
if TRI_Raster == '#' or not TRI_Raster:
    TRI_Raster = "C:\\Users\\Documents\\ArcGIS\\Default.gdb\\rastercalc1" # provide a default value if unspecified

# Local variables:
MIN = Input_raster
MAX = Input_raster

# Process: 3x3Max
arcpy.gp.FocalStatistics_sa(Input_raster, MAX, "Rectangle 3 3 CELL", "MAXIMUM", "DATA")

# Process: 3x3Min
arcpy.gp.FocalStatistics_sa(Input_raster, MIN, "Rectangle 3 3 CELL", "MINIMUM", "DATA")

# Process: Raster Calculator
arcpy.gp.RasterCalculator_sa("SquareRoot(Abs((Square(\"%MAX%\") - Square(\"%MIN%\"))))", TRI_Raster)

এটি প্রশ্নের মধ্যে বর্ণিত টিআরআই নয়। বাস্তবে, এটি একেবারে "অভদ্রতা" পরিমাপ করার কথা ভাবা যায় না, কারণ আপনি যখন কেবল উল্লম্ব ডেটামটি স্থানান্তর করেন তখন এটি পরিবর্তিত হয়। উদাহরণস্বরূপ, মানগুলির সাথে আপনার একটি 3x3 পাড়ার টিআরআই (1,2, ..., 9) স্কয়ার্ট (9 ^ 2-1 ^ 2) = 8.9 হবে, তবে মানগুলিতে 100 যোগ করে (যা কেবলমাত্র ড্যাটামকে ছাড়াই পরিবর্তন করে) পৃষ্ঠের আকৃতিটি পরিবর্তন করা) sqrt দেয় (109 ^ 2-101 ^ 2) = 41।
হোবার

0

এটি টপোগ্রাফিক পজিশন সূচকটির মতো দেখতে খুব সুন্দর লাগে, এমন একটি প্রক্রিয়া যা আমি আমার এক প্রকল্পের জন্য সম্প্রতি ব্যবহার করেছি। ইএসআরআই সমর্থন পৃষ্ঠায় একটি আর্কস্ক্রিপ্ট রয়েছে, ইএসআরআই রিসোর্স সেন্টার পৃষ্ঠায় একটি টোগোগ্রাফি টুলবক্স এবং জেনেস এন্টারপ্রাইজ পৃষ্ঠাতে প্রক্রিয়া সম্পর্কিত আরও কিছু তথ্য রয়েছে ।


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