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)

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

1 个答案: