用模型构建器批量处理比较慢,因此可以利用Arcpy进行处理
# -*- coding: utf-8 -*-
import arcpy
# 多个shp文件所在文件夹
z_rasters = "H:/wd30/ww/sci/global30/z2001-2020/"
# 工作环境
arcpy.env.workspace = z_rasters
# 提取文件夹中的tif栅格文件
rasters = arcpy.ListFiles("*.tif")
# 输出raster所在路径
output_raster = "H:/wd30/ww/sci/global30/project/"
# 循环转换投影
for raster in rasters:
# 获取名称
catchment_name = raster[:-4]
out_raster_name = output_raster + catchment_name + ".tif"
# 投影坐标系WKID代码
out_coor_system = 32639
# 这个“30 30”是重采样的大小
arcpy.ProjectRaster_management (raster, out_raster_name, out_coor_system, "NEAREST", "30 30")
print "finished"
当数据需要进行跨地理坐标系进行投影转换时,例如我想从GCS_WGS_1894投影到CGCS2000_3_Degree_GK_Zone_39的坐标系,那么要首先进行一个地理坐标系的转换,这里可以直接进行定义坐标系,定义地理坐标系以后再进行投影转换。
# -*- coding: utf-8 -*-
import arcpy
# 多个shp文件所在文件夹
z_rasters = "H:/wd30/ww/sci/global30/z2001-2020/"
# 工作环境
arcpy.env.workspace = z_rasters
# 提取文件夹中的tif栅格文件
rasters = arcpy.ListFiles("*.tif")
# 输出raster所在路径
output_raster = "H:/wd30/ww/sci/global30/project/"
# 循环转换投影
for raster in rasters:
# 获取名称
catchment_name = raster[:-4]
out_raster_name = output_raster + catchment_name + ".tif"
# 这里的4490和4527是坐标系的WKID码
arcpy.DefineProjection_management(raster, 4490)
arcpy.ProjectRaster_management (raster, out_raster_name, 4527, "NEAREST", "30 30")
print "finished"