17-3.m

matlab图像处理源码

I = imread('cameraman.tif');
I = im2double(I);
%图像存储类型转换
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)
title('原始图像')
figure;
imshow(I2)
title('压缩后的图像')