点密度（Point Density） (Spatial Analyst)
使用法
-
[検索範囲] パラメータの値が大きいほど、単純化された密度のラスタが作成されます。値が小さいほど、詳細を示すラスタが作成されます。
-
近傍内にあるポイントだけが、密度の計算時に考慮されます。特定セルの近傍内にポイントがない場合、そのセルには NoData が割り当てられます。
-
面積単位の係数単位がポイント間の距離に比べて小さい場合、出力ラスタの値が小さくなる可能性があります。大きな値を得るには、大きな単位の面積単位の係数（平方キロメートル対平方メートルなど）を使用します。
-
出力ラスタの値は、必ず浮動小数になります。
構文
| パラメータ | 説明 | データ タイプ |
in_point_features |
密度を計算する入力ポイント フィーチャ。 | Feature Layer |
population_field |
各ポイントの人口値を示すフィールド。[Population フィールド] は、連続サーフェスの計算に使用する個数や数量です。 population フィールドの値は、整数値または浮動小数点値です。 フィールドのオプションとデフォルトの動作は以下のとおりです。
| Field |
cell_size （オプション） | 出力ラスタ データセットのセル サイズ。 環境設定で特定の値が指定されている場合、その値を使用します。設定されていない場合、セル サイズは、出力空間参照において、出力の範囲の幅または高さ（どちらか短い方）を 250 で割った値になります。 | Analysis Cell Size |
neighborhood （オプション） |
密度値の計算に使用する各セル周辺のエリアの形状を指定します。 これは、近傍解析クラスです。 近傍解析クラスは次の 4 種類です。NbrAnnulus、NbrCircle、NbrRectangle、NbrWedge。 クラスの形式と説明:
デフォルトは [NbrCircle] です。半径は、出力空間参照において、出力範囲の幅または高さ（どちらか短い方）を 30 で割った値になります。 | Neighborhood |
area_unit_scale_factor （オプション） |
出力密度値の望ましい面積単位。 デフォルトの単位は、出力空間参照の投影法の距離単位に基づいて選択されます。密度の出力を変換する場合は、これを適切な単位に変更できます。線密度の値は長さと面積の両方の単位を変換します。 たとえば、入力単位がメートルの場合、デフォルトの出力の面積密度単位はポイント フィーチャでは平方キロメートル、ポリライン フィーチャではキロメートル/平方キロメートルになります。 入力フィーチャの単位に基づくデフォルトの密度単位は以下のとおりです。
| String |
リターン
| 名前 | 説明 | データ タイプ |
| out_raster |
出力ポイント密度ラスタ。 常に浮動小数点ラスタです。 | Raster |
コードのサンプル
次の例では、ポイント シェープ ファイルから密度ラスタを計算しています。
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
pdensOut = PointDensity("rec_sites.shp", "NONE", 60, NbrCircle(2500, "MAP"))
pdensOut.save("C:/sapyexamples/output/pointdensity")
次の例では、ポイント シェープ ファイルから密度ラスタを計算しています。
# Name: PointDensity_Ex_02.py
# Description: Calculates a magnitude per unit area from point
# features that fall within a neighborhood around each cell.
# 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
inFeatures = "rec_sites.shp"
populationField = "NONE"
cellSize = 60
# Create the Neighborhood Object
radius = 2500
myNbrCirc = NbrCircle(radius, "MAP")
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute PointDensity
outPdens = PointDensity(inFeatures, populationField, cellSize,
myNbrCirc, "SQUARE_KILOMETERS")
# Save the output
outPdens.save("C:/sapyexamples/output/outpdens")