基于Matlab的小波分析在图像处理中的应用资料讲解.doc

《基于Matlab的小波分析在图像处理中的应用资料讲解.doc》由会员分享，可在线阅读，更多相关《基于Matlab的小波分析在图像处理中的应用资料讲解.doc（59页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、Good is good, but better carries it.精益求精，善益求善。基于Matlab的小波分析在图像处理中的应用-基于Matlab的小波分析在图像处理中的应用摘要:本文先介绍了小波分析得基本理论，包括连续小波变换、离散小波变换和小波包分析。小波变换具有时频局部化的特点,因此不但能对图像提供较精确的时域定位,也能提供较精确的频域定位。经过小波变换的图像具有频谱划、方向选择、多分辨率分析和天然塔式数据结构特点。基于小波变换这些特性，讨论了MATLAB语言环境下图像压缩，图像去噪，图像融合，图像分解，图像增强的基本方法。关键词:小波分析；图像压缩；图像去噪；图像融合；图像分解

2、；图像增强1引言小波分析诞生于20世纪80年代,被认为是调和分析即现代Fourier分析发展的一个崭新阶段。众多高新技术以数学为基础,而小波分析被誉为“数学显微镜”,这就决定了它在高科技研究领域重要的地位。目前,它在模式识别、图像处理、语音处理、故障诊断、地球物理勘探、分形理论、空气动力学与流体力学上的应用都得到了广泛深入的研究,甚至在金融、证券、股票等社会科学方面都有小波分析的应用研究。在传统的傅立叶分析中，信号完全是在频域展开的，不包含任何时频的信息，这对于某些应用来说是很恰当的，因为信号的频率的信息对其是非常重要的。但其丢弃的时域信息可能对某些应用同样非常重要，所以人们对傅立叶分析进行了

3、推广，提出了很多能表征时域和频域信息的信号分析方法，如短时傅立叶变换，Gabor变换，时频分析，小波变换等。其中短时傅立叶变换是在傅立叶分析基础上引入时域信息的最初尝试，其基本假定在于在一定的时间窗内信号是平稳的，那么通过分割时间窗，在每个时间窗内把信号展开到频域就可以获得局部的频域信息，但是它的时域区分度只能依赖于大小不变的时间窗，对某些瞬态信号来说还是粒度太大。换言之，短时傅立叶分析只能在一个分辨率上进行。所以对很多应用来说不够精确，存在很大的缺陷。而小波分析则克服了短时傅立叶变换在单分辨率上的缺陷，具有多分辨率分析的特点，在时域和频域都有表征信号局部信息的能力，时间窗和频率窗都可以根据信

4、号的具体形态动态调整，在一般情况下，在低频部分（信号较平稳）可以采用较低的时间分辨率，而提高频率的分辨率，在高频情况下（频率变化不大）可以用较低的频率分辨率来换取精确的时间定位。本文介绍了小波变换的基本理论，并介绍了一些常用的小波函数，它们的主要性质包括紧支集长度、滤波器长度、对称性、消失矩等，都做了简要的说明。然后研究了小波分析在图像处理中的应用，包括图像压缩，图像去噪，图像融合，图像分解，图像增强等。2小波分析的基本理论2.1连续小波变换定义：设，其傅立叶变换为，当满足允许条件（完全重构条件或恒等分辨条件）350)c(i)=2*c(i);elsec(i)=0.5*c(i);endend%下

5、面对处理后的系数进行重构xx=waverec2(c,s,sym4);%画出重构后的图像subplot(122);image(xx);colormap(map);title(增强图像);axissquare输出结果如图所示：图10小波分析用于图像增强本文将主要讨论图像增强中的钝化和锐化两种方法，钝化操作主要是提出图像中的低频成分，抑制尖锐的快速变化成分，锐化操作正好相反，将图像中尖锐的部分进可能得提取出来，用于检测和识别等领域。下面将以例子说明这两种方法在Matlab中的实现，并对于基于傅立叶变换的传统频域方法同小波方法做一下比较。4.3.2图像钝化图像钝化在时域中的处理相对简单，只需要对图像作

6、用一个平滑滤波器，使得图像中的每个点与其相邻点做平滑即可，这里不做详细介绍，我们来介绍一下基于傅立叶变换的频域处理方法。下面我们以chess信号为例，通过两种方法对图像钝化的结果做一下比较。loadchess分别保存用DCT方法和小波方法的变换系数blur1=X;blur2=X;对原图像做二维离散余弦变换ff1=dct2(X);对变换结果在频域做BUTTERWORTH滤波fori=1:256forj=1:256ff1(i,j)=ff1(i,j)/(1+(i*j+j*j)/8192)2);endend重建变换后的图像blur1=idct2(ff1);对图像做2层的二维小波分解c,l=wavede

7、c2(X,2,db3);csize=size(c);对低频系数进行放大处理，并抑制高频系数fori=1:csize(2);if(c(i)300)c(i)=c(i)*2;elsec(i)=c(i)/2;endend通过处理后的小波系数重建图像blur2=waverec2(c,l,db3);显示三幅图像subplot(221);image(wcodemat(X,192);colormap(gray(256);title(原始图像,fontsize,18);subplot(223);image(wcodemat(blur1,192);colormap(gray(256);title(采用DCT方法钝

8、化图像,fontsize,18);subplot(224);image(wcodemat(blur2,192);colormap(gray(256);title(采用小波方法钝化图像,fontsize,18);图11图像钝化从图中可以看出，采用DCT在频域做滤波的方法得到钝化结果更为平滑，这是因为其分辨率最高，而小波方法得到的结果在很多地方有不连续的现象，因为我们对系数做放大或抑制在阈值两侧有间断，而且分解层数很低，没有完全分离出频域的信息。而且我们在做系数放大或抑制的时候，采用的标准是根据系数绝对值的大小，没有完全体现出其位置信息，但是在小波系数中，我们很容易在处理系数的过程中加入位置信息。

9、4.3.3图像锐化与图像钝化所做的工作相反，图像锐化的任务是突出高频信息，抑制低频信息，从快速变化的成分中分离出标识系统特性或区分子系统边界的成分，以便于进一步的识别、分割等操作。在时域（空域）中，锐化的方法不外乎是作用掩码或做差分，同钝化的道理一样，无论是掩码和差分都很难识别点之间的关联信息，我们下面的例子同样是在频域完成的，用传统的傅立叶分析方法（这里采用的是DCT变换）得到的频域系数。loadchess;分别保存用DCT方法和小波方法的变换系数blur1=X;blur2=X;对原图像做二维离散余弦变换ff1=dct2(X);对变换结果在频域做BUTTERWORTH滤波fori=1:256

10、forj=1:256ff1(i,j)=ff1(i,j)/(1+(32768/(i*i+j*j)2);endend重建变换后的图像blur1=idct2(ff1);对图像做2层的二维小波分解c,l=wavedec2(X,2,db3);csize=size(c);对低频系数进行放大处理，并抑制高频系数fori=1:csize(2);if(abs(c(i)300)c(i)=c(i)*2;elsec(i)=c(i)/2;endend通过处理后的小波系数重建图像blur2=waverec2(c,l,db3);显示三幅图像subplot(221);image(wcodemat(X,192);colormap(gray(256);title(原始图像,fontsize,18);subplot(223);image(wcodemat(blur1,192);colormap(gray(256);title(采用DCT方法锐化图像