library(lattice)
p1 = xyplot(1:10 ~ 1:10)
p2 = xyplot(10:1 ~ 1:10)
(
p = latticeCombineGrid(
list(p1, p2)
)
)
if (requireNamespace("raster", quietly = TRUE)) {
# load data
# Use a probability map assuming high potential for city expansion is just
# resulting from proximity to current urban area:
prd = raster::raster(system.file("extdata/probability.rst", package = "Orcs"))
# observed city growth between 1990 and 2006
obs = raster::raster(system.file("extdata/citygrowth.tif", package = "Orcs"))
# masking current urban area since these pixels have no potential for change
msk = raster::raster(system.file("extdata/citymask.tif", package = "Orcs"))
# create data list
dat <- list(prd, obs, msk)
# create list of lattice plots
plist <- lapply(dat, raster::spplot, scales = list(draw = TRUE))
# # draw individually
# plist[[1]]
# plist[[2]]
# plist[[3]]
# combine to grid, using c(1, 3) layout
p = latticeCombineGrid(plist, layout = c(1, 3))
print(p)
}
Run the code above in your browser using DataLab