基于Python的批量重采样

需求

现有一些tif格式的MODIS的EVI数据，空间分辨率250m，希望重采样为1km。

简介

Resample_management (in_raster, out_raster, {cell_size}, {resampling_type})

• in_raster:想要更改空间分辨率的栅格数据集。
• out_raster:要创建的数据集的名称、位置和格式。
• cell_size:指定像元大小
  • 使用单个数字指定方形像元大小
  • 使用两个数字（以空格分隔）指定 X 和 Y 像元大小
  • 使用栅格数据集（从其导入方形像元大小）的路径
• resampling_type：重采样方法
  • NEAREST 最邻近法
  • BILINEAR 双线性插值
  • CUBIC三次卷积
  • MAJORITY众数重采样

代码

import arcpy
arcpy.env.workspace = "I:\\MODIS\\EVI\\All_AlbersProjection"
rasterlist = arcpy.ListRasters("*", "tif")   
for raster in rasterlist:
    print str(raster)
    out = "I:\\MODIS\\EVI\\2006_2017\\"+"evi"+raster[26:34]+".tif"  
    arcpy.Resample_management(raster, out, "1000", "NEAREST")
print("OK!")