用高斯分布增加泊松噪声

时间:2013-10-27 22:08:47

标签: python gaussian poisson fits

我设法为我的.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)

1 个答案:

答案 0 :(得分:0)

对于高斯分布(正态分布),请使用np.random.normal

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

在此处详细信息:NumPy Normal以及更常见的所有类型:NumPy Random