1、小波变换及其在 图像处理中的典型应用,赵丹培 图像处理中心 2013年9月,目 录,8.1 从傅里叶变换到小波变换的时频分析法 8.2 小波变换分类 8.3 小波变换的多分辨分析特性 8.4 尺度函数与小波 8.5 小波变换的快速实现 8.6 图像的多分辨分解与重建 8.7 小波变换在图像边缘检测中的应用 8.8 小波变换在图像去噪中的应用 8.9 小波变换在图像融合中的应用,Fourier变换一直是信号处理领域中应用最广泛、效果最好的一种分析手段，是时域到频域互相转化的工具，从物理意义上讲，傅里叶变换的实质是把对原函数的研究转化为对其傅里叶变换的研究。但是傅里叶变换只能提供信号在整个时间域上
2、的频率，不能提供信号在某个局部时间段上的频率信息。,8.1 从傅里叶变换到小波变换的 时频分析法,8.1.1 傅里叶变换,8.1.1 傅里叶变换,傅里叶变换：对于时域的常量函数，在频域将表现为冲击函数，表明具有很好的频域局部化性质。,傅里叶变换,反傅里叶变换,8.1.1 傅里叶变换,x=sin(2*pi*50*t)+sin(2*pi*300*t);%产生50HZ和300HZ的信号 f=x+3.5*randn(1,length(t);%在信号中加入白噪声,时间,由于傅立叶变换无法作局部分析，为此，人们提出了短时傅里叶变换（STFT）的概念，即窗口傅里叶变换。短时傅里叶变换将整个时间域分割成一些小
3、的等时间间隔，然后在每个时间段上用傅里叶分析，它在一定程度上包含了时间频率信息，但由于时间间隔不能调整，因而难以检测持续时间很短、频率很高的脉冲信号的发生时刻。,8.1.2 短时傅里叶变换,8.1.2 短时傅里叶变换,基本思想是：把信号划分成许多小的时间间隔，用傅立叶变换分析每一个时间间隔，以便确定该时间间隔存在的频率。 STFT的处理方法是对信号施加一个滑动窗(反映滑动窗的位置)后，再作傅立叶变换。即：,8.1.2 短时傅里叶变换,8.1.2 短时傅里叶变换,短时傅里叶变换的分析特点 (a)频率变化的影响 (b) 基本分析单元的特点,小波起源： 1984年Morlet提出;1985年Meye
4、r构造出小波;1988年，Daubechies证明了离散小波的存在;1989年，Mallat提出多分辨分析和二进小波变换的快速算法;1989年Coifman、 Meyer引入小波包;1990年崔锦泰等构造出样条单正交小波基;1994年Sweldens提出二代小波提升格式小波(Lifting Scheme)。 小波定义： “小”是指在时域具有紧支集或近似紧支集，“波”是指具有正负交替的波动性，直流分量为0。 小波概念：是定义在有限间隔而且其平均值为零的一种函数。,8.1.3 小波变换,持续宽度相同,振荡波,波与小波的差异：,用镜头观察目标 (待分析信号)。代表镜头所起的作用(如滤波或卷积)。相当
5、于使镜头相对于目标平行移动。的作用相当于镜头向目标推进或远离。,小波变换的粗略解释,8.1.4 小波变换的时频分析,尺度a较大,距离远 视野宽,概貌观察,尺度a较小,距离近 视野窄,细节观察,分析 频率低,分析 频率高,由粗到精,多分辨分析,品质因数保持不变,小波变换的时频分析特点：,小波变换的分析特点 (a) 尺度a不同时时域的变化 (b)尺度a不同时频域的变化,小波变换的多分辨分析特性：,不同a值下小波分析区间的变化,不同a值下分析小波频率范围的变化,频窗,时窗,小波变换的时频局部特性：,8.1.5 连续小波变换,尺度因子 的作用是将基本小波 做伸缩， 越大 越宽。,小波的位移与伸缩,设
6、，当 满足允许条件时：,8.1.5 连续小波变换,称 为一个“基小波”或“母小波”。 小波变换的含义是： 把基本小波(母小波)的函数 作位移后，再在不同尺度下与待 分析信号作内积，就可以得到一个小波序列。,连续情况时，小波序列为： (基本小波的位移与尺度伸缩)其中 为尺度参量， 为平移参量。 离散的情况，小波序列为 ：,根据容许条件要求，当=0时，为使被积函数是有效值，必须有 ，所以可得到上式的等价条件为：此式表明 中不含直流，只含有交流，即具有震荡性，故称为“波”，为了使 具有局部性，即在有限的区间之外很快衰减为零，还必须加上一个衰减条件：,衰减条件要求小波具有局部性，这种局部性称为“小”，
7、所以称 为小波。 对于任意的函数 的连续小波变换定义为：逆变换为： 是尺度因子， 反映位移。,线性 设： 平移不变性 若 ，则 伸缩共变性如果 的CWT是 则 的CWT是 冗余性(自相似性) 由连续小波变换恢复原信号的重构公式不是唯一的,8.1.6 连续小波的性质,目 录,8.1 从傅里叶变换到小波变换的时频分析法 8.2 小波变换分类 8.3 小波变换的多分辨分析特性 8.4 尺度函数与小波 8.5 小波变换的快速实现 8.6 图像的多分辨分解与重建 8.7 小波变换在图像边缘检测中的应用 8.8 小波变换在图像去噪中的应用 8.9 小波变换在图像融合中的应用,8.2 小波变换分类,小波函数
8、中 三个变量均为连续变量，称为连续小波。可以对 三个变量施加不同的离散化条件，并相应地对小波及小波变换进行分类。其中，最重要的两种分类：离散小波及离散小波变换二进小波及二进小波变换,8.2.1 离散小波变换,如果设定 ，则 对于任意函数 ，定义相应的离散小波变换为：如果这时 构成空间 的一组规范正交基，对于任一函数 的反演式为一展开式：,8.2.2 二进小波及二进小波变换,在连续小波变换中，令参数 而参数 仍取连续值，则有二进小波：这时， 的二进小波变换定义为,目 录,8.1 从傅里叶变换到小波变换的时频分析法 8.2 小波变换分类 8.3 小波变换的多分辨分析特性 8.4 尺度函数与小波 8
9、.5 小波变换的快速实现 8.6 图像的多分辨分解与重建 8.7 小波变换在图像边缘检测中的应用 8.8 小波变换在图像去噪中的应用 8.9 小波变换在图像融合中的应用,多分辨分析是小波分析中最重要的概念之一，它将一个函数表示为一个低频成分与不同分辨率下的高频成分，并且多分辨分析能提供一种构造小波的统一框架，提供函数分解与重构的快速算法。由理想滤波器引入多分辨率分析的概念：,8.3 小波变换的多分辨分析特性,多分辨分析定义：,空间 中的一系列闭子空间 ，称为 的多分辨率分析或逼近，若下列条件满足： 单调性： ，对任意 逼近性： 伸缩性： 平移不变性： Riesz基：存在 ，使 构成 的Ries
10、z基，即是线性无关的，且存在常数 与 ，满足 使得对任意的 ，总存在序列 使得 且 ，称 为尺度函数，并称 生成 的一个多分辨分析 。,是一个无限维向量空间，称为平方可积空间，将 用它的子空间 ， 表示，其中 称为尺度空间， 称为小波空间。尺度空间的递归嵌套关系：小波空间 是 和 之间的差，即 ，它捕捉由 逼近 时丢失的信息。推出：,多分辨率的空间关系图,目 录,8.1 从傅里叶变换到小波变换的时频分析法 8.2 小波变换分类 8.3 小波变换的多分辨分析特性 8.4 尺度函数与小波 8.5 小波变换的实现 8.6 图像的多分辨分解与重建 8.7 小波变换在图像边缘检测中的应用 8.8 小波变
11、换在图像去噪中的应用 8.9 小波变换在图像融合中的应用,两尺度方程,若 是尺度函数，它生成 的多分辨分析 ，则必然存在系数序列 ，使得以下尺度关系成立:这就是两尺度方程，必须满足下列条件：定义函数 为尺度函数，若其经过整数平移 和 尺度 上的伸缩，得到一个尺度和位移均可变化的函数集 合：,和 的基本性质是两尺度差分方程：两尺度方程的频域表示为：,目 录,8.1 从傅里叶变换到小波变换的时频分析法 8.2 小波变换分类 8.3 小波变换的多分辨分析特性 8.4 尺度函数与小波 8.5 小波变换的快速实现 8.6 图像的多分辨分解与重建 8.7 小波变换在图像边缘检测中的应用 8.8 小波变换在
12、图像去噪中的应用 8.9 小波变换在图像融合中的应用,8.5.1 Mallat算法与塔式分解,系数分解的快速算法：,Mallat小波快速分解算法的流程图,系数重构的快速算法：,Mallat小波快速重构算法的流程图,目 录,8.1 从傅里叶变换到小波变换的时频分析法 8.2 小波变换分类 8.3 小波变换的多分辨分析特性 8.4 尺度函数与小波 8.5 小波变换的快速实现 8.6 图像的多分辨分解与重建 8.7 小波变换在图像边缘检测中的应用 8.8 小波变换在图像去噪中的应用 8.9 小波变换在图像融合中的应用,8.6.1 二维小波变换的实现,假定二维尺度函数可分离，则有 其中 、 是两个一维
13、尺度函数。若 是相应的小波，那么下列三个二维基本小波：与 一起就建立了二维小波变换的基础。,8.6.2 图像小波变换的正变换,正变换图像小波分解的正变换可以依据二维小波变换按如下方式扩展，在变换的每一层次，图像都被分解为4个四分之一大小的图像。,8.6.3 图像小波变换的逆变换,逆变换 在每一层(如最后一层)都通过在每一列的左边插入一列零来增频采样前一层的4个阵列(即4个分解图像)； 接着用重构低通滤波器h和重构高通滤波器g来卷积各行，再成对地把这几个的阵列加起来； 然后通过在每行上面再插入一行零来将刚才所得两个阵列(图像)的大小增频采样为NN； 再用h和g与这两个阵列的每列进行卷积。这两个阵
14、列的和就是这一层次重建的结果。,对于二维图像信号，在每一层分解中，由原始图像信号与一个小波基函数的内积后再经过在x和y方向的二倍间隔抽样而生成四个分解图像信号。对于第一个层次(j=1)可写成：,8.6.4 二维小波变换的Mallat算法,将上式内积改写成卷积形式，则得到离散小波变换的Mallat算法的通用公式：,二维小波变换Mallat算法的通用公式：,8.6.5 二维Mallat多分辨率分解与重构,图像的Mallat快速塔式分解实验,8.6.6 多孔算法,多孔算法的分解实验,目 录,8.1 从傅里叶变换到小波变换的时频分析法 8.2 小波变换分类 8.3 小波变换的多分辨分析特性 8.4 尺
15、度函数与小波 8.5 小波变换的实现 8.6 图像的多分辨分解与重建 8.7 小波变换在图像边缘检测中的应用 8.8 小波变换在图像去噪中的应用 8.9 小波变换在图像融合中的应用,边缘像素实质上是局部图像范围内灰度的急剧变化点(奇异点)，图像边缘就是二维图像中奇异点的集合。边缘点在频域表现为高频信号，而图像噪声也多为高频信号，这使得两者难以区分。边缘检测的目的就是既要将高频信号从图像中分离出来，又要区分边缘与噪声，准确地标定边缘的位置。,8.7 小波变换在图像边缘检测中的应用,传统的边缘检测方法 为什么要用小波来进行边缘检测？ 边缘检测的分类 局部模极大值点 Canny准则,小波多尺度局部模
16、极大值边缘检测的原理 小波多分辨率边缘检测的具体实现 小波函数的选取 自适应阈值的选取 利用边缘信息进行目标定位 仿真实验,8.7.1小波多尺度局部模极大值边缘检测的原理,假设 是二维平滑函数，且满足 可把它沿 两个方向上的导数作为基本小波：对于一幅图像 ，其小波变换为：,整个图像的二进小波变换即矢量：模值为 ：相角为 ：,8.7.2 小波多分辨率边缘检测的具体实现,搜寻模极大值：,噪声的滤除：,(1)阈值法硬阈值，软阈值，自适应阈值； (2)多分辨率分割利用模极大值在各个尺度的传播特性去除噪声；利用小波变换尺度间相关性去除噪声；平移不变量的小波去噪方法；,边缘跟踪算法的四个约束条件： a)“
17、方向不变性”原则； b)角度平滑条件； c)幅值最邻近条件； d)“互认”原则。,边缘点的链接：,三种小波多尺度局部模极大值 边缘检测方法的比较,方法一：小波变换模极大值用于边缘检测的原始方法 具体步骤： 利用多孔算法对原图像 进行保持图像大小不变的小波变换，生成水平分量 和垂直分量 ； 计算出梯度向量的模值 计算出梯度向量的相角按照相角量化方法划分为8个方向，求出不同方向的局部模极大值点； 对各尺度的边缘图像进行阈值处理； 链化模极大值点，去除长度小于一定阈值的边缘链，就得到各 个尺度上的边缘图像。,方法二：小波变换模极大值用于边缘检测的简化方法,具体步骤： 对所给图像的每一行进行小波变换，
18、求出 ； 对所给图像的每一列进行小波变换，求出 ； 计算出梯度向量的模值 计算出梯度向量的相角 将相角 按8方向量化，求出 在不同方向的局部模极大值点； 对各尺度的边缘图像进行阈值处理； 链化模极大值点，去除长度小于一定阈值的边缘链，就得到各个尺度上的边缘图像。,方法三：小波变换模极大值的多尺度综合方法,具体步骤： 求出各尺度的模图像簇 和相角图像簇 ； 对各尺度的边缘图像进行阈值处理； 将相角 按8方向量化，求出 在不同方向的局部模极大值点； 由粗到细的边缘链接：对经过阈值处理的最粗分辨率级上的模极大值开始，链接模极大值图像中模值相近，相角相似的非零像素点，删除长度小于链长阈值的边缘链，得到
19、大尺度下单像素宽的图像边缘 ；,针对尺度 的每一个边缘像素，搜索 尺度下以这点为中心的3X3邻域，将该邻域内出现的所有可能边缘点均标记为候选边缘点，得到 尺度下的候选边缘点图像 非候选边缘点标记为零； 将 尺度下的候选边缘点图像 中模值相近，相角相似的非零像素点链接，删除长度小于阈值的边缘链，得到 尺度下单像素宽的图像边缘 ； 重复步骤，直到 为止，边缘图像即为综合后形成的边缘图像，也就是我们最终所要得到的边缘图像。,8.7.3 小波函数的选取,在实际边缘检测中，小波函数的选取直接关系到边缘检测的结果，小波变换相当于对图像进行带通滤波，在一定程度上减少了噪声对图像的影响，但同时也去掉了一些模糊
20、边缘。这就要求寻找一种具有好的去噪特性同时又能精确地提取边缘的小波函数，而且这种小波在满足Canny准则的同时应具有限紧支撑特性。在关于边缘提取的文献中，主要采用的是正交小波，双正交小波和B样条小波作为边缘提取的小波函数。通常选取高斯函数作为平滑函数。三次B样条函数已很接近高斯函数，能满足大多数的应用，而它的紧支性优于高斯函数。在有噪声环境下进行边缘检测时，需要平衡噪声抑制能力与边缘定位精度，此时三次B样条函数是最优的。,其中， 为阈值初值， 为比例系数，N为采样点 数。根据 和 的大小来决定窗口内均值对阈值 的作用。,8.7.4 自适应阈值选取,选择一个窗口在小波系数上滑动，窗口大小可以 取
21、3232或1616，将尺度下小于阈值 的梯度 值置为0，自适应阈值为：,8.7.5 利用边缘信息进行目标定位,对经过综合后得出的图像边缘，根据形心公式就可以计算出图像的形心坐标，判断出目标在视场中的位置，实现目标定位。,8.7.6 实验结果（一）,五个尺度的模极大值提取,第一个尺度的噪声很多,多尺度链接后的模极大值,实验二:阈值的选取对边缘检测结果影响的实验,用固定阈值在一阶尺度提取的边缘,用自适应阈值在一阶尺度提取的边缘,用固定阈值在二阶尺度提取的边缘,用自适应阈值在二阶尺度提取的边缘,实验三：三种模极大值边缘检测方法的性能比较,实验四：基于三次B样条小波的边缘检测实验,由于小波变换具有多分
22、辨分析特性和时频局部化能力，在边缘检测、去噪和图像增强等方面都具有很强的优势； 更适合用来检测受噪声污染严重的模糊图像和低对比度图像，尤其对微弱目标，它首先能抑制噪声、增强对比度，然后利用多尺度的模极大值方法有效检测出目标边缘，从而实现3对比度下的目标精确定位； 这种算法要进行多尺度运算，所以计算量很大。,小 结,目 录,8.1 从傅里叶变换到小波变换的时频分析法 8.2 小波变换分类 8.3 小波变换的多分辨分析特性 8.4 尺度函数与小波 8.5 小波变换的实现 8.6 图像的多分辨分解与重建 8.7 小波变换在图像边缘检测中的应用 8.8 小波变换在图像去噪与增强中的应用 8.9 小波变
23、换在图像融合中的应用,8.8 小波变换在图像去噪与增强中的应用,高频系数置零的线性去噪方法 小波系数硬阈值去噪方法 小波系数软阈值去噪方法 小波系数自适应阈值去噪方法 基于小波模极大值的去噪方法 基于信号奇异性的去噪方法 基于小波系数相关性的去噪方法 基于小波变换的图像增强方法,8.8.1 高频系数置零去噪,1.对噪声图像进行二维离散小波分解，分解层数一般取2或3 层，分解过程如下图。 2.对每一层的高频系数LH，HL，HH，置零。 3.对小波系数进行重构。,这是最简单的利用小波变换性质的去噪方法。,8.8.2 基于小波硬阈值的去噪方法,1.首先将图像信号求小波变换。 2.除了最粗尺度信号外，
24、将各细节信号作阈值处理，当某位置小波变换值大于阈值时，保留原值，否则置零。即,3.利用小波变换重构，求出信号的滤波值。,8.8.3 基于小波软阈值的去噪方法,1.首先对图像信号进行小波变换，得出带有噪声的小波系数。 2.将各细节信号作阈值处理，当某位置小波变换值大于阈值时，作下面运算，下式中sgn(x)代表符号函数,否则置零。 3.利用小波变换重构，求出信号的滤波值。,软硬阈值滤波器,8.8.4 非线性软阈值去噪方法,其中，噪声方差的估计为 ， MAD 为图像中位值，n为信号采样点数。,阈值选取,1.2.3.4.,对含有高斯噪声的Lena图像利用硬阈值法、软阈值法和自适应阈值法去噪后实验结果,
25、基于小波变换阈值方法去噪的不足, 只适合高斯白噪声的去噪，对椒盐噪声效果不明显。原因在于阈值法保留的是大于阈值的小波系数，而椒盐噪声在图像上表现为或者是灰度值特别大的白像素，或者特别小的黑像素，椒盐噪声的小波系数都很大，所以不能用阈值分离出来。 阈值方法的去噪效果依赖于信噪比的大小，它特别适合信噪比高的图像去噪。 在图像信号不连续点处会有伪吉普斯现象。 阈值方法的关键在于阈值的选取，而选择一种普适性很好的阈值选取方法是很困难的，事实上，人们已经证明在均方误差意义上阈值方法能得到原始图像信号的最优估计，然而在实际应用中还是需要根据具体的情况和经验来对一些阈值进行改进。,8.8.5 基于多尺度模极
26、大值的小波去噪方法,根据图像边缘点与噪声点具有不同的奇异性，则它们小波变换后的模极大值在不同尺度下的传播特性也不同这一特性进行图像去噪。当信号 在 处Lipschitz指数为 时， 反映了函数在点 处的奇异性大小。是当信号在 处的奇异性 时，表明小波系数模极大值将随尺度 的增大而增大；当 时，则随 的增大而减小。通常的图像具有时域上的相关性，因而 ；而噪声由于在该点不可导， 。因此，边缘点的模值和噪声点的模值随尺度的变化具有不同的规律，所以可以利用这一特性将噪声分离出来。,基于多尺度模极大值的小波去噪方法 的具体实现,通过分析图像和噪声在小波域中对应的系数模极大值在不同尺度上的分布情况，来研究
27、图像和噪声的突变点在不同尺度上的传递特性。在从低到高的分解尺度中，图像突变点对应的小波系数极大值具有传递性，而噪声突变点不具有这种传递性。根据这一性质，在确定出各尺度小波系数极大值的基础上，由粗到精的跟踪不同尺度上的小波系数极大值，并依据其不同尺度间的传递性，识别信噪属性，剔除噪声部分对应的小波系数极大值，从而抑制噪声，提高图像质量。,模极大值方法去噪的过程,1.小波系数极大值的确定 2.图像极大值跟踪 3.噪声极大值滤除,模极大值多尺度链接方法对高斯噪声 图像去噪的实验结果,对含有高斯噪声（ 0.01）的Lena图像利用多尺度模极大值链接方法去噪后图像,模极大值多尺度链接方法对椒盐噪声 图像
28、去噪的实验结果,对含有椒盐噪声（ 0.005）的图像利用多尺度模极大值链接方法去噪后图像,含有椒盐噪声的图像、均值滤波去噪图像、中值滤波去噪图像、硬阈值去噪图像、软阈值去噪图像和模极大值奇异性去噪图像的实验结果对比,基于多尺度模极大值奇异性对椒盐噪声的去噪实验结果,模极大值多尺度链接方法对混和噪声图像的去噪实验,含有高斯和椒盐混合噪声图像、硬阈值去噪图像、软阈值去噪图像、自适应阈值去噪图像、中值滤波去噪图像和模极大值多尺度链接去噪图像,8.8.6 基于小波系数相关性的去噪方法,信号与噪声的小波变换在各尺度下的不同传播特性表明，信号的小波变换在各尺度间有较强的相关性，而且在边缘处具有很强的相关性
29、；而噪声的小波变换在各尺度间确没有明显的相关性，而且噪声的小波变换主要集中在小尺度各层次中。根据信号与噪声的小波变换在不同尺度间的上述不同特点，可以通过将相邻尺度的小波系数直接相乘来增强信号，抑制噪声。由于噪声主要分布在小尺度上，所以这种现象在小尺度上非常明显。,首先定义 为尺度 上 点处的 相关系数。为使相关系数与小波系数具有可比性，定义 规范化相关系数：其中和 分别表示对应于尺度的小波系数与相关系数的能 量。显然，在尺度 下，小波系数与规范化相关系数具有相同 的能量，这为它们之间提供了可比性。,基于小波系数相关法去噪的实验结果（一）,（ 0.01）Lena图像和基于小波系数相关法去噪后的图
30、像,小 结,对小波系数使用阈值的去噪方法速度快，实现简单，所以一直被广泛使用在高斯噪声去除上，其中的硬阈值方法容易产生伪吉普斯效应，软阈值方法克服了这个缺点，图像比较平滑，去噪效果也比硬阈值好，而自适应阈值使用的是局部阈值，比软、硬阈值都好。 模极大值多尺度链接和小波系数相关去噪两种方法虽然速度都比较慢，但是它们都比阈值方法效果好，对各种噪声效果都比较稳定，特别是对混杂噪声的去除方面。但是在椒盐噪声图像的去噪方面，只有模极大值链接的方法有效果，而且比传统的中值滤波方法要保留更多的边缘信息。在对同时带有高斯噪声和椒盐噪声的图像去噪时，模极大值链接的方法去噪效果非常突出，而且稳定。,8.8.7 基
31、于小波的图像增强方法,小波多分辨分析由于它能多尺度多角度提取信号特征，并在不同尺度上让噪声和信号明显地区分开来，所以它在图像去噪和增强方面有很大优势。基于小波多分辨分析的图像增强，就是突出图像的边缘细节，尽可能的消除负面因素，从而达到增强图像的目的。因此，基于小波分析的图像增强可以转化为两步。第一步是区分噪声和图像的边缘细节；第二步是根据需要对图像的边缘细节适度增强。通过一个合适的增益函数完成上面所述的两个步骤，使去噪和增强一次完成。在增强图像的同时应考虑到噪声问题，所以对小波分解后高频部分的处理很重要。为此，设计合适的增益函数是关键的一步。,简单的线性图像增强的方法,第一类是简单的增益函数，
32、即对每个方向的高频采用不同的常数 ； 第二类是比较复杂的增益函数，有分段线性函数和非线性函数。根据处理高频增益函数的不同方法，将图像增强分为线性增强和非线性增强 。,线性图像增强方法的实现,线性增强是对小波分解的高频采取的常数增益，具体的算法是： a、对图像进行多级小波变换，得到 b、计算图像的原始噪声水平； C、计算图像在各个尺度上不同方向的噪声水平；d、去噪并增强。认为幅值小于 的系数为噪声，给予抑制，而幅值大于 的系数给予一定的增益。这样得到 e、对增强后的小波系数，进行小波逆变换得到增强后的图像。,基于边缘提取的去噪和图像增强相结合的实现,1.对图像进行小波变换 2.计算小波变换的多尺
33、度模和幅角； 3.用局部模极大值方法检测出边缘细节点； 4.利用合适的非线性增益函数增强边缘点，同时去除高频噪声点； 5.对增强后的小波系数进行小波逆变换得到增强后的图像。,基于小波变换的图像增强方法的实验结果,对Lena图像利用小波系数线性增强前后的效果对比,目 录,8.1 从傅里叶变换到小波变换的时频分析法 8.2 小波变换分类 8.3 小波变换的多分辨分析特性 8.4 尺度函数与小波 8.5 小波变换的实现 8.6 图像的多分辨分解与重建 8.7 小波变换在图像边缘检测中的应用 8.8 小波变换在图像去噪中的应用 8.9 小波变换在图像融合中的应用,8.9 小波变换在图像融合中的应用,利
34、用小波进行图像融合的基本思想是： 1、先对多源图像进行二维小波分解； 2、然后在小波变换域内通过比较各图像的细节信息，在不同尺度上实现图像融合，提取出重要的小波系数；(融合法则) 3、最后进行小波反变换，得到数据融合之后的图像。,8.9.1 基于小波的图像融合方法,8.9.2 基于小波图像融合的过程,基于小波变换的图像融合过程,8.9.3 像素级融合的基本原理,基于像素的融合规则,融合算子的确定低频融合算子高频融合算子,8.9.4 区域融合的基本原理,基于区域的融合规则,融合算子的确定低频融合算子高频融合算子,8.9.5 实现中涉及的的一些关键问题,1、小波基的选取 ： 通常在设计滤波器时，考
35、虑的主要方面有： 紧支撑性(Compact Support)：h(n)的取值范围是有限的，在该取值范围内，h(n)的值不为零，在该范围以外，h(n)的值为零。具有紧支撑的小波基收敛速度快。 规则性(Regularity)：即平滑性(Smoothness)或连续性(Continuity)。在图像融合中，若滤波器的规则性较好，则进行重构时的视觉效果较好。 对称性(Symmetry)：也称线性相位(Linear Phase)。h(n)在时域的对称，在频域上就表现为具有线性相位。,8.9.5 实现中涉及的的一些关键问题,2、边界延拓问题 ：实际的图像处理中，由于图像都是有限尺寸的，在把滤波器应用于边界
36、时，存在如何处理边界的问题。比较常用的方法有零延拓、边界重复延拓、周期延拓和对称周期延拓和双对称延拓等。主要目的是要降低边界不连续性所产生的在边界上变换系数衰减慢的问题。,8.9.5 实现中涉及的的一些关键问题,3、分解层数的选择 ： 小波分解层数J 的大小与所选用小波函数的光滑度有关，J 选择小，小波函数光滑度就差，但运算量小；J 选择大，小波函数光滑度就好，融合的频率范围越丰富，融合结果的细节也就越丰富。但是它的运算量大，而且分解的层数越多，顶层融合损失的信息量越大，均是小波反变换不能恢复的损失，因此基于小波分解的层数不宜过高。 实验证明选择J = 3 或J = 5 从融合图像质量来看几乎
37、没有差别，因此，我们选择J = 3 即可以满足要求。,8.9.6 基于小波的图像融合实验,实验一,蝴蝶眼睛序列显微图像,三种融合方法的融合效果比较,(a) 改进的拉普拉斯融合方法 (b) 基于小波像素级的融合方法,实验一,实验三,8.9.6 基于小波的图像融合实验,融合效果比较,实验三,(a)改进的拉普拉斯融合方法 (b) 小波区域融合方法,参考书目,小波分析及其应用(孙延奎) 基于MATLAB的系统分析与设计小波分析 小波变换的工程分析与应用(杨福生) 小波分析算法与应用(程正兴) A wavelet tour of signal processing Ten Lectures on Wavelets,
