MATLAB实现Harris角点检测与图像配准

function im=picMatched(pic1,newLoc1,pic2,newLoc2) % 功能：获取拼接之后的图片 % 输入： % pic1、pic2：待拼接图片 % newLoc1、newLoc2：变换控制点矩阵 % 输出： % im：拼接后的图片 if length(size(pic1))==2 pic1=cat(3,pic1,pic1,pic1); end if length(size(pic2))==2 pic2=cat(3,pic2,pic2,pic2); end SZ=2000; X1=newLoc1(:,2); Y1=newLoc1(:,1); X2=newLoc2(:,2); Y2=newLoc2(:,1); sel=randperm(length(newLoc1),3); x=X2(sel)'; y=Y2(sel)'; X=X1(sel)'; Y=Y1(sel)'; U=[x;y;ones(1,3)]; V=[X;Y;ones(1,3)]; T=V/U; cntrX=SZ/2; cntrY=SZ/2; im=zeros(SZ,SZ,3); for i=1:size(pic2,1) for j=1:size(pic2,2) tmp=T*[j;i;1]; nx=round(tmp(1))+cntrX; ny=round(tmp(2))+cntrY; if nx>=1 &&nx<=SZ &&ny>=1 &&ny<=SZ im(ny,nx,:)=pic2(i,j,:); end end end im=imresize(im,1,'bicubic'); tpic1=zeros(SZ,SZ,3); tpic1(1+cntrY:size(pic1,1)+cntrY,1+cntrX:size(pic1,2)+cntrX,:)=pic1; re=rgb2gray(uint8(im))-rgb2gray(uint8(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); im=uint8(im); if length(size(pic1))==2 im=rgb2gray(im); end