ポイントで抽出（Extract by Points） (Spatial Analyst)

サマリ

座標ポイントに基づいてラスタのセルを抽出します。

使用法

• 選択されていないセル位置には NoData 値が代入されます。

• 入力としてマルチバンド ラスタを指定すると、新しいマルチバンド ラスタが出力として作成されます。入力マルチバンド ラスタの各バンドは個別に解析されます。

  デフォルトの出力形式は Esri Grid スタックです。Esri Grid スタックには、数字で始まる名前、スペースを含む名前、9 文字を超える名前は使用できません。

• 入力が、4 バンド以上のマルチバンド ラスタから作成したレイヤである場合、抽出操作ではレイヤによって読み込まれた（シンボル表示された）バンドだけが考慮されます。そのため、出力マルチバンド ラスタには、入力レイヤの表示で使用されるバンドに該当する 3 つのバンドしか出力されません。

• 入力ラスタが整数の場合、出力ラスタも整数になります。入力のいずれかが浮動小数点値である場合、出力は浮動小数点値になります。

構文

ExtractByPoints (in_raster, points, {extraction_area})

パラメータ	説明	データ タイプ
in_raster	セルを抽出する入力ラスタ。	Raster Layer
points [point,...]	ポイント クラス オブジェクトの Python リストは、値をラスタから抽出する位置を表します。 ポイント オブジェクトは、X、Y 座標の組のリストで指定されます。オブジェクトの形式: [point(x1,y1), point(x2,y2),...] ポイントは、in_raster と同じマップ単位です。	Point
extraction_area （オプション）	指定したポイント位置に基づいてセルを抽出するか（内部）、ポイント位置の外のセルを抽出するか（外部）識別します。 INSIDE — 選択したポイント内にあるセルを出力ラスタに書き込むことを指定するキーワード。外部にあるすべてのセルは、出力ラスタでは値が NoData になります。 OUTSIDE — 入力ポイントの外部のセルを選択し、出力ラスタに書き込むことを指定するキーワード。	String

リターン

名前	説明	データ タイプ
out_raster	入力ラスタから抽出したセル値が格納される出力ラスタ。	Raster

コードのサンプル

ExtractByPoints（ポイントで抽出）の例 1（Python ウィンドウ）

次の例では、指定したポイント座標に基づいてラスタからセルを抽出しています。

import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
pointList = [arcpy.Point(743050, 4321275), arcpy.Point(743100, 4321200),
             arcpy.Point(734500,4322000)]
outPointExtract = ExtractByPoints("soil", pointList,"INSIDE")
outPointExtract.save("c:/sapyexamples/output/pntextract")

ExtractByPoints（ポイントで抽出）の例 2（スタンドアロン スクリプト）

次の例では、指定したポイント座標に基づいてラスタからセルを抽出しています。

# Name: ExtractByPoints_Ex_02.py
# Description: Extracts the cells of a raster based on a set of points.
# 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
inRaster = "soil"
pointList = [arcpy.Point(743050, 4321275), 
             arcpy.Point(743100, 4321200), 
             arcpy.Point(743500, 4322000),
             arcpy.Point(742900, 4321800)]

# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")

# Execute ExtractByPoints
outPointExtract = ExtractByPoints("soil", pointList,"INSIDE")

# Save the output 
outPointExtract.save("c:/sapyexamples/output/pntext")

環境

セル サイズ, 現在のワークスペース, マスク, 出力データの座標系, 範囲, テンポラリ ワークスペース, スナップ対象ラスタ

ライセンス情報

ArcView: 必須 Spatial Analyst

ArcEditor: 必須 Spatial Analyst

ArcInfo: 必須 Spatial Analyst