phasekorel.m

图像对齐（inage registration）工具箱

function [km,dx,dy]=phasekorel(obra,obrb);


sz=max([size(obra,1),size(obra,2),size(obrb,1),size(obrb,2)]);
szz=ceil(sz/2)*2;

   fa=fft2(obra,szz,szz);
   fb=fft2(obrb,szz,szz);


ww=(fa).*(conj(fb));
www=ww./abs(ww);

km=abs(ifft2(www));

[val,dx,dy]=maximum(km);