使用Python的GDAL模块打开同一投影但不同尺寸的图像

Question

我不太确定这个地方是否适合问这个问题，不过，

我正在使用Python 2.7中的gdal.Open()模块来打开Landsat的GeoTIFF图像。这些图像使用的是相同的UTM地图投影，但它们的大小不一样。

我该如何调整这些图像的大小，使它们变得一样呢？我想要逐个像素进行比较。

Answer 1

你需要把所有的栅格图像调整成相同的形状和大小。你可以使用gdalwarp这个工具来处理每个栅格文件，比如：

gdalwarp -te xmin ymin xmax ymax -tr xres yres -r resampling_method srcfile dstfile

你可以选择用-ts width height来替代-tr（目标分辨率），这样可以指定行和列的数量。

选项-te（目标范围）和-tr/-ts（目标分辨率/大小）会强制每个源栅格图像都有相同的形状，这样你就可以进行逐像素的比较了。此外，你还需要根据你的分析目的选择一个合适的-r 重采样方法。

Answer 2

gdal_translate这个工具可以很简单地用来改变图像的大小。你只需要使用-outsize这个参数，它需要两个整数值，分别代表图像的宽度和高度，或者你也可以用两个百分比值来缩放图像。点击这里了解更多

gdal_translate -outsize newxsize newysize imageFile outputFile