function im=imRegist(im1,im2,T)
sz=3*max(length(im1),length(im2));
dim=length(size(im1));
if dim==3
[row1,col1,~]=size(im1);
[row2,col2,~]=size(im2);
im=zeros(sz,sz,3);
else
[row1,col1]=size(im1);
[row2,col2]=size(im2);
im=zeros(sz,sz);
end
cX=sz/3;
cY=sz/3;
if dim==3
im(1+cY:row1+cY,1+cX:col1+cX,:)=im1;
else
im(1+cY:row1+cY,1+cX:col1+cX)=im1;
end
T=T^(-1);
for i=1:size(im,1)
for j=1:size(im,2)
xy1=[j-cX;i-cY;1];
xy2=round(T*xy1);
nx=xy2(1);
ny=xy2(2);
if nx>=1&& nx<=col2 && ny>=1 && ny<=row2
if i<=cY || i>=cY+row1 || j<=cX ||j>=cX+col1
if dim==3
im(i,j,:)=im2(ny,nx,:);
else
im(i,j)=im2(ny,nx);
end
end
end
end
end
im=imCrop(im);
im=uint8(im);
function im=imCrop(pic)
if length(size(pic))==3
gray=rgb2gray(pic);
else
gray=pic;
end
SZ=length(gray);
k=1;
while k<SZ
if any(any(gray(k,:)))
break
end
k=k+1;
end
ceil=k;
k=SZ;
while k>0
if any(any(gray(k,:)))
break
end
k=k-1;
end
bottom=k;
k=1;
while k<SZ
if any(any(gray(:,k)))
break
end
k=k+1;
end
left=k;
k=SZ;
while k>0
if any(any(gray(:,k)))
break
end
k=k-1;
end
right=k;
if length(size(pic))==3
im=pic(ceil:bottom,left:right,:);
else
im=pic(ceil:bottom,left:right);
end
Matlab实现SIFT特征检测+图像拼接
1星 需积分: 48 92 浏览量
2021-10-23
20:27:39
上传
评论 26
收藏 3.7MB ZIP 举报 
sift特征检测+图像拼接.zip （19个子文件） 
main.m 654B sift.m 1KB
test3.png 782KB appendingImages.m 453B getTransMat.m 1KB test4.png 863KB optIndex.m 326B tmp.pgm 580KB imRegist.m 1KB drawRotAglHist.m 193B getRotAgl.m 100B
siftWin32.exe 92KB img2gray.m 136B featureMatch.m 872B tmp.key 1.26MB drawPoints.m 195B test1.png 638KB linePoints.m 451B test2.png 671KB
评论2