因为前期对栅格数据进行裁剪处理的时候,一个用了别人的批量裁剪工具包,一个自己又构建了一个模型,所以两个模型最后裁剪下来的栅格数据行数先等,列数差“1”。看到网上说可以用列数小的对大的进行掩膜处理,而同学刚好又做过,还分享了代码!所以在此记录一下,方便以后使用。
如有侵权,请私信删除
#ArcPy栅格裁剪:对齐多个栅格图像的范围、统一行数与列数
import arcpy
tif_file_path = r"E:\mask\2002"
result_file_path = r"E:\tuha_25_19902020\1/"
snap_file_name = r"E:\19900101.tif"
arcpy.env.workspace = tif_file_path
arcpy.env.snapRaster = snap_file_name
tif_file_list = arcpy.ListRasters("*", "tif")
for tif_file in tif_file_list:
key_name = tif_file.split(".tif")[0] + ".tif"
arcpy.Clip_management(tif_file,
"#",
result_file_path + key_name,
snap_file_name,
"#",
"#",
"MAINTAIN_EXTENT")
print('finish')