function?[mumask]=kmeans（imak）
ima=imread（‘b.jpg‘）;
figure;
imshow（ima）;
k=3;
ima=double（ima）;
copy=ima;?????????
ima=ima（:）;???????
s=length（ima）;
%建立灰度直方圖
m=max（ima）+1;
h=zeros（1m）;
hc=zeros（1m）;??????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????
for?i=1:s
??if（ima（i）>0）?h（ima（i））=h（ima（i））+1;end;
end
ind=find（h）;
hl=length（ind）;
%初始化聚類中心
mu=（1:k）*m/（k+

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????917??2011-12-21?11:11??kmeans.m

-----------?---------??----------?-----??----

??????????????????917????????????????????1