matlab编写的中值滤波器

%中值滤波 clc,clear; f=imread('2.bmp'); subplot(221),imshow(f); f1=imnoise(f,'gaussian',0.002,0.0008); subplot(222),imshow(f1); 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(223),imshow(Y); %均值滤波 k=floor(3*3/2)+1; [M,N]=size(X); uint8 Z=zeros(M,N); funBox=zeros(3,3); temp=zeros(3*3); for i=1:M-3 for j=1:N-3 funBox=X(i:i+3,j:j+3); temp=funBox(:); tempSort=sort(temp); Z(i,j)=tempSort(k); end; end; subplot(224),imshow(Z);