How to plot 10 maps as one figure in R?

Question

I am using this code to plot 10 maps as one figure in r? i want aplabeltical orders (a) (b) etc below the maps and years 2009-2019 above the figure. I have given Greens as result colour of ndvi in R. Its not prominant. Which colour can i use to make it more visible?

library(raster)
library(rgdal)
library(rasterVis)
library(lattice)
library(RColorBrewer)

r1 <- raster("E:/GB extent(2009-2019)/rasterstack.tif", band=1)

r2 <- raster("E:/GB extent(2009-2019)/rasterstack.tif", band=2)
r3 <- raster("E:/GB extent(2009-2019)/rasterstack.tif", band=3)
r4 <- raster("E:/GB extent(2009-2019)/rasterstack.tif", band=4)
r5 <- raster("E:/GB extent(2009-2019)/rasterstack.tif", band=5)
r6 <- raster("E:/GB extent(2009-2019)/rasterstack.tif", band=6)
r7 <- raster("E:/GB extent(2009-2019)/rasterstack.tif", band=7)
r8 <- raster("E:/GB extent(2009-2019)/rasterstack.tif", band=8)
r9 <- raster("E:/GB extent(2009-2019)/rasterstack.tif", band=9)
r10 <- raster("E:/GB extent(2009-2019)/rasterstack.tif", band=10)

boudary <- readOGR(dsn= "E:/Shapefile/New folder/GB_Bondary.shp", layer =  "GB_Bondary")





raster_stack <- stack(r1,r2,r3,r4,r5,r6,r7,r8,r9,r10)

#levelplot(raster_stack, layout=c(5,2)


colr <- colorRampPalette((brewer.pal(3, 'Greens')))

levelplot(raster_stack,
          layout=c(5,2),
          margin=FALSE, 
          #main=list("title",col="Black",fontface='bold'),
          colorkey=list(
            space='right',                   
            labels=list(at=c(0.5,0.6,0.7,0.8,0.9,1.0), font=5,fontface='bold'),
            axis.line=list(col='black'),
            width=0.75
          ),    
          par.settings=list(
            strip.border=list(col='transparent'),
            strip.background=list(col='transparent'),
            axis.line=list(col='transparent')
          ),
          scales=list(draw=FALSE),
          col.regions=colr,                   
          at=seq(0.5,1.0, len=101),
          names.attr=(c('2009','2010','2011',
                        '2013','2014','2015','2016'
                        ,'2017','2018','2019')),font=list(face="bold"),cex=.4)+
  layer(sp.polygons(boudary, lwd=0.5))

###########################################################################

Answer 1

等待大神答复