目录
- 一、算法原理
- 1、原理概述
- 2、参考文献
- 二、代码实现
- 三、结果展示
本文由CSDN点云侠原创,原文链接。如果你不是在点云侠的博客中看到该文章,那么此处便是不要脸的爬虫与GPT。
一、算法原理
1、原理概述
相关系数的计算公式为:
2、参考文献
[1] 张宗营. SAR影像配准方法研究[D]. 山东科技大学, 2010.
二、代码实现
clc;
clear;
close all;%% -------------------------------读入两幅图像-----------------------------
imgA = imread('im1.png');
imgB = imread('im2.png');
%% --------------------------------精度转换--------------------------------
imgA = double(imgA);
imgB = double(imgB);
%% -------------------------------健壮性判断-------------------------------
[imgArow,imgAcol] = size(imgA);
[imgBrow,imgBcol] = size(imgB);if imgArow < 1||imgAcol < 1||imgBrow < 1||imgBcol < 1error('您的输入有误!维数不能小于1');
elseif imgArow ~= imgBrow||imgAcol ~= imgBcolerror('您输入德矩阵维数不相等!');
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);
%% -----------------------------结果可视化---------------------------------
figure;
subplot(121)
imagesc(imgA);
colorbar
title("图1");
hold on
subplot(122)
imagesc(imgB);
title('图2');
colorbar
hold on;