चलो थोड़ा (बस थोड़ा) बीजगणित करते हैं।
केंद्रीय वर्ग में x मान होने दें; x_i, i = 1, .., 8 को पड़ोसी वर्गों में मानों को अनुक्रमित करने दें; और आर स्थलाकृत बीहड़ता सूचकांक होने दें। यह नुस्खा कहता है कि r ^ 2 का योग (x_i - x) ^ 2 के बराबर है। दो चीजें जिन्हें हम आसानी से गणना कर सकते हैं वे हैं (i) पड़ोस में मानों का योग, s = Sum {x_i} + x के बराबर; और (ii) मानों के वर्गों का योग, t = सम {x_i ^ 2} + x ^ 2 के बराबर। (ये मूल ग्रिड और इसके वर्ग के लिए फोकल आँकड़े हैं।)
चौकों का विस्तार
r ^ 2 = Sum {(x_i - x) ^ 2}
= सम {x_i ^ 2 + x ^ 2 - 2 * x * x_i}
= Sum {x_i ^ 2} + 8 * x ^ 2 - 2 * x * Sum {x_i}
= [Sum {x_i ^ 2} + x ^ 2] + 7 * x ^ 2 - 2 * x * [Sum {x_i} + x - x]
= t + 7 * x ^ 2 - 2 * x * [सम {x_i} + x] + 2 * x ^ 2
= t + 9 * x ^ 2 - 2 * x * s ।
उदाहरण के लिए, एक पड़ोस पर विचार करें
1 2 3
4 5 6
7 8 9
यहाँ, x = 5, s = 1 + 2 + ... + 9 = 45, और t = 1 + 4 + 9 + ... + 81 + 285. फिर
(1-5) ^ 2 + (2-5) ^ 2 + ... + (9-5) ^ 2 = 16 + 9 + 4 + 1 + 1 + 4 + 9 + 16 = 60 = r ^ 2
और बीजीय तुल्यता कहते हैं
60 = आर ^ 2 = 285 + 9 * 5 ^ 2 -2 * 5 * 45 = 285 + 225 - 450 = 60, जो जांच करता है।
कार्यप्रवाह इसलिए है:
एक डीईएम दिया।
[डीईएम] की गणना एस = फोकल राशि (3 x 3 वर्ग पड़ोस में)।
DEM2 = [DEM] * [DEM] की गणना करें।
[DEM2] की गणना t = फोकल योग (3 x 3 वर्ग पड़ोस में)।
गणना करें r2 = [t] + 9 * [DEM2] - 2 * [DEM] * [s]।
वापसी r = Sqrt ([r2])।
यह टोटो में 9 ग्रिड संचालन के होते हैं , जो सभी तेज हैं। वे आसानी से रेखापुंज कैलकुलेटर (आर्कगिस 9.3 और पहले), कमांड लाइन (सभी संस्करण), और मॉडल बिल्डर (सभी संस्करण) में आसानी से किए जाते हैं।
BTW, यह एक "औसत ऊंचाई परिवर्तन" नहीं है (क्योंकि ऊंचाई में परिवर्तन सकारात्मक और नकारात्मक हो सकता है): यह एक मूल मतलब वर्ग उन्नयन परिवर्तन है। यह http://arcscripts.esri.com/details.asp?dbid=14156 पर वर्णित "स्थलाकृतिक स्थिति सूचकांक" के बराबर नहीं है , जो (प्रलेखन के अनुसार) x - (s - x) / 8 के बराबर है। ऊपर के उदाहरण में, TPI 5 के बराबर है - (45-5) / 8 = 0 जबकि TRI, जैसा कि हमने देखा, Sqrt (60) है।