DCT.zip

DCT.zip

%DCT
x=imread（‘frame11.png‘）; 
x=rgb2gray（x）; 
x=double（x）;%变成双精度
t=dctmtx（8）;%生成DCT变换矩阵
%亮度量化表，乘以一个系数，用来控制压缩率。系数越大，压缩率越高，得到的图像质量越差。
a1=1*[1611101624405161;
   1212141926586055;
   1413162440576956;
   1417222951878062;
   182237566810910377;
   243555648110411392;
   49647887103121120101;
   7292959811210010399];
y = blkproc（x [8 8] ‘P1 * x * P2‘ t t‘）; 
%量化
y10= blkproc（y [8 8] ‘round（x./P1）‘a1）; 
figure（1）;imshow（uint8（y10）[]）;title（‘DCT系数‘）; 
% figure（1）;imshow（log（abs（y10））[]）;title（‘DCT系数‘）; 
%掩膜矩阵用来压缩DCT系数，保留的AC系数越多，得到的图像质量越高。
mask1=[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 
0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0];%保留了左上3个 
mask2=[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];%保留了左上10 
mask3=[1 1 1 1 1 1 1 1 
       1 1 1 1 1 1 1 0 
       1 1 1 1 1 1 0 0 
       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];%保留了左上36 
mask4=[ 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1];%保留了全部 
%% 反量化与加掩膜
y1=blkproc（y10[8 8]‘P1.*P2.*x‘mask1a1）; 
y2=blkproc（y10[8 8]‘P1.*P2.*x‘mask2a1）; 
y3=blkproc（y10[8 8]‘P1.*P2.*x‘mask3a1）; 
y4=blkproc（y10[8 8]‘P1.*P2.*x‘mask4a1）; 

y11=uint8（blkproc（y1[8 8]‘P1*x*P2‘t‘t））; 
y22=uint8（blkproc（y2[8 8]‘P1*x*P2‘t‘t））; 
y33=uint8（blkproc（y3[8 8]‘P1*x*P2‘t‘t））; 
y44=uint8（blkproc（y4[8 8]‘P1*x*P2‘t‘t））; 

figure（2）;imshow（uint8（x））title（‘原图‘）; 
figure（3）;imshow（y11）title（‘还原保留左上3‘）; 
figure（4）;imshow（y22）title（‘还原保留左上10‘）; 
figure（5）;imshow（y33）title（‘还原保留左上36‘）; 
figure（6）;imshow（y44）title（‘还原全部‘）; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-04-22 16:39  DCT\
     文件        2055  2019-04-22 16:39  DCT\DCT.m
     文件      402536  2007-07-28 03:16  DCT\frame11.png