Организация растрового графика в R с использованием rasterVis и gridExtra

Мне трудно строить карты вместе в R, используя rasterVis и gridExtra. Вот как выглядит мой код (извините, конечно, не элегантно):

require(rasterVis)
require(gridExtra)

p1 <- levelplot(avit_bac12, margin=FALSE, xlab=NULL, ylab=NULL,  par.settings = RdBuTheme, main=NULL, scales=list(draw=FALSE))
p2 <- levelplot(avit_bac17, margin=FALSE, xlab=NULL, ylab=NULL,  par.settings = RdBuTheme, main=NULL, scales=list(draw=FALSE))
p3 <- levelplot(avit_eng, margin=FALSE, xlab=NULL, ylab=NULL,  par.settings = RdBuTheme, main=NULL, scales=list(draw=FALSE))
p4 <- levelplot(avit_sf, margin=FALSE, xlab=NULL, ylab=NULL,  par.settings = RdBuTheme, main=NULL, scales=list(draw=FALSE))
p5 <- levelplot(bac12_bac17, margin=FALSE, xlab=NULL, ylab=NULL,  par.settings = RdBuTheme, main=NULL, scales=list(draw=FALSE))
p6 <- levelplot(bac12_eng, margin=FALSE, xlab=NULL, ylab=NULL,  par.settings = RdBuTheme, main=NULL, scales=list(draw=FALSE))
p7 <- levelplot(bac12_sf, margin=FALSE, xlab=NULL, ylab=NULL,  par.settings = RdBuTheme, main=NULL, scales=list(draw=FALSE))
p8 <- levelplot(bac17_eng, margin=FALSE, xlab=NULL, ylab=NULL,  par.settings = RdBuTheme, main=NULL, scales=list(draw=FALSE))
p9 <- levelplot(bac17_sf, margin=FALSE, xlab=NULL, ylab=NULL,  par.settings = RdBuTheme, main=NULL, scales=list(draw=FALSE))
p10 <- levelplot(eng_sf, margin=FALSE, xlab=NULL, ylab=NULL,  par.settings = RdBuTheme, main=NULL, scales=list(draw=FALSE))

p.avit <- levelplot(avit, margin=FALSE, xlab=NULL, ylab=NULL, par.settings = magmaTheme, main="Avitabile et al. (2016)", scales=list(draw=FALSE))
p.bac12 <- levelplot(bac12, margin=FALSE, xlab=NULL, ylab=NULL, par.settings = magmaTheme, main="Baccini et al. (2012)", scales=list(draw=FALSE))
p.bac17 <- levelplot(bac17, margin=FALSE, xlab=NULL, ylab=NULL, par.settings = magmaTheme, main="Baccini et al. (2017)", scales=list(draw=FALSE))
p.eng <- levelplot(eng, margin=FALSE, xlab=NULL, ylab=NULL, par.settings = magmaTheme, main="Englund et al. (2017)", scales=list(draw=FALSE))
p.sf <- levelplot(sf, margin=FALSE, xlab=NULL, ylab=NULL, par.settings = magmaTheme, main="Soares-Filho et al. (2017)", scales=list(draw=FALSE))

p.list <- list(p.avit, p.bac12, p.bac17, p.eng, p.sf,
               p1, p2, p3, p4, p5, p6, p7, p8, p9, p10)

lay <- rbind(c(NA,  2,  3,  4,  5),
             c(1,   6,  7,  8,  9),
             c(2,   NA, 10, 11, 12),
             c(3,   NA, NA, 13, 14),
             c(4,   NA, NA, NA, 15))

grid.arrange(grobs=p.list, layout_matrix=lay,
         vp = grid::viewport(width=0.4,height=0.3))

Я надеялся, что он вернет что-то вроде этого:

... но вместо этого размер карт искажается, и я получаю это:

Любые идеи о том, как решить эту проблему?

TWest 22.08.2019 источник

comment

Кстати, levelplot() просто считывает растры - TWest 22.08.2019

Ответы (1)

arrow_upward
1
arrow_downward

Проблема в том, что некоторые числа в layout_matrix повторяются. Повторяющиеся числа означают, что график должен охватывать ячейки сетки с одинаковыми номерами. Поскольку каждый из ваших графиков должен занимать одну ячейку в сетке, никакие числа в матрице макета не должны повторяться. Если график необходимо разместить в двух местах, передайте его grid.arrange дважды, но убедитесь, что layout_matrix содержит уникальные номера (и NA, если необходимо). Проще с примерами:

library(ggplot2)
library(gridExtra)
p <- ggplot(iris, aes(x=Sepal.Width, y=Sepal.Length, color=Species)) + geom_point(show.legend = FALSE)

lay1 <- rbind(c(NA, 1,2), 
              c(2,3, NA))
grid.arrange(p,p,p, layout_matrix=lay1) 
# adding a 4th plot in the above grid.arrange call results with error

lay2 <- rbind(c(NA, 1,2), 
              c(3,4, NA))
grid.arrange(p,p,p,p, layout_matrix=lay2)

lay3 <- rbind(c(1,1,NA,2),
              c(1,1,3, NA))
grid.arrange(p,p,p, layout_matrix=lay3)

Три примера должны давать следующие графики:

слой1:

слой2:

слой3:

teofil 22.08.2019

Организация растрового графика в R с использованием rasterVis и gridExtra

Ответы (1)

Похожие вопросы