数字图像处理实验八图像变换压缩编码

《数字图像处理实验八图像变换压缩编码》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、实验八图像变换压缩编码实验目的：1，掌握变换压缩编码的原理与方法2，了解有损压缩图像质量评价的指标二，实验条件1，MATLAB软件2,典型的灰度，彩色图像。三，实验原理1，去除数据冗余度可以有效地压缩数据2,图像编码压缩的主耍技术指标：压缩比，客观评价SNR,主观评价四，实验内容1,以DCT变换为例，采用标准亮度量化表，得有损压缩的变换值，计算并显示其直方图2,采用标准量化表将上述变换值反量化回去后用IDCT变换回图像，比较并计算压缩前后图像的差异(MSE,PSNR等)3,根据由1计算得到的直方图试判断采用Huffman编码的平均码长。五，实验步骤1，以DCT变换为例，采用

2、标准亮度量化表，得有损压缩的变换值，计算并显示其直方图(1)代码如下所示I=imread(,a.jpg,);J=rgb2gray(I);J=imresize(J,[512512]);%imshow(J);A=zeros(8,8);B=zeros(512,512);lighttable=...%标准亮度量化表[1611101624405161,1212141926586055,1413162440576956,1417222951878062,182237566810910377,243555648110411392，49647887103121120101，72929598

3、11210010399];fori=l:64%用i,j来控制子块forj=1:64form=(8*(i-l)+l):8*iforn=(8*(j-l)+l):8*jA(m-(8*(i-l)+l)+l，n-(8*(j-l)+l)+l)=J(m，n);%提取8*8子块endendCl=dct2(A);%对子块进行dct2变换C2=round(Cl./lighttable);%将Cl进行量化，并取整forp=l:8%将子块按序存进512*512的矩阵中forq=l:8B(p+8*(i-l)，q+8*(H))=C2(p，q);endendendendH=zeros(1,256*2-1

4、);fori=l:mforj=l:nH(B(i，j)+255)=H(B(i，j)+255)+l;%统计各值概率endendX=[-255:255];plot(X,H);%绘制坐标为-255—255的直方图(2)直方图如下所示x101002003002,采用标准量化表将上述变换值反量化回去后用IDCT变换回图像，比较并计算压缩前后图像的差异(MSE,PSNR等)(1)代码如下所示fori=1:64%用i,j来控制子块forj=1:64form=(8*(i-l)+l):8*iforn=(8*(j-!)+l):8*jA(m-(8*(i-1)+1)+1,n-(8*(j-1)+1)+

5、1)=B(m，n);endendC3=A.*lighttable;°/。将子块反量化回去C4=idct2(C3);%对C3进行IDCT变换forp=l:8forq=l:8D(p+8*(i-l),q+8*(j-l))=C4(p，q);endendendendD=uint8(D);subplot(1，2,1)，imshow(J);subplot(l,2,2),imshow(D);%计算J和D均方根误差e=double(J)-double(D);fm,n]=size(e);MSE=sqrt(sum(sum(e.A2))/(m*n));k=8;%k为阁像是表示地个像素点所用的二进制

6、位数，即位深。fmax=2.八k-1;a=fmax.八2;b=sum(sum(e.A2));PSNR=10*log(m*n氺a/b);(2)显示压缩前后图像(3)显示MSE和PSNR的值MSE=4.4339PSNR=81.043,根据由1计算得到的直方图试判断采用Huffman编码的平均码长。答：由1的直方图可以确定一个大的分布值为-3—3。即编码需要的比特数位3位。其中有一位来表示正负。六，讨论与分析1，有可能增强压缩效率的部分在哪里？在取直方图的时候，综合整幅图像的数据。使得值分布得更集中，可以更有效地去除数据冗余度从而压缩数据。