%???首先利用PCA對(duì)其進(jìn)行主成分變化???然后利用K-means方法對(duì)其進(jìn)行分類
clc;
clear?all;
close?all;
%讀取影像
I=imread（‘12.tif‘）;
subplot（121）;
imshow（I）;
%生成矩陣
I=double（I）;
[n1n2m]=size（I）;
%至0
P=zeros（n1n2m）;
%重新分配灰度值（三個(gè)波段等間隔分配0~255）
for?i=1:m
????P（::i）=（（I（::i）-min（min（I（::i））））/（max（max（I（::i）））-min（min（I（::i）））））.*255;
end
%計(jì)算像素點(diǎn)總個(gè)數(shù)
n=n1*n2;
%將Q的值重新分配成3波段
p=reshape（Pnm）;
%進(jìn)行主成分分析得到新的主元
[COEFFSCORElatent]=princomp（p）;
pcaData1=SCORE（:1）;
%設(shè)定分的類別個(gè)數(shù)
k=5;
%進(jìn)行k-means聚類分析
[idxbestCbestsumDbestDbest]=kmeans（pcaData1k）;
n=length（idxbest）;
%的到新的聚類后在進(jìn)行重排列
o=reshape（idxbestn1n2）;
%設(shè)置每一類的元素
cm=[255?0?0;128?128?128;0?255?0;0?0?255;0?64?64];
cpm=reshape（cm（o:）n1n23）;
subplot（122）;
imshow（cpm）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-06-26?22:45??非監(jiān)督分類\
?????文件??????234649??2018-06-26?14:05??非監(jiān)督分類\12.tif
?????文件?????????848??2018-06-26?21:11??非監(jiān)督分類\multiSpectral.m