www.gusucode.com > peizhun配准matlab程序源码 > peizhun配准matlab程序源码/GLPF.m
function [I,J] = GLPF(handles) I = handles.I; J = handles.J; I = double(I); J = double(J); PQ = 2*size(I); M = PQ(1); N = PQ(2); u = 0:(M-1); v = 0:(N-1); idx = find(u>M/2); u(idx) = u(idx)-M; idy = find(v>N/2); v(idy) = v(idy)-N; [V,U] = meshgrid(v,u); D0 = 0.05*PQ(2); H = exp(-(U.^2+V.^2)/(2*(D0^2))); FI = fft2(I,size(H,1),size(H,2)); FJ = fft2(J,size(H,1),size(H,2)); gi = real(ifft2(H.*FI)); gj = real(ifft2(H.*FJ)); gi = gi(1:size(I,1),1:size(I,2)); gj = gj(1:size(I,1),1:size(I,2)); I = unit8(gi); J = unit8(gj);