栅格重采样非常慢

Question

) 我想对 Sentinel-2 图像重新采样，以便两个波段的分辨率适合。我为此使用 raster::resample(raster1, raster2, method="ngb")，但执行时间非常长。它是运行至少 15 分钟，到目前为止我没有获得任何输出。我不确定这是否会在某个时候完成。

我使用的是 Sentinel-2 的 08 波段（10 m 分辨率）和 12 波段（20 m 分辨率）。

使用 QGIS 对同一个栅格进行重新采样需要一分钟左右的时间。我想问一下这是否是此功能的常见问题，是否有任何解决方法或其他 function/package 我可以用来解决这个问题。

谢谢 :-)

Answer 1

您可以尝试包 gdalUtils.

中的函数 gdalwarp

这里有一个关于如何使用 gdalwarp 进行重采样的例子：

这里有一个基于gdalwarp的重采样函数的例子：

如果您安装了 GRASS GIS，您还可以查看具有 fasterProjectRaster 功能的 fasterRaster 软件包。

Answer 2

也许 terra 包更快，虽然我没有为此任务做过速度测试。

相同的语法和参数：

resampled_raster <- terra::resample(raster1, raster2, method = "ngb")

Resampling of raster very slow