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