资源描述:
《小波变换在图像处理中的运用及其matlab实现 - 副本》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、file:///C/DocumentsandSettings/Administrator/桌面/代码/第20章/代码20-1.txtclear%清理工作空间loadwbarb;%装入图像figure;%新建窗口image(X);%显示图像colormap(map)%设置色彩索引图title('原始图像');%设置图像标题axissquare%设置显示比例disp('压缩前图像X的大小');%显示文字whos('X')%显示图像属性%对图像用小波进行层小波分解[c,s]=wavedec2(X,2,'bior3.7');%提取小波分解结构中的一层的低频系数和高频系数cal=appcoef2(
2、c,s,'bior3.7',1);%水平方向ch1=detcoef2('h',c,s,1);%垂直方向cv1=detcoef2('v',c,s,1);%斜线方向cd1=detcoef2('d',c,s,1);%各频率成份重构a1=wrcoef2('a',c,s,'bior3.7',1);h1=wrcoef2('h',c,s,'bior3.7',1);v1=wrcoef2('v',c,s,'bior3.7',1);d1=wrcoef2('d',c,s,'bior3.7',1);c1=[a1,h1;v1,d1];%显示分频信息figure;%新建窗口image(c1);%显示图像colorm
3、ap(jet)%设置色彩索引图axissquare;%设置显示比例title('分解后低频和高频信息');%设置图像标题ca1=appcoef2(c,s,'bior3.7',1);ca1=wcodemat(ca1,440,'mat',0);%改变图像高度并显示ca1=0.5ca1;figure;%新建窗口image(ca1);%显示图像colormap(map);%设置色彩索引图axissquare;%设置显示比例title('第一次压缩图像');%设置图像标题disp('第一次压缩图像的大小为:');%显示文字whos('ca1')%显示图像属性%保留小波分解第二层低频信息进行压缩ca
4、2=appcoef2(c,s,'bior3.7',2);%首先对第二层信息进行量化编码ca2=wcodemat(ca2,440,'mat',0);file:///C/DocumentsandSettings/Administrator/桌面/代码/第20章/代码20-1.txt(第1/2页)[2011-5-1614:36:27]file:///C/DocumentsandSettings/Administrator/桌面/代码/第20章/代码20-1.txt%改变图像高度并显示ca2=0.25ca2;figure;%新建窗口image(ca2);%显示图像colormap(map);%设
5、置色彩索引图axissquare;%设置显示比例title('第二次压缩图像');%设置图像标题disp('第二次压缩图像的大小为:');%显示文字whos('ca2')%显示图像属性file:///C/DocumentsandSettings/Administrator/桌面/代码/第20章/代码20-1.txt(第2/2页)[2011-5-1614:36:27]file:///C/DocumentsandSettings/Administrator/桌面/代码/第20章/代码20-2.txtclear;%清理工作空间loadwbarb;%装入图形信号figure;%新建窗口image
6、(X);%显示图像colormap(map);%设置色彩索引图axissquare;%设置显示比例title('原始图像');%设置图像标题disp('压缩前图像的大小');%显示文字whos('X');%显示图像属性axissquare;%设置显示比例%对图像用db3小波进行二层小波分解[c,s]=wavedec2(X,5,'db3');[thr,sorh,keepapp]=ddencmp('cmp','wv',X);[Xcomp,cxc,lxc,perf0,perfl2]=wdencmp('gbl',c,s,'db3',5,thr,sorh,keepapp);figure;%新建窗口
7、image(Xcomp);%显示图像colormap(map);%设置色彩索引图axissquare;%设置显示比例title('压缩后的图像');%设置图像标题disp('压缩后图像的大小');%显示文字whos('Xcomp')%显示图像属性disp('小波分解系数中值为0的系数个数百分比');%显示文字disp(perf0);%显示数值disp('压缩后剩余能量百分比');%显示文字disp(perfl2);%显示数值file