Python+ArcPy 对栅格数据进行批量重投影

在对遥感图像进行投影时，常常遇到投影信息不一致的问题，因此需要对数据进行投影转换。

采用ArcGIS中Data Management --> Projection and Transformation --> Raster --> Project Raster可以对栅格文件进行逐一投影，但是当文件数量较大时，该操作并不是很方便。因此这里我们采用ArcPy将栅格数据批量转换到与shp文件相同的投影。

import os
import arcpy


arcpy.env.workspace = 'E:/data/'


def GetRaster(file_path, out_path):
    n=0
    files = os.listdir(file_path)
    for file in files:
        type = file.split(".")
        if type[-1] == "tif":
            print file
            raster = os.path.join(file_path,file)
            out_raster = os.path.join(out_path,file)
            projecttype = "Albers_Province.prj"      ## 目标投影信息的文件路径,可以是shapefile对应的prj文件
            arcpy.ProjectRaster_management(raster,out_raster, \
                                           projecttype, "NEAREST","#" , "", "","#")
            n=n+1
    print str(n)+" rasters are processed !!!"



if __name__ ==  "__main__":
    data_path = "./data/"    ## 栅格数据路径
    out_path = "./out/"      ## 输出文件路径
    GetRaster(data_path, out_path)