s083.m

简单的数字图像处理程序

% 灰度图像层叠中值滤波  s083
I83=imread ('2.jpg');
I830=rgb2gray(I83);
I831=imnoise (I830,'salt & pepper',0.05);
[m,n]=size(I830);
I832=zeros(m,n);
for k=1:255
    d=I831;
    for i=1:m
       for j=1:n
          if d(i,j)>=k
             I832(i,j)=1;
          else
             I832(i,j)=0;
          end
          I832=im2double(I832);
       end
    end
    I832=I832+medfilt2(I832);
end
subplot(1,3,1);imshow(I830),title('原图');
subplot(1,3,2);imshow(I831),title('噪声');
subplot(1,3,3);imshow(im2uint8(I832)),title('滤波')