www.pudn.com > 领域平均法与均值滤波.rar > med_filtering.m

function med_filtering()  
X=imread('d:\1\LENA128.bmp','bmp') 
F=double(X); 
X=im2double(X); 
figure(10); 
imshow(im2uint8(X)); 
noise_x=imnoise(X,'salt & pepper',0.02); 
figure(15); 
imshow(im2uint8(noise_x)); 
[m,n]=size(X);  
output_image=X; 
for i=2:m-1 
    for j=2:n-1 
        x=[X(i-1,j-1),X(i-1,j),X(i-1,j+1),X(i,j-1),X(i,j),X(i,j+1),X(i+1,j-1),X(i+1,j),X(i+1,j+1)]; 
     output_image(i,j)=median(x);  
     end 
 end 
figure(100); 
imshow(im2uint8(output_image));