python - 如何使用 GDAL 更改栅格的分辨率？

Question

我正在寻找更改 GDAL 栅格数据集分辨率的最佳方法。

例如，我有一个像素大小为 (30, -30) 的栅格，我想将像素大小更改为 (5, -5)，将给定像素的所有值内插到输出栅格中。

因此，对于输入栅格的每个像素，我希望输出栅格中有 36 个像素都共享相同的值。

如果我运行gdalwarp -tr 5 -5 inputRaster.tif outputRaster.tif，我会得到我正在寻找的结果，所以我假设我应该能够使用一些 GDAL 函数来复制这个功能。

如果可能的话，我宁愿避免使用对 python 的 Subprocess 类的调用。

score 7 · Accepted Answer

您需要重新投影栅格。例如，从交互式 Python shell：

from osgeo import gdal
help(gdal.ReprojectImage)

测试套件中提供了一个 Python 示例。

score 4 · Accepted Answer

4

使用 gdal.Warp 函数：

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

于 2019-02-28T19:56:13.287 回答

2 回答 2