使用python gdal模块以相同的投影但不同的大小打开图像

Question

我不确定这是否是正确的网站，但

我在python 2.7上使用gdal.Open（）模块打开Landsat GeoTIFF图像。它们具有相同的UTM地图投影，但图像尺寸不同。

如何将图像尺寸固定为相同？我想直接进行逐像素比较。

Answer 1

gdal_translate可以简单地用于使用-outsize参数更改图像的大小，该参数分别采用两个整数值作为xsize和ysize，或两个百分比值来缩放图像。 see here

gdal_translate -outsize newxsize newysize imageFile outputFile

Answer 2

您需要将所有栅格重新投影到相同的栅格形状/大小。您可以使用gdalwarp对每个光栅文件执行此操作，例如：

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

您可以选择将-tr（目标分辨率）替换为-ts width height以指定行数/列数。

-te（目标范围）和-tr / -ts（目标分辨率/尺寸）选项会强制每个源栅格具有相同的栅格形状叠加，因此您可以正确进行逐像素比较。您还需要选择合适的-r resampling method，具体取决于您的分析尝试实现的目标。