Jak zmienić rozdzielczość rastra przy użyciu GDAL?

Question

Szukam najlepszego sposobu na zmianę rozdzielczości zestawu danych rastrowych GDAL.

Na przykład mam raster, który ma rozmiar piksela (30, -30), i chciałbym zmienić rozmiar piksela na (5, -5), interpolując wszystkie wartości dla danego piksela w wyjściowy raster.

Tak więc dla każdego piksela rastra wejściowego chciałbym mieć 36 pikseli w wyjściowym rastrze, z których wszystkie mają tę samą wartość.

Jeśli uruchomię gdalwarp -tr 5 -5 inputRaster.tif outputRaster.tif, otrzymam dokładnie taki wynik, którego szukam, a więc przypuszczam, że powinienem móc powielić tę funkcjonalność za pomocą funkcji GDAL.

Wolałbym, jeśli to możliwe, unikać używania wywołania podprocesora Pythona.

Answer 1

Musisz ponownie wykonać rastr. Na przykład, z interaktywnym Python powłoki:

from osgeo import gdal 
help(gdal.ReprojectImage) 

Pythona przykład podano w test suite.

Pełniejsza dokumentacja dotyczy funkcji C++ GDALReprojectImage.