求取图像相关系数matlab

%两幅影像的行列数必须一致才能运行成功
%否则请用envi对影像重采样，再导出成tif格式即可
%如果是多光谱影像请将各波段单独进行运算
%读入两幅图像
imgA=imread（‘读取影像a.tif‘）;
imgB=imread（‘读取影像b.tif‘）;
%精度转换
imgA=double（imgA）;
imgB=double（imgB）;
%CC = XGXS（imgAimgB）;

%fprintf（‘%8.5f\n‘CC）;


%健壮性判断
[imgArowimgAcol]=size（imgA）;
[imgBrowimgBcol]=size（imgB）;
% [imgArowimgAcol];
% [imgBrowimgBcol];
if imgArow<1||imgAcol<1||imgBrow<1||imgBcol<1
    error（‘您的输入有误！维数不能小于1‘）;
elseif imgArow~=imgBrow||imgAcol~=imgBcol
    error（‘您输入德矩阵维数不相等！‘）;
end

%求两幅图像的均值差
imgA=imgA-mean2（imgA）;
imgB=imgB-mean2（imgB）;

%求求两幅图像的相关系数
CC=sum（sum（imgA.*imgB））./（sqrt（sum（sum（imgA.^2）））.*sqrt（sum（sum（imgB.^2））））;

fprintf（‘%8.5f\n‘CC）;
 

    
        tle>鏈皢瀵硅薄寮曠敤璁剧疆鍒板璞＄殑瀹炰緥銆?/title>
        ta name=“viewport“