kmean之matlab

clc;
close all;
n=2;
a=imread('3.jpg');
a=rgb2hsv(a);
M=size(a,1);N=size(a,2);
for i=1:100
b=double(a);b=b(:);
b=reshape(b,M*N,3);
[A,B]=kmeans(b,n);
for j=1:n
index=find(A==j);
b(index,:)=ones(length(index),1)*[rand rand rand];
end
b=reshape(b,M,N,3);
imshow(mat2gray(b));pause(0.1);
end

——————————————————–

clear all;
clc;
I_rgb=imread('dog.jpg');
figure();imshow(I_rgb);title('原始图像');
%去噪
filter=ones(