基于Python的批量重采样
2019-06-26 本文已影响31人
荔枝猪
需求
现有一些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!")