function [xtran,ytran,delta]=computedelta(im1,im2)
% pre-align two images in frequecy domain
% author: bugzhao  bugzhao@sdu.edu.cn
% May 09th,2005
% Please refer to "Mosaics of Scenes with Moving Objectc" for this method
% 粗估算兩幅圖像之間的平移差，返回水平方向平移差到xtran，返回垂直方向平移差到ytran
% 用在圖像配準中的粗匹配中，對噪聲、圖像中物體發生運動等極為魯棒！
if(size(im1,3)==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,[]);
[ytran,xtran]=find(delta==maxvalue);
if(max(max(delta))==1) %完全相關，即表明兩幅圖像完全一樣
    xtran=0;
    ytran=0;
else if(ytran==size(im1,1))
    ytran=0;
    else if (xtran==size(im1,2))
        xtran=0;
        end
    end
end
fprintf('圖像2相對于圖像1的位移量(單位-像素):\n 水平方向(x方向) :  %d \n 豎直方向(y方向) :  %d \n',xtran,ytran);
end