Блог

Я ищу лучший способ изменить разрешение набора растровых данных GDAL.

Например, у меня есть растр с размером пикселя (30, -30), и я хотел бы изменить размер пикселя на (5, -5), интерполируя все значения для данного пикселя в выходной растр.

Итак, для каждого пикселя входного растра я хотел бы иметь 36 пикселей в выходном растре, которые имеют одинаковое значение.

Если я запущу gdalwarp -tr 5 -5 inputRaster.tif outputRaster.tif , я получу именно тот результат, который я ищу, и поэтому я бы предположил, что я смогу воспроизвести эту функциональность с помощью некоторой функции GDAL.

Я бы предпочел избегать использования вызова класса подпроцесса python, если это возможно.

Вам нужно перепроектировать растр. Например, из интерактивной оболочки Python:

 from osgeo import gdal
help(gdal.ReprojectImage)
  

Пример Python приведен в наборе тестов.

Более полная документация предоставляется для функции C GDALReprojectImage .

Используйте gdal.Функция деформации:

 gdal.Warp('outputRaster.tif', 'inputRaster.tif', xRes=5, yRes=5)