Desktop Help 10.0 - 投影栅格 (数据管理)

摘要

将栅格数据集从一个投影变换到另一个投影。

用法

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

语法

ProjectRaster_management (in_raster, out_raster, out_coor_system, {resampling_type}, {cell_size}, {geographic_transform}, {Registration_Point}, {in_coor_system})

参数	说明	数据类型
in_raster	输入栅格数据集。	Raster Layer
out_raster	要创建的输出栅格数据集。以文件格式存储栅格数据集时，需要指定文件扩展名，具体如下： .bil - ESRI BIL .bip - ESRI BIP .bmp - BMP .bsq - ESRI BSQ .dat - ENVI DAT .gif - GIF .img - ERDAS IMAGINE 文件 .jpg - JPEG .jp2 - JPEG 2000 .png - PNG .tif - TIFF 无扩展名 - ESRI GRID 以地理数据库形式存储栅格数据集时，不应向栅格数据集的名称添加文件扩展名。将栅格数据集存储到 JPEG 文件、JPEG 2000 文件、TIFF 文件或地理数据库时，可以指定压缩类型和压缩质量。	Raster Dataset
out_coor_system	投影输入栅格的坐标系。默认值将基于“输出坐标系”环境设置进行设定。此参数的有效值为扩展名为“.prj”的文件（ArcGIS 附带的 prj 文件可在“C:\Program Files\ArcGIS\Coordinate Systems”中找到）。现有要素类、要素数据集、栅格目录（基本上包含了与坐标系相关的所有内容）。坐标系的字符串表示。要生成此类较长的字符串，可向模型构建器添加一个坐标系变量，根据需要设置该变量的值，然后将模型导出到 Python 脚本。	Coordinate System
resampling_type （可选）	要使用的重采样算法。默认设置为 NEAREST。 NEAREST —最邻近分配法 BILINEAR —双线性插值法 CUBIC —三次卷积插值法 MAJORITY —多数重采样法 NEAREST 和 MAJORITY 选项用于分类数据，如土地利用分类。NEAREST 选项是默认设置，因为它是最快的插值法，同时也因为它不会更改像元值。请勿对连续数据（如高程表面）使用 NEAREST 或 MAJORITY。 BILINEAR 选项和 CUBIC 选项最适用于连续数据。不推荐对分类数据使用 BILINEAR 或者 CUBIC，因为像元值可能被更改。	String
cell_size （可选）	新栅格数据集的像元大小。默认像元大小为所选栅格数据集的像元大小。	Analysis Cell Size
geographic_transform （可选）	在两个地理坐标系或基准面之间实现变换的方法。当输入和输出坐标系的基准面相同时，地理（坐标）变换为可选参数。如果输入和输出基准面不同，则必须指定地理（坐标）变换。有关各个受支持的地理（基准面）变换的详细信息，请参阅位于 <安装位置>\ArcGIS\Desktop10.0\Documentation 中的 geographic_transformations.pdf。	String
Registration_Point （可选）	对齐像素时使用的 x 和 y 坐标（位于输出空间中）。注册点的工作原理与捕捉栅格的概念类似。通过注册点可指定原点以便对输出像元进行定位，而不是将输出捕捉到现有栅格像元对齐方式。所有输出像元与该点之间必须间隔一个像元大小以上。该点的坐标不必位于一角，也不必落入栅格数据集中。捕捉栅格环境设置将优先于注册点参数。因此，如果您要设置注册点，请确保尚未设置捕捉栅格。	Point
in_coor_system （可选）	输入栅格数据集的坐标系。	Coordinate System

代码示例

ProjectRaster 示例 1（Python 窗口）

这是 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 示例 2（独立脚本）

这是 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()

环境

压缩, 当前工作空间, 输出配置关键字, 输出坐标系, 范围, 金字塔, 栅格统计, 临时工作空间, 捕捉栅格, 切片大小

许可信息

ArcView: 是

ArcEditor: 是

ArcInfo: 是

7/10/2012

投影栅格 (数据管理)

摘要

用法

语法

代码示例

环境

相关主题

许可信息