%Kurt's 2D discrete cosine transform; uses fft call. Uses my own kdct
%command.
function Y = kdct2d(y)

Y0 = kdct(y)';
Y = kdct(Y0)';