0

我设法将泊松噪声添加到我的 .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)
4

1 回答 1

0

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

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

在此处查看更多信息:NumPy Normal,更一般地适用于所有类型:NumPy Random

于 2013-10-27T22:26:06.860 回答