原创力文档- 1、本文档共4页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
数字图像处理实验八图像变换压缩编码
实验实验目的：实验条件
MATLAB软件
2，典型的灰度，彩色图像。
实验原理实验内容实验步骤I=imread(a.jpg);
J=rgb2gray(I);
J=imresize(J,[512 512]);
%imshow(J);
A=zeros(8,8);
B=zeros(512,512);
lighttable=... %标准亮度量化表
[16 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];
for i=1:64 %用i，j来控制子块
for j=1:64
for m=(8*(i-1)+1):8*i
for n=(8*(j-1)+1):8*j
A(m-(8*(i-1)+1)+1,n-(8*(j-1)+1)+1)=J(m,n);%提取8*8子块
end
end
C1=dct2(A); %对子块进行dct2变换
C2=round(C1./lighttable); %将C1进行量化，并取整
for p=1:8 %将子块按序存进512*512的矩阵中
for q=1:8
B(p+8*(i-1),q+8*(j-1))=C2(p,q);
end
end
end
end
H=zeros(1,256*2-1);
for i=1:m
for j=1:n
H(B(i,j)+255)=H(B(i,j)+255)+1; %统计各值概率
end
end
X=[-255:255];
plot(X,H); %绘制坐标为-255—255的直方图
（2）直方图如下所示
采用标准量化表将上述变换值反量化回去后用IDCT变换回图像，比较并计算压缩前后图像的差异（MSE，PSNR等）
（1）代码如下所示
for i=1:64 %用i，j来控制子块
for j=1:64
for m=(8*(i-1)+1):8*i
for n=(8*(j-1)+1):8*j
A(m-(8*(i-1)+1)+1,n-(8*(j-1)+1)+1)=B(m,n);
end
end
C3=A.*lighttable; %将子块反量化回去
C4=idct2(C3); %对C3进行IDCT变换
for p=1:8
for q=1:8
D(p+8*(i-1),q+8*(j-1))=C4(p,q);
end
end
end
end
D=uint8(D);
subplot(1,2,1),imshow(J);
subplot(1,2,2),imshow(D);
%计算J和D均方根误差
e = double(J) - double(D);
[m, n] = size(e);
MSE = sqrt(sum(sum(e.^2))/(m*n));
k = 8;
%k为图像是表示地个像素点所用的二进制位数，即位深。
fmax = 2.^k - 1;
a = fmax.^2;
b = sum(sum(e.^2));
PSNR = 10*log(m*n*a/b);
（2）显示压缩前后图像
（3）显示MSE和PSNR的值
MSE=4.4339
PSNR=81.04
根据由1计算得到的直方图试判断采用Huffman编码的平均码长。
答：由1的直方图可以确定一个大的分布值为-3—3。即编码需要的比特数位3位。其中有一位来表示正负。
讨论与分析
文档评论（0）