%兩幅影像的行列數必須一致才能運行成功
%否則請用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“