基于lbg矢量量化图像压缩编码实验

基于lbg矢量量化图像压缩编码实验

ID:20411346

大小:790.00 KB

页数:11页

时间:2018-10-09

基于lbg矢量量化图像压缩编码实验_第1页
基于lbg矢量量化图像压缩编码实验_第2页
基于lbg矢量量化图像压缩编码实验_第3页
基于lbg矢量量化图像压缩编码实验_第4页
基于lbg矢量量化图像压缩编码实验_第5页
资源描述:

《基于lbg矢量量化图像压缩编码实验》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、基于LBG的矢量量化图像压缩编码实验一、实验原理矢量量化:要想得到好的性能编码,仅采用标量量化是不可能的。当把多个信源符号联合起来形成多维矢量,再对矢量进行标量量化时自由度将更大,同样的失真下,量化基数可进一步减少,码率可进一步压缩。这种量化叫矢量量化。应用:在航天、军事、气象、医学、多媒体等领域中经常需要大量存储和传输各种静态图像和视频图像。为了提高传输效率和减少存储空间,必须采取有效的压缩编码算法消除图像中所包含的各种冗余信息并在给定的失真条件下使用尽量少的比特数来描述图像。矢量量化(VQ)作为一种有效的有损压

2、缩技术,其突出优点是压缩比大以及解码算法简单,因此它已经成为图像压缩编码的重要技术之一。矢量量化压缩技术的应用领域非常广阔,如军事部门和气象部门的卫星(或航天飞机)遥感照片的压缩编码和实时传输、雷达图像和军用地图的存储与传输、数字电视和DVD的视频压缩、医学图像的压缩与存储、网络化测试数据的压缩和传输、语音编码、图像识别和语音识别等等。LGB算法:一种有效和直观的矢量量化码书设计算法——LBG算法(也叫GLA算法)是由Linde、Buzo和Gray于1980年首先提出来的。该算法基于最佳矢量量化器设计的最佳划分和最

3、佳码书这两个必要条件,且是Lloyd算法在矢量空间的推广,其特点为物理概念清晰、算法理论严密及算法实现容易。设训练矢量集为,待产生的码书为,其中,,,则码书设计过程就是需求把训练矢量集分成个子集的一种最佳聚类方案,而子集的质心矢量作为码字。假设平方误差测度用来表征训练矢量和码字之间的失真,即:则码书设计的准则可用下列数学形式表达:最小化约束条件,其中为矩阵,其元素满足:矩阵可看作训练矢量的聚类结果。根据,可计算码字:其中代表子集中训练矢量的数目,或者说是矩阵第行中非零元素的数目。针对训练矢量集为,其LBG算法的具体

4、步骤如下:步骤1:给定初始码书,令迭代次数,平均失真,给定相对误差门限。步骤2:用码书中的各码字作为质心,根据最佳划分原则把训练矢量集划分为个胞腔,满足步骤3:计算平均失真判断相对误差是否满足若满足,则停止算法,码书就是所求的码书。否则,转步骤4。步骤4:根据最佳码书条件,计算各胞腔的质心,即由这个新质心形成新码书,置,转步骤2。二、实验目的采用矢量量化算法(LBG)获得图像压缩所需要的码书,通过码书实现图像压缩编码。三、实验内容对给定的一幅图片四、实验步骤(1)对训练图片,采用LBG算法获取最佳码书设计;(2)采

5、用熵编码实现图像索引编号的压缩。五、程序代码clearall;data=imread('cameraman.tif');%调入原始图像data=double(data)/255;%归一化[m,n]=size(data);%求出图像的行数和列数figure(1)subplot(1,2,1);imshow(data);%显示原始图像title('原始图像')subplot(1,2,2);imhist(data);title('直方图')siz_word=4;%设置码字的大小siz_book=512;%设置码书的大小da

6、ta1=zeros(m*n,1);fori=1:mforj=1:ndata1((i-1)*n+j)=data(i,j);endendM1=floor(m*n/siz_word);r=mod(m*n,siz_word);ifr>0M1=M1+1;enddata2=zeros(M1,siz_word);l=1;A=zeros(siz_word,1);r=1;fori=1:m*nA(r)=data1(i);ifr==siz_worddata2(l,:)=A;l=l+1;r=1;elser=r+1;endendcode_b

7、ook=zeros(siz_book,siz_word);%LBG算法开始%初始化码书l=1;r=1;A=zeros(siz_word,1);fori=1:siz_book*siz_wordA(r)=data1(i);ifr==siz_wordcode_book(l,:)=A;l=l+1;r=1;elser=r+1;endendMIU=zeros(M1,siz_book);t=1;whilet==1fori=1:M1B=zeros(siz_word,1);B=data2(i,:);A=zeros(siz_word,

8、1);A=code_book(1,:);tep=0.0;forl=1:siz_wordtep=tep+(A(l)-B(l))^2;endr=1;forj=2:siz_bookA=code_book(j,:);temp=sum((A-B).^2);iftemp

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。