r - 栅格的重采样非常慢

Question

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

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

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

谢谢：-）

score 1 · Accepted Answer

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

相同的语法和参数：

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

score 0 · Accepted Answer

您可以尝试使用gdalwarppackage中的功能gdalUtils。

这里有一个关于如何gdalwarp用于重采样的示例：

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

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

2 回答 2