在R中合并多个栅格

我一直在尝试找到一种省时的方法来在 R 中合并多个光栅图像。这些是来自乞力马扎罗南部地区的相邻 ASTER 场景，我的目标是将它们放在一起以获得一张大图像。

这是我到目前为止得到的（对象“ast14dmo”代表 RasterLayer 对象列表）：

# Loop through single ASTER scenes
for (i in seq(ast14dmo.sd)) {
  if (i == 1) {
    # Merge current with subsequent scene
    ast14dmo.sd.mrg <- merge(ast14dmo.sd[[i]], ast14dmo.sd[[i+1]], tolerance = 1)
  } else if (i > 1 && i < length(ast14dmo.sd)) {
    tmp.mrg <- merge(ast14dmo.sd[[i]], ast14dmo.sd[[i+1]], tolerance = 1)
    ast14dmo.sd.mrg <- merge(ast14dmo.sd.mrg, tmp.mrg, tolerance = 1)
  } else {
    # Save merged image
    writeRaster(ast14dmo.sd.mrg, paste(path.mrg, "/AST14DMO_sd_", z, "m_mrg", sep = ""), format = "GTiff", overwrite = TRUE)
  }
}

正如您肯定猜到的那样，该代码有效。然而，考虑到每个栅格对象大约 70 mb 大，合并需要相当长的时间。我还尝试了Reduce和do.call，但是失败了，因为我无法传递绕过光栅文件不同来源的参数“tolerance”。

有人知道如何加快速度吗？

使用“terra”，您首先创建一个 SpatRasterCollection

library(terra)
r1 <- rast(xmax=-150, ymin=60, ncols=30, nrows=30)
values(r1) <- 1:ncell(r1)
r2 <- rast(xmin=-100, xmax=-50, ymax=50, ymin=30)
res(r2) <- c(xres(r1), yres(r1))
values(r2) <- 1:ncell(r2)

s <- sprc(r1, r2)
m <- merge(s)

使用已弃用的“raster”包。

您可以使用do.call

ast14dmo.sd$tolerance <- 1
ast14dmo.sd$filename <- paste(path.mrg, "/AST14DMO_sd_", z, "m_mrg.tif", sep = "")
ast14dmo.sd$overwrite <- TRUE
mm <- do.call(merge, ast14dmo.sd)

这里有一些数据，来自示例raster::merge

r1 <- raster(xmx=-150, ymn=60, ncols=30, nrows=30)
r1[] <- 1:ncell(r1)
r2 <- raster(xmn=-100, xmx=-50, ymx=50, ymn=30)
res(r2) <- c(xres(r1), yres(r1))
r2[] <- 1:ncell(r2)

x <- list(r1, r2)
names(x) <- c("x", "y")
x$filename <- 'test.tif'
x$overwrite <- TRUE
m <- do.call(merge, x)