DCT.zip_DCT变换编码_dct图像压缩

clear all; % close I0=imread('3.jpg'); I=rgb2gray(I0); subplot(2,2,1),imshow(I); xlabel('a) 原始图像'); disp('原始图像大小：') whos('I') I=double(I); %图像类型转换，将图像矩阵转换成双精度类型 T=dctmtx(8); %DCT矩阵 mask=[1 0 0 0 0 0 0 0; %0AC系数遮罩矩阵 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 0; 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 0]; % mask=[1 1 1 1 0 0 0 0; %9AC系数遮罩矩阵 % 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]; % m=1*[1 11 10 16 24 40 51 61; %量化矩阵 % 12 12 14 19 26 58 60 55; % 14 13 16 24 40 57 69 56; % 14 17 22 29 51 87 80 62; % 18 22 37 56 68 109 103 77; % 24 35 55 64 81 104 113 92; % 49 64 78 87 103 121 120 101; % 72 92 95 98 112 100 103 99]; B=blkproc(I,[8 8],'P1*x*P2',T,T'); %DCT变换 % LB=blkproc(B,[8 8],'round(x./P1)',m); %量化 MB=blkproc(B,[8 8],'x.*P1',mask); %对DCT系数进行遮罩处理 I2=blkproc(MB,[8 8],'P1*x*P2',T',T); %反DCT变换IDCT subplot(2,2,2),imshow(uint8(I-I2)); xlabel('b) 保留0个交流系数'); disp('压缩后的图像大小：') whos('I2') mean(mean(abs(I-I2))) mask=[1 1 0 0 0 0 0 0; %3AC系数遮罩矩阵 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; 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 0]; MB=blkproc(B,[8 8],'x.*P1',mask); %对DCT系数进行遮罩处理 I2=blkproc(MB,[8 8],'P1*x*P2',T',T); %反DCT变换IDCT subplot(2,2,3),imshow(uint8(I-I2)); xlabel('c) 保留2个交流系数'); mean(mean(abs(I-I2))) mask=[1 1 1 1 1 0 0 0; %5AC系数遮罩矩阵 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]; MB=blkproc(B,[8 8],'x.*P1',mask); %对DCT系数进行遮罩处理 I2=blkproc(MB,[8 8],'P1*x*P2',T',T); %反DCT变换IDCT subplot(2,2,4),imshow(uint8(I-I2)); xlabel('d) 保留5个交流系数'); mean(mean(abs(I-I2)))