2024年5月3日发(作者：)

基于Harris角点的图像拼接技术（MATLAB源代码）

clc

clear all

% 读入图片

pic1=imread('');

pic2=imread('');

% Harris角点检测

points1=myHarris(pic1);

points2=myHarris(pic2);

% 画出Harris角点

figure(1)

drawHarrisCorner(pic1,points1,pic2,points2);

% 角点特征描述

des1=myHarrisCornerDescription(pic1,points1);

des2=myHarrisCornerDescription(pic2,points2);

% 角点粗匹配

matchs=myMatch(des1,des2);

% 获取各自出匹配角点位置

matchedPoints1=points1(matchs(:,1),:);

matchedPoints2=points2(matchs(:,2),:);

% 粗匹配角点连线

figure(2)

drawLinedCorner(pic1,matchedPoints1,pic2,matchedPoints2);

% 角点精匹配

[newLoc1,newLoc2]=pointsSelect(matchedPoints1,matchedPoints2);

% 精匹配角点连线

figure(3)

drawLinedCorner(pic1,newLoc1,pic2,newLoc2);

% 图像拼接

im=picMatched(pic1,newLoc1,pic2,newLoc2);

% 显示拼接图像

figure(4)

imshow(im);

set(gcf,'Color','w');

function points=myHarris(pic)

% 功能：寻找Harris角点

% 输入：RGB图像或gray图像

% 输出：角点所在的行、纵的N×2矩阵

if length(size(pic))==3

pic=rgb2gray(pic);

end

pic=double(pic);

hx=[-1 0 1];

Ix=filter2(hx,pic);

hy=[-1;0;1];

Iy=filter2(hy,pic);

Ix2=Ix.*Ix;

Iy2=Iy.*Iy;

Ixy=Ix.*Iy;

h=fspecial('gaussian',[7 7],2);

Ix2=filter2(h,Ix2);

Iy2=filter2(h,Iy2);

Ixy=filter2(h,Ixy);

[heigth,width]=size(pic);

alpha=0.06;

R=zeros(heigth,width);

for i=1:heigth

for j=1:width

M=[Ix2(i,j) Ixy(i,j);Ixy(i,j) Iy2(i,j)];

R(i,j)=det(M)-alpha*(trace(M)^2);

end

end

Rmax=max(max(R));

pMap=zeros(heigth,width);

for i=2:heigth-1

for j=2:width-1

if R(i,j)>0.01*Rmax

tm=R(i-1:i+1,j-1:j+1);

tm(2,2)=0;

if R(i,j)>tm

pMap(i,j)=1;

end

end

end

end

[row,col]=find(pMap==1);

points=[row,col];

function drawHarrisCorner(pic1,points1,pic2,points2)

% 功能：画出Harris角点的连接