コマンドラインArcInfo Workstationにアクセスせずに、ArcGIS Desktopで地形の耐久性指数を計算する方法を知っている人はいますか？
「地形的耐久性指数（TRI）は、デジタル標高グリッドの隣接セル間の標高差の量を表すために、Riley、et al。（1999）によって開発された測定です。このプロセスは、基本的に中央セルからの標高値の差を計算します。そして、8つの標高差値のそれぞれを2乗してすべて正の値にし、平方を平均します。次に、この平均の平方根を取得することにより、地形的耐久性指数が導出され、平均標高変化に対応します。グリッド上の任意の点とその周囲の領域との間。」- ジェフリー・エヴァンスによるamlアークスクリプトから
回答:
ArcGISの外を見ることをお勧めします）無料のgdalソフトウェアを使用して非常に簡単：http ://www.gdal.org/gdaldem.html
gdaldem TRI input_dem output_TRI_map
または、saga gisでそれを好む場合：http : //www.saga-gis.org/saga_modules_doc/ta_morphometry/ta_morphometry_16.html
少し（ほんの少し）代数をしましょう。
xを中央の正方形の値とします。x_i、i = 1、..、8で、隣接する正方形の値にインデックスを付けます。rを地形の耐久性指数とします。このレシピでは、r ^ 2は（x_i-x）^ 2の合計に等しいと言います。簡単に計算できる2つのことは、（i）s = Sum {x_i} + xに等しい近傍の値の合計です。（ii）t = Sum {x_i ^ 2} + x ^ 2に等しい値の二乗和。（これらは元のグリッドとその正方形の焦点統計です。）
正方形を展開すると
r ^ 2 = Sum {（x_i-x）^ 2}
= Sum {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 * [Sum {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 = r ^ 2 = 285 + 9 * 5 ^ 2 -2 * 5 * 45 = 285 + 225-450 = 60、チェックします。
したがって、ワークフローは次のとおりです。
DEMが与えられた。
s = [DEM]の焦点和（3 x 3平方近傍）を計算します。
DEM2 = [DEM] * [DEM]を計算します。
t = [DEM2]の焦点和（3 x 3平方近傍）を計算します。
r2 = [t] + 9 * [DEM2]-2 * [DEM] * [s]を計算します。
r = Sqrt（[r2]）を返します。
これはtotoの9つのグリッド操作で構成され、すべてが高速です。これらは、ラスター計算機（ArcGIS 9.3以前）、コマンドライン（すべてのバージョン）、およびモデルビルダー（すべてのバージョン）で簡単に実行できます。
ところで、これは「平均標高変化」ではありません（標高の変化は正と負になる可能性があるためです）。これは二乗平均平方根の標高の変化です。それはしないで説明した「地形位置インデックス」に等しいhttp://arcscripts.esri.com/details.asp?dbid=14156（S - - X）/ 8×に等しい（マニュアルに従って）。上記の例では、TPIは5-（45-5）/ 8 = 0ですが、TRIはSqrt（60）です。
Riley et al。、（1999）TRIは、合計平方偏差の平方根です。これは、スケールなしの分散に非常に近いです。RileyのTRIの実装が必要な場合は、@ whuberで概説されている方法に従ってください（@ user3338736が提供する方法は、ウィンドウ内のメトリックを最大化し、セルのバリエーションごとのセルを表しません）。
Geomorphometry＆Gradient Metrics ArcGIS Toolboxには、指定されたウィンドウの分散であるTRIのバリエーションがあります。これはより柔軟で正当なものだと思います。また、しわや解離など、その他の表面構成のメトリックもあります。
-編集：以下の情報は正しくありません。正しいプロセスを説明するwhuberの投稿をご覧ください。
TRI（Riley 1999）とTPI（Jenness 2002）は似ていますが、異なります。
ArcGIS 10.xを使用してTRIおよびTPIを計算するには...
ステップ1： Focal Statisticsツールを使用して、DEMから2つの新しいラスターデータセットを作成します。
ラスター1 "MAX"）近傍：長方形、高さ：3、幅：3、単位：セル、統計タイプ：最大
ラスター2 "MIN"）近傍：長方形、高さ：3、幅：3、単位：セル、統計タイプ：最小
ステップ2： Raster Calculatorを使用して、作成した2つのラスターデータセットに対して次の機能を実行します。
TRIの場合：SquareRoot（Abs（（Square（ "％MAX％"）-Square（ "％MIN％"））））
TPIの場合：（ "％Input DEM％"-"％MIN％"）/（ "％MAX％"-"％MIN％"）
以下は、TRI用に作成したモデルからエクスポートされたサンプルPythonコードです。
# -*- 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)