如何在 R 中使用 levelplot 在绘图中添加子图？

Question

我想在我的主关卡中添加一个子图，类似于this image。我想要两张地图的纬度/经度轴，但我还没有设法找到一种方法来使用 levelplot 做到这一点。另外，我见过this question，但没有答案。 我对 subplot 和新的 par() 感到很困惑，所以我真的需要一些帮助。

library("raster")
library ("rasterVis")
library("viridis")

r <- raster(system.file("external/test.grd", package="raster"))#subplot
r2<-r/2 #main plot
ckey <- list(labels=list(cex=1.5), height=1,space='top')
levelplot(r2,colorkey=ckey,margin=FALSE,xlab="", ylab="",col.regions=viridis,at=seq(100,500, len=20),
          scales=list(x=list(cex=1.5),y=list(cex=1.5))) #from this point, I do not know how to add a subplot

Answer 1

gridExtra 包中的grid.arrange 可以为您做到这一点

library(raster)
library(rasterVis)
library(viridis)
library(gridExtra)

r <- raster(system.file("external/test.grd", package = "raster")) # subplot
r2 <- r / 2 # main plot
ckey <- list(labels = list(cex = 1.5), height = 1, space = "top")

lv1 <- levelplot(r2,
  colorkey = ckey, margin = FALSE, xlab = "", ylab = "", col.regions = viridis, at = seq(100, 500, len = 20),
  scales = list(x = list(cex = 1.5), y = list(cex = 1.5))) # from this point, I do not know how to add a subplot

lv2 <- levelplot(r2,
  colorkey = ckey, margin = FALSE, xlab = "", ylab = "", col.regions = viridis, at = seq(100, 500, len = 20),
  scales = list(x = list(cex = 1.5), y = list(cex = 1.5)))

lay <- rbind(
  c(1, 1, NA, NA, NA),
  c(1, 1, 2, 2, 2),
  c(NA, NA, 2, 2, 2),
  c(NA, NA, 2, 2, 2)
)

grid.arrange(lv1, lv2, layout_matrix = lay,
             top = "Figure 1.", 
             bottom = "bottom\nlabel")

^{由reprex package (v0.3.0) 于 2019 年 11 月 22 日创建}