function img=imaging(us,ur)

% function img=imaging(us,ur)
%
% us,ur --> nz,nx,nt
% 
nz=size(us,1);
nx=size(us,2);
nf=size(us,3);

for i=1:nx;
    for j=1:nz;
        us(j,i,:)=fftshift(fft(us(j,i,:)));
        ur(j,i,:)=fftshift(fft(ur(j,i,1:nf)));
    end;
end;

img=zeros(nz,nx);
for i=1:nf;
    img=img+conj(us(:,:,i)).*ur(:,:,i);
end;
img=real(img);