1

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

不正确的布局

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

这里的解决方案很接近，但不处理奇数个地块。如何在levelplot命令中实现所需的布局？

4

0 回答 0