I read the 22 raster files and made the stack of all and then cropped them using shapefile (shp_Bhote). I divided the 22 rasters into four groups according to different time raster. I estimated the frequency of the values in raster for each group and made a four stack from them. And then I plotted the four different raster stack using level plot and layout function. However, the title of each plot on the graph comes automatically like layer.1, layer.2, layer.3, and layer.4. I want to remove this title from my plot, however, I can not figure out the issue. Help will be highly appreciated. I am sorry for the bad question.
Here is the code that I used for the plot.
List of tif files
library(raster)
library(sp)
library(rgdal)
library(lattice)
library(RColorBrewer)
library(rasterVis)
list_tif=list.files(getwd(),pattern=".tif")
a=stack(list_tif)
a_subset = subset(a, 1:8)
b_subset = subset(a,9:16)
c_subset = subset(a, 17:22)
shp_Bhote <- readOGR(dsn = ".", layer = "Bhote_Kosi")
Bhote = crop(a,extent(shp_Bhote), snap = "out")
Bhote_a = crop(a_subset,extent(shp_Bhote), snap = "out")
Bhote_b = crop(b_subset,extent(shp_Bhote), snap = "out")
Bhote_c = crop(c_subset,extent(shp_Bhote), snap = "out")
count_lake_Bhote=sum(Bhote)
count_lake_Bhote_a=sum(Bhote_a,na.rm=T)
count_lake_Bhote_b=sum(Bhote_b,na.rm=T)
count_lake_Bhote_c=sum(Bhote_c,na.rm=T)
frequency_Bhote=count_lake_Bhote/nlayers(Bhote)*100
frequency_Bhote_a=count_lake_Bhote_a/nlayers(Bhote_a)*100
frequency_Bhote_b=count_lake_Bhote_b/nlayers(Bhote_b)*100
frequency_Bhote_c=count_lake_Bhote_c/nlayers(Bhote_c)*100
#####remove NA Values
frequency_Bhote[frequency_Bhote==0]<-NA
frequency_Bhote_a[frequency_Bhote_a==0]<-NA
frequency_Bhote_b[frequency_Bhote_b==0]<-NA
frequency_Bhote_c[frequency_Bhote_c==0]<-NA
Bhote_stack <- stack(frequency_Bhote_a, frequency_Bhote_b, frequency_Bhote_c, frequency_Bhote)
breaks <- seq(0, 100, by=1)
cols<-colorRampPalette(c("cyan","blue","pink","red"))(length(breaks)-1)
levelplot(Bhote_stack,colorkey = TRUE, layout=c(4, 1),scales=list(x=list(draw=FALSE),y=list(draw=FALSE)), at=breaks, col.regions=cols, margin = F))+
layer(sp.polygons(shp_Bhote))