泊松分布噪声图像
泊松分布的概率函數(shù)為：
泊松分布的參數(shù)λ是單位時(shí)間(或單位面積)內(nèi)隨機(jī)事件的平均發(fā)生率。 泊松分布適合于描述單位時(shí)間內(nèi)隨機(jī)事件發(fā)生的次數(shù)。
clc,clear,close all
warning off
feature jit off
im = imread('coloredChips.png');
Z1 = imnoise_poission(im,size(im,1),size(im,2),0.3);
Z1 = im2uint8(Z1); % 類型轉(zhuǎn)換
figure('color',[1,1,1]),
im(:,:,1) = im(:,:,1) + Z1; % R
im(:,:,2) = im(:,:,2) + Z1; % G
im(:,:,3) = im(:,:,3) + Z1; % B
subplot(121); imshow(im);title('加泊松分布噪聲圖像')
subplot(122); imhist(Z1); title('加泊松分布噪聲圖像直方圖')
function R = imnoise_poission(im,M, N, lamda) % input: % 泊松poission分布，噪聲的類型； % M,N：輸出噪聲圖像矩陣的大小 % a,b：各種噪聲的分布參數(shù) % output: % R: 輸出的噪聲圖像矩陣，數(shù)據(jù)類型為double型 % 設(shè)定默認(rèn)值 if nargin < 4lamda = 0.5; end% 產(chǎn)生poission分布噪聲for i = 1:Mfor j=1:Nb=1;c = double( floor(im(i,j)/30)+1 );for k =1:cb=b*k;endR(i,j) = exp(-lamda).*lamda.^(c)./b;endend end
function R = imnoise_poission(im,M, N, lamda) % input: % 泊松poission分布，噪聲的類型； % M,N：輸出噪聲圖像矩陣的大小 % a,b：各種噪聲的分布參數(shù) % output: % R: 輸出的噪聲圖像矩陣，數(shù)據(jù)類型為double型 % 設(shè)定默認(rèn)值 if nargin < 4lamda = 0.5; end% 產(chǎn)生poission分布噪聲for i = 1:Mfor j=1:Nb=1;c = double( floor(im(i,j)/30)+1 );for k =1:cb=b*k;endR(i,j) = exp(-lamda).*lamda.^(c)./b;endend end
總結(jié)
