Question

我需要在Python 3中编写一个函数，它返回一个矩形场（例如100x100点）上的位置数组（x，y），这些场根据同质空间泊松过程分散。

到目前为止，我已经使用Python代码找到了这个资源，但不幸的是，我无法为Python 3找到/安装scipy：

http://connor-johnson.com/2014/02/25/spatial-point-processes/

它帮助我理解泊松点过程究竟是什么以及它是如何工作的。

我一直在玩numpy.random.poisson一段时间了，但我很难解释它的回归。

http://docs.scipy.org/doc/numpy/reference/generated/numpy.random.poisson.html

>>> import numpy as np
>>> np.random.poisson(1, (1, 5, 5))
array([[[0, 2, 0, 1, 0],
        [3, 2, 0, 2, 1],
        [0, 1, 3, 3, 2],
        [0, 1, 2, 0, 2],
        [1, 2, 1, 0, 3]]])

我认为该命令的作用是创建一个5x5字段 =（1,5,5）并在该字段上以 lambda = 1 的速率散布对象。结果矩阵中显示的数字是物体位于该特定位置的概率。

根据同质空间泊松过程，如何在该5x5场上散布10个物体？我的第一个猜测是遍历整个数组并在每个位置上插入一个“3”的对象，然后在每个其他位置上插入一个“2”，依此类推，但我不确定实际的概率我应该用来确定是否应该插入一个对象。

根据以下资源，我可以通过简单地将速率和对象计数（10 * 1 = 10）相乘并使用该值作为我的lambda来模拟10个以1的速率散布在字段上的对象，即< / p>

>>> np.random.poisson(10, (1, 5, 5))
array([[[12, 12, 10, 16, 16],
        [ 8,  6,  8, 12,  9],
        [12,  4, 10,  3,  8],
        [15, 10, 10, 15,  7],
        [ 8, 13, 12,  9,  7]]])

但是，我不知道如何让事情变得更容易。我只是以这种方式增加对象出现的速率。

Poisson point process in matlab

总结一下，我的主要问题是：如何使用numpy.random.poisson(lam, size)对分散在二维字段n上的dx*dy个对象进行建模？

Answer 1

我似乎以错误的方式看待了这个问题。经过更多的离线研究后，我发现实际上足以创建一个代表对象数量的随机泊松值 n = np.random.poisson(100)并在0和1之间创建相同数量的随机值

x = np.random.rand(n)
y = np.random.rand(n)

现在我只需要将x和y值的两个数组连接到（x，y）元组的数组。那些是我正在寻找的随机位置。我可以将每个x和y值乘以我的字段的边长，例如100，将值缩放到我想要显示的100x100字段。

我认为＆＃34;随机性＆＃34;那些位置应该由随机泊松过程确定，但似乎只需要数位置需要由它确定，而不是实际的位置值。

Answer 2

完全正确。您绝对不需要SciPy，但是当我第一次在Python中模拟Poisson点过程时，我也使用了SciPy。我在这篇文章中为原始代码提供了模拟过程的详细信息：

https://hpaulkeeler.com/poisson-point-process-simulation/

我只是在最近的代码中使用NumPy：

import numpy as np; #NumPy package for arrays, random number generation, etc
import matplotlib.pyplot as plt #for plotting

#Simulation window parameters
xMin=0;xMax=1;
yMin=0;yMax=1;
xDelta=xMax-xMin;yDelta=yMax-yMin; #rectangle dimensions
areaTotal=xDelta*yDelta;

#Point process parameters
lambda0=100; #intensity (ie mean density) of the Poisson process

#Simulate a Poisson point process
numbPoints = np.random.poisson(lambda0*areaTotal);#Poisson number of points
xx = xDelta*np.random.uniform(0,1,numbPoints)+xMin;#x coordinates of Poisson points
yy = yDelta*np.random.uniform(0,1,numbPoints)+yMin;#y coordinates of Poisson points

代码也可以在这里找到：

https://github.com/hpaulkeeler/posts/tree/master/PoissonRectangle

我还在那里上传了更多的Python（以及MATLAB和Julia）代码，用于模拟多个点过程，包括各种形状的Poisson点过程和聚类点过程。

https://github.com/hpaulkeeler/posts

Python 3中的泊松点过程，带有numpy，没有scipy

2 个答案: