我有这段代码,它是更大代码的一部分。
newimg = pf.getdata(name)
print newimg
newhdr = pf.getheader(name)
normalNoise = np.random.normal(0,poisson)
test = newimg + normalNoise
pf是“pyfits”,“name”是我读过的.fits文件的名称。我想把这个.fits文件添加到其中,以便噪声量随机分布在所有像素中。我希望能够为更大的标准偏差添加更多的噪音;换句话说,随着我增加变量“泊松”的值,图像应该变得越来越不可辨别。我该怎么做呢?