function [U V]=dftuv(M,N)
%DFTUV Computes meshgrid frequency matrices
%
u=0:(M-1);
v=0:(N-1);
% Compute indices
idx=find(u>M/2);
u(idx)=M-u(idx);   
idy=find(v>N/2);
v(idy)=N-v(idy);
[V U]=meshgrid(v,u);