1

我正在研究R 3.0.1并为集群泊松过程进行​​了模拟,R 通常有一个默认区域,基本上是一个框,在下一张图片中你可以看到我的模拟:

在此处输入图像描述

到目前为止一切都很好,问题是我想做的是模拟相同distribution但使用地理区域,但我不知道如何更改参数以便使用地理坐标拥有不同的区域。例如:

在此处输入图像描述

总而言之,基本上我想做的是弄清楚如何将这个区域更改为更大的区域,以便进行相同的模拟但使用新区域。这是我尝试过的代码:

library(spatstat)

sim1 = rpoispp(100)
plot(sim1)
4

1 回答 1

6

你可以试试这个:

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)

最终地图

于 2013-06-27T06:11:32.093 回答