《图像拼接matlab12055.pdf》由会员分享，可在线阅读，更多相关《图像拼接matlab12055.pdf（14页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、-.z.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);%角点粗匹配 mat
2、chs=myMatch(des1,des2);%获取各自出匹配角点位置 matchedPoints1=points1(matchs(:,1),:);matchedPoints2=points2(matchs(:,2),:);%粗匹配角点连线 figure(2)drawLinedCorner(pic1,matchedPoints1,pic2,matchedPoints2);-.z.%角点精匹配 newLoc1,newLoc2=pointsSelect(matchedPoints1,matchedPoints2);%精匹配角点连线 figure(3)drawLinedCorner(pic1,new
3、Loc1,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);h*=-1 0 1;I*=filter2(h*,pic);hy=-1;0;1;-.z.Iy=filt
4、er2(hy,pic);I*2=I*.*I*;Iy2=Iy.*Iy;I*y=I*.*Iy;h=fspecial(gaussian,7 7,2);I*2=filter2(h,I*2);Iy2=filter2(h,Iy2);I*y=filter2(h,I*y);heigth,width=size(pic);alpha=0.06;R=zeros(heigth,width);for i=1:heigth for j=1:width M=I*2(i,j)I*y(i,j);I*y(i,j)Iy2(i,j);R(i,j)=det(M)-alpha*(trace(M)2);end end Rma*=ma*(m
5、a*(R);pMap=zeros(heigth,width);for i=2:heigth-1 for j=2:width-1 if R(i,j)0.01*Rma*-.z.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 角点的连接%输入：%pic1、pic2：待拼接的图像%points1、point
6、s2：Harris角点位置*1=points1(:,2);Y1=points1(:,1);*2=points2(:,2);Y2=points2(:,1);dif=size(pic1,2);imshowpair(pic1,pic2,montage);hold on -.z.plot(*1,Y1,b*);plot(*2+dif,Y2,b*);set(gcf,Color,w);function des=myHarrisCornerDescription(pic,points)%功能：Harris角点特征描述%输入：%pic：原图像%points：角点位置%输出：%des：8N的角点特征描述矩阵 if
7、 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)-pc;des(4,k)=pic(p(1)+1,p(2)-pc;-.z.des(5,k)=pic(p(1)+1,p(2)+1)-pc;des(6,k)=pic(p(1),p(2)+1)
8、-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：匹配的特征点对应关系 len1=length(des1);len2=length(des2);match1=zeros(len1,2);cor1=zeros(1,len2);for i=1:len1 d1=des1(:,i);for j
9、=1:len2 d2=des2(:,j);cor1(j)=(d1*d2)/sqrt(d1*d1)*(d2*d2);end -.z.,ind*=ma*(cor1);match1(i,:)=i,ind*;end match2=zeros(len2,2);cor2=zeros(1,len1);for i=1:len2 d2=des2(:,i);for j=1:len1 d1=des1(:,j);cor2(j)=(d1*d2)/sqrt(d1*d1)*(d2*d2);end ,ind*=ma*(cor2);match2(i,:)=ind*,i;end matchs=;for i=1:length(ma
10、tch1)for j=1:length(match2)if match1(i,:)=match2(j,:)matchs=matchs;match1(i,:);end end end -.z.function drawLinedCorner(pic1,loc1,pic2,loc2)%功能：画出匹配角点的连接%输入：%pic1、pic2：待拼接的图像%loc1、loc2：匹配角点位置*1=loc1(:,2);Y1=loc1(:,1);*2=loc2(:,2);Y2=loc2(:,1);dif=size(pic1,2);imshowpair(pic1,pic2,montage);hold on fo
11、r k=1:length(*1)plot(*1(k),Y1(k),b*);plot(*2(k)+dif,Y2(k),b*);line(*1(k),*2(k)+dif,Y1(k),Y2(k),Color,r);end set(gcf,Color,w);function newLoc1,newLoc2=pointsSelect(loc1,loc2)%功能：筛选匹配特征点对，获取高精度的控制点%输入：%loc1、loc2：粗匹配特征点位置-.z.%输出：%newLoc1、newLoc2：精匹配控制点位置 slope=(loc2(:,1)-loc1(:,1)./(loc2(:,2)-loc1(:,2)
12、;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&n*=1&ny=SZ im(ny,n*,:)=pic2(i,j,:);end end end im=imresize(im,1,bicubic);tpic1=zeros(SZ,SZ,3);tpic1(1+trY:size(pic1,1)+trY,1+tr*:size(pic1,2)+tr*,:)=pic1;re=rgb2gray(uint8(im)-rgb2gray(uin
13、t8(tpic1);for k=1:3 ta=im(:,:,k);tb=tpic1(:,:,k);ta(re=0)=tb(re=0);im(:,:,k)=ta;end clear ta tb re tpic1 im=getPicture(im,SZ);-.z.im=uint8(im);if length(size(pic1)=2 im=rgb2gray(im);end function im=getPicture(pic,SZ)%功能：获取图像有用区域%输入：%pic：拼接图像%SZ：预定图像尺寸%输出：%im：有用区域图像 if length(size(pic)=2 pic=cat(3,pic,pic,pic);end k=1;while k0 if any(any(pic(k,:,:)break end k=k-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,:);