Czy ktoś wie, jak obliczyć Indeks odporności topograficznej w ArcGIS Desktop bez dostępu do linii poleceń ArcInfo Workstation?
„Wskaźnik wytrzymałości topograficznej (TRI) to pomiar opracowany przez Riley i wsp. (1999) w celu wyrażenia wielkości różnicy wysokości między sąsiednimi komórkami cyfrowej siatki wysokości. Proces ten zasadniczo oblicza różnicę wartości wysokości z komórki środkowej i ośmioma komórkami bezpośrednio go otaczającymi. Następnie obciąża każdą z ośmiu wartości różnicy wysokości, aby wszystkie były dodatnie, i uśrednia kwadraty. Wskaźnik chropowatości topograficznej jest następnie uzyskiwany na podstawie pierwiastka kwadratowego z tej średniej i odpowiada średniej zmianie wysokości między dowolnym punktem na siatce a otaczającym obszarem. ” - z arcriptu amla autorstwa Jeffreya Evansa
źródło
Odpowiedzi:
Polecam szukać poza ArcGIS) Bardzo łatwe w użyciu darmowego oprogramowania gdal: http://www.gdal.org/gdaldem.html
Lub jeśli wolisz w saga gis: http://www.saga-gis.org/saga_modules_doc/ta_morphometry/ta_morphometry_16.html
źródło
Zróbmy małą (tylko trochę) algebrę.
Niech x będzie wartością w centralnym kwadracie; niech x_i, i = 1, .., 8 indeksuje wartości w sąsiednich kwadratach; i niech r będzie topograficznym wskaźnikiem odporności. Ten przepis mówi, że r ^ 2 równa się sumie (x_i - x) ^ 2. Dwie rzeczy, które możemy łatwo obliczyć, to (i) suma wartości w sąsiedztwie, równa s = Suma {x_i} + x; oraz (ii) suma kwadratów wartości, równa t = Suma {x_i ^ 2} + x ^ 2. (Są to ogniskowe statystyki oryginalnej siatki i jej kwadratu).
Rozszerzanie kwadratów daje
r ^ 2 = Suma {(x_i - x) ^ 2}
= Suma {x_i ^ 2 + x ^ 2 - 2 * x * x_i}
= Suma {x_i ^ 2} + 8 * x ^ 2 - 2 * x * Suma {x_i}
= [Suma {x_i ^ 2} + x ^ 2] + 7 * x ^ 2 - 2 * x * [Suma {x_i} + x - x]
= t + 7 * x ^ 2 - 2 * x * [Suma {x_i} + x] + 2 * x ^ 2
= t + 9 * x ^ 2 - 2 * x * s .
Na przykład rozważ okolicę
Tutaj x = 5, s = 1 + 2 + ... + 9 = 45, a t = 1 + 4 + 9 + ... + 81 = 285. Następnie
(1-5) ^ 2 + (2-5) ^ 2 + ... + (9-5) ^ 2 = 16 + 9 + 4 + 1 + 1 + 4 + 9 + 16 = 60 = r ^ 2
a algebraiczna równoważność mówi
60 = r ^ 2 = 285 + 9 * 5 ^ 2 -2 * 5 * 45 = 285 + 225 - 450 = 60, co sprawdza.
Przepływ pracy jest zatem następujący:
Biorąc pod uwagę DEM.
Oblicz s = suma ogniskowa (ponad 3 x 3 dzielnice kwadratowe) [DEM].
Oblicz DEM2 = [DEM] * [DEM].
Oblicz t = suma ogniskowa (ponad 3 x 3 dzielnice kwadratowe) z [DEM2].
Oblicz r2 = [t] + 9 * [DEM2] - 2 * [DEM] * [s].
Zwraca r = Sqrt ([r2]).
Składa się z 9 operacji siatki w całości , z których wszystkie są szybkie. Są one łatwo przeprowadzane w kalkulatorze rastrowym (ArcGIS 9.3 i wcześniejszych), w wierszu poleceń (wszystkie wersje) i w Konstruktorze modeli (wszystkie wersje).
BTW, nie jest to „średnia zmiana wysokości” (ponieważ zmiany wysokości mogą być dodatnie i ujemne): jest to podstawowa zmiana średniej wysokości kwadratu. To jest , nie są równe z „topograficznych wskaźnika położenia” opisanej w http://arcscripts.esri.com/details.asp?dbid=14156 , który (zgodnie z dokumentacją) jest równa x - y - (x) / 8. W powyższym przykładzie TPI wynosi 5 - (45-5) / 8 = 0, podczas gdy TRI, jak widzieliśmy, to Sqrt (60).
źródło
Riley i wsp., (1999) TRI jest pierwiastkiem kwadratowym z sumowanych kwadratowych odchyleń. Jest to bardzo zbliżone do nieskalowanej wariancji. Jeśli chcesz zaimplementować TRI Riley, postępuj zgodnie z metodologią przedstawioną przez @whuber (metodologia podana przez @ user3338736 uogólniła metrykę na maksimum w oknie i nie reprezentuje komórki według wariantu komórki).
Mam odmianę TRI w naszym Geomorfometrii i Gradient Metrics ArcGIS Toolbox, która jest wariantem określonego okna. Uważam to za bardziej elastyczne i uzasadnione. Istnieją również inne parametry konfiguracji powierzchni, w tym chropowatość i rozwarstwienie.
źródło
-Edytuj: poniższe informacje są nieprawidłowe. Proszę zobaczyć post whuber wyjaśniający prawidłowy proces .....
TRI (Riley 1999) i TPI (Jenness 2002) są podobne, ale różne.
Aby obliczyć TRI i TPI za pomocą ArcGIS 10.x ...
Oto przykładowy kod Pythona wyeksportowany z modelu zbudowanego dla TRI ....
źródło
Brzmi to bardzo podobnie do wskaźnika pozycji topograficznej, procesu, którego ostatnio użyłem w jednym z moich projektów. Na stronie wsparcia ESRI znajduje się ArcScript , przybornik Topografia na stronie Centrum zasobów ESRI oraz kilka dodatkowych informacji na temat procesu na stronie Jenness Enterprises .
źródło