读图

mix imread( romems.jpg )

pan imread( rome_pan.jpg )

显示原多光谱图

subplot(221)

imshow(mix),title( 多光谱 )

subplot(222)

imshow(pan),title( 全色图像 )

预处理

mix double(mix)/255

pan double(pan)/255

求相关矩阵

[r c bands] size(mix)

pixels r*c

mix reshape(mix, [pixels bands])

correlation (mix *mix)/pixels

求特征向量与特征值-PCA conversion code