Layout 使用levelplot和rasterVis的奇数图

我正在使用

rasterVis

包和

levelplot

函数在一个绘图上绘制七个光栅

以下是我想要的布局：

我使用以下代码实现了这一点：

png("E:/all_files/production/plots/for_final/hist_vis_split2.png", 
width=6, height=9, units="in", res=72)
layout(matrix(c(1,1,2,3,4,5,6,7), 4, 2, byrow = T))
print(maurer+layer(sp.lines(sr)), split=c(2,1,3,4), more=TRUE)
print(nar.ccsm+layer(sp.lines(sr)), split=c(1,2,3,4), more=TRUE)
print(nar.gfdl+layer(sp.lines(sr)), split=c(3,2,3,4), more=TRUE)
print(bcsd.ccsm+layer(sp.lines(sr)), split=c(1,3,3,4), more=TRUE)
print(bcsd.gfdl+layer(sp.lines(sr)), split=c(3,3,3,4), more=TRUE)
print(bcca.ccsm+layer(sp.lines(sr)), split=c(1,4,3,4), more=TRUE)
print(bcca.gfdl+layer(sp.lines(sr)), split=c(3,4,3,4))
dev.off()

除了笨重和缺乏控制之外，这张图片还缺少一个主标题和一个通用的

colorkey

我更喜欢使用

levelplot

命令。下面是我用一个

colorkey

在一个

levelplot

中绘制所有七个光栅的代码。遗憾的是，此布局不正确：

---

crop.stack可能值得签出
gridExtra:：grid.arrange
：可能值得签出
gridExtra:：grid.arrange
：
crop.stack <- stack(maurer, bcsd.ccsm.crop, bcsd.gfdl.crop, bcca.ccsm.crop, 
                bcca.gfdl.crop, nar.ccsm.crop, nar.gfdl.crop)
plot.names <- c("Maurer", "BCSD CCSM", "BCSD GFDL", "BCCA CCSM", 
            "BCCA GFDL", "NARCCAP CCSM", "NARCCAP GFDL")

png("E:/all_files/production/plots/for_final/hist_vis.png", 
    width=6, height=9, units="in", res=72)
hist <- levelplot(crop.stack, main="Historical Average Production Days", 
          col.regions=cols, 
          names.attr=plot.names, 
          scales=list(draw=F), 
          layout=c(2,4))
hist + layer(sp.lines(sr))
dev.off()