图像压缩综合实验.docVIP

华 北 电 力 大 学 实 验 报 告 院/系：电子与通信工程系 课程名称：数字图像处理 日期：2010-11-2 专业 通信工程 学号 实验时间：14：10-16:10 班级 通信0703 姓名 实验名称 图像压缩综合实验 操作成绩 报告 成绩 教师签名 实验 环境 MATLAB 实验任务 要求 （1）进一步熟悉DCT的概念和原理； （2）掌握对灰度和彩色图像作离散余弦变换和反变换的方法； （3）能选择适当的量化间隔、阈值和区域进行编码； （4）综合利用变换、量化、Huffman编码等知识，实现JPEG压缩编码。 实 验 原 理 与 方 案 JPEG压缩原理： 基于DCT的JPEG压缩过程主要是对图像的空间冗余进行压缩，其基本过程为：首先进行DCT正变换，再对DCT系数进行量化，并对量化后的直流系数和交流系数分别进行差分编码和行程编码，最后进行熵编码。 编码的简化框图如图3.1所示，解码为其逆过程。 量化表 Huffman码表 图3.1 JPEG压缩编码器框图 2-DCT变换原理： JPEG将源数据图像分成8×8大小的子块，然后进行DCT变换，2-DCT变换公式如下： 其中： f（x,y）—输入/输出图像取样值（基准系统的取值为[-128，127]）； C(u,v)—DCT系数（基准系统中C（u,v）的取值范围为[-1023，1023]）； C（0，0）代表DC系数，其余63个为AC系数。 实 验 程 序 DCT变换、量化、区域编码、阈值编码参考程序 num2=0; num3=0; I = imread(cameraman.tif); I = im2double(I)*255; T = dctmtx(8);% 计算离散余弦变换矩阵 figure,imshow(I/255), Q=10; B = blkproc(I,[8 8],P1*x*P2,T,T);%对图像进行不同的快处理 mask = [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 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); B2=round(B2/Q);%round功能为取整 for i=1:256 for j=1:256 if(B2(i,j)==0) num3=num3+1; end end end B3=B2*Q; I2 = blkproc(B3,[8 8],P1*x*P2,T,T); figure,imshow(I2/255) for i=1:256 for j=1:256 if(abs(B(i,j))40) B(i,j)=0; num2=num2+1; end end end I3 = blkproc(B,[8 8],P1*x*P2,T,T); figure,imshow(I3/255) a2=num2, a3=num3 JPEG压缩程序： clear all clc disp(Enter ...) disp(0 (default) - load a 64 x 64 image, or ) disp(1 - Use text book 8 x 8 data (example 8.28)) 实 验 程 序 chos=input(Enter your choice: ); if isempty(chos) chos=0; end if chos==0, load lena.mat f=x; imshow(mat2gray(f)) clear x else f=[139 144 149 153 155 155 155 155 144 151 153 156 159 156 156 156 150 155 160 163 158 156 156 156 159 161 162 160 160 159 159 159 159 160 161 162 162 155 155 155 161 161 161 161 160