几何分布噪声图像
clc,clear,close all
warning off
feature jit off
im = imread('coloredChips.png');
Z0 = imnoise_geometry(im,size(im,1),size(im,2),0.9);
Z1 = im2uint8(Z0); % 類型轉(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_geometry(im,M, N, b)
% input:
% 幾何geometry分布，噪聲的類型；
% M,N：輸出噪聲圖像矩陣的大小
% a,b：各種噪聲的分布參數(shù)
% output:
% R: 輸出的噪聲圖像矩陣，數(shù)據(jù)類型為double型
% 設(shè)定默認(rèn)值
if nargin < 3b = 0.5;
end% 產(chǎn)生幾何分布噪聲for i = 1:Mfor j=1:Na = double( floor(im(i,j)/30)+1 );R(i,j) = b .* (1-b).^(a-1);endend
end
總結(jié)
