function dcty = dct2D(img)

dcty = zeros(8,8);
tem  = 1;
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 = 1;
    end;
end;