以某数据集为例:
P=load('PaviaU.mat');
[m,n,band]=size(P.paviaU);
ReshapedP=reshape(P.paviaU,m*n,band);
[coeff,score,latent,tsquared] = pca(ReshapedP);
percent=cumsum(latent)./sum(latent);
disp(percent);
feature_after_PCA=score(:,1:3);
RES=reshape(feature_after_PCA,m,n,3);
subplot(131);imshow(RES(:,:,1));
subplot(132);imshow(RES(:,:,2));
subplot(133);imshow(RES(:,:,3));