% 真彩圖像層疊中值濾波  s084
clear;
I83=imread('2.jpg');
I830=imnoise(I83,'salt & pepper',0.05);
I831=I830(:,:,1);
I832=I830(:,:,2);
I833=I830(:,:,3);
[m,n]=size(I831);
I834=zeros(m,n);
I835=zeros(m,n);
I836=zeros(m,n);
for k=1:255
    d1=I831;
    d2=I832;
    d3=I833;
    for i=1:m
       for j=1:n
            if d1(i,j)>=k
             d1(i,j)=1;
               else
              d1(i,j)=0;
            end
            if d2(i,j)>=k
             d2(i,j)=1;
               else
              d2(i,j)=0;
            end
            if d3(i,j)>=k
              d3(i,j)=1;
               else
              d3(i,j)=0;
            end
       end
    end
    d1=im2double(d1);
    d2=im2double(d2);
    d3=im2double(d3);
    d1=medfilt2(d1);
    d2=medfilt2(d2);
    d3=medfilt2(d3);
    I834=I834+d1;
    I835=I835+d2;
    I836=I836+d3;
end
subplot(1,3,1);imshow(I83),title('原圖');
subplot(1,3,2);imshow(I830),title('噪聲');
subplot(1,3,3);imshow(cat(3,I834,I835,I836)),title('濾波')