複数の抽出値 → ポイント（Extract Multi Values to Points） (Spatial Analyst)
サマリ
1 つ以上のラスタから、ポイント フィーチャクラスで指定された位置にあるセル値を抽出し、その値をポイント フィーチャクラスの属性テーブルに記録します。
使用法
-
入力には、任意のラスタの組み合わせ（シングル バンドまたはマルチバンド）を指定できます。
-
各入力ラスタのセル値が抽出され、各入力ラスタのセル値を含む新しいフィールドが入力ポイント フィーチャクラスに追加されます。
-
デフォルトでは、出力フィールド名は、入力ラスタの名前から作成されます。または、ラスタ値を格納する各フィールドに一意の名前を付けることができます。
-
入力がマルチバンド ラスタの場合、すべてのバンドのフィールドが追加され、出力フィールドの名前にバンド数を表す「b1_、b2_、…bn」という接頭辞が追加されます。
-
内挿オプションは、ラスタから値を取得する方法を決定します。デフォルト オプションでは、サンプリング中のセルの中央にある値を使用します。内挿オプションでは、共一次内挿法を使用してセル中心の値が内挿されます。
-
シェープファイル形式のフィールドには、最大 10 文字という制限があります。デフォルトでは、入力シェープファイルに追加された出力フィールドは切り捨てられ、一意になります。これにより、名前が長くて非常に似ている場合は、入力ラスタ間の区別が難しくなる可能性があります。この場合、フィーチャをファイル ジオデータベースに変換することをお勧めします。
-
値ラスタの NoData のセルは、NULL 値になります。シェープファイルの場合、NULL はサポートされていないため、値は 0（ゼロ）になります。
注意:このツールを使用すると、入力データが変更されます。詳しい説明および不適切なデータ変更を防ぐための方法については、「出力を伴わないツール」をご参照ください。
構文
| パラメータ | 説明 | データ タイプ |
in_point_features |
ラスタ値に追加する入力ポイント フィーチャ。 | Feature Layer |
in_rasters [[Raster, {Output Field Name}],...] |
入力ポイント フィーチャの位置に基づいて抽出する入力ラスタの値。 オプションで、ラスタ値を格納するフィールドの名前を指定できます。デフォルトでは、入力ラスタのデータセット名に基づいて、一意のフィールド名が作成されます。 | Value Table |
bilinear_interpolate_values （オプション） |
内挿を使用するかどうかを指定します。
| Boolean |
コードのサンプル
複数のラスタからポイント シェープファイル フィーチャクラスの属性にセル値を抽出します。
import arcpy
from arcpy.sa import *
from arcpy import env
env.workspace = "c:/sapyexamples/data"
ExtractMultiValuesToPoints("observers.shp", [["elevation", "ELEV"],
["costraster", "COST"], ["flowdir", "DIR"]], "NONE")
内挿を使用して、複数のラスタからポイント シェープファイル フィーチャクラスの属性にセル値を抽出します。
# Name: ExtractMultiValuesToPoints_Ex_02.py
# Description: Extracts the cells of multiple rasters as attributes in
# an output point feature class. This example takes a multiband IMG
# and two GRID files as input.
# 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 = "poi.shp"
inRasterList = [["doqq.img", "doqqval"], ["redstd", "focalstd"],
["redmin", "focalmin"]]
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute ExtractValuesToPoints
ExtractMultiValuesToPoints(inPointFeatures, inRasterList, "BILINEAR")