Rotate registration testing FFT Matlab code

T =[
1 0 1 2 3 4 5 6 6 5 4 3 2 1 0 1
0 1 2 3 4 5 6 7 7 6 5 4 3 2 1 0
1 2 3 4 5 6 7 8 8 7 6 5 4 3 2 1
2 3 4 5 6 7 8 9 9 8 7 6 5 4 3 2
3 4 5 6 7 8 9 10 10 9 8 7 6 5 4 3
4 5 6 7 8 9 10 9 9 10 9 8 7 6 5 4
5 6 7 8 9 10 9 8 8 9 10 9 8 7 6 5
6 7 8 9 10 9 8 7 7 8 9 10 9 8 7 6
6 7 8 9 10 9 8 7 7 8 9 10 9 8 7 6
5 6 7 8 9 10 9 8 8 9 10 9 8 7 6 5
4 5 6 7 8 9 10 9 9 10 9 8 7 6 5 4
3 4 5 6 7 8 9 10 10 9 8 7 6 5 4 3
2 3 4 5 6 7 8 9 9 8 7 6 5 4 3 2
1 2 3 4 5 6 7 8 8 7 6 5 4 3 2 1
0 1 2 3 4 5 6 7 7 6 5 4 3 2 1 0
1 0 1 2 3 4 5 6 6 5 4 3 2 1 0 1];

Im =[
1 0 1 2 3 4 5 6 6 5 4 3 2 1 0 1 2 3 4 5 0 0 0 0 0 0 0 0 0 0 0 0
0 1 2 3 4 5 6 7 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6 5 4 3 2
1 2 3 4 5 6 7 8 8 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6 5 4 3
2 3 4 5 6 7 8 9 9 8 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6 5 4
3 4 5 6 7 8 9 10 10 9 8 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6 5
4 5 6 7 8 9 10 9 9 10 9 8 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6
5 6 7 8 9 10 9 8 8 9 10 9 8 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7
6 7 8 9 10 9 8 7 7 8 9 10 9 8 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8
6 7 8 9 10 9 8 7 7 8 9 10 9 8 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8
5 6 7 8 9 10 9 8 8 9 10 9 8 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7
4 5 6 7 8 9 10 9 9 10 9 8 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6
3 4 5 6 7 8 9 10 10 9 8 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6 5
2 3 4 5 6 7 8 9 9 8 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6 5 4
1 2 3 4 5 6 7 8 8 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6 5 4 3
0 1 2 3 4 5 6 7 7 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6 5 4 3 2
1 0 1 2 3 4 5 6 6 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6 5 4 3 2 1
2 1 0 1 2 3 4 5 5 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6 5 4 3 2 1 0
3 2 1 0 1 2 3 4 4 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6 5 4 3 2 1 0 1
4 3 2 1 0 1 2 3 3 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6 5 4 3 2 1 0 1 2
5 4 3 2 1 0 1 2 2 1 0 1 2 3 4 5 6 7 8 9 8 7 6 5 4 3 2 1 0 1 2 0
6 5 4 3 2 1 0 1 1 0 1 2 3 4 5 6 7 8 9 8 7 6 5 4 3 2 1 0 1 2 0 0
7 6 5 4 3 2 1 0 0 1 2 3 4 5 6 7 8 9 8 7 6 5 4 3 2 1 0 1 2 0 0 0
8 7 6 5 4 3 2 1 1 2 3 4 5 6 7 8 9 8 7 6 5 4 3 2 1 0 1 2 0 0 0 0
9 8 7 6 5 4 3 2 2 3 4 5 6 7 8 9 8 7 6 5 4 3 2 1 0 1 2 0 0 0 0 0
8 9 8 7 6 5 4 3 3 4 5 6 7 8 9 8 7 6 5 4 3 2 1 0 1 2 0 0 0 0 0 0
7 8 9 8 7 6 5 4 4 5 6 7 8 9 8 7 6 5 4 3 2 1 0 1 2 0 0 0 0 0 0 0
6 7 8 9 8 7 6 5 5 6 7 8 9 8 7 6 5 4 3 2 1 0 1 2 0 0 0 0 0 0 0 0
5 6 7 8 9 8 7 6 6 7 8 9 8 7 6 5 4 3 2 1 0 1 2 0 0 0 0 0 0 0 0 0
4 5 6 7 8 9 8 7 7 8 9 8 7 6 5 4 3 2 1 0 1 2 0 0 0 0 0 0 0 0 0 0
3 4 5 6 7 8 9 8 8 9 8 7 6 5 4 3 2 1 0 1 2 0 0 0 0 0 0 0 0 0 0 0
2 3 4 5 6 7 8 9 9 8 7 6 5 4 3 2 1 0 1 2 0 0 0 0 0 0 0 0 0 0 0 0
1 2 3 4 5 6 7 8 8 7 6 5 4 3 2 1 0 1 2 0 0 0 0 0 0 0 0 0 0 0 0 0];
%Im=zeros(32,32);
%Im(1:16,1:16)=T;
TT=imrotate(T,45);
temp=zeros(32,32);
temp(4:28,4:28)=TT;
TT=temp;
IF=fftshift(fft2(Im));
TF=fftshift(fft2(TT));
pIF=imgpolarcoord(IF,15,360);
pTF=imgpolarcoord(TF,15,360);
pIF=abs(pIF);
pTF=abs(pTF);
fIF=fft2(pIF);
fTF=fft2(pTF);
tt=abs(fIF.*fTF);
zr=(tt==0);
GetMaxPos(zr)

Ra=fIF.*conj(fTF)./tt;
re=ifft2(Ra);

GetMaxPos(abs(re))

结果不对,我猜也许是笛卡尔坐标转极坐标时,命令不好. 应该是园形转换.