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