基于python实现tif批量转投影

主要运用的包：arcpy

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

 1 import os
 2 import sys
 3 import re
 4 sys.path.append(r"D:\arcgis\Desktop10.2\arcpy")
 5 import arcpy
 6 
 7 arcpy.env.workspace = 'F:\\paper_graduate\\ERA5_4326\\'
 8 
 9 def GetRaster(file_path, out_path):
10     n = 0
11     files = os.listdir(file_path)
12     for file in files:
13         type = file.split(".")
14         if type[1] == "tif":
15             # print file
16             pattern = '.tif$'
17             file_out = re.split(pattern, file)[0]+ "_T.tif"
18             raster = os.path.join(file_path, file)
19             out_raster = os.path.join(out_path, file_out)
20 
21             dataset = "F:\\Analysis_data\\2021tif_MPv4\\MOD09A1.A2021057.MP.tif"
22             mask = arcpy.Describe(dataset).spatialReference
23             cell_size = 25000
24             arcpy.ProjectRaster_management(raster, out_raster, mask, "BILINEAR", cell_size, \
25                                            "#", "#", "#")
26             n = n+1
27     print str(n)+" rasters are processed !!!"
28 
29 
30 if __name__ ==  "__main__":
31     data_path = "F:\\paper_graduate\\ERA5_4326\\"
32     out_path = "F:\\paper_graduate\\ERA5_LM"
33     GetRaster(data_path, out_path)

PS：

ProjectRaster_management各个参数的含义见：https://webhelp.esri.com/arcgisdesktop/9.3/index.cfm?TopicName=project_raster_%28data_management%29
posted @ 2022-05-08 10:39  搞研究不秃头  阅读(491)  评论(0)    收藏  举报