function dcty = dct2D(img)

dcty = zeros(8,8);
tem  = 0;
for u=1:8
    for v=1:8
        for x=1:8
            for y=1:8
                tem=tem + (img(x,y)*cos((2*x-1)*(u-1)*pi/16)*cos((2*y-1)*(v-1)*pi/16));
            end;
        end;
        
        if (u==1 & v~=1)
            dcty(u,v) = tem / (4*sqrt(2));
        
        elseif (v==1 & u~=1)
            dcty(u,v) = tem / (4*sqrt(2));
        
        elseif (u==1 & v==1)
            dcty(u,v) = tem / 8;
        
        else
            dcty(u,v) = tem / 4;
        end;
            tem = 0;
    end;
end;