function?[?Result?]?=?K_means（?Kthreshold?）
%?此函數用于K-means圖像分類
%?K為待分類圖像

[RowColDim]=size（K）;
if?Dim>1?????????%如果大于1，就是多波段圖像，需要灰度化為單波段圖像
????A=rgb2gray（K）;?
else
????A=K;?????????%單波段圖像
end
B=double（A）;
Result=zeros（RowCol）;

n=input（‘請選擇類別總數：‘）;
for?i=1:n
????Cen（i）=min（min（B））+（max（max（B））-min（min（B）））*（i/（n+1））;
end

delta=Cen;
while?delta>threshold
????sum=zeros（1n）;
????count=zeros（1n）;
????avg=zeros（1n）;
????for?i=1:Row
????????for?j=1:Col
????????????Min=256;
????????????for?k=1:n
????????????????if?abs（Cen（k