function X=invertTwoDimDCT(Y,C)
%     C = dctMat(size(Y,1));
    X = C'*Y*C;