在R中的同一图中用两个样本模拟雷场

Question

我试图通过在同一图中绘制两个泊松分布样本来模拟雷区，一个样本的强度更高，面积小于另一个样本。这是雷区，另一个只是噪点（石头，洞，金属）被视为点。我不能让R用轴绘制相同单位的点。无论我做什么，这些点都会跨越整个情节，即使我只希望X点覆盖四分之一的情节。我的R代码如下：

library(spatstat)
Y = rpoispp(c(5),win=owin(c(0,10),c(0,10)))
X = rpoispp(c(10),win=owin(c(0,5),c(0,5)))

如果你能帮助我，请告诉我。

Answer 1

我的猜测是你做的事情如下：

> plot(Y)
> plot(X)

绘制积分。

这个问题是类ppp（这是rpoispp函数返回的）的plot函数的默认行为是创建一个只包含其点的新图。因此，第二个plot调用基本上会删除第一个绘图，并在不同缩放的窗口中绘制自己的点。您可以通过为第二个图设置选项add=TRUE来覆盖此行为。所以代码

> plot(Y)
> plot(X, add=TRUE, cols="red")

应该为您提供类似：

的内容

查看文档（help(plot.ppp)）以获取更多解释和其他选项来美化情节。