ลองทำพีชคณิตนิดหน่อย
ให้ x เป็นค่าในจตุรัสกลาง ให้ x_i, i = 1, .. , 8 จัดทำดัชนีค่าในพื้นที่ใกล้เคียง และปล่อยให้ r เป็นดัชนีความทนทานของภูมิประเทศ สูตรนี้บอกว่า r ^ 2 เท่ากับผลรวมของ (x_i - x) ^ 2 สองสิ่งที่เราสามารถคำนวณได้อย่างง่ายดายคือ (i) ผลรวมของค่าในย่านใกล้เคียงเท่ากับ s = Sum {x_i} + x; และ (ii) ผลรวมของกำลังสองของค่าเท่ากับ t = Sum {x_i ^ 2} + x ^ 2 (นี่คือสถิติโฟกัสสำหรับกริดดั้งเดิมและสแควร์)
การขยายกำลังสองให้
r ^ 2 = ผลรวม {(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 * [รวม {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 + 1 + 4 + 9 + 16 = 60 = r ^ 2
และความเท่าเทียมกันของพีชคณิตบอกว่า
60 = r ^ 2 = 285 + 9 * 5 ^ 2 -2 * 5 * 45 = 285 + 225 - 450 = 60 ซึ่งตรวจสอบ
เวิร์กโฟลว์จึงเป็น:
ได้รับ DEM
คำนวณ s = ผลรวมของจุดโฟกัส (พื้นที่ละแวกใกล้เคียง 3 x 3 ตาราง) ของ [DEM]
คำนวณ DEM2 = [DEM] * [DEM]
คำนวณ t = ผลรวมโฟกัส (พื้นที่ใกล้เคียง 3 x 3 ตร.) ของ [DEM2]
คำนวณ r2 = [t] + 9 * [DEM2] - 2 * [DEM] * [s]
ส่งคืน r = Sqrt ([r2])
ประกอบด้วยการทำงานของกริด 9 รายการในโตโต้ซึ่งทั้งหมดนั้นรวดเร็ว พวกมันจะถูกนำไปใช้ในเครื่องคิดเลขแบบแรสเตอร์ (ArcGIS 9.3 และรุ่นก่อนหน้า) บรรทัดคำสั่ง (ทุกรุ่น) และตัวสร้างแบบจำลอง (ทุกรุ่น)
BTW นี่ไม่ใช่ "การเปลี่ยนแปลงระดับความสูงโดยเฉลี่ย" (เนื่องจากการเปลี่ยนแปลงระดับความสูงสามารถเป็นค่าบวกและค่าลบ): มันเป็นค่าเฉลี่ยของการเปลี่ยนแปลงระดับความสูงแบบจัตุรัส มันเป็นเรื่องที่ไม่เท่ากับ "ดัชนีภูมิประเทศตำแหน่ง" อธิบายไว้ที่http://arcscripts.esri.com/details.asp?dbid=14156ซึ่ง (ตามเอกสาร) เท่ากับ x - (s - x) / 8 ในตัวอย่างข้างต้น TPI เท่ากับ 5 - (45-5) / 8 = 0 ในขณะที่ TRI ที่เราเห็นคือ Sqrt (60)