ゾーン統計ツールの仕組み
Spatial Analyst ライセンスで利用可能。
Image Analyst ライセンスで利用可能。
ゾーン統計の動作は、別のデータセットによって定義されたゾーン内のラスター (値ラスター) のセル値に対して統計を計算します。 ゾーン別に統計を計算するツールには、ゾーン特性、ゾーン統計、ゾーン統計をテーブルに出力 があります。
ゾーン統計 ツールは、一度に 1 つの統計を計算して、ラスター出力を作成します。 この値は、そのゾーンに対応するセルのラスター出力のセル値になります。 ゾーン フィーチャがオーバーラップするゾーンを含んでいる場合、出力ラスター内のセルが 1 つの値のみを表すことができるため、1 つのゾーンに対してのみ統計が計算されます。
ゾーン統計をテーブルに出力 は、事前に定義されたサブセットを使用して 1 つ以上の統計を計算するか、すべての統計を計算して、テーブル出力を作成します。 ゾーン統計 と同様に、得られる統計は、ゾーンごとの単一の値です。 出力テーブルには、1 つのゾーンにつき 1 つのレコードが存在し、統計値が、事前に定義されたフィールドに報告されます。 ゾーン入力がフィーチャであり、オーバーラップするゾーンを含んでいる場合、すべてのゾーンに対して統計が計算され、出力がゾーンごとの個別のレコードに報告されます。
ゾーン特性 ツールは、複数の値ラスターの統計を計算して、テーブル出力を作成します。 出力テーブルには、1 つのゾーンにつき 1 つのレコードが存在し、統計値が、値ラスター名と選択した統計の種類から派生した事前に定義されたフィールドに報告されます。 ゾーン統計をテーブルに出力 ツールと同様に、このツールはオーバーラップするゾーンを処理し、出力が、ゾーンごとの個別のレコードに報告されます。
入力ゾーン レイヤーは、ラスターまたはフィーチャのいずれかであることができるゾーンの形状、値、および位置を定義します。 ゾーン演算時に、最初にフィーチャ データがラスターに変換されます。 ラスター データでは、ゾーンは、セルが隣接しているかどうかにかかわらず、すべて同じ値を含んでいるセルです。 各ゾーンは、一意の ID を持つ必要があり、ラスターである場合は、整数データ型を含む必要があります。 ゾーン入力内の一意の値の任意の整数フィールドまたは文字列フィールドを指定して、ゾーンを定義できます。
入力値ラスターには、各ゾーンの出力統計の計算に使用される値が含まれます。 整数データ型でも浮動小数点データ型でもかまいません。
下の図では、入力された値の平均がゾーンごとに識別されています:
値ラスター内のセルがラスター ゾーンとして識別される方法
統計を計算するために、ツールは、最初に、各ゾーンに含まれるすべてのセルについて、値ラスターからセル値を抽出します。 値ラスターにおけるゾーン内のセルのこの識別は、ゾーンを値ラスターの上にオーバーレイすることによって実行されます。 ゾーンおよび入力値が、両方とも同じセル サイズのラスターであり、セル同士が位置揃えされている場合、ゾーンのセルにオーバーレイする値ラスターのセル値が抽出され、統計が計算されます。
ゾーン ラスターのセル サイズ、セル配置、または空間参照のいずれかが、値ラスターと異なっている場合、ゾーン ラスターと値ラスターの間で、セルを互いに完全にオーバーレイできません。 その場合、このツールは、1 つまたは両方のラスターを調整し、セルの完全なオーバーレイを達成します。 この調整は、出力セル サイズ、出力空間参照、およびスナップ対象ラスターが内部で値ラスターによって決まるように行われます。 セル サイズは同じだが、セルが位置揃えされていない場合、値ラスターがスナップ対象ラスターとして内部で使用されます。 いずれの場合も、ゾーン演算を実行する前に内部でリサンプリングが開始されます。
ジオプロセシング環境設定で、セル サイズ、スナップ対象ラスター、出力座標系、またはこれらの組み合わせが指定された場合、これらの設定を遵守して作成された解析ウィンドウで、ゾーン演算が実行されます。 詳細については、「Spatial Analyst における解析ウィンドウの決定方法」をご参照ください。
値ラスター内のセルがフィーチャ ゾーンとして識別される方法
ゾーン演算は、基本的に、1 つがゾーン、もう 1 つが値である 2 つのラスターに対して実行されるラスター解析です。 ゾーンがフィーチャによって定義される場合、フィーチャからラスターへの内部変換が発生します。 ポリゴン ゾーンの内部変換では、ポリゴン → ラスター ツールのセルの中心メソッドを使用して、値ラスターのセル サイズおよびスナップ対象ラスターで入力をラスター化します。 これは、ラスター化のグリッドのどのセルの中心もフィーチャ ゾーンに含まれない場合、出力においてゾーンが欠落するという予期しない結果を引き起こすことがあります。 この結果は、内部のゾーン ラスターのセルの面積より小さいゾーンで発生する可能性があり、さらに大きいゾーンでも発生することがあります。
下の例では、図 (1) が、入力フィーチャ ゾーン、入力値ラスター、およびセルの中心を表しています。 入力フィーチャは 3 つのゾーン (黄色い形状) を含んでおり、次のことが当てはまります:
zone1 が個々のセルより大きい。
zone2 および zone3 が 1 つのセルより小さい。
セルの中心が zone2 から外れているが、zone3 に含まれている。
図 (2) のゾーン ラスター化プロセス中に、セルの中心が zone1 および zone2 に含まれていないため、zone3 のみがラスター化され、他の 2 つのゾーンが基本的に非表示になる。
出力でゾーンが非表示になるのを防ぐには、各ゾーンが、値ラスターの 1 つ以上のセルの中心を必ず含むようにします。 これを実行する 1 つの方法は、環境内で小さなセル サイズを指定して、作成するセルの中心を増やすことです。 デフォルトでは、解析セル サイズは値ラスターのセル サイズになります。 解析環境で、値ラスターのセル サイズより小さいセル サイズを指定した場合、上の図 (3) に示されているように、より多くのゾーンを捕捉できるようになります。 小さなセル サイズを指定すると、大規模な出力ラスターが生成されることに注意してください。 追加の詳細が入力値ラスターに実際に存在しないため、高解像度の出力は、結果が高品質に見えることに必ずしもつながりません。
値ラスターの同じセル サイズ、セルの位置揃え、および空間参照を使用してフィーチャ ゾーンがラスター ゾーンに変換された後に、ゾーンを値ラスターにオーバーレイすることによって、ゾーン内の値ラスターからのセルの抽出が実行されます。
ジオプロセシング環境設定で、セル サイズ、スナップ対象ラスター、出力座標系、またはこれらの組み合わせが指定された場合、これらの設定によって定義された解析ウィンドウで、フィーチャからラスターへの内部変換を含むゾーン演算が実行されます。 詳細については、「Spatial Analyst における解析ウィンドウの決定方法」をご参照ください。
算術統計および円周統計の計算
すべてのセル値を加算して求められた値をセル数で除算して平均を計算する方法は、標高などのデータに有効です。 一方で、傾斜方向などの周期的変化量 (0 ～ 360 度のコンパス方向) や 1 日の時間数 (0 ～ 24 時間) を示すデータの場合は、算術平均を計算しても、正確な出力が生成されません。これは、最小値と最大値が同じ数量を表しているためです。 このようなデータの場合は、円周統計を計算する必要があります。
たとえば、0 度と 360 度という 2 つのセル値の平均を計算する場合、算術平均では、計算結果が 180 度になります。 0 度と 360 度は同じコンパス方向を表しているため、この値は正確ではありません。 円形平均を計算すると、正確な統計が取得され、計算結果が 0 度になります。
円周統計の計算を指定するには、円周統計の計算 (Python では circular_calculation = "CIRCULAR") パラメーターをオンにします。 円周統計を計算する際には、周期的データを表す最小値と最大値に注意してください。 最小値は 0 と見なされます。 最大値は、円形折り返し値 (circular_wrap_value in Python) パラメーターとして指定できます。 このパラメーターのデフォルト値は 360 です。
データのタイプに応じて、正確な出力を生成するために、統計計算のタイプと円周統計に適した円形折り返し値を選択します。 使用できる円周統計オプションは、平均値、最頻値、最少頻値、標準偏差、および 種類 です。
多次元ラスターを使用してゾーン統計を計算
多次元ラスター データは、複数の回数および複数の深度または高さで、データを表します。 このタイプのデータは、一般に、大気科学、海洋学、および地球科学で使用され、プラットフォームを監視することによって観測されるか、衛星によって捕捉されるか、データがさまざまな統計手法を使用して処理、集約、または内挿される数値シミュレーション モデルから生成されます。 多次元ラスターの詳細については、「多次元ラスター データ」をご参照ください。
ゾーン統計 ツールおよび ゾーン統計をテーブルに出力 ツールは、多次元ゾーン データおよび値ラスター データを入力としてサポートします。 多次元として処理 パラメーターをオンにした場合 (Python では process_as_multidimensional パラメーターの ALL_SLICES)、ゾーン統計は多次元ラスターのすべてのスライスに対して計算されます。 多次元として処理 パラメーターをオフにした場合 (Python では CURRENT_SLICES)、現在のスライスのみが処理されます。
多次元データに対するゾーン統計解析の例としては、次のようなものが挙げられます:
気象学者は、特定の期間のハリケーンの軌跡に沿って、ハリケーンの移動および降水分布に関する知見を得る必要があります。 気象学者は、ゾーン統計 ツールで多次元処理を使用して、時間と共に変化したハリケーン ゾーンに関して、タイム スライスごとの平均降雨量を知ることができます。
環境保護活動家は、特定の河川流域に関して、過去 30 年間の毎日の最大降雨量データから、極端な出来事の分布を確認する必要があります。 パーセンタイル値のリストでのパーセンタイルの統計の種類と共に ゾーン統計をテーブルに出力 ツールを使用して、多次元として処理する場合の時系列データに関して、毎日の最大降雨量データの分布を見つけることができます。
サポートされている多次元ラスター データ タイプには、多次元ラスター レイヤー、多次元モザイク、イメージ サービス、および Esri の CRF (クラウド ラスター形式) が含まれます。
ArcGIS Pro で多次元ラスター レイヤーを追加するには、マップ タブの レイヤー グループで、データの追加 にある 多次元ラスター レイヤー オプションをクリックします。 または、多次元ラスター レイヤーの作成 ツールを使用して、ゾーン演算に適した変数を選択し、多次元ラスター レイヤーを生成します。
ゾーン統計の多次元出力
ゾーン統計 ツールが入力を多次元として処理するように指定した場合、ツールは、多次元ラスター出力を作成します。 ゾーン演算は、ゾーン ラスターのスライスと値ラスターの現在の変数のスライスの間で、スライスごとに発生します。 計算済みの統計値が多次元変数に格納され、この多次元変数の名前は、値ラスターの変数名と計算されている統計を組み合わせて作成されます。 出力変数の次元数およびスライス数は、ゾーン ラスター入力および値ラスター入力の特定の性質によって決まります。
プロパティ ウィンドウで、ラスター出力の多次元情報を調べることができます。
ArcPy で Raster オブジェクトの mdinfo プロパティを使用して、次元、次元値の数、および変数内のスライスの総数の詳細を確認することもできます。
ゾーン統計をテーブルに出力 ツールで、データを多次元として処理するように指定した場合、すべてゾーンおよびスライスに対して計算された統計を含むフラットなテーブル出力が生成されます。 このテーブルは、変数名、ディメンション名、ゾーンごとに計算された統計を示すための追加フィールドを含みます。
多次元処理が、ゾーン ラスターと値ラスターの間でスライスごとに発生するため、ゾーン統計 ツールからの出力多次元ラスター内のスライスの数、および ゾーン統計をテーブルに出力 ツールからの出力テーブル内のレコードの数は、入力ラスターのタイプおよび入力ラスター内のスライスの数によって決まります。 以下のサブセクションで、各例について説明します。
同じ次元を含む多次元ゾーンおよび値ラスター
海洋のさまざまな深さで、対応する深さでのさまざまな温度範囲にわたって最大塩分濃度を知るには、温度ゾーンを表す多次元ゾーンおよび塩分濃度を表す多次元値ラスターを使用してゾーン統計を実行する必要があります。 ゾーン演算は、値ラスターからの対応するスライスを使用して、ゾーンのスライスごとに実行されます。 出力多次元ラスターは、値ラスターと同じ数のスライスを含みます。
下の図では、ゾーン ラスターと値ラスターの両方における変数が、同じ 3 つの次元 x、y、d、およびディメンション値での同じ数のスライス d0、d1、d2 を含んでいます。 出力多次元ラスターにおける変数も、同じ 3 つの次元 x、y、d、およびディメンション値での同じ数のスライス d0、d1、d2 を含んでいます。
ゾーン統計をテーブルに出力 ツールの出力内のレコードの総数は、各スライス内のゾーン数を加算することによって決定されます。 深さ d0、d1、および d2 でのゾーン数がそれぞれ 5、4、および 3 である場合、レコードの総数は 12 (5 + 4 + 3 = 12) になります。
異なる次元を含む多次元ゾーンおよび値ラスター
ROV (遠隔操作探査機) などのアセットを配置するための適切な位置およびタイム ウィンドウは、異なる時間での ROV の可能性のある位置を表す多次元ゾーン、および異なる深さと時間での海流を表す HYCOM (Hybrid Coordinate Ocean Model) モデルの出力などの多次元値ラスターを使用してゾーン統計を実行することによって、決定できます。
ゾーン演算は、値ラスターからの各スライスを使用して、ゾーン ラスターからのスライスごとに実行されます。 出力多次元ラスター内のスライスの数は、ゾーン ラスター内のスライスの数に値ラスター内のスライスの数を掛けることによって決定されます。
下の図では、ゾーン ラスターにおける変数が、3 つの次元 x、y、d、およびディメンション値での 3 つのスライス d0、d1、d2 を含んでいます。 値ラスターにおける変数が、3 つの次元 x、y、t、およびディメンション値での 2 つのスライス t0 および t1 を含んでいます。 出力多次元ラスター内の変数も、4 つの次元 x、y、d、t を含みます。
ゾーン統計 ツールの出力内のスライスの総数は、ゾーン ラスター内の深さの数を値ラスター内の時間ステップの数で乗算することによって決定され、この場合は、6 (3 つの深さ x 2 つの時間ステップ = 6) になります。 ゾーン統計をテーブルに出力 ツールの出力内のレコードの総数は、各スライス内のゾーン数を乗算することによって決定されます。 ゾーン数が 5 である場合、この場合のレコードの総数は 30 (5 つのゾーン x 3 つの深さ x 2 つの時間ステップ =30) になります。
多次元値ラスターのみ
その年の毎日の各郡内の最高気温を知るには、毎日の気温を表す多次元値ラスターおよび郡を表すゾーン ラスターを使用してゾーン統計を実行する必要があります。 ゾーン演算は、同じゾーン ラスターを使用して、値ラスターからのスライスごとに実行されます。 出力多次元ラスターは、値ラスターと同じ数のスライスを含みます。
下の図では、値ラスターにおける変数が、3 つの次元 x、y、t、およびディメンション値での 3 つのスライス t0、t1、t2 を含んでいます。 出力多次元ラスターにおける変数も、同じ 3 つの次元 x、y、t、およびディメンション値での同じ数のスライス t0、t1、t2 を含んでいます。
ゾーン統計をテーブルに出力 ツールの出力内のレコードの総数は、ゾーン数および値ラスター内のスライス数を乗算することによって決定されます。 ゾーン数が 5 である場合、レコードの総数は 15 (5 x 3 =15) になります。
多次元値ゾーン ラスターのみ
生態学的景観計画のために、時間と共に変化する各氾濫原ゾーン カテゴリー内の 10 年間の最大降水量の平均を知るには、氾濫原ゾーンを表す多次元ゾーン ラスターおよび 10 年間の最大降水量を表す値ラスターを使用してゾーン統計を実行する必要があります。 ゾーン演算は、同じ値ラスターを使用して、ゾーン ラスターからのスライスごとに実行されます。 出力多次元ラスターは、ゾーン ラスターと同じ数のスライスを含みます。
下の図では、ゾーン ラスターにおける変数が、3 つの次元 x、y、t、およびディメンション値での 3 つのスライス t0、t1、t2 を含んでいます。 出力多次元ラスターにおける変数も、同じ 3 つの次元 x、y、t、およびディメンション値での同じ数のスライス t0、t1、t2 を含んでいます。
ゾーン統計をテーブルに出力 ツールの出力内のレコードの総数は、ゾーン数をゾーン ラスター内のスライス数で乗算することによって決定されます。 ゾーン数が 5 である場合、レコードの総数は 15 (5 x 3 =15) になります。
統計
ゾーン統計を計算するために使用可能な統計の種類を、詳細および入力例に対する各オプションの結果を示す図と共に、以下に示します。
最頻値
各ゾーン内の最も頻出する値が、そのゾーン内のすべてのセルに割り当てられます。
最頻値が複数ある場合は、最小値が採用され、その値がゾーン内のすべてのセルに割り当てられます。
例:
OutRas = ZonalStatistics(ZoneRas, "VALUE", ValRas, "MAJORITY")最頻値の数
- 各ゾーン内の最も頻出する値の数が、そのゾーン内のすべてのセルに割り当てられます。
最頻値の割合
- 各ゾーン内の最も頻出する値の数の割合が、そのゾーン内のすべてのセルに割り当てられます。
最大
- 各ゾーン内の最大値が、そのゾーン内のすべてのセルに割り当てられます。
例:
OutRas = ZonalStatistics(ZoneRas, "VALUE", ValRas, "MAXIMUM")平均
各ゾーン内の値の平均値がそのゾーン内のすべての出力セルに割り当てられます。
算術平均の式は次のとおりです:
ここでは:
x̄ は平均値
xi は観測値
N は観測数
円形平均の式は次のとおりです:
ここでは:
x̄ は円形平均の値
xi は観測値
N は観測数
Σsin xi と Σcos xi がどちらも 0 に相当する不完全なケースでは、特別な値として -1 が使用されます。この値は、円形平均が明確に定義されていないことを意味します。
例:
OutRas = ZonalStatistics(ZoneRas, "VALUE", ValRas, "MEAN")中央値
各ゾーン内の値の中央値がそのゾーン内のすべての出力セルに割り当てられます。
この統計の種類の値は、Hyndman および Fan (1996 年) の手法 Q1 を使用して計算されます。 2 つの格納された値がターゲットの中央値に等しく近い場合、2 つの値のうちの小さいが使用されます。
中央値を計算するために、ゾーン内のすべてのセルがランク付けされます。 ゾーン内に n 個のセルがあり、n が奇数の場合、ゾーン内の各セルに中央 ((n+1)/2) 番目の値が書き込まれます。 偶数のセルがある場合は、(n/2) 番目の値が出力されます。
例:
OutRas = ZonalStatistics(ZoneRas, "VALUE", ValRas, "MEDIAN")最小
- 各ゾーン内の最小値がそのゾーン内のすべてのセルに割り当てられます。
例:
OutRas = ZonalStatistics(ZoneRas, "VALUE", ValRas, "MINIMUM")最少頻値
各ゾーン内の最も頻出しない値が、そのゾーン内のすべてのセルに割り当てられます。
最少頻値が複数ある場合は、小さい方の値が採用され、その値でゾーンが埋め尽くされます。
例:
OutRas = ZonalStatistics(ZoneRas, "VALUE", ValRas, "MINORITY")最少頻値の数
- 各ゾーン内の最も頻出しない値の数が、そのゾーン内のすべてのセルに割り当てられます。
最少頻値の割合
- 各ゾーン内の最も頻出しない値の数の割合が、そのゾーン内のすべてのセルに割り当てられます。
パーセンタイル
各ゾーン内の値のパーセンタイルが、そのゾーン内のすべての出力セルに割り当てられます。
この統計の種類の値は、Hyndman および Fan (1996 年) の手法 Q1 を使用して計算されます。 2 つの格納された値がターゲットの中央値に等しく近い場合、2 つの値のうちの小さいが使用されます。
パーセンタイルを計算するために、式 R = P/100 x (n - 1) + 1 を使用して値ラスター内のすべてのセルがランク付けされます。ここで、P は目的のパーセンタイル、n はセルの数です。
例:
OutRas = ZonalStatistics(ZoneRas, "VALUE", ValRas, "PERCENTILE")範囲
各ゾーン内の最大値と最小値の差が、そのゾーン内のすべてのセルに割り当てられます。
範囲は次のように定義されます:
ゾーンの範囲 = ゾーンの最大値 - ゾーンの最小値
例:
OutRas = ZonalStatistics(ZoneRas, "VALUE", ValRas, "RANGE")標準偏差
各ゾーン内の値の標準偏差が、そのゾーン内のすべてのセルに割り当てられます。
算術標準偏差の式は次のとおりです:
ここでは:
σ は標準偏差
xi は観測値
x̄ は平均値
N は観測数
注意:
標準偏差は、サンプルから推定される (N - 1 法) のではなく、母集団全体で計算されます (N 法)。 つまり、標準偏差の計算は、Microsoft Excel の STDEV の方法ではなく、STDEVP の方法と同じです。
円周標準偏差の式は次のとおりです:
ここでは:
σ は円周標準偏差
R̄ は
度単位の n 個の角度の例では、角度 a1、a2、～ an が集計され、それぞれの角度が単位ベクトル
で表されます。この単位ベクトルは、対応する観測の方向を示します。
例:
OutRas = ZonalStatistics(ZoneRas, "VALUE", ValRas, "STD")合計
各ゾーン内のすべてのセル値の合計が、そのゾーン内のすべてのセルに割り当てられます。
出力ラスターのデータ タイプは浮動小数点数です。 これは、合計の値は非常に大きくなる傾向があり、整数値で表すことができない場合があるためです。
たとえば、1 つのゾーンに 2,500 個のロウと 2,500 個のカラムのセルがあり、各セルの値が 1,000 である場合、そのゾーンの合計値は、2,500 x 2,500 x 1,000 = 62.5 億になります。 整数値の出力が必要で、その範囲が ± 21.47 億以内の場合は、Int ツールを適用できます。
例:
OutRas = ZonalStatistics(ZoneRas, "VALUE", ValRas, "SUM")種類
- 各ゾーン内の個別値の数が、そのゾーン内のすべてのセルに割り当てられます。
例:
OutRas = ZonalStatistics(ZoneRas, "VALUE", ValRas, "VARIETY")出力データ タイプ
出力データ タイプ (整数または浮動小数) は、実行するゾーン計算と入力値ラスターのタイプの両方によって決定されます。 次の表は、出力可能なラスターのデータ タイプを示しています:
|
統計 |
値の入力タイプ |
出力 |
|---|---|---|
|
最頻値 |
Integer* |
Integer |
|
最頻値の数 |
Integer* |
Integer |
|
最頻値の割合 |
Integer* |
Float |
|
最大 |
整数、浮動小数点数 |
値の入力と同じ |
|
平均 |
整数、浮動小数点数 |
Float |
|
中央値 |
整数、浮動小数点数 |
Integer |
|
最小 |
整数、浮動小数点数 |
値の入力と同じ |
|
最少頻値 |
Integer* |
Integer |
|
最少頻値の数 |
Integer* |
Integer |
|
最少頻値の割合 |
Integer* |
Float |
|
パーセンタイル |
整数、浮動小数点数 |
Integer |
|
範囲 |
整数、浮動小数点数 |
値の入力と同じ |
|
標準偏差 |
整数、浮動小数点数 |
Float |
|
合計 |
整数、浮動小数点数 |
Float |
|
種類 |
Integer* |
Integer |
注意:
*\ 整数のみがサポートされています。
ゾーン データセット内のいずれかのセル位置が NoData である場合、出力では、その位置に NoData が割り当てられます。
参考文献
Rob J. Hyndman and Yanan Fan (1996) "Sample Quantiles in Statistical Packages" The American Statistician, Vol. 50, No. 4 (Nov., 1996), pp. 361-365.