2013-10-27 63 views
0

我設法將泊松噪聲添加到我的.fits圖像,但我需要添加分佈像高斯分佈的平均/中值(mu_0)爲0和分佈越來越寬的噪聲(西格瑪)。我找不到用這種方式添加噪音的語法,那麼有人可以通過它來指導我嗎?目前,我添加的泊松噪聲是均勻分佈的,這不是我想要的;我需要高斯隨機性。添加具有高斯分佈的泊松噪聲

這裏是代碼的相關位:

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) 

回答

0

高斯分佈(正態分佈)使用np.random.normal

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

查看更多在這裏:NumPy Normal一般多爲所有類型:NumPy Random

+0

我想這一點,但我得到了以下錯誤: ---> 44測試= im0.data + poissonNoise 45打印np.var(測試) 46 im0.data =測試 ValueError:操作數無法與形狀一起廣播(130,130)(2) – vdogsandman

+0

忘記轉換爲浮點數。如果這不起作用,你發送給'normal'的確切參數是什麼? –