数字图像处理---中值/最大值/最小值滤波的实现以及测试对不同噪声的结果。

数字图像处理的滤波器，基于空域的处理。

代码实现如下，有问题请留言，如果有帮助的话，评论支持一下。

function zy=filter_zy(image,option,num)
%image是输入的图片，option是使用哪种滤波器的选项，num指滤波器模板的大小。
image=im2double(image);
[r,c]=size(image);
jieguo=zeros(r,c);
n=floor(num/2);
a=padarray(image,[n,n],'symmetric');

[M,N]=size(a);
    for i=1+n:M-n
        for k=1+n:N-n
            b=a(i-n:i+n,k-n:k+n);
            b=reshape(b,1,numel(b));
            jieguo(i-n,k-n)=median(b,2);
        end
    end

if option==0%最小值滤波
    for i=1+n:M-n
        for k=1+n:N-n
            b=a(i-n:i+n,k-n:k+n);
            %b=reshape(b,1,numel(b));
            jieguo(i-n,k-n)=min(min(b));
        end
    end
end
if option==9%最大值滤波
    for i=1+n:M-n
        for k=1+n:N-n
            b=a(i-n:i+n,k-n:k+n);
            %b=reshape(b,1,numel(b));
            jieguo(i-n,k-n)=max(max(b));
        end
    end
end
zy=double(jieguo);

代码测试对不同噪声的效果：

ii=imread('cameraman.tif');
%ii = imnoise(ii, 'gaussian');
ii= imnoise(ii,'salt & pepper',0.01);
subplot(2,4,1);imshow(ii),title('添加椒盐噪声后');
%a=fspecial('average',15);
 %Y=filter2(a,ii)/255;  
 %用生成的滤波器进行滤波,并归一化
 Y=ordfilt2(ii,5,ones(3,3));
 yy=filter_zy(ii,5,3);
 yyMin=filter_zy(ii,0,3);
 yyMax=filter_zy(ii,9,3);
 
 Ymin=ordfilt2(ii,1,ones(3,3));
Ymax=ordfilt2(ii,9,ones(3,3));
subplot(2,4,2);imshow(Y),title('系统的的中值滤波：'); 
subplot(2,4,3);imshow(Ymin),title('系统的的最小值滤波：');
subplot(2,4,4);imshow(Ymax),title('系统的的最大值滤波：');%显示滤波后的图象
subplot(2,4,6);imshow(yy),title('自编的中值滤波:');
subplot(2,4,7);imshow(yyMin),title('自编的最小值滤波：');
subplot(2,4,8);imshow(yyMax),title('自编的最大值滤波：');

如果有帮助的话，还请评论一下，鼓励一下作者0.0