# NOT RUN {
# generate random images
img <- matrix(1:16, ncol = 4, byrow = TRUE)
r <- raster(img)
# assign the limit of the data in the raster stack
r2 <- genLimitRasterRange(r, mn = 4, mx = 10)
r3 <- genLimitRasterRange(r, mn = 4, mx = 10, rm.values = TRUE)
# plot limited data
spplot(stack(r, r2 ,r3))
# }
Run the code above in your browser using DataLab