#legend
Вопрос:
Я хочу сделать многопанельный график с несколькими растрами с категориальными значениями. Не могли бы вы помочь мне создать многопанельный сюжет с общей легендой? Мои растры категоричны со значениями в диапазоне от 0 до 11. Например, в одном растре у меня есть значения 0,1,2 и 10, а в другом-0, 5 и 11. Вот мой код. Я не смог сгенерировать воспроизводимый пример 🙁
library(raster)
library(rasterVis)
#read the files
r_cur <- raster(file.path(mydircur,"Pab_cur_best.tif"))
r_45 <- raster(file.path(mydir45,"Pab_45_best.tif"))
r_85 <- raster(file.path(mydir85,"Pab_85_best.tif"))
#Convert the rasters to factors and add attributes
r1 <- as.factor(r_cur)
r2<- as.factor(r_45)
r3<- as.factor(r_85)
r4<- as.factor(clus)
## Add a column "Seed" to the Raster Attribute Table
rat1 <- levels(r1)[[1]]
rat1[["Seed"]] <- c("Local", "CL1", "CL2", "CL3", "CL4", "CL5", "CL6", "CL7", "CL8", "CL9", "CL11")
levels(r1) <- rat1
rat2 <- levels(r2)[[1]]
rat2[["Seed"]] <- c("Local", "CL1", "CL2", "CL3", "CL4", "CL5", "CL6", "CL7", "CL9", "CL11")
levels(r2) <- rat2
rat3 <- levels(r3)[[1]]
rat3[["Seed"]] <- c("Local", "CL1", "CL2", "CL3", "CL4", "CL5", "CL6", "CL7", "CL8", "CL9", "CL11")
levels(r3) <- rat3
#Specify colours for each plot seperately
# then I make the multiplot
p1 <- levelplot(r1, col.regions=cl_cur,margin=FALSE, xlab=NULL, ylab=NULL,scales=list(draw=FALSE)) layer(sp.polygons(EU))
p2 <- levelplot(r2, col.regions=cl_45,margin=FALSE, xlab=NULL, ylab=NULL,scales=list(draw=FALSE)) layer(sp.polygons(EU))
p3 <- levelplot(r3, col.regions=cl_85,margin=FALSE, xlab=NULL, ylab=NULL,scales=list(draw=FALSE)) layer(sp.polygons(EU))
# then I make the multiplot
grid.arrange(p1, p2,p3,p4, ncol=4)
g <- arrangeGrob(p1, p2,p3,p4, ncol=4)
Большое спасибо