添加具有高斯分布的泊松噪声答案

【问题标题】：Adding Poisson noise with a Gaussian distribution添加具有高斯分布的泊松噪声
【发布时间】：2013-10-27 22:08:47
【问题描述】：

我设法将泊松噪声添加到我的 .fits 图像中，但我需要添加像高斯分布的噪声，平均/中值 (mu_0) 为 0，分布范围越来越广 (sigma)。我找不到以这种方式控制添加噪音的语法，所以有人可以帮我看看吗？目前，我添加的泊松噪声是均匀分布的，这不是我想要的；我需要高斯随机性。

以下是相关代码：

    im = pf.open(name)
    isinstance(im,list)
    im0 = im[0]
    print im0.data.shape
    print np.var(im0.data)
    poissonNoise = np.random.poisson(poisson, im0.data.shape).astype(float)
    test = im0.data + poissonNoise
    print np.var(test)
    im0.data = test
    stringee = 'POISSON'
    pf.writeto(stringee+poisson+name, data=test, clobber=True, header=im0.header)
    check = pf.open(stringee+poisson+name)
    np.var(check[0].data)

【问题讨论】：

标签： python gaussian poisson fits

【解决方案1】：

对于高斯分布（正态分布），使用np.random.normal。

normalNoise = np.random.normal(center, scale, shape).astype(float)

在此处查看更多信息：NumPy Normal，更一般地查看所有类型：NumPy Random

【讨论】：

我试过这个，但是我得到了以下错误：---> 44 test = im0.data + poissonNoise 45 print np.var(test) 46 im0.data = test ValueError: operands could not与形状一起广播 (130,130) (2)
忘记转换为浮点数。如果这不起作用，您发送给normal 的确切参数是什么？