my_dct.m

来自「清华大学胡广书教授< 数字信号处理理论算法实现>>第二版一书的」· M 代码 · 共 25 行

function C=my_DCT(A,k)
s=size(A);
for m=1:s(1)
    if rem(m,8)==0
        for n=1:s(2)
            if rem(n,8)==0
                B((m-8+1):m,(n-8+1):n)=dct2(A((m-8+1):m,(n-8+1):n));
                B((m-k+1):m,(n-8+1):n)=0;
                B((m-8+1):m,(n-k+1):n)=0;
            end
        end
    end
end
%
for m=1:128
    if rem(m,8)==0
        for n=1:128
            if rem(n,8)==0
                C((m-8+1):m,(n-8+1):n)=idct2(B((m-8+1):m,(n-8+1):n));
            end
        end
    end
end
% 
C=uint8(C);