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

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