%%%%%%%%%%%%%%%%%%%%%%%%%%%%%圖像的讀入%%%%%%%%%%%%%%%%%%%%%%%%
[filenamepathname]=uigetfile（‘*.bmp‘‘select?a?picture‘）;
?if?（filename）
?????filename=strcat（pathnamefilename）;
?????filename=lower（filename）;
?else
?msgbox（‘please?select?a?image‘）;
?return;
?end
I=imread（filename）;
?if?（size（I3）==3）
?????I=rgb2gray（I）;
?end
%%%%%%%%%%%%%%%%%1..樣本點為所有的像素點%%%%%%%%%%
%%%%%%%%%%%%%%%?2...聚類中心點的確定%%%%%%%%%%%%%%
??????figureimshow（I[]）title（‘請選取預期分類數目個聚類中心點--->‘）;
??????[mn]=size（I）;
????????k=4;?????????????????????????%?k?用來記錄預期聚類的數目這里假定為2
????????c=k;?????????????????????????%?c來記錄實際分類時的數目
????????hold?on;
???????[pq]=getpts;?????????????????%獲取聚類中心，返回去類中心的坐標，p存放第一個坐標點，q存放第二個坐標點（數對）
??????J=zeros（13*c）;????????????????%？
??????for?i=1:c
???????f=round（q（i））;?????????????????%分為c類????J來記錄各類的中心像素值...（坐標和像素值）
???????d=round（p（i））;
???????J（i）=I（fd）;????????????????????
??????end
???????%disp（J）;
%%%%%%%%%%%%%%%%%3...將樣本點分到各個聚類中心去%%%%%%%%%%%%%
max_gen=10;?????