《提取图像纹理特征程序设计小波变换》由会员分享，可在线阅读，更多相关《提取图像纹理特征程序设计小波变换（9页珍藏版）》请在金锄头文库上搜索。
1、沈阳理工大学数字图象处理综合课程设计1 1. 设计目的小波变换可以使得信号的低频长时特性和高频短时特性同时得到处理，具有良好的局部化性质， 能有效地克服傅氏变换在处理非平稳复杂信号时存在的局限性，具有极强的自适应性，因此在图像处理中具有极好应用价值。与Fourier变换相比，小波变换是空间( 时间 ) 和频率的局部变换，因而能有效地从信号中提取信息。通过伸缩和平移等运算功能可对函数或信号进行多尺度的细化分析，解决了Fourier变换不能解决的许多困难问题。本设计主要分析了基于小波变换的图像分解和图像压缩技术，并运用Matlab 软件对图像进行分解，然后提取其中与原图像近似的低频信息，达到对图像
2、进行压缩的目的。分别作第一层分解和第二层分解，并比较图像压缩的效果。2. 设计要求利用小波变换的基本原理在MATLAB 环境下编写程序对静态图像处理并观察分析其处理效果，通过伸缩和平移等运算功能可对函数或信号进行多尺度的细化分析，解决了Fourier变换不能解决的许多困难问题。并运用Matlab软件对图像进行分解， 然后提取其中与原图像近似的低频信息，达到对图像进行压缩的目的。分别作第一层分解和第二层分解，并比较图像压缩的效果。沈阳理工大学数字图象处理综合课程设计2 3设计原理3.1 小波变换的分解算法小波是函数空间L2(R) 中满足下述条件的一个函数或者信号(x) 这里，R*=R-0 表示非
3、零实数全体。对于任意的函数或者信号)(xf，其小波变换定义为dxabxxf adxxxfba RRbafW)()( |1)()(),(),(因此，对任意的函数)(xf，它的小波变换时一个二元函数。另所谓多分辨分析是指设Vj；Zj 是L2(R) 上的一列闭子空间，其中的一函数，如果它们满足如下五个条件，即（1）单调性：VVjj1，Zj；（2）唯一性： 0VIj Zj；（3）稠密性：VYj ZjL2(R) ；（4）伸缩性：VVjjxfxf1)2()(，Zj；（5） Riesz 基存在性：存在Vt0)(， 使得)(2nxi；Zn 构成Vj的 Riesz基。称)(t为尺度函数。那么称 Vj；Zj ，)
4、(t是L2(R) 上的一个多分辨分析。若定义函数)()(222 ,nxxijnj，Znj,，则由多分辨分析的定义，容易得到一个重要结果，即函数族)()(222 ,nxxijnj；Zn 是空间Vj的标准正交基。 关于多分辨分析， 在这里以一个三层的分解进行说明，其小波分解树如图2.1 所示（ A表示低频， D表示高频）。*2| )(|RdxC沈阳理工大学数字图象处理综合课程设计3 SA1A2A3D2D3D1小波分解法从图中可以明显看出， 多分辨分析只是对低频部分进行进一步分解，而高频部分则不予考虑。分解具有关系DDDAS1233，另外强调一点，这里只是以一个层分解进行说明，如果要进行进一步分解，
5、则可以把低频部分分解成低频部分和高频部分，以下再分解， 依次类推。 在理解多分辨分析时， 必须牢牢把握一点， 即分解的最终目的是力求构造一个在频率上高度逼近空间的正交小波基， 这些频率分辨率不同的正交小波基相当于带宽各异的带通滤波器。从图 2.1 可以看出，多分辨分析只对低频空间进行进一步的分解，使频率的分辨率变得越来越高。3.2 小波变换的重构算法设pjl, 、qijl, （i=1,2,3 ）是由两个一元两尺度序列得到的二元尺度序列，即pjl,=ppjl21，pjl1,=qpjl21,qjl2,=pqjl21,pjl3,=qqjl21。则有重构算法为cmnk,;1=dqcpijlk jijm
6、lnjlkjmlnjl,;312,2,;2,2,小波重构的数据传递示意图如图所示沈阳理工大学数字图象处理综合课程设计4 222222),(mnc)(1 nq)(2 mq)(1 np)(2 mp)(1 nq)(1 npLLLHHLHH小波变换在图像压缩中的应用二维离散小波变换后的系数分布构成了信号),(yxf的二维正交小波分解系数 ( 如图所示 ), S3W33 W32 W31W13W23W12W22W11W21二维正交小波分解系数它们每一个都可被看做一幅图像,),(1mnfWj给出了),(yxf垂直方向的高频分量的小波分解系数，),(3mnfWj给出了),(yxf水平方向的高频分量的小波分解系
7、数，),(2mnfWj给出了),(yxf对角方向高频分量的小波分解系数，fSJ 给出了),(yxf的低频分量的小波分解系数。由此可见，若用SJ，Wj1，Wj2，Wj3分别表示),(mnfSj，),(1mnfWj，),(2mnfWj，),(3mnfWj经 21 亚抽样后的变换系数 (简称为子图像 ) ，则任一图像都可以分解为j=-J ， -1 之间的 3J+1 个离散子图像：SJ，Wj1，Wj2，Wj3其中SJ是原图像的一个近似，Wij(i=1 ，2，3；j=-J ， -1) 则是图像在不同方向、不同分辨率下的细节；ZZjjjjjjmnfWmnfWmnfWmnfS mn),(),(),(),(3
8、21 1,.,),(沈阳理工大学数字图象处理综合课程设计5 GH2fSJ 12GHHG2222W11W21W31SJ如果原图像有 N2个像素，则子图像SJ，Wj1，Wj2，Wj3分别有Ni2个像素，因而分解后总的像素数NT为NNNNjjiJT2212 344可见，分解后总的像素数不变。二维数字信号也即数字图像， 对它的处理时基于图像的数字化来实现的。图像的书画结果就是一个巨大数字矩阵，图像处理就在这个矩阵上完成。所以，可将二维数字信号dmn看做),(0mnfS，即dmn=),(0mnfS=Rdxdymynxyxfmnyxyxf 2),(),(),)(,( 0*),(并采用与一维情况类似的mal
9、lat算法。由于两次一维小波变换来实现一次二维小波变换， 所以先对该矩阵的行进行小波变换，再对列进行小波变换。 从信号滤波的角度实现二维小波变换的框图分别如图3.4 所示。沈阳理工大学数字图象处理综合课程设计6 4. 设计内容4.1 代码下面的实例是基于二维小波分析对图像进行压缩。一个图像作小波分解后， 可得到一系列不同分辨率的子图像, 不同分辨率的子图像对应的频率是不相同的。高分辨率 (即高频 ) 子图像上大部分点都接近于0，越是高频这种现象越明显。对一个图像来说， 表现一个图像最主要的部分是低频部分，所以一个最简单的压缩方法是利用小波分解，去掉图像的高频部分而只保留低频部分。clear a
10、ll a=imread(D:2.JPEG); % 读取图片figure; imshow(a) X=rgb2gray(a); subplot(221);image(X);colormap(gray(256); title(原始图像 ); axis square; c,s=wavedec2(X,2,bior3.7); ca1=appcoef2(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,
11、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; subplot(222);image(c1); axis square title(分解后低频和高频信息 ); ca1=appcoef2(c,s,bior3.7,1); ca1=wcodemat(ca1,440,mat,0); 沈阳理工大学数字图象处理综合课程设计7 ca1=0.5*ca1; subplot(223);image(ca1);colormap(gray(256); title(第一次压缩图像 ); axis
12、 square ca2=appcoef2(c,s,bior3.7,2); ca2=wcodemat(ca2,440,mat,0); ca2=0.25*ca2; subplot(224);image(ca2);colormap(gray(256); axis square; title(第二次压缩图像 ); 4.2. 仿真结果4.3 结果分析从图中可以看出，第一次压缩我们是提取原始图像中小波分解第一层的低频 信息，此时压缩效果较好，压缩比较小( 约为 1/3) ；第二次压缩是提取第一层分 解低频部分的低频部分 ( 即小波分解第二层的低频部分) ，其压缩比比较大 (1/12) ，压缩效果在视觉上也
13、基本过得去，它不需要经过其它处理即可获得较好 的压缩效果。沈阳理工大学数字图象处理综合课程设计8 结束语这次设计利用小波变换完成了对静态图像进行压缩的目的，基本上实现了设计的要求，在这里对老师的指导和帮助表示感谢。图像压缩是一个很有发展前途的研究领域，它的研究就是寻找高压缩比的方法且压缩后的图像要有合适的信噪比， 在压缩传输后还要恢复原信号， 且在压缩、传输、恢复的过程中， 还要求图像的失真度小。 而将小波分析引入图像压缩的研究范畴，当一个图像作小波分解后， 可得到一系列不同分辨率的子图像，不同分辨率的子图像对应的频率是不相同的. 高分辨率子图像上大部分点的数值都接近0，越高就越明显。而对于一
14、个图像来说，表现一个图像的最主要部分是低频部分。而且小波分析能使压缩比高、压缩速度快， 压缩后能保持信号与图像的特征基本不变。在数字图像处理中具有很强的使用价值。沈阳理工大学数字图象处理综合课程设计9 参考文献1 阮秋琦 . 数字图像处理学 M. 北京：电子工业出版社， 2001. 2 龚声蓉，刘纯平，王强等 . 数字图像处理与分析 M. 北京：清华大学出版社，2006. 3 贾永红 . 计算机图像处理与分析 M. 武汉：武汉大学出版社， 2001. 4 陈桂明 . 应用 MATLAB 语言处理数字信号与图像处理M. 北京： 科学出版社，2000. 5 夏德深，傅德胜 . 计算机图像处理及应用 M. 南京：东南大学出版社， 2004. 6 姚敏. 计算机图像处理 M. 北京：机械工业出版社， 2006. 7 容观澳 . 计算机图像处理 M. 北京：清华大学出版社， 2000. 8 吴健康 . 数字图像分析 M. 北京：邮电出版社， 1989.
