基于Harris的图像拼接技术(MATLAB源代码).pdfVIP

基于Harris 角点的图像拼接技术 （MATLAB 源代码） clc clear all % 读入图片 pic 1=imread (lena1.jpg ); pic2=imread (lena2.jpg ); % Harris 角点检测 points 1=myHarris (pic 1); points2=myHarris (pic2); % 画出Harris 角点 figure (1) drawHarrisCorner (pic 1,points 1,pic2,points2); % 角点特征描述 des1=myHarrisCornerDescription (pic 1,points 1); des2=myHarrisCornerDescription (pic2,points2); % 角点粗匹配 matchs=myMatch (des1,des2); % 获取各自出匹配角点位置 matchedPoints1=points 1 (matchs (:,1),:); matchedPoints2=points2 (matchs (:,2),:); % 粗匹配角点连线 figure (2) drawLinedCorner (pic 1,matchedPoints1,pic2,matchedPoints2); % 角点精匹配 [newLoc1,newLoc2]=pointsSelect (matchedPoints1,matchedPoints2); % 精匹配角点连线 figure (3) drawLinedCorner (pic 1,newLoc1,pic2,newLoc2); % 图像拼接 im=picMatched (pic 1,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 (pic 1,points 1,