使用matlab实现的基于兴趣点的图像合并

function [I,h]=ImgMatch(I1,I2,n,aspect) %对图像 I1 和图像 I2 进行拼接 % aspect=0 为 y 方向（水平）连接； aspect=1 为 x 方向（垂直）连接 %输入检查 error(nargchk(2,4,nargin)); if nargin==2 n=8; aspect=0; elseif nargin==3 aspect=0; end if ndims(I1)~=ndims(I2) error('Imgage I1 must be the same as Image I2'); end if ndims(I1)>3 error('I1 and I2 must be image'); end if ndims(I1)~=3&&ndims(I1)~=2 error('I1 and I2 must be RGB or gray image'); end %转换图像成灰度模式 if ndims(I1)==3 f1=rgb2gray(I1); f2=rgb2gray(I2); else f1=I1; f2=I2; end %兴趣点检测并进行相似度检测 [pos,avb]=LACC(f1,f2,n); if avb==0 % 2 幅图像没有相似度 error('match error'); end %检查 pos，去除明显错误的点 if aspect==0 %垂直拼接 for i=size(pos,1):-1:1 if pos(i,2)<pos(i,4) pos(i,:)=[]; end end else %水平拼接 for i=size(pos,1):-1:1 if pos(i,1)<pos(i,3) pos(i,:)=[]; end end end x=pos(1:end,1:2)'; x(3,1:size(x,2))=1; y=pos(1:end,3:4)'; %{ %使用用最小中值法去除错误点 [h,RX,RY,avb]=LMeds(x,y,3,52); if avb==0 error('Image match fail'); end pos=[RX(1:2,:);RY(:,:)]; pos=pos'; %} %{ %进行图像角度变换 %求 4 个角点 if aspect==0 % y 向连接 corner=[1 1 size(I2,1) size(I2,1);pos(1,4) size(I2,2) pos(1,4) size(I2,2)]; else % x 向连接 corner=[pos(1,3) pos(1,3) size(I2,1) size(I2,1);1 size(I2,2) 1 size(I2,2)]; end [I2,yy,Corner]=Affine(I2,h,y,corner); %仿射变换 %黑颜色消除 if aspect==0 minx=max(Corner(1,1),Corner(1,2)); maxx=min(Corner(1,3),Corner(1,4)); if (maxx-minx+1)/size(I2,1)>0.9 I2=I2(minx:maxx,:,:); yy(1,:)=yy(1,:)-minx+1; end else miny=max(Corner(2,1),Corner(2,3)); maxy=min(Corner(2,2),Corner(2,4)); if (maxy-miny+1)/size(I2,2)>0.9 I2=I2(:,miny:maxy,:); yy(2,:)=yy(2,:)-miny+1; end end %将图像 I2 重新转换成灰度图像 if ndims(I2)==3 f2=rgb2gray(I2); else f2=I2; end %求拼接点 f1x=x(1,1); f1y=x(2,1); f2x=yy(1,1); f2y=yy(2,1); %} %{' f1x=pos(1,1); f1y=pos(1,2); f2x=pos(1,3); f2y=pos(1,4); %} [m1,n1]=size(f1); [m2,n2]=size(f2); %图像拼接 if aspect==0 %水平向连接 %图像对齐 if m1-f1x>m2-f2x I1=I1(1:end-((m1-f1x)-(m2-f2x)),:,:); elseif m1-f1x<m2-f2x I2=I2(1:end-((m2-f2x)-(m1-f1x)),:,:); end if f1x>f2x I1=I1((f1x-f2x+1):end,:,:); elseif f1x<f2x I2=I2((f2x-f1x+1):end,:,:); end %连接 %{ ' I(:,1:f1y-f2y,:)=I1(:,1:f1y-f2y,:); %重叠段，用重叠区线性过度的方法 ymax=n1; ymin=f1y-f2y+1; if ymax>ymin %不是刚刚重叠 for y=ymin:ymax jq=(ymax-y)/(ymax-ymin); %计算权值 I(:,y,:)=jq*I1(:,y,:)+(1-jq)*I2(:,y-f1y+f2y,:); end end I(:,(n1+1):(n2-f2y+f1y),:)=I2(:,(n1-f1y+f2y+1):n2,:); %} %{ %直接连接 I(:,1:f1y,:)=I1(:,1:f1y,:); I(:,f1y+1:n2-f2y+f1y,:)=I2(:,f2y+1:n2,:); %} end if aspect==1 %垂直连接 %图像对齐 if n1-f1y>n2-f2y I1=I1(:,1:end-((n1-f1y)-(n2-f2y)),:); elseif n1-f1y<n2-f2y I2=I2(:,1:end-((n2-f2y)-(n1-f1y)),:); end if f1y>f2y I1=I1(:,(f1y-f2y+1):end,:); elseif f1y<f2y I2=I2(:,(f2y-f1y+1):end,:); end %连接 I(1:f1x-f2x,:,:)=I1(1:f1x-f2x,:,:); %重叠段，用重叠区线性过度的方法 xmax=m1; xmin=f1x-f2x+1; if xmax>xmin %不是刚刚重叠 for x=xmin:xmax jq=(xmax-x)/(xmax-xmin); %计算权值 I(x,:,:)=jq*I1(x,:,:)+(1-jq)*I2(x-f1x+f2x,:,:); end end I((m1+1):(m2-f2x+f1x),:,:)=I2((m1-f1x+f2x+1):m2,:,:); %} %{ %直接连接 I(1:f1x,:,:)=I1(1:f1y,:,:); I(f1x+1:n2-f2x+f1x,:,:)=I2(f2x+1:n2,:,:); %} end