close?all
clear?all
A=imread（‘1.jpg‘）;
figure;
imshow（A）;
title（‘原圖‘）;
[u?v?w]=size（A）;

cform=makecform（‘srgb2lab‘）;
lab_A=applycform（Acform）;?%?轉換到lab空間，對于顏色分割的吧?lab空間相互分量聯系性比較小?利于分割
ab?=?double（lab_A（::2:3））;
nrows?=?size（ab1）;
ncols?=?size（ab2）;
ab?=?reshape（abnrows*ncols2）;%矩陣轉換類型轉換為double型
nColors?=4;
%?采用k-means方法實現聚類，重復三次
[cluster_idx?cluster_center]?=?kmeans（abnColors‘distance‘‘sqEuclidean‘?‘Replicates‘3）;?
%對不同的類別進行標志
pixel_labels?=?reshape（cluster_idxnrowsncols）;
%?figure;
showImage_color（AnColorscluster_idx）;
%?imshow（pixel_labels[]）?title（‘聚類后圖片‘）;

%分類后的矩陣
segmented_images?=?cell（14）;
rgb_label?=?repmat（pixel_labels[1?1?3]）;
for?k?=?1:nColors
????color?=?A;
????color（rgb_label?~=?k）?=?0;
????segmented_images{k}?=?color;
end
figure;
subplot（221）imshow（segmented_images{1}）title（‘cluster?1‘）;
subplot（222）imshow（segmented_images{2}）title（‘cluster?2‘）;
subplot（223）imshow（segmented_images{3}）title（‘cluster?3‘）;
subplot（224）imshow（segmented_images{4}）title（‘cluster?4‘）;
%?imshow（segmented_images{1}）?title（‘cluster?1‘）;
%
figure;
imshow（segmented_images{2}）?title（‘cluster?2‘）;
%?f