function [r]=POC(f,g)
    F=fft2(f);
    G=fft2(g);
    conjG=conj(G);
    FconjG=F.*conjG;
    R=FconjG./(abs(FconjG));
    r=real(ifft2(R));
end