FFT的图像配准

一种基于FFT得快速图像匹配 matlab 带有图片 可运行

function [xtranytrandelta]=computedelta（im1im2）
% pre-align two images in frequecy domain
% author: bugzhao  bugzhao@sdu.edu.cn
% May 09th2005
% Please refer to “Mosaics of Scenes with Moving objectc“ for this method
% 粗估算两幅图像之间的平移差，返回水平方向平移差到xtran，返回垂直方向平移差到ytran
% 用在图像配准中的粗匹配中，对噪声、图像中物体发生运动等极为鲁棒！
if（size（im13）==3）
    im1=rgb2gray（im1）;
    im2=rgb2gray（im2）;
end
im1=double（im1）;
im2=double（im2）;
F1=fft2（im1）;
F2=fft2（im2）;
fz=F1.*conj（F2）;
fm=abs（F1.*F2）;
div=fz./fm;
delta=ifft2（div）;
maxvalue=max（max（delta））;
delta=fliplr（flipud（delta））;
imshow（delta[]）;
[ytranxtran]=find（delta==maxvalue）;
if（max（max（delta））==1） %完全相关，即表明两幅图像完全一样
    xtran=0;
    ytran=0;
else if（ytran==size（im11））
    ytran=0;
    else if （xtran==size（im12））
        xtran=0;
        end
    end
end
fprintf（‘图像2相对于图像1的位移量（单位-像素）:\n 水平方向（x方向） :  %d \n 竖直方向（y方向） :  %d \n‘xtranytran）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      22741  2005-05-09 03:38  FFT的图像配准\1g.jpg

     文件      20696  2005-05-09 05:04  FFT的图像配准\1g3.jpg

     文件      22082  2005-05-09 03:38  FFT的图像配准\2g.jpg

     文件       1023  2005-05-09 05:20  FFT的图像配准\computedelta.m

     文件        177  2005-05-09 05:33  FFT的图像配准\test.m

     文件       4862  2006-09-29 09:08  FFT的图像配准\中国图象图形网下载说明.html

     目录          0  2008-10-23 16:12  FFT的图像配准

----------- ---------  ---------- -----  ----

                71581                    7