ArcGIS中利用Arcpy批量对土地利用类型数据投影转换

用模型构建器批量处理比较慢，因此可以利用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"