matlab-均值濾波(原創(chuàng))

clearall;closeall;P=input('inputpicturedata:');d=size(P);if(d(3)>1)P=rgb2gray(P);endchoice=1;while(choice~=0)P=imnoise(P,'salt&pepper',0.02);subplot(221);imshow(P);title('加入椒鹽噪聲后的圖像');P=double(P);[mn]=size(P);g=zeros(m,n);m1=m+2;n1=n+2;PP=zeros(m1,n1);fori=2:m1-1forj=2:n1-1PP(i,j)=P(i-1,j-1);endendPP(1,2:n1-1)=P(2,1:n);PP(m1,2:n1-1)=P(m-1,1:n);PP(1:m1,1)=PP(1:m1,3);PP(1:m1,n1)=PP(1:m1,n1-2);disp('-----------------------------------------------');disp('主菜單');disp('TheNumberofNgeighbourhood-AverageingMethod');disp('4、4鄰域平均');disp('8、8鄰域平均，不考慮點(diǎn)本身亮度');disp('9、8鄰域平均，考慮點(diǎn)本身亮度');disp('----------------------------------------------');choice=input('Pleaseinputthenumberofmethod:(4,8,9),0表示結(jié)束\n')switch(choice)case4fori=2:m1-1forj=2:n1-1g(i,j)=(PP(i-1,j)+PP(i+1,j)+PP(i,j-1)+PP(i,j+1))/4;endendg=uint8(g);subplot(222);imshow(g);title('經(jīng)4鄰域平均去噪后的圖片');case8fori=2:m1-1forj=2:n1-1g(i,j)=(PP(i-1,j-1)+PP(i-1,j)+PP(i-1,j+1)+PP(i,j-1)+PP(i,j+1)+PP(i+1,j-1)+PP(i+1,j)+PP(i+1,j+1))/8;endendg=uint8(g);subplot(223);imshow(g);title('經(jīng)8鄰域平均,不考慮點(diǎn)本身亮度去噪后的圖片');case9fori=2:m1-1forj=2:n1-1g(i,j)=(PP(i-1,j-1)+PP(i-1,j)+PP(i-1,j+1)+PP(i,j-1)+PP(i,j)+PP(i,j+1)+PP(i+1,j-1)+PP(i+1,j)+PP(i+1,j+1))/9;endendg=uint8(g);subplot(224);imshow(g);title('經(jīng)8鄰域平均,考慮點(diǎn)本身亮度去噪后的圖片');