投影栅格 (Data Management)

许可等级:BasicStandardAdvanced

摘要

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

了解有关“投影栅格”工作原理的详细信息

用法

语法

ProjectRaster_management (in_raster, out_raster, out_coor_system, {resampling_type}, {cell_size}, {geographic_transform}, {Registration_Point}, {in_coor_system})
参数说明数据类型
in_raster

输入栅格数据集。

Mosaic Layer; 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.1\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()

环境

相关主题

许可信息

ArcGIS for Desktop Basic:是
ArcGIS for Desktop Standard:是
ArcGIS for Desktop Advanced:是
9/15/2013