求两幅多波段影像的相关系数，并且能实现影像的多波段读取.

clear clc %读取一幅多波段影像 im1 = multibandread('C:\multi1', [126 112 4], ... 'uint8', 0, 'bip', 'ieee-le'); %读取融合后的影像 im2 = multibandread('C:\fusion', [126 112 4], ... 'uint8', 0, 'bip', 'ieee-le') ; %提前生成一个4*4的o阵用于存放相关系数 A=zeros(4,4); %用一循环求取两两波段之间的相关系数 for i=1:4 for j=1:4 r=corr2(im1(:,:,i),im2(:,:,j)); A(i,j)=r; end end %输出相关系数矩阵A disp(A) %多波段影像的其中一个波段与融合