基于python实现tif批量转投影

主要运用的包：arcpy

配置方法：https://www.cnblogs.com/renyan1128/p/16225343.html

import os
import sys
import re
sys.path.append(r"D:\arcgis\Desktop10.2\arcpy")
import arcpy

arcpy.env.workspace = 'F:\\paper_graduate\\ERA5_4326\\'

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
            pattern = '.tif$'
            file_out = re.split(pattern, file)[0]+ "_T.tif"
            raster = os.path.join(file_path, file)
            out_raster = os.path.join(out_path, file_out)

            dataset = "F:\\Analysis_data\\2021tif_MPv4\\MOD09A1.A2021057.MP.tif"
            mask = arcpy.Describe(dataset).spatialReference
            cell_size = 25000
            arcpy.ProjectRaster_management(raster, out_raster, mask, "BILINEAR", cell_size, \
                                           "#", "#", "#")
            n = n+1
    print str(n)+" rasters are processed !!!"


if __name__ ==  "__main__":
    data_path = "F:\\paper_graduate\\ERA5_4326\\"
    out_path = "F:\\paper_graduate\\ERA5_LM"
    GetRaster(data_path, out_path)

PS：

ProjectRaster_management各个参数的含义见：https://webhelp.esri.com/arcgisdesktop/9.3/index.cfm?TopicName=project_raster_%28data_management%29