forwarddct.m

来自「实现图象,可以直接运行,里面有示例图象,十分便捷的实现图象压缩」· M 代码 · 共 26 行

function JQ=forwardDCT(Im,array)
Im=Im-128;
T = dctmtx(8);
y = blkproc(Im,[8 8],'P1*x*P2',T,T');
%y=blkproc(y,[8 8],'blk_thres_coding', 10);
JQ = blkproc(y,[8 8],'round(x ./ P1)',array);





% for i=1:8:256
%     for j=1:8:256
%         Jdct(i:i+7,j:j+7)= (T * Im(i:i+7,j:j+7)) * transpose(T);
%     end
% end
% 
% 
% for i=1:8:256
%     for j=1:8:256
%         JQ(i:i+7,j:j+7)= floor(Jdct(i:i+7,j:j+7)./array)+0.5;
%     end
% end