




已阅读5页,还剩3页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
实验三、图像压缩编码技术一、实验目的1、理解有损压缩和无损压缩的概念;2、理解图像压缩的主要原则和目的;3、了解几种常用的图像压缩编码方式;4、利用MATLAB程序进行图像压缩编码。二、实验原理1、图像压缩原理图像压缩主要目的是为了节省存储空间,增加传输速度。图像压缩的理想标准是信息丢失最少,压缩比例最大。不损失图像质量的压缩称为无损压缩,无损压缩不可能达到很高的压缩比;损失图像质量的压缩称为有损压缩,高的压缩比是以牺牲图像质量为代价的。压缩的实现方法是对图像重新进行编码,希望用更少的数据表示图像。信息的冗余量有许多种,如空间冗余,时间冗余,结构冗余,知识冗余,视觉冗余等,数据压缩实质上是减少这些冗余量。高效编码的主要方法是尽可能去除图像中的冗余成分,从而以最小的码元包含最大的图像信息。编码压缩方法有许多种,从不同的角度出发有不同的分类方法,从信息论角度出发可分为两大类。(1)冗余度压缩方法,也称无损压缩、信息保持编码或熵编码。具体说就是解码图像和压缩编码前的图像严格相同,没有失真,从数学上讲是一种可逆运算。(2)信息量压缩方法,也称有损压缩、失真度编码或烟压缩编码。也就是说解码图像和原始图像是有差别的,允许有一定的失真。应用在多媒体中的图像压缩编码方法,从压缩编码算法原理上可以分为以下3类:(1)无损压缩编码种类哈夫曼(Huffman)编码,算术编码,行程(RLE)编码,Lempel zev编码。(2)有损压缩编码种类预测编码,DPCM,运动补偿;频率域方法:正交变换编码(如DCT),子带编码;空间域方法:统计分块编码;模型方法:分形编码,模型基编码;基于重要性:滤波,子采样,比特分配,向量量化;(3)混合编码。有JBIG,H.261,JPEG,MPEG等技术标准。本实验主要利用MATLAB程序进行赫夫曼(Huffman)编码和行程编码(Run Length Encoding, RLE)。三、实验内容1、实现基本JPEG的压缩和编码分三个步骤:(1)首先通过DCT变换去除数据冗余;(2)使用量化表对DCT系数进行量化;(3)对量化后的系数进行Huffman编码。四、实验步骤1打开计算机,启动MATLAB程序;2选择一幅图像,并进行赫夫曼和行程编码压缩处理;3 将原图像在Photoshop软件中打开,分别以不同的位图文件格式进行“另保存”,比较它们的数据量。4记录和整理实验报告五、实验仪器1计算机;2 MATLAB、Photoshop等程序;3移动式存储器(软盘、U盘等)。4记录用的笔、纸。六、实验报告内容1叙述实验过程;2提交实验的原始图像和结果图像。clearload woman;X=imread(circuit.tif,tif);data=uint8(X);zipped,info=huffencode(data);unzipped=huffdecode(zipped,info,data);subplot(121);imshow(data);subplot(122);imshow(unzipped);%erms=compare(data(:),unzipped(:)cr=info.ratiowhos data unzipped zipped%huffencodevectorHuffman% function zipped,info=huffencode(vector)%unit8%info%info.pad%info.huffcodesHuffman%info.rows%info.cols%info.length%info.maxcodelenif isa(vector,uint8) error(input argument must be a uint8 vector);endm,n=size(vector);vector=vector(:);f=frequency(vector); %symbols=find(f=0);f=f(symbols);f,sortindex=sort(f); %symbols=symbols(sortindex);len=length(symbols);symbols_index=num2cell(1:len);codeword_tmp=cell(len,1);while length(f)1 %Huffman index1=symbols_index1; index2=symbols_index2; codeword_tmp(index1)=addnode(codeword_tmp(index1),uint8(0); codeword_tmp(index2)=addnode(codeword_tmp(index2),uint8(1); f=sum(f(1:2) f(3:end); symbols_index=index1,index2 symbols_index(3:end); f,sortindex=sort(f); symbols_index=symbols_index(sortindex);endcodeword=cell(256,1);codeword(symbols)=codeword_tmp;len=0;for index=1:length(vector) % len=len+length(codeworddouble(vector(index)+1);endstring=repmat(uint8(0),1,len);pointer=1;for index=1:length(vector) % code=codeworddouble(vector(index)+1; len=length(code); string(pointer+(0:len-1)=code; pointer=pointer+len;endlen=length(string);pad=8-mod(len,8); %8pad0if pad0 string=string uint8(zeros(1,pad);endcodeword=codeword(symbols);codelen=zeros(size(codeword);weights=2.(0:23);maxcodelen=0;for index=1:length(codeword) len=length(codewordindex); if lenmaxcodelen maxcodelen=len; end if len0 code=sum(weights(codewordindex=1); code=bitset(code,len+1); codewordindex=code; codelen(index)=len; endendcodeword=codeword:;%cols=length(string)/8;string=reshape(string,8,cols);weights=2.(0:7);zipped=uint8(weights*double(string);%huffcodes=sparse(1,1);for index=1:nnz(codeword) huffcodes(codeword(index),1)=symbols(index);end%info.pad=pad;info.huffcodes=huffcodes;info.ratio=cols./length(vector);info.length=length(vector);info.maxcodelen=maxcodelen;info.rows=m;info.cols=n;%huffdecodevectorHuffman%function vector=huffdecode(zipped,info,image)ifisa(zipped,uint8) error(input argument must be a uint8 vector);end%01len=length(zipped);string=repmat(uint8(0),1,len.*8);bitindex=1:8;for index=1:lenstring(bitindex+8.*(index-1)=uint8(bitget(zipped(index),bitindex);endstring=logical(string(:);len=length(string);%weights=2.(0:51);vector=repmat(uint8(0),1,info.length);vectorindex=1;codeindex=1;code=0;for index=1:len code=bitset(code,codeindex,string(index); codeindex=codeindex+1; byte=decode(bitset(code,codeindex),info); if byte0 vector(vectorindex)=byte-1; codeindex=1; code=0; vectorindex=vectorindex+1; endend%vector=reshape(vector,info.rows,info.cols);%addnodefunction f=frequency(vector)ifisa(vector,uint8) error(input argument must be a uint8 vector);endf=repmat(0,1,256);len=length(vector);for index=0:255 f(index+1)=sum(vector=uint8(index);endf=f./len;%decodefunction codeword_new=addnode(codeword_old,item)codeword_new=cell(size(codeword_old);for index=1:length(codeword_old) codeword_newindex=item codeword_oldindex;end%frequencyfunction byte=decode(code,info)byte=info.huffcodes(code);clear allI=imread(circuit.tif);zipped,info=RLEncode(I);unzipped=RLEdecode(zipped,info);subplot(121);imshow(I);subplot(122);imshow(unzipped);%erms=compare(I(:), unzipped(:)cr=info.ratiowhos I unzipped zipped function unzipped=RLEdecode(zip,info)zip=uint8(zip); m,n=size(zip);unzipped=;for i=1:m section=repmat(zip(i,1),1,double(zip(i,2); unzipped=unzipped section;endunzipped=reshape(unzipped,info.rows ,info.cols);functionzipped,info=RLEncode(vector)m,n=size(vector);vector=vector(:);vector=uint8(vector(:);l=length(vector);c=vector(1);e(
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2025年电厂电工基础试题及答案
- 美术创作指导机器人企业制定与实施新质生产力项目商业计划书
- 有机肥料制作创新创业项目商业计划书
- 肾脏保护儿童茶企业制定与实施新质生产力项目商业计划书
- 美白护肤技术合作行业跨境出海项目商业计划书
- 软包锂离子电池创新创业项目商业计划书
- 网络文学原创与出版创新创业项目商业计划书
- 郑州工会知识题库及答案
- 创业培训讲师考试题及答案
- 潜意识模拟测试题及答案
- 注销代理协议书
- 2025-2026学年高一上学期第一次月考化学试卷(福建)
- LED道路照明灯具采购安装合同(含安全防护与智能化升级)
- 《童年》整本书解读与教学设计
- 卡林巴琴课件
- 甘肃工装装修施工方案
- 2025江苏省药品监督管理局审评中心招聘3人考试参考题库及答案解析
- AI技术在遗传学课程教学中的创新路径
- 设备健康管理模型-第2篇-洞察及研究
- 监狱罪犯教育改造课件
- 新版2026统编版小学道德与法治三年级上册 第4课《 科技力量大》第1课时 科技改变生活和科技改变观念 教学课件
评论
0/150
提交评论