均值滤波函数

%============================================================================================================
%author:王同乐
%date:2016/10/30
img=imread('lena.png');
img=rgb2gray(img);
img=im2double(img)%输入图像类型为uint8,将其转换为double类型进行运算
gau=imnoise(img,'gaussian');
subplot(2,2,1);imshow(gau);title('添加高斯噪声的图像');%imshow()函数在显示图像时会自动将类型转换为uint8
subplot(2,2,2);imshow(median_filter(gau,3));title('滤除高斯噪声后的图像');
salt=imnoise(img,'salt & pepper');
subplot(2,2,3);imshow(salt);title('添加椒盐噪声的图像');
subplot(2,2,4);imshow(median_filter(salt,3));title('滤除椒盐噪声后的图像');

%================================================================================================================
%均值滤波开始
subplot(2,2,1);imshow(gau);title('添加高斯噪声的图像');%imshow()函数在显示图像时会自动将类型转换为uint8
subplot(2,2,2);imshow(avg_filter(gau,3));title('均值滤波除去高斯噪声的图像');%imshow()函数在显示图像时会自动将类型转换为uint8
subplot(2,2,3);imshow(salt);title('添加椒盐噪声的图像');
subplot(2,2,4);imshow(avg_filter(salt,3));title('均值滤除椒盐噪声后的图像');

function d=avg_filter(x,n)
filter=ones(n);
d=x;
[width,height]=size(x);
for ii=1:width-(n-1)
    for jj=1:height-(n-1)
    tmp=(x(ii:ii+n-1,jj:jj+n-1).*filter);%注意这里必须用x,x为原来的矩阵,而d已经为经过部分处理的矩阵
    avg=sum(sum(tmp))/(n*n);
    d(ii+(n-1)/2,jj+(n-1)/2)=avg;
    end
end 
%均值滤波函数

效果图：