Использование модуля Python GDAL для открытия изображений в одной проекции, но разных размеров

Question

Я не уверен, что это правильный сайт для этого, но,

Я использую модуль gdal.Open () в python 2.7, чтобы открывать изображения Landsat GeoTIFF. Они имеют одинаковую проекцию карты UTM, но разные размеры изображений.

Как мне исправить размеры изображения, чтобы они были одинаковыми? Я хотел бы сделать прямое попиксельное сравнение.

Answer 1

Вам необходимо перепроецировать все растры в одну и ту же форму / размер. Вы можете сделать это с gdalwarp для каждого растрового файла, например ::

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

При желании можно заменить -tr (целевое разрешение) на -ts width height, чтобы указать количество строк / столбцов.

Параметры -te (целевые экстенты) и -tr / -ts (целевое разрешение / размер) заставляют каждый исходный растр иметь одинаковое наложение формы растра, поэтому вы можете корректно выполнять поэлементное пиксельное сравнение. Вам также необходимо выбрать подходящий -r метод повторной выборки , в зависимости от того, чего пытается достичь ваш анализ.

Answer 2

gdal_translate можно использовать просто для изменения размера изображения с помощью параметра -outsize, который принимает два целочисленных значения как xsize и ysize соответственно или два процентных значения для масштабирования изображения. см. Здесь

gdal_translate -outsize newxsize newysize imageFile outputFile