kmeans算法matlab实现

% K-mean方法聚类
clear all;

%%%%%%%%%%%%%%%%%%%读取7幅遥感图像数据%%%%%%%%%%%%%%%%%%
for i=1:7
    I(i,:,:)=imread(['L',num2str(i),'.bmp']);
end
I=double(I);
[t,M,N]=size(I);

%%%%%%%%%%%%%%%%%%%将图像分为三类，分别根据图像选取初始分类中心%%%%%%%%%%%%%%%%%%

%等灰度值法确定初始分类中心
% for i=1:7
% m(i)=sum(sum(I(i,:,:)))/(256*256);
% p(i)=sqrt(sum(sum((I(i,:,:)-m(i)).^2))/(256*256-1));
% end
% N1=(m-p)';
% N2=(m)';
% N3=(m+p)';

% %读取相应点的坐标，确定中心
% for i=1:7
%     I(i,:,:)=imread(['L',num2str(i),'.bmp']);
%     figure
% imshow(['L',num2str(i),'.bmp']);
% [x,y]=ginput(3);
% % % 第一类水域
% N1(i)=I(i,x(1),y(1));
% % %第二类居民区