1、自定义函数都在下面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,po
2、ints2); % 角点粗匹配 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
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 图% 输出：角点所在的行、纵的N2 矩阵 if length(size(pic)=3 pic=rgb2gray(pic); end pic=double(pic); hx=-1 0 1
4、; 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,
5、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,poi
6、nts2) % 功能：画出 Harris 角点的连接% 输入： % pic1、pic2 ：待拼接的图像 % points1、points2：Harris角点位置 X1=points1(:,2); Y1=points1(:,1); X2=points2(:,2); Y2=points2(:,1); dif=size(pic1,2); imshowpair(pic1,pic2,montage); hold on plot(X1,Y1,b*); plot(X2+dif,Y2,b*); set(gcf,Color,w); function des=myHarrisCornerDescription(p
7、ic,points) % 功能：Harris角点特征描述 % 输入： % pic：原图像 % points：角点位置 % 输出： % des：8N的角点特征描述矩阵 if length(size(pic)=3 pic=rgb2gray(pic); end len=length(points); des=zeros(8,len); for k=1:len p=points(k,:); pc=pic(p(1),p(2); des(1,k)=pic(p(1)-1,p(2)-1)-pc; des(2,k)=pic(p(1),p(2)-1)-pc; des(3,k)=pic(p(1)+1,p(2)-1)
8、-pc; des(4,k)=pic(p(1)+1,p(2)-pc; des(5,k)=pic(p(1)+1,p(2)+1)-pc; des(6,k)=pic(p(1),p(2)+1)-pc; des(7,k)=pic(p(1)-1,p(2)+1)-pc; des(8,k)=pic(p(1)-1,p(2)-pc; des(:,k)=des(:,k)/sum(des(:,k); end function matchs=myMatch(des1,des2)% 功能：特征点双向匹配 % 输入： % des1、des2：特征点描述信息构成的矩阵 % 输出： % matchs：匹配的特征点对应关系 len
9、1=length(des1); len2=length(des2); match1=zeros(len1,2); cor1=zeros(1,len2); for i=1:len1 d1=des1(:,i); for j=1:len2 d2=des2(:,j); cor1(j)=(d1*d2)/sqrt(d1*d1)*(d2*d2); end ,indx=max(cor1); match1(i,:)=i,indx; end match2=zeros(len2,2); cor2=zeros(1,len1); for i=1:len2 d2=des2(:,i); for j=1:len1 d1=de
10、s1(:,j); cor2(j)=(d1*d2)/sqrt(d1*d1)*(d2*d2); end ,indx=max(cor2); match2(i,:)=indx,i; end matchs=; for i=1:length(match1) for j=1:length(match2) if match1(i,:)=match2(j,:) matchs=matchs;match1(i,:); end end end function drawLinedCorner(pic1,loc1,pic2,loc2)% 功能：画出匹配角点的连接 % 输入： % pic1、pic2 ：待拼接的图像 %
11、loc1、loc2 ：匹配角点位置 X1=loc1(:,2); Y1=loc1(:,1); X2=loc2(:,2); Y2=loc2(:,1); dif=size(pic1,2); imshowpair(pic1,pic2,montage); hold on for k=1:length(X1) plot(X1(k),Y1(k),b*); plot(X2(k)+dif,Y2(k),b*); line(X1(k),X2(k)+dif,Y1(k),Y2(k),Color,r); end set(gcf,Color,w); function newLoc1,newLoc2=pointsSelect
12、(loc1,loc2) % 功能：筛选匹配特征点对，获取高精度的控制点 % 输入： % loc1、loc2 ：粗匹配特征点位置 % 输出： % newLoc1、 newLoc2：精匹配控制点位置 slope=(loc2(:,1)-loc1(:,1)./(loc2(:,2)-loc1(:,2); for k=1:3 slope=slope-mean(slope); len=length(slope); t=sort(abs(slope); thresh=t(round(0.5*len); ind=abs(slope)=1 end bottom=k; %下边界 k=1; while k0 if any(any(pic(:,k,:) break end k=k-1; end right=k; %右边界 %获取图像 im=pic(ceil:bottom,left:right,:);
