基于harris角点和RANSAC算法的图像拼接matlab仿真

1.算法运行效果图预览

2.算法运行软件版本

MATLAB2022a

3.算法理论概述

       Harris角点检测是一种局部特征检测方法，它寻找图像中具有显著局部曲率变化的位置，即边缘转折点或角点。主要通过计算图像窗口内的自相关矩阵M，并对其特征值进行评估。Harris响应函数H由自相关矩阵M的两个主特征值 λ1​ 和λ2​ 计算得到：

当H值较大时，窗口内像素的变化足够大，表明可能存在角点。

       在图像拼接过程中，RANSAC用于在两幅图像的重叠区域找到正确的对应点对，即使数据中存在大量噪声和异常点。

1.随机选择一组候选点对作为基础模型（通常是仿射或透视变换模型），计算此模型参数。

2.应用模型参数A预测所有剩余点对是否符合模型，统计一致样本数。

3.重复步骤1和2一定次数（迭代次数T），选取一致样本数最多的模型作为最优模型。

4.设定阈值（如变换残差阈值t），确定最终的内点集合（即那些变换误差小于阈值的所有点对）。

5.使用内点集合重新估计变换参数，以提高精度。

综上所述，图像拼接流程：

在每幅图像中检测Harris角点，并提取特征描述符。

利用匹配算法（如SIFT、SURF等）在重叠区域找到对应的角点对。

应用RANSAC算法找出最优变换模型。

根据最优变换模型对一幅图像进行几何校正，使两幅图像的重叠部分对齐。

最后，对齐后的图像通过融合算法（如加权平均、高斯金字塔融合等）拼接成全景图像。

4.部分核心程序

I1_harris   = func_harris(img1,3,para_harris,scales);% 使用自定义函数计算 img1 的角点响应
I2_harris   = func_harris(img2,3,para_harris,scales);% 使用自定义函数计算 img2 的角点响应
.........................................................................
% RANSAC 方法参数设定
% RANSAC 迭代次数
Miter = 500;
% 内点距离阈值
lvls  = 0.003;
% 最小内点数量要求
nums  = 10;
% RANSAC
[ransac_points,~]=func_ransac(harris_p1,harris_p2,mpoint,Miter,lvls,nums);
 
% 存储 RANSAC 匹配得到的正确匹配点坐标
Lens         =length(ransac_points);
X1_r=zeros(1,Lens);
Y1_r=zeros(1,Lens);
X2_r=zeros(1,Lens);
Y2_r=zeros(1,Lens);
 
for i=1:Lens
    p1=ransac_points(i,1);
    p2=ransac_points(i,2);
    X1_r(i)=x1(p1);
    Y1_r(i)=y1(p1);
    X2_r(i)=x2(p2);
    Y2_r(i)=y2(p2);
end
 
figure
subplot(1,2,1);
imshow(img1);
title('RANSAC匹配点')
hold on;
plot(Y1_r,X1_r,'ro');
 
subplot(1,2,2);
imshow(img2);
title('RANSAC匹配点')
hold on;
plot(Y2_r,X2_r,'go');
% 合成一幅图像展示 RANSAC 匹配结果
figure
img_match=[img1,img2];
imshow(img_match);
title('匹配结果')
hold on;
plot(Y1_r,X1_r,'rx','LineWidth',1,'MarkerSize',8);
plot(Y2_r+size(img1,2),X2_r,'gx','LineWidth',1,'MarkerSize',8);
 
Xm1_ransac=X1_r; 
Ym1_ransac=Y1_r;
match1=zeros(Lens,2);
match1(:,1)=Xm1_ransac; 
match1(:,2)=Ym1_ransac;
 
Xm2_ransac=X2_r;
Ym2_ransac=Y2_r+size(img1,2);
match2=zeros(Lens,2);
match2(:,1)=Xm2_ransac; 
match2(:,2)=Ym2_ransac;
% 绘制匹配线段连接匹配点
for i=1:Lens
    hold on;
    plot([match1(i,2) match2(i,2)], [match1(i,1) match2(i,1)],'LineWidth',1)
end
% 计算 RANSAC 方法得到的仿射变换矩阵
H_ransac = func_affine(X2_r,Y2_r,X1_r,Y1_r);
 
 
%生成一张新的全景图
[I1_ransac,I2_ransac]=func_trans(img1,img2,H_ransac);
 
 
 
figure
I_ransac= I1_ransac+I2_ransac ;
imshow(I_ransac);
title('RANSAC拼接结果')