r - 如何生成带有图案的空间点

Question

我正在做一些工作，我需要在多边形上生成 a) 随机空间点 b) 非随机空间点，即 b) 点概率取决于例如东西梯度或与某个点源的距离或者是其他东西

spsample()对于 a) 我可以使用包中的命令在多边形上生成随机点sp，如下所示：

# Load a spatial polygon from maptools package
library(maptools)
nc <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1],        proj4string=CRS("+proj=longlat +datum=NAD27"))
plot(nc)

library(sp)
pts <- spsample(nc, 100, type="random")
plot(nc)
points(pts, pch=19, col="red")

这正是我想要的 a)。但是，这是否可以针对 b) 进行修改，例如，东方的点比西方更有可能？（虽然仍然能够指定我想要 100 分？）

Answer 1

我只知道如何使用spatstat包来做到这一点。使用函数rpoint您可以使用坐标 x,y 的任何函数来定义点的不均匀密度。在这里，我将函数定义为在该区域的西端具有值 0，并以斜率 100 向东线性增长：

library(maptools)
library(spatstat)
nc <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1],
                    proj4string=CRS("+proj=longlat +datum=NAD27"))
nc <- as.owin(nc)
west0 <- nc$xrange[1]
f <- function(x, y, ...){ 100 * (x - west0) }
pts <- rpoint(1000, f, win = nc)
plot(pts)