%*****************************************?
%Copyright?（c）?2015?Jingshuang?Hu?
?
%@filename:demo.cpp?
%@datetime:2015.08.20?
%@author:HJS?
%@e-mail:eleftheria@163.com?
%@blog:http://blog.csdn.net/hujingshuang?
%*****************************************
%%
%相位相關法
clear?all;
close?all;
%%
img=rgb2gray（imread（‘lena.jpg‘））;%載入圖像
%%
%參數設置，可更改
dl=64;%模板大小64x64
x0=200;y0=200;%起點，最小（11），不超過模板尺寸-dl
dx0=20;dy0=30;%偏移（不超過模板尺寸）
%%
f1=img（x0:x0+dl-1y0:y0+dl-1）;%截取一個圖當做模板
f2=img（x0+dx0:x0+dx0+dl-1y0+dy0:y0+dy0+dl-1）;%截取一個圖當做粗匹配找到的子圖
[m?n]=size（f1）;%m行n列
F1=fft2（f1mn）;%傅里葉變換
F2=fft2（f2mn）;
%FC1=conj（F1）;%共軛
FC2=conj（F2）;%共軛
delta=abs（ifft2（（F1.*FC2）./abs（F1.*FC2）））;
figure;surf（1:dl1:dldelta）;shading?interp;
max_delat=max（max（delta））;
[dx?dy]=find（delta==max_delat）;
dx=dx-1;dy=dy-1;
dx
dy

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????229598??2014-03-28?15:42??相位相關法\lena.JPG

?????文件????????907??2015-08-20?20:04??相位相關法\PC.m

?????目錄??????????0??2015-08-20?20:00??相位相關法

-----------?---------??----------?-----??----

???????????????230505????????????????????3