伽玛分布噪声图像
clc,clear,close all
warning off
feature jit off
im = imread('coloredChips.png');
Z1 = imnoise_gamma(size(im,1),size(im,2),1,0.1);
Z1 = im2uint8(Z1); % 類型轉換
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_gamma(M, N, a,b)
% input:
% 伽馬gamma分布，噪聲的類型；
% M,N：輸出噪聲圖像矩陣的大小
% a,b：各種噪聲的分布參數
% output:
% R: 輸出的噪聲圖像矩陣，數據類型為double型
% 設定默認值% 產生伽瑪分布噪聲if nargin <= 3a = 2; b = .5;endc=0.1;R = zeros(M, N); for i = 1:Mfor j=1:Nx= rand;if x>cR(i,j) = b.^a./(gam(a)) .*(x-c).^(a-1).*exp(-b.*(x-c));elseR(i,j)=0;endendend
endfunction T = gam(z)
umax= 10;
T=0;
for i =0:1:umaxT =T + i^(z-1)*exp(-i);
end
end
其中 ? , 稱為伽瑪函數，伽瑪函數是階乘在實數上的泛化，滿足性質 ? 。
伽瑪分布期望是 ? ，方差是 ? ，中數是
伽瑪分布具有加成性，即如果 ? ，那么
《新程序員》：云原生和全面數字化實踐50位技術專家共同創作，文字、視頻、音頻交互閱讀
《新程序員》：云原生和全面數字化實踐50位技術專家共同創作，文字、視頻、音頻交互閱讀
總結
