library(raster)
library(ggplot2)
library(rasterVis)
library(gstat)
# dummy grids, with spatial autocorrelation, normalized
xy <- expand.grid(1:50, 1:50)
names(xy) <- c('x', 'y')
set.seed(2) # 2
gdummy <- gstat(formula = z ~ 1, locations = ~x + y, dummy = TRUE, beta = 1,
model = vgm(psill = 0.35, range = 30, model = 'Sph'),
nmax = 20)
yy <- predict(gdummy, newdata = xy, nsim = 4)
gridded(yy) <- ~x + y
yy <- raster(yy)
yy <- focal(yy, w = matrix(1, 3, 3), mean, na.rm = TRUE, pad = TRUE)
png("inst/slides/figures/random-raster.png", height = 5, width = 5, res = 300,
units = "in", bg = "transparent")
lattice.options(layout.heights = list(bottom.padding = list(x = 0),
top.padding = list(x = 0)),
layout.widths = list(left.padding = list(x = 0),
right.padding = list(x = 0)))
levelplot(yy, scales = list(draw = FALSE), axes = FALSE,
colorkey = list(axis.line = list(col = "white"),
axis.text = list(col = "white")))
dev.off()