kmeans-圖像分割代碼

1、he=imread('02.png');%讀入圖像imshow(he),title('H&Eimage');text(size(he,2),size(he,1)+15,.'ImagecourtesyofAlanPartin,JohnsHopkinsUniversity',.'FontSize',7,'HorizontalAlignment','right');cform=makecform('srgb2lab');%色彩空間轉(zhuǎn)換lab_he=applycform(he,cfo

2、rm);ab=double(lab_he(:,:,2:3);%數(shù)據(jù)類型轉(zhuǎn)換nrows=size(ab,1);%求矩陣尺寸ncols=size(ab,2);%求矩陣尺寸ab=reshape(ab,nrows*ncols,2);%矩陣形狀變換nColors=3;%重復(fù)聚類3次，以避免局部最小值cluster_idxcluster_center=kmeans(ab,nColors,'distance','sqEuclidean',.'Replicates',3);pixel_labels=reshape(cluster_idx,nrows,ncols)

3、;%矩陣形狀改變imshow(pixel_labels,);%顯示圖像title('imagelabeledbyclusterindex');%設(shè)置圖像標(biāo)題segmented_images=cell(1,3);%細(xì)胞型數(shù)組rgb_label=repmat(pixel_labels,113);%矩陣平鋪fork=1:nColorscolor=he;color(rgb_label=k)=0;segmented_imagesk=color;endimshow(segmented_images1),%顯示處理后的圖像title('objectsincluster1');

4、%設(shè)置圖像標(biāo)題imshow(segmented_images2),%顯示處理后的圖像title('objectsincluster2');%設(shè)置圖像標(biāo)題imshow(segmented_images3),%顯示處理后的圖像title('objectsincluster3');%設(shè)置圖像標(biāo)題mean_cluster_value=mean(cluster_center,2);tmp,idx=sort(mean_cluster_value);blue_cluster_num=idx(1);L=lab_he(:,:,1);blue_idx=find(pixel_labels=blue_cluster_num);L_blue=L(blue_idx);is_light_blue=im2bw(L_blue,graythresh(L_blue);%圖像黑白轉(zhuǎn)換nuclei_labels=repmat(uint8(0),nrowsncols);%矩陣平鋪nuclei_labels(blue_idx(is_light_blue=false)=1;nuclei_labels=repmat(nuclei_labels,113);%矩陣平鋪blue_nuclei=he;顯示處理后的圖像blue