untitled3.m

来自「简单的DCT的MATLAB图象压缩实验」· M 代码 · 共 15 行

I=imread('cameraman.tif');
I=double(I)/255;
T=dctmtx(8);
B=blkproc(I,[8 8],'P1*x*P2',T,T);
mask=[1 1 1 1 1 0 0 0
      1 1 1 1 0 0 0 0
      1 1 1 0 0 0 0 0
      1 1 0 0 0 0 0 0
      1 0 0 0 0 0 0 0
      0 0 0 0 0 0 0 0
      0 0 0 0 0 0 0 0
      0 0 0 0 0 0 0 0];
  B2=blkproc(B,[8 8],'P1.*x',mask);
  I2=blkproc(B2,[8 8],'P1*x*P2',T,T);
  imshow(I),figure,imshow(I2)