ラスタの投影変換(Project Raster) (データ管理)
使用法
-
ラスタ データの投影方法は、座標系によって定義されます。
-
このツールは、誤差を 1/2 ピクセル未満に抑えることができます。
-
先に存在する空間参照を選択して別のデータセットからインポートすることも、新規に作成することもできます。
-
必要に応じて、すべてのデータが同じ投影に含まれるように、座標系を変更できます。
-
このツールが出力できるのは、矩形のセル サイズのみです。
-
出力は、BIL、BIP、BMP、BSQ、DAT、GIF、Esri Grid、IMG、JPEG、JPEG 2000、PNG、TIFF、または任意のジオデータベース ラスタ データセットに保存することができます。
-
ラスタ データセットを JPEG ファイル、JPEG 2000 ファイル、またはジオデータベースに格納するときに、環境設定で [圧縮] タイプと [圧縮品質] を指定できます。
-
共一次内挿近似法を使用して、ラスタ データセットを新しい空間参照に投影変換します。これにより、ピクセルを粗いメッシュ グリッドに投影変換し、ピクセル間に共一次内挿法を使用します。
-
NEAREST オプションは、最近隣内挿法を実行する、4 つの内挿法の中で最も高速の選択肢です。この方法は、セル値を変更しないため、主に土地利用分類などのカテゴリ データに使用されます。NEAREST を標高サーフェスなどの連続データに使用することは推奨されません。
-
BILINEAR オプションは、共一次内挿法を使用し、周囲の最も近い 4 つのセルの加重距離平均に基づいて、セルの新しい値を決定します。CUBIC オプションは、三次たたみ込み内挿法を使用し、周囲のポイントの間を通る滑らかなカーブを適合して、セルの新しい値を決定します。これらは連続データに対する最も適切な選択ですが、スムージングが発生する場合があります。三次たたみ込み内挿法は、入力ラスタの範囲外の値を含んだ出力ラスタを生成することがあります。望ましくない異なるセル値が生じる可能性があるため、BILINEAR または CUBIC をカテゴリ データに使用することは推奨されません。
-
ラスタ データセットのセルは矩形で、マップ座標空間内では等面積ですが、地表面上にセルで表される形状と面積はラスタ全体で一定になることはありません。なぜなら、マップ投影では形状と面積の両方を同時には保持できないためです。セルで表される面積は、ラスタ全体では異なります。そのため、出力ラスタ内のセル サイズ、行数および列数が変動することがあります。
-
球面(緯度 - 経度)座標と平面座標系間に投影変換を行う場合と、適当なセル サイズが不明なとき以外は、常に出力セル サイズを指定してください。
-
出力ラスタのデフォルトのセル サイズは、出力ラスタの中心位置の投影セル サイズに基づいて計算されます。これは通常、厳密な縮尺の中央子午線と緯度の交差部であり、最小歪み領域でもあります。入力ラスタの境界を投影変換すると、最小範囲と最大範囲に基づいて出力ラスタのサイズが提示されます。セルの値を計算するために、各セルを元の入力座標系に投影変換します。
-
入力座標系と出力座標系に同じ測地基準系が使用されている場合、地理座標系変換はオプション パラメータです。入力測地基準系と出力測地基準系が異なる場合は、地理座標系変換の指定が必須になります。
-
登録ポイントを使用して、出力セルを固定する原点に指定することができます。すべての出力セルは、このポイントを基点とするセル サイズの間隔になります。このポイントはコーナー座標である必要はなく、ラスタ データセット内に収まる必要もありません。[環境設定] で [スナップ対象ラスタ] が設定してある場合、登録ポイントは無視されます。
-
デフォルト回転楕円面は CLARKE 1866 です。NEWZEALAND_GRID など、投影変換に固有のものでない場合、SPHEROID サブコマンドを使用して別の値を指定してください。
-
両方が設定されている場合、優先的に使用されるのは、登録ポイントでなくスナップ ラスタ設定のほうです。
構文
パラメータ | 説明 | データ タイプ |
in_raster |
入力ラスタ データセット。 | Mosaic Layer; Raster Layer |
out_raster |
作成される出力ラスタ データセット。 ラスタ データセットをファイル形式で格納する場合は、ファイル拡張子を指定する必要があります。
ジオデータベースにラスタ データセットを格納する場合、ラスタ データセットの名前にファイル拡張子は付けません。 ラスタ データセットを JPEG ファイル、JPEG 2000 ファイル、TIFF ファイル、またはジオデータベースに格納するときに、圧縮タイプと圧縮品質を指定できます。 | Raster Dataset |
out_coor_system |
入力ラスタの投影先となる座標系。デフォルト値は、[出力データの座標系] 環境設定に基づいて設定されます。 このパラメータの有効値は、次のとおりです。
| Coordinate System |
resampling_type (オプション) |
使用するリサンプリング アルゴリズム。デフォルトは NEAREST です。
NEAREST オプションと MAJORITY オプションは、土地利用分類などのカテゴリ データに使用されます。デフォルトは NEAREST オプションです。これは、このオプションが最速であり、セル値を変更しないからです。NEAREST または MAJORITY を標高サーフェスなどの連続データに使用しないでください。 BILINEAR オプションと CUBIC オプションは、連続データに最も適しています。セル値が変更される可能性があるため、BILINEAR または CUBIC をカテゴリ データに使用することは推奨されません。 | String |
cell_size (オプション) |
新規ラスタ データセットのセル サイズ。 デフォルトのセル サイズは、選択したラスタ データセットのセル サイズです。 | Analysis Cell Size |
geographic_transform (オプション) |
2 つの地理座標系または測地基準系の間で使用される変換方式。 入力座標系と出力座標系に同じ測地基準系が使用されている場合、地理座標系変換はオプションです。入力測地基準系と出力測地基準系が異なる場合は、地理座標系変換の指定が必須になります。 サポートされている各地理座標系(測地基準系)変換については、<インストール ディレクトリ>\ArcGIS\Desktop10.1\Documentation にある geographic_transformations.pdf をご参照ください。 | String |
Registration_Point (オプション) |
(出力領域内の)ピクセル配置に使用される X 座標と Y 座標。 登録ポイントの動作は、スナップ ラスタの概念と同様です。既存のラスタセル配置に出力をスナップするかわりに、登録ポイントを出力セルを固定する原点を指定することができます。すべての出力セルは、このポイントを基点とするセル サイズの間隔になります。このポイントはコーナー座標である必要はなく、ラスタ データセット内に収まる必要もありません。 [スナップ対象ラスタ] 環境設定は、[登録ポイント] パラメータよりも優先的に使用されます。この理由から、登録ポイントを設定する必要がある場合は [スナップ対象ラスタ] を設定しないようにしてください。 | Point |
in_coor_system (オプション) |
入力ラスタ データセットの座標系。 | Coordinate System |
コードのサンプル
以下は、ProjectRaster(ラスタの投影変換)ツールを実行する Python サンプルです。
import arcpy
from arcpy import env
arcpy.ProjectRaster_management("c:/data/image.tif", "c:/output/reproject.tif",\
"World_Mercator.prj", "BILINEAR", "5",\
"NAD_1983_To_WGS_1984_5", "#", "#")
以下は、ProjectRaster(ラスタの投影変換)ツールを実行する Python スクリプト サンプルです。
##====================================
##Project Raster
##Usage: ProjectRaster_management in_raster out_raster out_coor_system {NEAREST | BILINEAR
## | CUBIC | MAJORITY} {cell_size} {geographic_transform;
## geographic_transform...} {Registration_Point} {in_coor_system}
try:
import arcpy
arcpy.env.workspace = r"C:/Workspace"
##Reproject a TIFF image with Datumn transfer
arcpy.ProjectRaster_management("image.tif", "reproject.tif", "World_Mercator.prj",\
"BILINEAR", "5", "NAD_1983_To_WGS_1984_5", "#", "#")
##Reproject a TIFF image that does not have a spatial reference
##Set snapping point to the top left of the original image
snapping_pnt = "1942602 304176"
arcpy.ProjectRaster_management("nosr.tif", "project.tif", "World_Mercator.prj", "BILINEAR",\
"5", "NAD_1983_To_WGS_1984_6", snapping_pnt,\
"NAD_1983_StatePlane_Washington_North.prj")
except:
print "Project Raster example failed."
print arcpy.GetMessages()