%?K-mean方法聚類(lèi)
clear?all;

%%%%%%%%%%%%%%%%%%%讀取7幅遙感圖像數(shù)據(jù)%%%%%%%%%%%%%%%%%%
for?i=1:7
????I（i::）=imread（[‘L‘num2str（i）‘.bmp‘]）;
end
I=double（I）;
[tMN]=size（I）;

%%%%%%%%%%%%%%%%%%%將圖像分為三類(lèi)，分別根據(jù)圖像選取初始分類(lèi)中心:三種方法%%%%%%%%%%%%%%%%%%

%等灰度值法確定初始分類(lè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）‘;

%?%讀取相應(yīng)點(diǎn)的坐標(biāo)，確定中心
%?for?i=1:7
%?????I（i::）=imread（[‘L‘num2str（i）‘.bmp‘]）;
%?????figure
%?imshow（[‘L‘num2str（i）‘.bmp‘]）;
%?[xy]=ginput（3）;
%?%?%?第一類(lèi)水域
%?N1（i）=I（ix（1）y（1））;
%?%?%第二類(lèi)居民區(qū)
%?N2（i）=I（ix（2）y（2））;
%?%?%第三類(lèi)
%?N3（i）=I（ix（3）y（3））;
%?close（figure）;
%?end
%?N1=N1‘;
%?N2=N2‘;
%?N3=N3‘;

%直接讀取圖像灰度值確定中心
N1=[42.3?63.06?61.7?21.6?15.6?48.1?9.6]‘;
N2=[10

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件??????66616??2007-03-05?08:43??遙感圖像分析\L7.bmp

?????文件???????2478??2011-01-20?16:04??遙感圖像分析\Kmeans.m

?????文件??????66616??2007-03-05?08:43??遙感圖像分析\L1.bmp

?????文件??????66616??2007-03-05?08:43??遙感圖像分析\L2.bmp

?????文件??????66616??2007-03-05?08:43??遙感圖像分析\L3.bmp

?????文件??????66616??2007-03-05?08:43??遙感圖像分析\L4.bmp

?????文件??????66616??2007-03-05?08:43??遙感圖像分析\L5.bmp

?????文件??????66616??2007-03-05?08:43??遙感圖像分析\L6.bmp

?????文件??????49664??2011-01-20?16:01??遙感圖像分析\遙感圖像分割.doc

?????目錄??????????0??2011-01-20?16:02??遙感圖像分析

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

???????????????518454????????????????????10