图像无损压缩程序设计 霍夫曼编码.docx
《图像无损压缩程序设计 霍夫曼编码.docx》由会员分享，可在线阅读，更多相关《图像无损压缩程序设计 霍夫曼编码.docx（14页珍藏版）》请在冰豆网上搜索。
图像无损压缩程序设计霍夫曼编码
成绩评定表
学生姓名
班级学号
1203030119
专业
电子信息工程
课程设计题目
图像的无损压缩程序设计——霍夫曼编码
评
语
组长签字：
成绩
日期
2015年7月20日
课程设计任务书
学院
信息科学与工程
专业
电子信息工程
学生姓名
班级学号
1203030119
课程设计题目
图像的无损压缩程序设计——霍夫曼编码
实践教学要求与任务：
本课题通过MATLAB编写适当的函数，对一个随机信源进行哈夫曼编码，得出码字，平均码长和编码效率。
从而理解信源编码的基本思想与目的以及哈夫曼编码方法的基本过程与特点，并且提高综合运用所学理论知识独立分析和解决问题的能力。
工作计划与进度安排：
2015年7月8日—11日：
熟悉编程环境，查阅相关资料。
2015年7月11日—12日：
图像的霍夫曼编码程序设计。
2015年7月12日—13日：
编码、调试、实验与分析。
2015年7月13日—15日：
撰写课程设计报告。
2015年7月15日—19日：
准备答辩。
指导教师：
2015年7月2日
专业负责人：
2015年7月2日
学院教学副院长：
2015年7月2日
摘要
哈夫曼编码（HuffmanCoding）是一种编码方式，以哈夫曼树—即最优二叉树，带权路径长度最小的二叉树，经常应用于数据压缩。
在计算机信息处理中，“哈夫曼编码”是一种一致性编码法（又称"熵编码法"），用于数据的无损耗压缩。
这一术语是指使用一张特殊的编码表将源字符（例如某文件中的一个符号）进行编码。
这张编码表的特殊之处在于，它是根据每一个源字符出现的估算概率而建立起来的（出现概率高的字符使用较短的编码，反之出现概率低的则使用较长的编码，这便使编码之后的字符串的平均期望长度降低，从而达到无损压缩数据的目的）。
本课题通过MATLAB编写适当的函数，对一个随机信源进行哈夫曼编码，得出码字，平均码长和编码效率。
从而理解信源编码的基本思想与目的以及哈夫曼编码方法的基本过程与特点，并且提高综合运用所学理论知识独立分析和解决问题的能力。
关键字：
哈夫曼；信源编码；MATLAB
1设计目的及相关知识
1.1设计目的
1）了解霍夫曼编码的原理。
2）理解图像的霍夫曼编码原理，了解其应用，掌握图像的霍夫曼编码的方法。
3）对图像编码程序设计进行较深入的认识，对知识牢固掌握。
4）掌握图像霍夫曼编码的整个过程及其中的注意事项。
5）了解图像无损压缩的目的及好处。
1.2图像的霍夫曼编码概念
所谓霍夫曼编码的具体方法：
先按出现的概率大小排队，把两个最小的概率相加，作为新的概率和剩余的概率重新排队，再把最小的两个概率相加，再重新排队，直到最后变成1。
每次相加时都将“0”和“1”赋与相加的两个概率，读出时由该符号开始一直走到最后的“1”，将路线上所遇到的“0”和“1”按最低位到最高位的顺序排好，就是该符号的霍夫曼编码
1.3Matlab图像处理通用函数
colorbar 显示彩色条
语法：
colorbar\colorbar（'vert'）\colorbar（'horiz'）\colorbar（h）\h=colorbar（...）\colorbar（...,'peer',axes_handle）
getimage从坐标轴取得图像数据
语法：
A=getimage（h）\[x,y,A]=getimage（h）\[...,A,flag]=getimage（h）\[...]=getimage
imshow显示图像
语法：
imshow（I,n）\imshow（I,[lowhigh]）\imshow（BW）\imshow（X,map）\imshow（RGB）\imshow（...,display_option）\imshow（x,y,A,...）\imshowfilename\h=imshow（...）
montage在矩形框中同时显示多幅图像
语法：
montage（I）\montage（BW）\montage（X,map）\montage（RGB）\h=montage（...）
immovie创建多帧索引图的电影动画
语法：
mov=immovie（X,map）\mov=immovie（RGB）
subimage在一副图中显示多个图像
语法：
subimage（X,map）\subimage（I）\subimage（BW）\ subimage（RGB）\subimage（x,y,...）\subimage（...）
truesize调整图像显示尺寸
语法：
truesize（fig,[mrowsmcols]）\truesize（fig）
warp将图像显示到纹理映射表面
语法：
warp（X,map）\warp（I,n）\warp（z,...）warp（x,y,z,...）\ h=warp（...）
zoom缩放图像
语法：
zoomon\zoomoff\zoomout\zoomreset\zoom\zoomxon\zoomyon\zoom（factor）\zoom（fig,option）
2课程设计分析
2.1图像的霍夫曼编码概述
赫夫曼（Huffman）编码是1952年提出的，是一种比较经典的信息无损熵编码，该编码依据变长最佳编码定理，应用Huffman算法而产生。
Huffman编码是一种基于统计的无损编码。
根据变长最佳编码定理，Huffman编码步骤如下：
（1）将信源符号xi按其出现的概率，由大到小顺序排列。
（2）将两个最小的概率的信源符号进行组合相加，并重复这一步骤，始终将较大的概率分支放在上部，直到只剩下一个信源符号且概率达到1.0为止；
（3）对每对组合的上边一个指定为1，下边一个指定为0（或相反：
对上边一个指定为0，下边一个指定为1）；
（4）画出由每个信源符号到概率1.0处的路径，记下沿路径的1和0；
（5）对于每个信源符号都写出1、0序列，则从右到左就得到非等长的Huffman码。
Huffman编码的特点是：
（1）Huffman编码构造程序是明确的，但编出的码不是唯一的，其原因之一是两个概率分配码字“0”和“1”是任意选择的（大概率为“0”，小概率为“1”，或者反之）。
第二原因是在排序过程中两个概率相等，谁前谁后也是随机的。
这样编出的码字就不是唯一的。
（2）Huffman编码结果，码字不等长，平均码字最短，效率最高，但码字长短不一，实时硬件实现很复杂（特别是译码），而且在抗误码能力方面也比较差。
（3）Huffman编码的信源概率是2的负幂时，效率达100%，但是对等概率分布的信源，产生定长码，效率最低，因此编码效率与信源符号概率分布相关，故Huffman编码依赖于信源统计特性，编码前必须有信源这方面的先验知识，这往往限制了霍夫曼编码的应用。
（4）Huffman编码只能用近似的整数位来表示单个符号，而不是理想的小数，这也是Huffman编码无法达到最理想的压缩效果的原因。
2.2图像的霍夫曼编码举例
假设一个文件中出现了8种符号S0，S1，S2，S3，S4，S5，S6，S7，那么每种符号要编码，至少需要3比特。
假设编码成000，001，010，011，100，101，110，111那么符号序列S0S1S7S0S1S6S2S2S3S4S5S0S0S1编码后变成000001111000001110010010011100101000000001，共用了42比特。
我们发现S0，S1，S2这三个符号出现的频率比较大，其它符号出现的频率比较小，如果我们采用一种编码方案使得S0，S1，S2的码字短，其它符号的码字长，这样就能够减少占用的比特数。
例如，我们采用这样的编码方案：
S0到S7的码字分别01，11，101，0000，0001，0010，0011，100，那么上述符号序列变成011110001110011101101000000010010010111，共用了39比特，尽管有些码字如S3，S4，S5，S6变长了（由3位变成4位），但使用频繁的几个码字如S0，S1变短了，所以实现了压缩。
可由下面的步骤得到霍夫曼码的码表
（1）首先把信源中的消息出现的频率从小到大排列。
（2）每一次选出频率最小的两个值，作为二叉树的两个叶子节点，将和作为它们的根节点，这两个叶子节点不再参与比较，新的根节点参与比较。
（3）重复
（2），直到最后得到和为1的根节点。
（4）将形成的二叉树的左节点标0，右节点标1。
把从最上面的根节点到最下面的叶子节点途中遇到的0，1序列串起来，就得到了各个符号的编码。
上面的例子用Huffman编码的过程如图下图所示，其中圆圈中的数字是新节点产生的顺序。
图2-1 Huffman编码的二叉树示意图
信源的各个消息从S0到S7的出现概率分别为4/14，3/14，2/14，1/14，1/14，1/14，1/14，1/14。
计算编码效率为98.5%，编码的冗余只有1.5%，可见霍夫曼编码效率很高。
产生Huffman编码需要对原始数据扫描两遍。
第一遍扫描要精确地统计出原始数据中，每个值出现的频率，第二遍是建立Huffman树并进行编码。
由于需要建立二叉树并遍历二叉树生成编码，因此数据压缩和还原速度都较慢，但简单有效，因而得到广泛的应用。
3仿真
主程序：
%以下为主程序mainp.m
clc
clear
closeall;
%定义HufData/Len为全局变量的结构体
globalHufData;
globalLen
disp（'计算机正在准备输出霍夫曼编码结果,请耐心等待……'）;
%原始码字的灰度
a=imread（'kids.tif'）;
%分区画出原始图像和灰度直方图
figure;
subplot（1,2,1）
imshow（a）;
%取消坐标轴和边框
axisoff
boxoff
title（'MATLAB自带图像','fontsize',13）;
subplot（1,2,2）;
axisoff
boxoff
imhist（a）;
title（'图像灰度直方图','fontsize',13）;
%图像的灰度统计
GrayStatistics=imhist（a）;
GrayStatistics=GrayStatistics';
GrayRatioo=GrayStatistics/sum（GrayStatistics）;
GrayRatioNO=find（GrayRatioo~=0）;
Len=length（GrayRatioNO）;
%初始化灰度集，防止系统随即赋予其垃圾值
GrayRatio=ones（1,Len）;
fori=1:
Len
GrayRatio（i）=GrayRatioo（i）;
end
GrayRatio=abs（sort（-GrayRatio））;
%将图像灰度概率赋予结构体
fori=1:
Len
HufData（i）.value=GrayRatio（i）;
end
%霍夫曼编码/霍夫曼编码
HuffmanCode（Len）;
%输出码字
zippedHuffman=1;
fori=1:
Len
tmpData=HufData（i）.code;
str='';
forj=1:
length（tmpData）
str=strcat（str,num2str（tmpData（j）））;
zippedHuffman=zippedHuffman+1;
end
disp（strcat（'a',num2str（i）,'=',str））
end
i;
%计算计算机一共输出多少个霍夫曼编码/霍夫曼编码
zippedHuffman;
%计算在删去0灰度级压缩之前的原始图像字节容量
unzipped_delete=i*8;
%计算压缩比率
ratio_delete=zippedHuffman/unzipped_delete;
%计算图像的压缩比率
ad=num2str（ratio_delete*100）;
str2=strcat（ad,'%'）;
disp（strcat（'霍夫曼编码压缩比率','=',str2））
4结果及分析
结果:
图4-1输出原图像与该图像像灰度直方图
计算机正在准备输出霍夫曼编码结果,请耐心等待……
a1=110
a2=11110
a3=11101
a4=01100
a5=01010
a6=01000
a7=00101
a8=00011
a9=111111
a10=111001
a11=101111
a12=101100
a13=101011
a14=101010
a15=101001
a16=100111
a17=100110
a18=100100
a19=100011
a20=100010
a21=100001
a22=100000
a23=011111
a24=011110
a25=011011
a26=011010
a27=010111
a28=010110
a29=010011
a30=001111
a31=001101
a32=001100
a33=001001
a34=001000
a35=000101
a36=000011
a37=000010
a38=000001
a39=000000
a40=1111101
a41=1111100
a42=1110001
a43=1110000
a44=1011101
a45=1011100
a46=1011011
a47=1010001
a48=1010000
a49=1001011
a50=1001010
a51=0111011
a52=0111010
a53=0111001
a54=0111000
a55=0100101
a56=0100100
a57=0011101
a58=0011100
a59=0001001
a60=0001000
a61=10110101
a62=101101001
a63=101101000
霍夫曼编码压缩比率=78.9683%
分析：
从输出灰度直方图可得出该图像的量化值主要集中在低灰度级处，通过输出可以看到该灰度级对应的霍夫曼编码，并且输出了该图像的压缩效率。
明显可得出霍夫曼编码大大的节省了空间，可以明显的减少发送时间。
5附录
子程序：
%子程序：
霍夫曼编码/霍夫曼编码函数HuffmanCode.m
functionHuffmanCode（OriginSize）
globalHufData;
globalLen
fori=1:
Len
%%霍夫曼编码树左边纪录为1
HufData（i）.left=1;
%%霍夫曼编码树右边纪录为0
HufData（i）.right=0;
%%输出码初始化为0
HufData（i）.code=[];
%%排序列表初始化
SortList（i）.symbol=i;
SortList（i）.value=HufData（i）.value;
end
%初始化原始消息数目
newsymbol=OriginSize;
forn=OriginSize:
-1:
2
%将N个消息进行排序
SortList=sortdata（SortList,n）;
%将最后两个出现概率最小的消息合成一个消息
newsymbol=newsymbol+1;
HufData（newsymbol）.value=SortList（n-1）.value+SortList（n）.value;
HufData（newsymbol）.left=SortList（n-1）.symbol;
HufData（newsymbol）.right=SortList（n）.symbol;
%将消息添加到列队的最后，为N-1个消息重新排序作好准备
SortList（n-1）.symbol=newsymbol;
SortList（n-1）.value=HufData（newsymbol）.value;
end
%遍历霍夫曼树，获得霍夫曼编码/霍夫曼编码
visit（newsymbol,Len,[]）;
end
%子程序：
冒泡排序法函数sortdata.m
functionreData=sortdata（SortList,n）
%根据消息概率进行排序
fork=n:
-1:
2
forj=1:
k-1
min=SortList（j）.value;
sbl=SortList（j）.symbol;
if（minSortList（j）.value=SortList（j+1）.value;
SortList（j+1）.value=min;
SortList（j）.symbol=SortList（j+1）.symbol;
SortList（j+1）.symbol=sbl;
end
end
end
reData=SortList;
end
%子程序：
遍历霍夫曼编码/霍夫曼编码树搜索函数visit.m
functionvisit（node,n,ocode）
globalHufData
ifnode<=n
%如果没有霍夫曼编码/霍夫曼编码树的子接点直接输出原始码，这里为空码（[]）HufData（node）.code=ocode;
else
if（HufData（node）.left>0）
%遍历左分支接点输出1，这里采用子函数嵌套调用
ocode1=[ocode1];
visit（HufData（node）.left,n,ocode1）;
end
if（HufData（node）.right>0）
%遍历右分支接点输出0，这里采用子函数嵌套调用
ocode2=[ocode0];
visit（HufData（node）.right,n,ocode2）;
end
end
end
结束语
数字图像处理应用很广泛，在科学研究、工农业生产、军事、公安、医疗卫生、教育等领域都具有广泛的应用，所以近来有了长足的发展。
MATLAB既是一种直观、高效的计算机语言,同时又是一个科学计算平台。
它为数据分析和数据可视化、算法和应用程序开发提供了最核心的数学和高级图形工具。
根据它提供的许多个数学和工程函数，工程技术人员和科学工作者可以在它的集成环境中交互或编程以完成各自的计算。
MATLAB中集成了功能强大的图像处理工具箱。
本次设计使我的知识面有了很大的提高，专业知识有了进一步的强化，对图像编码程序设计有了较深入的认识。
由于对知识掌握的不够牢固、全面，在设计中遇到了挺多问题，书中的程序有一些错误，通过查找资料和同学的帮助，才勉强得到结果，总体设计还算好，不过因为时间比较仓促，还有很多瑕疵，这也为我进一步学习这方面的知识提供了帮助。
这次课程设计给我们提供了很好的实践机会，有利于磨练我们的耐心，去认真的去做好一件事，对以后的学习机工作有很大的帮助！
参考文献
[1]刘刚.MATLAB数字图像处理[M].北京:
机械工业出版社,2010:
34-45.
[2]王家文.MATLAB6.5图形图像处理[M].上海:
国防工业出版社,2009:
6-14.
[3]王晓丹.MATLAB的系统分析与设计[M].西安:
西安电子科技大学出版社,2000:
168-220.
[4]余成波.数字图像处理及MATLAB实现[M].重庆:
重庆大学出版社,2003:
34-54.
[5]郝文化.MATLAB图形图像处理应用教程[M].河北:
中国水利水电出版社,2004:
66-78.
[6]苏金明.MATLAB图形图像[M].四川:
成都电子工业出版社,2005:
87-90.
升级会员 