中值滤波3*3滤波窗口

程序如下：

>> i=imread('Winter.jpg'); %读入图像

subplot(2,2,1);

imshow(i); %显示原始图像

title('original'); %设置图像标题

j = imnoise(i,'salt & pepper',0.02); %加均值为0，方差为0.02的椒盐噪声

subplot(2,2,2);

imshow(j); %显示处理后的图像

title('加入椒盐噪声');

k = imnoise(i,'gaussian',0.02); %加均值为0，方差为0.02的加入高斯噪声

l = imnoise(k,'salt & pepper',0.02); %加均值为0，方差为0.02的椒盐噪声

subplot(2,2,3);

imshow(l);

title('加入椒盐噪声和高斯噪声');

r=l(:,:,1);

g=l(:,:,2);

b=l(:,:,3);

[m,n]=size(r);

for i=1:m

for j=1:n

if(i==1|i==m|j==1|j==n)

G(i,j)=r(i,j);

else

hao(1)=r(i-1,j-1);

hao(2)=r(i-1,j);

hao(3)=r(i-1,j+1);

hao(4)=r(i,j-1);

hao(5)=r(i,j);

hao(6)=r(i,j+1);

hao(7)=r(i+1,j-1);

hao(8)=r(i+1,j);

hao(9)=r(i+1,j+1);

hao=sort(hao);

G(i,j)=hao(5);

end

end

end

>> for i=1:m

for j=1:n

if(i==1|i==m|j==1|j==n)

K(i,j)=g(i,j);

else

hao(1)=g(i-1,j-1);

hao(2)=g(i-1,j);

hao(3)=g(i-1,j+1);

hao(4)=g(i,j-1);

hao(5)=g(i,j);

thao(6)=g(i,j+1);

hao(7)=g(i+1,j-1);

hao(8)=g(i+1,j);

hao(9)=g(i+1,j+1);

hao=sort(hao);

K(i,j)=hao(5);

end

end

end

>> for i=1:m

for j=1:n

if(i==1|i==m|j==1|j==n)

L(i,j)=b(i,j);

else

hao(1)=b(i-1,j-1);

hao(2)=b(i-1,j);

hao(3)=b(i-1,j+1);

hao(4)=b(i,j-1);

hao(5)=b(i,j);

hao(6)=b(i,j+1);

hao(7)=b(i+1,j-1);

hao(8)=b(i+1,j);

hao(9)=b(i+1,j+1);

hao=sort(hao);

L(i,j)=hao(5);

end

end

end

>> l(:,:,1)=G;

>> l(:,:,2)=K;

>> l(:,:,3)=L;

>> subplot(2,2,4);

>> imshow(l);

用5*5的滤波窗口进行中值滤波处理：

程序如下:

>> i=imread('Winter.jpg'); %读入图像

subplot(2,2,1);

imshow(i); %显示原始图像

title('original');

j = imnoise(i,'salt & pepper',0.02); %加均值为0，方差为0.02的椒盐噪声。

subplot(2,2,2);

imshow(j); %显示处理后的图像

title('加入椒盐噪声');

k = imnoise(i,'gaussian',0.02); %加均值为0，方差为0.02的加入高斯噪声

l = imnoise(k,'salt & pepper',0.02); %加均值为0，方差为0.02的椒盐噪声

subplot(2,2,3);

imshow(l);

title('加入椒盐噪声和高斯噪声');

r=l(:,:,1);

g=l(:,:,2);

b=l(:,:,3);

[a,b]=size(r);

for i=1:a

for j=1:b

if(i==1|i==2|i==a-1|i==a|j==1|j==2|j==b-1|j==b)

P(i,j)=r(i,j);

else

hao(1)=r(i-2,j-2);

hao(2)=r(i-2,j-1);

hao(3)=r(i-2,j);

thao(4)=r(i-2,j+1);

hao(5)=r(i-2,j+2);

hao(6)=r(i-1,j-2);

hao(7)=r(i-1,j-1);

hao(8)=r(i-1,j);

hao(9)=r(i-1,j+1);

hao(10)=r(i-1,j+2);

hao(11)=r(i,j-2);

hao(12)=r(i,j-1);

hao(13)=r(i,j);

hao(14)=r(i,j