r - 如何用图案生成空间点

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

对于a)我可以使用spsample()包中的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分？)

最佳答案

我只知道如何使用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)

关于r - 如何用图案生成空间点，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/33874879/