在对遥感图像进行投影时,常常遇到投影信息不一致的问题,因此需要对数据进行投影转换。
采用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(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)