example8_10.m

来自「matlab图像处理源代码 大家多多分享」· M 代码 · 共 21 行

%用原图像与压缩后的图像相减
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 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
  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)
M=I2-I;
figure, imshow(mat2gray(M))