图像压缩编码程序设计

*******************实践教学*******************兰州理工大学计算机与通信学院2023年秋季学期图像处理综合训练题目：图像压缩编码专业班级：姓名：学号：指导教师：成绩：目录TOC\o"1-2"\u摘要1一、前言2二、算法分析与描述3三、详细设计过程5四、调试过程中出现的问题及相应解决方法7五、程序运行截图及其说明8六、简单操作手册11设计总结14参考资料15致谢16附录17摘要本实验是在根据图像压缩根本原理的根底上，独立设计图像的有损压缩编码算法，并在MATLAB编程环境中编程实现。其目的是理解有损压缩的概念，特别是通过K-L变换或DCT变换，掌握压缩编码的根本原理、算法，并通过编程给予实现，对原始图像数据和解压后的图像数据进行失真估计，了解压缩性能及其重建效果的影响、作用，给出适当的理论分析。关键词：Matlab；图像压缩；解码；行程编码。一、前言人类正在实践的21世纪是一个信息时代。今天的社会，信息技术已经全面效劳于社会生产和生活的方方面面，人们所做工作的相当一局部就是对信息的处理和传输，其中图像是人类获取信息、表达信息和传递信息的重要手段，是人类感知和认识世界的根底。有关研究说明，日常生活中人们所接受的各种信息中图像信息占总信息量的80%左右，从这一角度看，“百闻不如一见〞正是图像处理重要性的形象表达和经验总结。因此，数字图像处理技术无论是对于21世纪的科学理论研究，还是工程应用都将具有重要的影响。国内外许多有识之士指出，从某种意义上讲，数字图像处理是实现智能计算机、智能机器人或多媒体信系统的根底，未来计算机及智能机器人的开展与进步将在一定程度上依赖于机器视觉信息处理理论和技术的突破。虽然相对于经典学科，数字图像处理还很年轻，但日趋成熟的数字图像处理技术已经在很多方面得到深入而广泛的应用，一定程度上改变了人类的生活，给人们的日常生活、学习、工作带来极大的方便。例如，Internet上的视频播送、现代卫星或遥感照片的合成和处理、工业产品的自动检测、各种医学影像和图像的处理、远程医疗诊断及手术以及视频会议、视频等都采用了图像处理技术实现实时信息交互。如果说1964年美国喷气推进实验室首次处理了由太空船“徘徊者七号〞发回的月球照片，拉开了数字图像处理技术进入普遍应用的序幕，那么CT的创造、应用及诺贝尔奖的获得，那么使得数字图像处理技术大放异彩，并迅速进入了广泛应用阶段。目前，数字图像处理技术已在工程科学、计算机科学、信息科学、遥感、遥测、采矿、地质勘探、工业探伤、自动控制、机器人、军事、公安、生物学、医学、视频、多媒体、统计学甚至社会科学等领域得到了广泛应用，并显示出了更加诱人的前景，成为了包括计算机科学与技术、信息科学、航空航天和生物医学工程等在内的多学科的研究重点和热点。这些学科的研究成果又促使了图像处理技术向更高水平开展，数字图像处理技术正是在这种应用的迫切需要和自身的不断开展之中逐步完善的新兴学科。未来，图像处理技术的开展及应用与经济建设联系之紧密、影响之深远是不可估量的。算法分析与描述编码是方法建立在图像统计特性的根底上的。例如，在通信中的文件大多是二值图像，即每个像素的灰度值只有0和1两种取值。将一行中颜色值相同的相邻象素用一个计数值和该颜色值来代替。例如aaabccccccddeee可以表示为3a1b6c2d3e，，即有3个a，1个b，6个c，2个d，3个e。如果一幅图象是由很多块颜色相同的大面积区域组成，那么采用行程编码的压缩效率是惊人的。然而，该算法也导致了一个致命弱点，如果图象中每两个相邻点的颜色都不同，用这种算法不但不能压缩，反而数据量增加一倍。因此对有大面积色块的图像用行程编码效果比拟好。行程编码的可行性讨论：行程编码的压缩方法对于自然图片来说是不太可行的，因为自然图片像素点错综复杂，同色像素连续性差，如果硬要用行程编码方法来编码就适得其反，图像体积不但没减少，反而加倍。鉴于计算机桌面图，图像的色块大，同色像素点连续较多，所以行程编码对于计算机桌面图像来说是一种较好的编码方法。2.1图像的调入：functionyc%行程编码算法%例如aaabccccccddeee才可以表示为3a1b6c2d3e%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%读图I=imread('Lena.bmp');[mnl]=size(I);fid=fopen('yc.txt','w');%yc.txt是行程编码算法的灰度级及其相应的编码表%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%2.2行程编码算法：sum=0;fork=1:lfori=1:mnum=0;J=[];value=I(i,1,k);forj=2:nifI(i,j,k)==valuenum=num+1;%统计相邻像素灰度级相等的个数ifj==nJ=[J,num,value];endelseJ=[J,num,value];%J的形式是先是灰度的个数及该灰度的值value=I(i,j,k);num=1;endendcol(i,k)=size(J,2);%记录Y中每行行程行程编码数sum=sum+col(i,k);Y(i,1:col(i,k),k)=J;%将I中每一行的行程编码J存入Y的相应行中endend%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%输出相关数据[m1,n1,l1]=size(Y);disp('ԭͼÏñ´óС:')whos('I');disp('ѹËõͼÏñ´óС:')whos('Y');disp('ͼÏñµÄѹËõ±È:');disp(m*n*l/sum);2.3图像显示和保存:subplot(1,2,1),imshow(I),title('原图像')subplot(1,2,2),imshow(Y),title('行程编码解码后的图像')save('Y')%存储，以便解码用save('col')fclose(fid);三、详细设计过程3.1算法流程图读入图像读入图像将行程编码写入yc.txt中行程算法计算压缩比3.2由I=imread('Lena.bmp')把图像调入等待处理，再由行程编码算法进行计算和压缩。实现具体的压缩功能。最后输出原图像和压缩后的像和输出压缩比例。3.2行程编码算法sum=0;fork=1:lfori=1:mnum=0;J=[];value=I(i,1,k);forj=2:nifI(i,j,k)==valuenum=num+1;%统计相邻像素灰度级相等的个数ifj==nJ=[J,num,value];endelseJ=[J,num,value];%J的形式是先是灰度的个数及该灰度的值value=I(i,j,k);num=1;endendcol(i,k)=size(J,2);%记录Y中每行行程行程编码数sum=sum+col(i,k);Y(i,1:col(i,k),k)=J;%将I中每一行的行程编码J存入Y的相应行中endend%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%输出相关数据[m1,n1,l1]=size(Y);disp('原图像大小:')whos('I');disp('压缩图像大小:')whos('Y');disp('图像的压缩比:');disp(m*n*l/sum);四、调试过程中出现的问题及相应解决方法首先，我们的代码不完全是自己写出来的，所以刚开始调式时出现很多错误。比方；代码中的函数matlab软件不认识，不能实现函数的功能。还有函数变量定义不符合，但更多的是调入图像显示不出结果。经过研究和修改，功夫不负有心人，大局部错误都改正了。可是很纠结的是只输出图像压缩的数据和压缩比例，没有输出原图像和压缩后的图像。最后经过同学的帮助和自己的努力，发现代码缺少输出图像语句。最终顺利完成了本次的课程设计。五、程序运行截图及其说明图5.1图像的大小和压缩比例图5.2图像原图和压缩后图像图5.3行程编码表六、简单操作手册首先我们完成代码局部。经过查找资料和网上查询，完整的把代码的功能能在matlab软件上是先出来。文档局部比拟繁琐，需要结合代码的功能一步一步经过说明实现图像的行程压缩编码。在经过matlab软件运行输出图像，把图像截图放在文档局部加以说明。界面制作；图6.1画出两个axes1、axes2图像显示框图6.2把对应的程序添加到界面图6.3显示出原图像设计总结涉及到图像压缩的算法很多,我所用的是行程编码算法。这种算法是无损压缩，利用统计的方法，最终建立了数据与编码的关系，即利用了编码冗余实现压缩。在这两周的图像处理的综合训练中,通过对图象处理的设计过程,我加深了对图像压缩的理解,对DCT变换和行程编码和哈夫曼有所掌握,对课本中所学的各种图像知识的进一步理解和掌握,学会了如何把学到的知识用于解决实际问题,锻炼了自己动手的能力。在这次的综合训练中我认识到了自己的缺乏。在以后的学习中我会更加注意各个方面的能力的协调开展。在课程设计时遇到了很多的问题，在老师的帮助，和对各种资料的查阅中，将问题解决，培养了我自主动手，独立研究的能力，为今后在学习工作中能更好的开展打下了坚实的根底。通过对此次图象处理综合训练的实际操作及分析，加深对图像的理解及增强实际动手能力，锻炼了我们分析与编写软件代码的能力及团结协作的能力。两周的课程设计很短暂，但其间的内容是很充实的，在其中我学习到了很多平时书本中无法学到的东西，积累了经验，锻炼了自己分析问题，解决问题的能力，并学会了如何将所学的各科知识融会，组织，来配合学习，为以后的开展打下坚实的根底。参考文献张志涌，杨祖樱.MATLAB教程.北京：北京航空航天大学出版社，2006,9-201.刘卫国.MATLAB程序设计教程[M].北京：中国水利水电出版社，2005,101-150章毓晋.图像工程〔上册〕——图像分析.第二版.北京：清华大学出版社，2005，7-30章毓晋.图像工程〔中册〕——图像分析.第二版.北京：清华大学出版社，2005，30-87黄爱民，安向京，骆力.数字图像处理与分析根底.北京:中国水利水电出版社，2005,7-50致谢本课程设计的选题、材料来源、程序设计和制作和内容的安排、总体篇章的写作等都得到了指导老师徐志刚的真诚指导和关注，以至于我的课程设计能够按时顺利的完成。因此，特别感谢指导老师徐志刚，他以热心和耐心成全了我的课程设计，为我制作的全过程提出了许多珍贵的意见，特别是在编程方面给我很大程度的指导，辅助我尝试了各种函数的应用，使我掌握了许多珍贵的知识。所以真的很感谢徐志刚老师的指导！另外，我也很感谢班里的很多同学给我提出的意见和想法，对我的课程设计提供了很多的帮助！最后，我希望同学们能获得更多的收获!附录functionyc%%行程编码算法%例如aaabccccccddeee才可以表示为3a1b6c2d3e%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%读图I=imread('Lena.bmp');[mnl]=size(I);fid=fopen('yc.txt','w');%yc.txt是行程编码算法的灰度级及其相应的编码表%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%行程编码算法sum=0;fork=1:lfori=1:mnum=0;J=[];value=I(i,1,k);forj=2:nifI(i,j,k)==valuenum=num+1;%统计相邻像素灰度级相等的个数ifj==nJ=[J,num,value];endelseJ=[J,num,value];%J的形式是先是灰度的个数及该灰度的值value=I(i,j,k);num=1;endendcol(i,k)=size(J,2);%记录Y中每行行程行程编码数sum=sum+col(i,k);Y(i,1:col(i,k),k)=J;%将I中每一行的行程编码J存入Y的相应行中endend%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%输出相关数据[m1,n1,l1]=size(Y);disp('ԭͼÏñ´óС:')whos('I');disp('ѹËõͼÏñ´óС:')whos('Y');disp('ͼÏñµÄѹËõ±È:');disp(m*n*l/sum);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%将编码写入yc.txt中fork=1:l1fori=1:m1forj=1:col(i,k)fprintf(fid,'%d',Y(i,j,k));fwrite(fid,'');endendfwrite(fid,'');endsubplot(1,2,1),imshow(I),title('原图像)%¶Á³öԭͼsubplot(1,2,2),imshow(Y),title('行程编码解码后的图像')save('Y')%存储，以便解码用save('col')fclose(fid);functionvarargout=yc(varargin)%YCMATLABcodeforyc.fig%YC,byitself,createsanewYCorraisestheexisting%singleton*.%%H=YCreturnsthehandletoanewYCorthehandleto%theexistingsingleton*.%%YC('CALLBACK',hObject,eventData,handles,...)callsthelocal%functionnamedCALLBACKinYC.Mwiththegiveninputarguments.%%YC('Property','Value',...)createsanewYCorraisesthe%existingsingleton*.Startingfromtheleft,propertyvaluepairsare%appliedtotheGUIbeforeyc_OpeningFcngetscalled.An%unrecognizedpropertynameorinvalidvaluemakespropertyapplication%stop.Allinputsarepassedtoyc_OpeningFcnviavarargin.%%*SeeGUIOptionsonGUIDE'sToolsmenu.Choose"GUIallowsonlyone%instancetorun(singleton)".%%Seealso:GUIDE,GUIDATA,GUIHANDLES%Edittheabovetexttomodifytheresponsetohelpyc%LastModifiedbyGUIDEv2.511-Jan-202311:24:40%Begininitializationcode-DONOTEDITgui_Singleton=1;gui_State=struct('gui_Name',mfilename,...'gui_Singleton',gui_Singleton,...'gui_OpeningFcn',@yc_OpeningFcn,...'gui_OutputFcn',@yc_OutputFcn,...'gui_LayoutFcn',[],...'gui_Callback',[]);ifnargin&&ischar(varargin{1})gui_State.gui_Callback=str2func(varargin{1});endifnargout[varargout{1:nargout}]=gui_mainfcn(gui_State,varargin{:});elsegui_mainfcn(gui_State,varargin{:});end%Endinitializationcode-DONOTEDIT%Executesjustbeforeycismadevisible.functionyc_OpeningFcn(hObject,eventdata,handles,varargin)%Thisfunctionhasnooutputargs,seeOutputFcn.%hObjecthandletofigure%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%varargincommandlineargumentstoyc(seeVARARGIN)%Choosedefaultcommandlineoutputforychandles.output=hObject;%Updatehandlesstructureguidata(hObject,handles);%UIWAITmakesycwaitforuserresponse(seeUIRESUME)%uiwait(handles.figure_yc);%Outputsfromthisfunctionarereturnedtothecommandline.functionvarargout=yc_OutputFcn(hObject,eventdata,handles)%varargoutcellarrayforreturningoutputargs(seeVARARGOUT);%hObjecthandletofigure%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%Getdefaultcommandlineoutputfromhandlesstructurevarargout{1}=handles.output;%functionm_file_Callback(hObject,eventdata,handles)%hObjecthandletom_file(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)%functionm_file_open_Callback(hObject,eventdata,handles)%hObjecthandletom_file_open(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)[filename,pathname]=uigetfile(...{'*.bmp;*.jpg;*.png;*.jpeg','ImageFiles(*.bmp,*.jpg,*.png,*.jpeg)';...'*.*','AllFiles(*.*)'},...'Pickanimage');axes(handles.axes1);%用axes命令设定当前操作的坐标轴是axes_srcfpath=[pathnamefilename];%将文件名和目录名组合成一个完整的路径img_src=imread(fpath);imshow(img_src);%用imread读入图片，并用imshow在axes_src上显示setappdata(handles.figure_yc,'img_src',img_src);title('原图');%functionm_file_save_Callback(hObject,eventdata,handles)%hObjecthandletom_file_save(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)[filename,pathname]