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