Spatial Analyst のライセンスで利用可能。
概要
頂点を指定したポリゴンに基づいて、ラスターのセルを抽出します。
図
使用法
入力ラスターからの追加属性 (存在する場合) は、そのまま出力ラスター属性テーブルに引き継がれます。記録されるプロパティに応じて、一部の属性値は再計算が必要な場合があります。
[入力ラスター] (Python では in_raster) としてマルチバンド ラスターが指定された場合、すべてのバンドが使用されます。
マルチバンド ラスターから選択したバンドを処理するには、まず [コンポジット バンド (Composite Bands)] ツールを使用して、特定のバンドで構成された新しいラスター データセットを作成します。その後、結果を [入力ラスター] (Python では in_raster) として使用します。
デフォルトの出力形式は、ジオデータベース ラスターです。出力形式として Esri Grid スタックを指定した場合、スタック名を数字で開始すること、名前にスペースを使用すること、名前の長さが 9 文字を超えることはできません。
一連の X,Y 座標の代わりに、フィーチャクラスのポリゴンに基づいて抽出するには、[マスクで抽出 (Extract By Mask)] ツールを使用します。
ポリゴン オブジェクトは、1 つのパートまたは複数のパート (ポリゴン クラス) から構成されます。後者の場合、1 つのポリゴン輪郭を描けるように、すべてのポリゴン パートが連続している必要があります。不連続な複数のパートを含むポリゴン フィーチャに基づいて抽出するには、[マスクで抽出 (Extract By Mask)] ツールを使用します。
セルがポリゴンの内部か外部かを判定するには、セルの中心が使用されます。中心がポリゴンのアーク内にある場合、そのセルの一部がポリゴンの外にあっても、そのセル全体が内部にあると見なされます。
ポリゴンの頂点の上限は、1,000 です。ポリゴンの頂点は、時計回りの順番で入力する必要があります。最初と最後の頂点が一致していない場合、ポリゴンを閉じることができません。これは、複数のポリゴンを使用する場合に特に重要です。その場合、各ポリゴンの最後のポイントが最初の頂点と同じでないときに、最初の頂点に直接接続され、ポリゴンが自動的に閉じられる可能性があります。ただし、これにより想定とは異なる結果が生じる可能性があるため、この手法を利用する場合は注意してください。
ポリゴンのアークは他のアークと交差してもかまいませんが、複雑なポリゴンは推奨しません。
選択されていないセル位置には NoData 値が代入されます。
入力ラスターが整数の場合、出力ラスターも整数になります。入力のいずれかが浮動小数点値である場合、出力は浮動小数点値になります。
このツールに適用されるジオプロセシング環境の詳細については、「解析環境と Spatial Analyst」をご参照ください。
構文
ExtractByPolygon(in_raster, polygon, {extraction_area})| パラメーター | 説明 | データ タイプ |
in_raster | セルを抽出する入力ラスター。 | Raster Layer |
polygon [point,...] | 抽出する入力ラスターのエリアを定義するポリゴン。 それぞれのポリゴン パーツは、Point クラスによって定義される頂点のリストです。ポイントは、入力ラスターと同じマップ単位で、X,Y 座標の組として指定されます。 単一ポリゴンのオブジェクトの形式は次のとおりです。
必要に応じて、Polygon クラスを使用してポリゴン パートのリストを定義することで、複数のポリゴンのグループを指定できます。これを実行するには、すべてのパートが連続していて、1 つの境界で囲まれるようにする必要があるので注意してください。この場合のオブジェクトの形式は、次のような連続リストになります。
すべての場合で、ポリゴンを閉じるため、各ポリゴン パートの最後の座標は、最初と同じ座標にする必要があります。 | Point |
extraction_area (オプション) | 入力ポリゴンの内部と外部のどちらのセルを抽出するかを識別します。
| String |
戻り値
| 名前 | 説明 | データ タイプ |
| out_raster | 入力ラスターから抽出したセル値が格納される出力ラスター。 | Raster |
コードのサンプル
次の例では、定義されたポリゴン座標に基づいて、ラスターからセルを抽出しています。
import arcpy
from arcpy import env
from arcpy.sa import *
polyPoints = [arcpy.Point(743050, 4321275), arcpy.Point(743100, 4321200),
arcpy.Point(743500, 4322000),arcpy.Point(742900, 4321800)]
env.workspace = "C:/sapyexamples/data"
extPolygonOut = ExtractByPolygon("soil", polyPoints, "INSIDE")
extPolygonOut.save("c:/sapyexamples/output/extpoly")次の例では、定義されたポリゴン座標に基づいて、ラスターからセルを抽出しています。
# Name: ExtractByPolgyon_Ex_02.py
# Description: Extracts the cells of a raster based on a polygon.
# 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"
polyPoints = [arcpy.Point(743050, 4321275), arcpy.Point(743100, 4321200),
arcpy.Point(743500, 4322000),arcpy.Point(742900, 4321800)]
# Execute ExtractByPolygon
extPolygonOut = ExtractByPolygon(inRaster, polyPoints, "INSIDE")
# Save the output
extPolygonOut.save("c:/sapyexamples/output/extpoly02")環境
ライセンス情報
- Basic: 次のものが必要 Spatial Analyst
- Standard: 次のものが必要 Spatial Analyst
- Advanced: 次のものが必要 Spatial Analyst