如何用R模拟任意区域的空间泊松过程？

我正在研究

R3.0.1

，我对聚集泊松过程进行了模拟，R通常有一个默认区域，基本上是一个方框，在下一张图片中，您可以看到我的模拟：

到目前为止一切正常，问题是我想做的是模拟相同的

分布

，但使用地理区域，但我不知道如何更改参数，以便使用地理坐标获得不同的区域。例如：

总而言之，我想做的基本上是找出如何将这个区域改变为一个更大的区域，以便用新的区域进行相同的模拟。以下是我尝试的代码：

library(spatstat)

sim1 = rpoispp(100)
plot(sim1)

您可以尝试以下方法：

require(spatstat)
require(maps)

lambda = 0.1 # intensity of the process
lon = c(-100,-70) # domain's longitude
lat = c(-40,10)   # domain's latitude

sim = rpoispp(lambda, win=c(lon,lat))

# do the plot
par(mar=c(1,1,1,1))
map("world", xlim=lon, ylim=lat, fill=TRUE)
map.axes() # add axes
plot(sim, chars=19, cols="red",cex=0.5, add=TRUE)

# add other process

lon1 = c(-95,-85) # other area
lat1 = c(-5,5)
sim1 = rpoispp(5*lambda, win=c(lon1,lat1))
plot(sim1, chars=19, cols="blue",cex=0.5, add=TRUE)

您可以尝试以下方法：

require(spatstat)
require(maps)

lambda = 0.1 # intensity of the process
lon = c(-100,-70) # domain's longitude
lat = c(-40,10)   # domain's latitude

sim = rpoispp(lambda, win=c(lon,lat))

# do the plot
par(mar=c(1,1,1,1))
map("world", xlim=lon, ylim=lat, fill=TRUE)
map.axes() # add axes
plot(sim, chars=19, cols="red",cex=0.5, add=TRUE)

# add other process

lon1 = c(-95,-85) # other area
lat1 = c(-5,5)
sim1 = rpoispp(5*lambda, win=c(lon1,lat1))
plot(sim1, chars=19, cols="blue",cex=0.5, add=TRUE)

---

您的意思是，您想模拟在该地图中发生的过程吗？也就是说，不仅仅是一个简单的矩形区域？它可能是一个包含地图的矩形区域。Form reading

？RPOISP

我了解到有一个参数

win

，它接受类“owin”的对象，并且是“模拟图案的窗口”。你的意思是，您想模拟在该地图中发生的过程吗？也就是说，不仅仅是一个简单的矩形区域？嗯，它可能是一个包含地图的矩形区域。Form reading

？RPOISP

我了解到有一个参数

win

，它接受类“owin”的对象，并且是“模拟模式的窗口”。它工作得非常好。。。现在我有一个小问题，我可以把轴坐标设置为一个形状文件吗？我不知道，但是有一个shapefiles包。很高兴知道。这真是太好了。。。现在我有一个小问题，我可以把轴坐标设置为一个形状文件吗？我不知道，但是有一个shapefiles包。很高兴知道。