layout - levelplot と rasterVis を使用した奇数プロット

Question

rasterVisパッケージと関数を使用してlevelplot、1 つのプロットに 7 つのラスターをプロットしています。

これが私が目指しているレイアウトです: 正しいレイアウト

次のコードを使用してこれを達成しました。

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画像全体の 1 つのコモンが欠けています。

コマンドを使用することをお勧めしますlevelplot。以下は、7 つのラスターすべてlevelplotを 1 つの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コマンド内で目的のレイアウトを実現するにはどうすればよいですか?