如何使用R中的Terra软件包重新采样重叠的横梁？

Question

我有三个部分重叠的 DEM 栅格，它们具有不同的原点、分辨率（仅略有不同）和范围。我知道由于来源不同等原因，我需要使用 terra 的重新采样函数而不是合并或聚合，但我不确定如何启动一个空栅格以用于具有适当来源、分辨率和范围的重新采样，或者如何处理重叠区域。

是否有推荐的方法来选择现有栅格的原点和分辨率（例如中值与最小/最大值）？如何最好地处理重叠数据？有没有一种有效的方法来生成这样的栅格？

Answer 1

这是显示 的方式。

示例数据

library(terra)
x <- rast(xmin=0, xmax=10, ymin=0, ymax=10, res=1, vals=1)
y <- rast(xmin=9, xmax=19, ymin=0, ymax=10, res=0.9, vals=2)
z <- rast(xmin=0, xmax=10, ymin=8.1, ymax=18.1, res=1, vals=3)

您可能需要将输入栅格之一用作模板。假设您喜欢y。在这种情况下：

a <- list(x, z)
b <- lapply(a, \(i) {
     x <- extend(rast(y), i)
     resample(i, crop(x, i, "out"))
   })

b <- sprc(c(b, y))
m <- merge(b)
 

或创建一个新的模板栅格。在这种情况下，首先找到合并的范围

a <- list(x, y, z)
b <- sapply(a, \(i) ext(i) |> as.vector())
e <- ext(min(b[1,]), max(b[2,]), min(b[3,]), max(b[4,]))

# use the extent to create a raster with the desired spatial resolution 
r <- rast(e, res=1)

，现在如上所述：

gg <- lapply(a, \(i) resample(i, crop(r, i, "out")))
g <- merge(sprc(gg))

或者像这样，

src <- sprc(a)
ss <- impose(src, r)
s <- max(ss, na.rm=TRUE)

我想其中一些可以包裹在Terra方法中。

另请参见Mosaic作为Merge的替代方案。

至于选择最佳分辨率等，这取决于您的需求以及您的数据可能合理地支持的内容。但是，一个重要的一般考虑是，您要避免尽可能多地重新采样 - 因为它会恶化数据质量。