空域平滑 中值滤波代码

clear all;%清零所有变量
I=imread('pout.tif');%读入要分析的图
J=I;
[m,n]=size(I); 
T=40;
T3=zeros(m,n);
T4=zeros(m,n);
J=imnoise(I,'salt & pepper',0.01);%调用函数 添加椒盐噪声
subplot(4,3,10);
imshow(J);
xlabel('椒盐噪声后图像');
subplot(431);
imshow(I);
xlabel('原始图像');
I=im2double(uint8(I));%改成双精度型，进行数据转换
M=fft2(I);%进行二维傅里叶变换
N=fftshift(M);%将二维傅里叶变换的零频率移到频谱的中心
subplot(432);
imshow((log(1+abs(N))),[]);%绘制图像的对数幅度谱
xlabel('原图像傅里叶谱');
v1=real(120*M./abs(M));
subplot(433);
imshow(v1);
xlabel('原图像的相位谱');
for x=2:m-1                      
    for y=2:n-1  
        g(x,y)=round([0.125*J(x,y+1)+0.125*J(x,y-1)+0.125*J(x-1,y)+0.125*J(x+1,y)+0.125*J(x-1,y-1)+0.125*J(x-1,y+1)+0.125*J(x+1,y-1)+0.125*J(x+1,y+1)]);
       if (abs(J(x,y)- g(x,y))>T)   
          T3(x,y)=g(x,y);  
        else