读图 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) p