基于MATLAB实现对图像的矢量量化.docVIP

题目：基于MATLAB实现对图像的矢量量化 姓　 名： 学　 号： 系 别： 电子信息工程 专 业： 电子科学与技术 2014年 10月 11日基于MATLAB实现对图像的矢量量化 一、实验要求： 1.对于图像的数据可以只选用灰度图像； 2.实现对四个像素值的矢量量化，自己选择如何建立码表的方法及最佳匹配的方法； 3.比较原始图像和重构图像的效果； 二、实验内容： 1.实验素材： 在本次的对图像矢量量化的过程中，选取了琳娜（Linna）图像，256×256,12KB，如图2.1所示： 图2.1 2.码表的建立和最佳的匹配方法： 本次实验中选取的码表设计方法为LBG算法，该算法基于最佳矢量量化器设计的最佳划分和最佳码书这两个必要条件，且是Lloyd算法在矢量空间的推广，其特点为物理概念清晰、算法理论严密及算法实现容易。 设训练矢量集为，待产生的码书为，其中，，，则码书设计过程就是需求把训练矢量集分成个子集的一种最佳聚类方案，而子集的质心矢量作为码字。假设平方误差测度用来表征训练矢量和码字之间的失真，即： 则码书设计的准则可用下列数学形式表达： 最小化 约束条件 ， 其中为矩阵，其元素满足： 矩阵可看作训练矢量的聚类结果。根据，可计算码字： 其中 代表子集中训练矢量的数目，或者说是矩阵第行中非零元素的数目。 针对训练矢量集为，其LBG算法的具体步骤如下，流程图如图2.2所示： 图2.2 步骤1：给定初始码书，令迭代次数，平均失真，给定相对误差门限。 步骤2：用码书中的各码字作为质心，根据最佳划分原则把训练矢量集划分为个胞腔，满足 步骤3：计算平均失真 判断相对误差是否满足 若满足，则停止算法，码书就是所求的码书。否则，转步骤4。 步骤4：根据最佳码书条件，计算各胞腔的质心，即 由这个新质心形成新码书，置，转步骤2。 3.实验方法：在对同一幅琳娜（Linna）图像进行矢量量化时，分别采用在不同的码书条件下：512、256、1024观察比较相应的结果。 三、程序代码 clear all; I=imread(C:\linna.JPG); g=rgb2gray(I); imwrite(g,C:\Users\Johnson_Li\Desktop\Mr.Li.JPG,JPG) data=imread(Mr.Li.JPG); %调入原始图像 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; %设置码书的大小 data1=zeros(m*n,1); for i=1:m for j=1:n data1((i-1)*n+j)=data(i,j); end end M1=floor(m*n/siz_word); r=mod(m*n,siz_word); if r0 M1=M1+1; end data2=zeros(M1,siz_word); l=1; A=zeros(siz_word,1); r=1; for i=1:m*n A(r)=data1(i); if r==siz_word data2(l,:)=A; l=l+1; r=1; else r=r+1; end end code_book=zeros(siz_book,siz_word); %LBG算法开始 %初始化码书 l=1; r=1; A=zeros(siz_word,1); for i=1:siz_book*siz_word A(r)=data1(i); if r==siz_word