R/samplePts.R

In geospt: Geostatistical Analysis and Design of Optimal Spatial Sampling Networks

samplePts <- function(x, n, type, ...){
if(type=="hexagonal") {
pts.est <- spsample(x, 1.2*n, type, ...)
while (summary(pts.est)[[1]] != n) pts.est <- spsample(x, 1.2*n, type)
} else {
pts.est <- spsample(x, n, type, ...)
while (summary(pts.est)[[1]] != n) pts.est <- spsample(x, n, type)
}
return(pts.est)
}