数字图像去噪处理.docx

一．加噪 f=imread(mch.tif); figure(1);imshow(f); title(原始图像); g=imnoise(f,salt pepper,0.2); figure(2);imshow(g); title(椒盐噪声污染的图像); 二．图像去噪 均值滤波器 【1】均值滤波： 本例对一幅数字图像添加高斯噪声，并分别用3*3和7*7均值滤波模板进行噪声的滤除。所需程序代码如下： clear all %打开图像 I=imread(mch.tif); %创建均值滤波模板 H1=ones(3)/9; H2=ones(7)/49; %添加高斯噪声，均值为0，方差为0.02 J=imnoise(I,gaussian,0,0.02); %转换J为double数据类型 J=double(J); %均值滤波 G1=conv2(J,H1,same); G2=conv2(J,H2,same); %图像显示 subplot(2,2,1); imshow(I); title(原始图像); subplot(2,2,2); imshow(J,[]); title(添加高斯噪声图像); subplot(2,2,3); imshow(G1,[]); title(3*3均值滤波图像); subplot(2,2,4); imshow(G2,[]); title(7*7均值滤波图像); 【2】超限像素平滑法： clear all %打开图像 I=imread(pout.tif); [m n]=size(I);%求图像尺寸 T=50;%设定阈值 G=[]; %创建均值滤波模板 H1=ones(3)/9; %添加椒盐噪声 J=imnoise(I,salt pepper,0.05); %转换J为double数据类型 J=double(J); %均值滤波 G1=conv2(J,H1,same); %超限像素平滑 for i=1:m for j=1:n if abs(J(i,j)-G1(i,j))T G(i,j)=G1(i,j); else G(i,j)=J(i,j); end end end %图像显示 subplot(2,2,1); imshow(I); title(原始图像); subplot(2,2,2); imshow(J,[]); title(添加椒盐噪声图像); subplot(2,2,3); imshow(G1,[]); title(3*3均值滤波图像); subplot(2,2,4); imshow(G,[]); title(超限像素平滑图像); 2. 中值滤波 clear all %打开图像 I=imread(21.tif); %添加椒盐噪声 N1=imnoise(I,salt pepper,0.04); %添加高斯噪声 N2=imnoise(I,gaussian,0,0.02); %添加乘性噪声 N3=imnoise(I,speckle,0.02); G1=medfilt2(N1); G2=medfilt2(N2); G3=medfilt2(N3); %分别显示噪声图像和中值滤波图像 subplot(2,3,1); imshow(N1); title(添加椒盐噪声图像); subplot(2,3,2); imshow(N2); title(添加高斯噪声图像); subplot(2,3,3); imshow(N3); title(添加乘性噪声图像); subplot(2,3,4); imshow(G1); title(椒盐噪声中值滤波图像); subplot(2,3,5); imshow(G2); title(高斯噪声中值滤波图像); subplot(2,3,6); imshow(G3); title(乘性噪声中值滤波图像); 3. 维纳滤波 【1】 clear all I=imread(mch.tif); %添加高斯噪声 N1=imnoise(I,gaussian,0,0.00