(年整理)图像拼接基于harris检测matlab

图像拼接基于harris检测matlab代码

clc
clear all
% 读入图片
pic1=imread(’lena1.jpg’)；
pic2=imread（‘lena2。jpg’）;
% 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&