r - 将栅格中的值检索到具有不同范围和分辨率的另一个栅格中

Question

我有一个气候值的栅格层“w_wgs”。另一方面，我创建了一个空的栅格“网格”。我想将存储在“w_wgs”中的值分配到我的空栅格网格中。我尝试了 merge()、overlay()、over()，但这些都不起作用。

我怎样才能做到这一点 ？假设这两个栅格没有相同的范围和分辨率（显然也不相同数量的单元格）......但两者都具有相同的 CRS（WGS84）。

在这里，我的 2 个对象的描述：

> extent(w_wgs)
class       : Extent 
xmin        : -64.50344 
xmax        : 74.07016 
ymin        : 12.93039 
ymax        : 72.72534 
> res(w_wgs)
[1] 0.01320 0.00895

> res(grid)
[1] 0.08 0.08
> extent(grid)
class       : Extent 
xmin        : 5 
xmax        : 17.96 
ymin        : 40 
ymax        : 50 

Answer 1

我认为你需要使用resample

grid <- resample(w_wgs, grid)

一个应该与您的情况类似的小例子：

x <- matrix(1:100, nr = 10, nc = 10)
a <- raster(x)

x2 <- matrix(NA, nr = 3, nc = 3)
b <- raster(x2)
# Manually changing the extent and resolution of b
b@extent@xmax <- 0.5
b@extent@ymin <- 0.3

b <- resample(a, b)

par(mfrow = c(2, 1))
plot(a)
plot(extent(b), add = T)
plot(b)