Question

我正在使用rasterVis包和levelplot函数在单个地块上绘制七个栅格。

这是我要去的布局： correct layout

我使用以下代码实现了这个目标：

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命令。下面是我用于在单个levelplot中使用单个colorkey绘制所有七个栅格的代码。不幸的是，这没有正确的布局：

incorrect layout

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()

解决方案here已关闭，但不处理奇数个图。如何在levelplot命令中实现所需的布局？

奇数编号使用levelplot和rasterVis绘制

0 个答案: