Matlab使用PCA描绘多光谱图像

Matlab使用PCA描绘多光谱图像

```Matlabclear close allb1=imread('Fig1138(a)(WashingtonDC_Band1_564).tif');b2=imread('Fig1138(b)(WashingtonDC_Band2_564).tif');b3=imread('Fig1138(c)(WashingtonDC_Band3_564).tif');b4=imread('Fig1138(d)(WashingtonDC_Band4_564).tif');b5=imread('Fig1138(e)(WashingtonDC_Band5_564).tif');b6=imread('Fig1138(f)(WashingtonDC_Band6_564).tif');s= cat(3, b1,b2,b3,b4, b5, b6);[X, R]= imstack2vectors(s);P= princomp(X,6);g1= P.Y(:,1);g1= reshape(g1, 564, 564);figure,imshow(b1);figure,imshow(g1, []);d=diag(P.Cy);g2= P.Y(:,2);g2= reshape(g2, 564, 564);figure,imshow(b2);figure,imshow(g2, []);g3= P.Y(:,3);g3= reshape(g3, 564, 564);figure,imshow(b3);figure,imshow(g3, []);g4= P.Y(:,4);g4= reshape(g4, 564, 564);figure,imshow(b4);figure,imshow(g4, []);g5= P.Y(:,5);g5= reshape(g5, 564, 564);figure,imshow(b5);figure,imshow(g5, []);g6= P.Y(:,6);g6= reshape(g6, 564, 564);figure,imshow(b6);figure,imshow(g6, []);<div class="se-preview-section-delimiter"></div>




cat函数
Concatenation for heterogeneous arrays
不同数组/阵列的连接
Syntax
C = cat(dim,A,B)

## princomp函数 &covmatrix函数##
由于我使用的是冈萨雷斯的书籍，发现princomp函数与Matlab自身提供的函数不同，所以请将原来书本的代码复制到新的路径下，run即可。

由于原图是tif格式，本网站暂不支持tif格式的图像上传，故不在此发图。