编写用均值滤波去噪的matlab程序,用两种方法实现.(重谢)
方法一:filter2
clear all;
I=imread('lena.bmp');
%读入预处理图像
imshow(I)
%显示预处理图像
K1=filter2(fspecial('average',3),I)/255;
%进行3*3均值滤波
K2=filter2(fspecial('average',5),I)/255;
%进行5*5均值滤波
K3=filter2(fspecial('average',7),I)/255;
%进行7*7均值滤波
figure,imshow(K1)
figure,imshow(K2)
figure,imshow(K3)
方法二:双循环语句,移动平均法
%均值滤波
clc,clear;
f=imread('lena.bmp');
subplot(121),imshow(f),title('原图');
f1=imnoise(f,'gaussian',0.002,0.0008);
%subplot(222),imshow(f1),title('添加高斯噪声图');
k1=floor(3/2)+1;
k2=floor(3/2)+1;
X=f1;
[M,N]=size(X);
uint8 Y=zeros(M,N);
funBox=zeros(3,3);
for i=1:M-3
for j=1:N-3
funBox=X(i:i+3,j:j+3);
s=sum(funBox(:));
h=s/9;
Y(i+k1,j+k2)=h;
end;
end;
Y=Y/255;
subplot(122),imshow(Y),title('均值滤波');
实现图:
1、双循环语句,移动平均法。
<p>双循环语句,移动平均法</p> <p>%均值滤波</p>
<p>clc,clear;</p>
<p>f=imread('lena.bmp');</p>
<p>subplot(121),imshow(f),title('原图');</p>
<p>f1=imnoise(f,'gaussian',0.002,0.0008);</p>
<p>%subplot(222),imshow(f1),title('添加高斯噪声图');</p>
<p>k1=floor(3/2)+1;</p>
<p>k2=floor(3/2)+1;</p> <p>X=f1;</p>
<p>[M,N]=size(X);</p> <p>uint8 Y=zeros(M,N);</p>
<p>funBox=zeros(3,3);</p> <p>for i=1:M-3 </p>
<p> for j=1:N-3 </p> <p> funBox=X(i:i+3,j:j+3); </p>
<p> s=sum(funBox(:));</p> <p> h=s/9; </p>
<p> Y(i+k1,j+k2)=h; </p> <p> end;</p>
<p>end;</p> <p>Y=Y/255;</p>
<p>subplot(122),imshow(Y),title('均值滤波');</p>
<p>实现图:</p>
2、filter2。
<p>filter2</p>
<p>clear all;</p>
<p>I=imread('lena.bmp');</p>
<p>%读入预处理图像</p>
<p>imshow(I)</p> <p>%显示预处理图像</p>
<p>K1=filter2(fspecial('average',3),I)/255;</p> <p>%进行3*3均值滤波</p>
<p>K2=filter2(fspecial('average',5),I)/255;</p> <p>%进行5*5均值滤波</p>
<p>K3=filter2(fspecial('average',7),I)/255;</p> <p>%进行7*7均值滤波</p>
<p>figure,imshow(K1)</p> <p>figure,imshow(K2)</p>
<p>figure,imshow(K3)</p>