高斯去噪Matlab代码

I=imread('example.JPG');
%figure,imshow(I);
I_r=I(:,:,1);
I_g=I(:,:,2);
I_b=I(:,:,3);
I1 = rgb2gray(I);
II=I1;
[r,c]=find(II<180);
for i=1:size(r)
    II(r(i),c(i))=0;
end
I_r(~II) = 0;
I_g(~II) = 0;
I_b(~II) = 0;
I_1 = cat(3, I_r, I_g, I_b);%在彩图去对原图进行选择
figure,imshow(I_1);
I_r1=I_1(:,:,1);
I_g1=I_1(:,:,2);
I_b1=I_1(:,:,3);
%figure,imshow(I1);
I2 = imopen(I1, strel('diamond',3));
%figure,imshow(I2);
I2 = medfilt2(I2, [3 3]);
%figure,imshow(I2);
I2 = imclearborder(I2);
%figure,imshow(I2);
I2 = im2bw(I2, graythresh(I2));
%figure,imshow(I2);
I2 = imclose(I2, strel('disk',3));
%figure,imshow(I2);