# NOT RUN {
library(sp)
library(raster)
library(rgeos)
data(meuse)
coordinates(meuse) <- ~x+y
# create "known" locations
locs <- meuse[sample(1:nrow(meuse), 5),]
# systematic sample using extent polygon
e <- as(extent(meuse), "SpatialPolygons")
s <- background(e, p=1000, known=locs, d=300)
plot(s,pch=20)
points(locs, pch=20, col="red")
# systematic sample using irregular polygon
data(meuse.grid)
coordinates(meuse.grid) = c("x", "y")
gridded(meuse.grid) = TRUE
meuse.poly = gUnaryUnion(as(meuse.grid, "SpatialPolygons"))
s <- background(meuse.poly, p=1000, known=locs, d=200)
plot(s,pch=20)
plot(meuse.poly, add=TRUE)
points(locs, pch=20, col="red")
# random sample using irregular polygon
s <- background(meuse.poly, p=500, known=locs,
d=200, type="random")
plot(s,pch=20)
plot(meuse.poly, add=TRUE)
points(locs, pch=20, col="red")
# systematic sample using defined extent
extent(meuse)
s <- background(ext=c(178605, 181390, 329714, 333611),
p=1000, known=locs, d=300)
plot(s,pch=20)
points(locs, pch=20, col="red")
# }
Run the code above in your browser using DataLab