如何为 R 中的 SpatialPixelsDataFrame 类对象制作多图？

Question

我有一个类似于one provided in a StackOverflow response 的代码，使用的是 Ege Rubak 的这个答案，但是我正在对这一年进行分析，所以我有一个特定年份每个月的图像。

library(gstat)
library(sp)

lat <-  c(-23.49174, -23.49179, -23.49182, -23.49183, -23.49185, -23.49187)
long <- c(152.0718, 152.0718, 152.0717, 152.0717, 152.0717, 152.0717)
pH <- c(8.222411, 8.19931, 8.140428, 8.100752, 8.068141, 8.048852)
sample <- data.frame(lat, long, pH)

x.range <- range(sample$long)
y.range <- range(sample$lat)

x<-seq(x.range[1], x.range[2], length.out=20)
y<-seq(y.range[1], y.range[2], length.out=20)
grd<-expand.grid(x,y)

coordinates(sample) = ~long+lat
coordinates(grd) <- ~ Var1+Var2
gridded(grd) <- TRUE

proj4string(sample) <- CRS("+proj=longlat +datum=WGS84")
proj4string(grd) <- CRS("+proj=longlat +datum=WGS84")

dat.idw <- idw(formula=pH ~ 1, locations = sample, newdata = grd, idp = 2.0)
#> [inverse distance weighted interpolation]

plot(dat.idw, axes = T)

如下图所示，如何将 12 个图像合并到一个图中并执行类似于 ggplot facet wrap 的操作？

测试this answer中建议的代码：

library(gridExtra)
library(raster)
s <- stack(raster(dat.idw))
s2 <- stack(raster(dat.idw))

p1 <- levelplot(s)
p2 <- levelplot(s2)
grid.arrange(p1, p2, ncol=2)

# Error in UseMethod("levelplot") : 
#  no applicable method for 'levelplot' applied to an object of class 
# "c('RasterStack', 'Raster', 'RasterStackBrick', 'BasicRaster')"

Answer 1

您可以通过准备光栅堆栈对象并将其提供给函数来使用 rasterVis 包中的 levelplot 函数。

library(raster)
library(rasterVis)

r1 <- raster(dat.idw)
r2 <- raster(dat.idw)

# I am only using two rasters as an example
# You can stack 12 rasters for your raster stack 
s <- stack(list(r1, r2))

rasterVis::levelplot(s)

或使用sp 包中的spplot。

sp::spplot(s)

或tmap 包。

库（tmap）

tm_shape(s) +
  tm_raster()