Chúng ta hãy làm một ít (chỉ một chút) đại số.
Gọi x là giá trị trong hình vuông trung tâm; Đặt x_i, i = 1, .., 8 lập chỉ mục các giá trị trong các ô vuông lân cận; và gọi r là chỉ số độ chắc chắn địa hình. Công thức này cho biết r ^ 2 bằng tổng của (x_i - x) ^ 2. Hai thứ chúng ta có thể tính toán dễ dàng là (i) tổng các giá trị trong vùng lân cận, bằng s = Sum {x_i} + x; và (ii) tổng bình phương của các giá trị, bằng t = Sum {x_i ^ 2} + x ^ 2. (Đây là các thống kê đầu mối cho lưới ban đầu và cho hình vuông của nó.)
Mở rộng hình vuông cho
r ^ 2 = Tổng {(x_i - x) ^ 2}
= Tổng {x_i ^ 2 + x ^ 2 - 2 * x * x_i}
= Tổng {x_i ^ 2} + 8 * x ^ 2 - 2 * x * Tổng {x_i}
= [Tổng {x_i ^ 2} + x ^ 2] + 7 * x ^ 2 - 2 * x * [Tổng {x_i} + x - x]
= t + 7 * x ^ 2 - 2 * x * [Tổng {x_i} + x] + 2 * x ^ 2
= t + 9 * x ^ 2 - 2 * x * s .
Ví dụ, hãy xem xét một khu phố
1 2 3
4 5 6
7 8 9
Ở đây, x = 5, s = 1 + 2 + ... + 9 = 45 và t = 1 + 4 + 9 + ... + 81 = 285. Sau đó
(1-5) ^ 2 + (2-5) ^ 2 + ... + (9-5) ^ 2 = 16 + 9 + 4 + 1 + 1 + 4 + 9 + 16 = 60 = r ^ 2
và tương đương đại số nói
60 = r ^ 2 = 285 + 9 * 5 ^ 2 -2 * 5 * 45 = 285 + 225 - 450 = 60, sẽ kiểm tra.
Các quy trình làm việc do đó là:
Đưa ra một DEM.
Tính s = Tổng tiêu điểm (trên 3 x 3 vùng lân cận) của [DEM].
Tính toán DEM2 = [DEM] * [DEM].
Tính t = Tổng tiêu điểm (trên 3 x 3 vùng lân cận) của [DEM2].
Tính r2 = [t] + 9 * [DEM2] - 2 * [DEM] * [s].
Trở về r = Sqrt ([r2]).
Điều này bao gồm 9 hoạt động lưới trong toto , tất cả đều nhanh chóng. Chúng dễ dàng được thực hiện trong máy tính raster (ArcGIS 9.3 trở về trước), dòng lệnh (tất cả các phiên bản) và Model Builder (tất cả các phiên bản).
BTW, đây không phải là "thay đổi độ cao trung bình" (vì thay đổi độ cao có thể là dương và âm): đó là thay đổi độ cao bình phương gốc. Nó không bằng "chỉ số vị trí địa hình" được mô tả tại http://arcscripts.esri.com/details.asp?dbid=14156 , mà (theo tài liệu) bằng x - (s - x) / 8. Trong ví dụ trên, TPI bằng 5 - (45-5) / 8 = 0 trong khi TRI, như chúng ta đã thấy, là Sqrt (60).