t分布噪声图像
t分布
由于在實際工作中,往往σ是未知的,常用s作為σ的估計值,為了與u變換區別,稱為t變換,統計量t 值的分布稱為t分布。 假設X服從標準正態分布N(0,1),Y服從χ2(n)分布,那么Z=X/sqrt(Y/n)的分布稱為自由度為n的t分布,記為 Z~t(n)。 clc,clear,close all warning off feature jit off im = imread('coloredChips.png'); Z1 = imnoise_t(size(im,1),size(im,2),3); Z1 = uint8(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('加t分布噪聲圖像') subplot(122); imhist(Z1); title('加t分布噪聲圖像直方圖')function R = imnoise_t(M, N, a) % input: % t分布,噪聲的類型; % M,N:輸出噪聲圖像矩陣的大小 % a,b:各種噪聲的分布參數 % output: % R: 輸出的噪聲圖像矩陣,數據類型為double型 % 設定默認值 if nargin < 1a = 1; end% 產生F分布噪聲R1=zeros(M,N);R2=zeros(M,N);for i=1:aR1 = 5*randn(M, N).^2;R2 = 5*randn(M, N).^2;R = R1./sqrt(R2./a);end end
總結