GDAL根据行列数以及空间范围重采样栅格影像

近期有对栅格影像重采样的需要，即将源栅格依据目标栅格的分辨率和空间范围进行空间上的重采样，使像元完全匹配。网上教程多建议使用 ArcGIS 的重采样工具或者裁剪或者栅格计算器工具，使用过程中发现有点问题，特此记录。

1. 匹配空间分辨率

使用栅格计算器, Con(IsNull(原始栅格), 0, 原始栅格)，然后在环境中选择参考栅格的空间范围和像元大小，此时输出的目标栅格与参考栅格空间分辨率基本一致，但是行列数不一致，因此需要下一步统一行列数。经过裁剪还有重采样等工具尝试无果后，使用GDAL解决问题。

2. gdalwarp 指定栅格影像行列数以及空间范围进行重采样

使用 GDAL 附带的 gdalwarp 工具，

其命令行调用格式如下：

gdalwarp -ts <cols> <rows> -te xmin ymin xmax ymax -r bilinear input.tif output.tif

gdalwarp -ts 588 3613 -te 668061.548884 4741010.683808 669168.191065 4747169.293762 -r bilinear simu_gas_plumn_CH4_ppmm_HyTES_embed.tif simu_gas_plumn_CH4_ppmm_HyTES_embed_resample.tif

然后输出:

Creating output file that is 588P x 3613L.
Processing simu_gas_plumn_CH4_ppmm_HyTES_embed.tif [1/1] : 0Using internal nodata values (e.g. -3.40282e+38) for image simu_gas_plumn_CH4_ppmm_HyTES_embed.tif.
Copying nodata values from source simu_gas_plumn_CH4_ppmm_HyTES_embed.tif to destination simu_gas_plumn_CH4_ppmm_HyTES_embed_resample.tif.
...10...20...30...40...50...60...70...80...90...100 - done.

这样输出影像的行数就是 588，列数就是 3613。ok！