投影栅格 (Data Management)
用法
-
坐标系可定义栅格数据的投影方式。
-
此工具可保证投影过程中的误差低于半个像素。
-
您可以选择预先存在的空间参考，也可以选择从其他数据集导入或者自行创建一个新的空间参考。
-
您可能需要更改坐标系以便您的数据都使用相同的投影。
-
本工具仅能输出方形像元大小。
-
可将输出结果保存为 BIL、BIP、BMP、BSQ、DAT、GIF、Esri Grid、IMG、JPEG、JPEG 2000、PNG、TIFF 格式或任意地理数据库栅格数据集。
-
将栅格数据集存储到 JPEG 文件、JPEG 2000 文件或地理数据库时，可在“环境设置”中指定压缩类型和压缩质量。
-
通过双线性插值近似方法将栅格数据集投影到新的空间参考，该方法会将像素投影在粗糙的网格格网上并在像素之间应用双线性插值法。
-
NEAREST 选项（用于执行最邻近分配法）是四种插值法当中速度最快的插值方法。其主要用于分类数据（如土地利用分类），因为它不会更改像元值。不推荐对连续数据（如高程表面）使用 NEAREST。
-
BILINEAR 选项使用双线性插值法，可根据四个最近周围像元的加权平均距离确定像元的新值。CUBIC 选项使用三次卷积插值法，可通过拟合穿过周围点的平滑曲线确定新的像元值。这些对于连续数据来说是最适合的选择，但可能会引起平滑处理。请注意，三次卷积插值法可能会导致输出栅格包含输入栅格范围之外的值。不推荐对分类数据使用 BILINEAR 或者 CUBIC，因为不同的像元值可能会进行不必要的引入。
-
栅格数据集的各个像元为方形并且在地图坐标空间中显示为面积相等，但在整个栅格中，像元所代表的地球表面上的形状和面积不可能完全一致。这是因为地图投影无法同时保持形状和面积不变。在整个栅格中，各像元所代表的面积各不相同。因此，输出栅格中的像元大小以及行数和列数也可能发生变化。
-
请务必指定一个输出像元大小，除非您在球面（纬度-经度）坐标和平面坐标系之间进行投影，因为此时无法确定合适的像元大小。
-
默认的输出栅格像元大小由位于输出栅格中心的经过投影的像元大小所决定。这（通常）也是中央经线与真实比例纬线的交汇处以及形变最小的区域。对输入栅格的边界进行投影时，最小和最大范围值将决定输出栅格的大小。将某像元投影回输入坐标系中可确定该像元的值。
-
当输入和输出坐标系的基准面相同时，地理（坐标）变换为可选参数。如果输入和输出基准面不同，则必须指定地理（坐标）变换。
-
配准点可用于指定原点，以便对各输出像元进行定位。所有输出像元与该点之间必须间隔一个像元。该点的坐标不必位于一角，也不必落入栅格数据集中。如果在“环境设置”中设置了“捕捉栅格”，那么将忽略配准点。
-
如果 CLARKE 1866 与投影（例如 NEWZEALAND_GRID）不存在固有关系或者通过 SPHEROID 子命令指定了其他椭球体，则 CLARKE 1866 将作为默认椭球体。
-
如果两者都进行了设置，那么捕捉栅格的参数设置将优先于配准点。
语法
| 参数 | 说明 | 数据类型 |
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 (可选) |
在两个地理坐标系或基准面之间实现变换的方法。 当输入和输出坐标系的基准面相同时，地理（坐标）变换为可选参数。如果输入和输出基准面不同，则必须指定地理（坐标）变换。 有关各个受支持的地理（基准面）变换的详细信息，请参阅位于 <安装位置>\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()