example6_1.m

来自「《MATLAB 6.5辅助图像处理》 本书是基于MATLAB 6.5的图像处」· M 代码 · 共 16 行

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)