## Not run:
# library(sp)
# library(raster)
#
# data(meuse)
# coordinates(meuse) <- ~x+y
# proj4string(meuse) <- CRS("+init=epsg:28992")
#
# ## view different aspects of same data set
# m1 <- mapview(meuse, zcol = "soil", burst = TRUE)
# m2 <- mapview(meuse, zcol = "lead")
# m3 <- mapview(meuse, zcol = "landuse", map.types = "Esri.WorldImagery")
# m4 <- mapview(meuse, zcol = "dist.m")
#
# latticeView(m1, m2, m3, m4) # 4 panels
# sync(m1, m2, m3, m4) # 4 panels synchronised
# latticeView(m1, m2) # 2 panels, split vertical
# latticeView(m1, m2, ncol = 1) # 2 panels split horizontal
# sync(m1, m2, ncol = 1) # same but synchronised
# sync(m1, m2, m3, m4, sync = list(c(1, 2), c(3, 4))) # individual syncing
# sync(m1, m2, m3, m4, sync = list(c(1, 2, 4)))
#
# ## view all layers of raster stack
# map_list <- lapply(seq(nlayers(poppendorf)), function(i) {
# mapview(poppendorf[[i]], layer.name = names(poppendorf)[i])
# })
#
# latticeView(map_list, ncol = 5)
#
# ## view multiple data sets
# m1 <- mapview(meuse, zcol = "soil", burst = TRUE)
# m2 <- mapview(atlStorms2005, zcol = "Name")
# m3 <- mapview(poppendorf[[5]], use.layer.names = TRUE)
# m4 <- mapview(gadmCHE, color = "black")
#
# latticeView(m1, m2, m3, m4) # not synced
# sync(m1, m2, m3, m4) # synced
# sync(m1, m2, m3, m4, no.initial.sync = FALSE) # all maps zoomed to m4 extent
#
# ## End(Not run)
Run the code above in your browser using DataLab