ポイント統計（Point Statistics） (Spatial Analyst)
使用法
-
フィールドが整数の場合、使用できるオーバーレイ統計情報は、[MEAN]、[MAJORITY]、[MAXIMUM]、[MEDIAN]、[MINIMUM]、[MINORITY]、[RANGE]、[STD]、[SUM]、および [VARIETY] です。フィールドが浮動小数点の場合、使用できる統計情報は、[MEAN]、[MAXIMUM]、[MINIMUM]、[RANGE]、[STD]、[SUM] のみです。
-
[統計情報の種類（オプション）] が [MAJORITY]、[MAXIMUM]、[MEDIAN]、[MINIMUM]、[MINORITY]、[RANGE] および [SUM] の場合、ラスタの出力データ タイプは入力フィールド タイプと同じになります。[統計情報の種類（オプション）] が [MEAN] と [STD] の場合、出力ラスタは必ず浮動小数になります。[VARIETY] の場合、出力ラスタは必ず整数になります。
-
ラスタ セルの近傍内にポイントがない場合、[VARIETY] の統計情報は値 0 を割り当てます。他の統計情報については、NoData が割り当てられます。
構文
| パラメータ | 説明 | データ タイプ |
in_point_features |
各出力セルの近傍内の統計情報を計算する入力ポイント フィーチャ。 入力は、ポイントとマルチポイント フィーチャクラスが使用できます。 | Feature Layer |
field |
フィールドは、入力ポイント フィーチャの任意の数値フィールドにできます。 入力フィーチャが Z 値を含む場合、Shape フィールドを使用できます。 | Field |
cell_size (オプション) |
出力ラスタ データセットのセル サイズ 環境設定で特定の値が指定されている場合、その値を使用します。設定されていない場合は、出力空間参照において、入力フィーチャ データセットの範囲の幅または高さ（どちらか短い方）を 250 で割った値になります。 | Analysis Cell Size |
neighborhood (オプション) |
Neighborhood クラスは、統計情報の計算に使用する各入力ポイント周辺のエリアの形状を指定します。 利用可能な近傍のタイプは、NbrAnnulus、NbrCircle、NbrRectangle、NbrWedge です。 近傍の形式は以下のとおりです。
デフォルトの近傍は、幅と高さが 3 x 3 セルの正方形の NbrRectangle です。 | Neighborhood |
statistics_type (オプション) |
計算する統計情報の種類。 計算は、各出力ラスタ セルの近傍内にあるポイント入力の指定したフィールド値に対して実行されます。
| String |
戻り値
| 名前 | 説明 | データ タイプ |
| out_raster |
出力ポイント統計ラスタ。 | Raster |
コードのサンプル
次の例では、各出力ラスタ セルの円形の近傍内にある入力シェープファイル ポイント フィーチャの統計情報（合計値）を求めています。
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outPointStats = PointStatistics("ca_ozone_pts.shp", "OZONE", 500,
NbrCircle(10000, "MAP"), "SUM")
outPointStats.save("C:/sapyexamples/output/pointstatsout")
次の例では、各出力ラスタ セルの円形の近傍内にある入力シェープファイル ポイント フィーチャの統計情報（平均値）を求めています。
# Name: PointStatistics_Ex_02.py
# Description: Calculates a statistic on points over a specified
# neighborhood outputting a raster.
# Requirements: Spatial Analyst Extension
# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "C:/sapyexamples/data"
# Set local variables
inPointFeatures = "ca_ozone_pts.shp"
field = "OZONE"
cellSize = 500
neighborhood = NbrCircle(6000, "MAP")
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute PointStatistics
outPointStatistics = PointStatistics(inPointFeatures, field, cellSize,
neighborhood, "MEAN")
# Save the output
outPointStatistics.save("C:/sapyexamples/output/pointstatout")