图形图象

开发平台：

Matlab

phasec.m：源码内容

							function [row,col]=phasec(IL,IR)
Z=0;
A=fft2(double(IR));
B=fft2(double(IL));
AA=conj(A);
C1=B.*AA;
C2=norm(C1,1);
%取模
C3=C1/C2;
C=ifft2(C3);
%C＝exp(-j*2pi*(*x0+*y0))
[m,n]=size(C);
for i=1:m
    for j=1:n
    if C(i,j)>Z
        Z=C(i,j);
        row=i;
        col=j;
    end
end
end

						