图像算法五：【图像小波变换】多分辨率重构、Gabor滤波器、Haar小波

matlab设计：

与单纯运用某种自适应算法相比,基于小波分解的自适应滤波算法在收敛速度和稳定性上都有了很大的提高

小波变换：

示意：尺度越大，采用越大的时间窗，尺度越小，采用越短的时间窗，即尺度与频率成反比。在时频两域都具有表征信号局部特征的能力，其在低频部分具有较高的频率分辨率和较低的时间分辨率，在高频部分具有较高的时间分辨率和较低的频率分辨率，所以被誉为分析信号的数字显微镜。

小波变换与傅里叶变换对比的优点：短时傅立叶分析是把一个短时间的窗函数加在信号上，再对这一部分加窗的信号做傅立叶变换，当然，假定了在短时间内信号是平稳的。然后移动窗函数，对整个时间轴上的信号做短时傅立叶分析。但是，又有问题出现了，短时傅立叶分析一旦确定了窗函数，窗函数的形状便不再改变，于是短时傅立叶分析对于信号的时频分辨率也随之确定。

而我们希望稍微智能一点，即对于变化剧烈的信号，时间分辨率能高一点，频率分辨率可以降低，而对于变化平坦的信号，频率分辨率要高一些，时间分辨率可以低点。于是小波分析便产生了。

小波分析同样对信号在一短时间内做加窗（小波函数）分析，只是这个窗，即小波函数既可以在时间轴上移动，又可以伸缩，伸缩便意味着时频分辨率的改变。

小波变换的物理表现：分解的结果是产生长度减半的两个部分，一个是经低通滤波器产生的原始信号的平滑部分，另一个则是经高通滤波器产生的原始信号细节部分。

不同j所确定的频带是独立的，随j变化相互独立的频带覆盖了整个频率轴，分辨率j反应了频带的位置和带宽把信号分解到一系列相互独立的频带上分辨率j-1时的近似信号=f(x)在分辨率为j时的近似部分+细节部分

再对低频部分进行相似运算

3、小波变换在matlab中的使用

[C, S] = wavedec2(x, n, wname); % 对图像进行小波分解

小波变换的层数选择：通常分解层数过多，并且对所有的各层小波空间的系数都进行阈值处理会造成信号的信息丢失严重，去噪后的信噪比反而下降，同时导致运算量增大，使处理变慢.

分解层数过少则去噪效果不理想，信噪比提高不多，但不会出现信噪比下降的情况.

小波函数的选取：按常理选取了sym4

[C,S]的含义：C=[A(n)|H(n)|D(n)|………H(1)|V(1)|D(1)]

h ,v ,d 分别反映水平、垂直、对角线方向

S(1,:)=size of approximation coefficients(n)

部分程序运行：原始图像→加噪图像（同时加入高斯噪声和椒盐噪声）→小波分解后的图像显示

.M=imread('detfinger1.png');

%读取MATLAB中的名为detfinger的图像分辨率：512×512

subplot(1,2,1)

imshow(M); %显示原始图像

title('原始图像')

P1=imnoise(M,'gaussian',0.02);

%加入高斯躁声,方差为0.02的高斯噪声

P2=imnoise(P1,'salt &pepper',0.02);

%同时加入高斯噪声和椒盐躁声

subplot(1,2,2)

imshow(P2);%加入椒盐躁声后显示图像

title('加入高斯椒盐躁声后');

%%对加入高斯噪声的图像进行小波分解

P2=double(P2);

[CA,CH,CV,CD]= dwt2(P2,'db1 ');

y =[CA,CH;CV,CD];

y=uint8(y);

subplot(3,1,3);

imshow(y);

title('小波分解');

4、滤波器的特性

我的思路：先了解维纳滤波器的原理再进一步改进，理解自适应滤波器，进行自适应滤波器的设计，再进行自适应滤波器的组合设计，从而实现去噪功能

维纳滤波器：根据平稳随机信号的全部过去和当前的观察数据来估计信号的当前值，在最小均方误差的条件下得到系统的传递函数，参数是固定的，适用于平稳随机信号。

卡尔曼滤波器：根据当前时刻数据的观测值和前一时刻对该一时刻的预测值进行递推数据。它自动调节本身的冲击响应特性（自动调节数字滤波器的系数），以适应信号变化的特性，从而达到最优化滤波。它的参数是时变的，适用于非平稳随机信号。

这两种滤波器最优滤波的条件：噪声的统计特性先验已知

但实际应用中，常常无法得到这些统计特征的先验知识（统计特性是随时间变化的），因此实现不了最优滤波。

自适应滤波器：无法得到一些统计特性的先验知识（统计特性随时间变化）

而自适应滤波器在输入信号统计特性未知或者统计特性变化时，能够自动调节自身的参数，使其按照某种准则达到最优滤波。由于自适应滤波器具有这种特性，自提出以来，在实际工程中的众多领域得到了广泛应用。

5、如何设计和建立自适应滤波器

自适应滤波器工作原理是:系统能够按照某种算法自动调节权系数,使其实际输出和期望输出的均方误差达到最小值。

自适应滤波器的结构：有FIR和IIR两种。由于IIR滤波器存在稳定性的问题,因此一般采用FIR滤波器。由于FIR滤波器横向结构的算法具有容易实现和计算量少等优点,在对收敛速度不是很快的场合,多采用FIR作为自适应滤波器结构。

自适应滤波算法：

自适应滤波器的核心部分就是自适应滤波算法，基于不同的准则，最常用的自适应滤波算法是LMS算法和RLS算法，LMS算法以最小均方误差为准则,而RLS算法(递归最小二乘算法)以最小误差平方和为准则。

LMS算法由于其具有计算量小,稳定性好,易于实现等优点,被广泛使用。

LMS算法中均方误差表示为:

使均方误差达到最小值时，得到最佳权系数w*，应满足下列方程:在实际中，LMS 迭代算法是以最速下降法为原则进行的，用公式表示为:

式中，u 是自适应步长; ( n) 为迭代的梯度。

这里来估计均方误差 MSE 的梯度，并以此梯度估计值代替最速下降法中理论情况下的梯度真值，LMS( 最小均方) 算法进行梯度估计的方法是以误差信号的每一次迭代的瞬间时平方值代替其均方值，并以此来估计梯度。

在自适应滤波器的设计过程中，滤波器阶数 N和步长μ 等参数的选取对仿真结果的影响总结如下.

（1）自适应算法中步长 μ 的选择：μ 的大小影响算法的收敛速度。其中，为矩阵的最大特征值μ 取值小，收敛速度慢，计算工作量大，但滤波性能较好。μ 取值大，收敛速度快，计算工作量小，滤波性能较差。但 μ 大到一定值时，收敛速度变化不明显，且取值过大，会造成计算溢出.

(2)自适应滤波器阶数 N 的影响： N 不可以任意选取，需要根据经验加上实际仿真验证比较才能最终确定.。

当阶数 N 取值大时，迭代次数增加，收敛速度变快。但当阶数 N 大到一定程度，收敛速度变化不明显，且可能引起系数迭代过程不收敛.

自适应滤波器原理：

输入信号x(n)

、参考信号d(n)、误差信号e(n)

e(n)为y(n)与d(n)的误差信号，根据最小均方误差算法调节自适应滤波器算法的参数，来优化滤波器结构，从而实现尽可能地保留图像固有的信息。

LMS算法最核心的思想是用平方误差代替均方误差,基本的LMS算法为:

w(k,n+ 1)=w(k,n)+ 2μe(n)x(n-k)

其中,w(k,n)和w(k,n+ 1)分别为迭代前后的系数值；n和n+ 1为前后两个时刻;k= 0,1,

… ,N- 1,N为滤波器的阶数；μ为收敛因子；e(n)=d(n)-y(n)=d(n)-xT(n)w(n)=d(n)-wT(n)x(n)为误差信号；x(n-k)为输入信号；y(n)=xT(n)w(n)为输出信号。

6、matlab仿真

进行小波变换→ 进行滤波后→ 再进行图像重构→ 再小波反变换进行图像重构→先用维纳滤波器再用自适应滤波器替换

%%对加入高斯噪声的图像进行小波分解

J=double(J);

[CA,CH,CV,CD]= dwt2(J,'db1 ');

% [c,s]=wavedec2(j,1,'haar');

y =[CA,CH;CV,CD];

y=uint8(y);

subplot(2,2,3);

imshow(y);

title('小波分解');

%xx=wiener2(CA,[3 3]); %对加噪图像进行二维自适应维纳滤波 %滤波器窗口大小

yy=wiener2(CH,[3 3]); %对加噪图像进行二维自适应维纳滤波 %滤波器窗口大小

zz=wiener2(CV,[3 3]); %对加噪图像进行二维自适应维纳滤波 %滤波器窗口大小

qq=wiener2(CD,[3 3]); %对加噪图像进行二维自适应维纳滤波 %滤波器窗口大小

%%%图像重构

XX = idwt2(CA,yy,zz,qq,'db1');

XX=uint8(XX);

subplot(2,2,4);

imshow(XX);

title('重构以后的图像');

其他相关文章请点击：

图像算法一：【图像点运算】灰度直方图、伽马变换、灰度阈值变换、分段线性变换、直方图均衡化

图像算法二：【图像几何变换】平移、镜像、转置、缩放、旋转、插值

图像算法三：【图像增强--空间域】图像平滑、中值滤波、图像锐化

图像算法四：【图像增强--频率域】傅里叶变换、快速傅里叶变换、频域滤波、频域低通滤波、频域高通滤波

图像算法五：【图像小波变换】多分辨率重构、Gabor滤波器、Haar小波

图像算法六：【彩色图像处理】彩色模型、彩图处理

图像算法七：【形态学图像处理】二值运算、膨胀、腐蚀

图像算法八：【图像分割】边缘检测（Roberts，Sobel，canny）、霍夫变换（hough）、阈值分割、区域分割

图像算法九：【图像特征提取】特征降维、PCA人脸特征抽取、局部二进制

~~~~

图像算法五：【图像小波变换】多分辨率重构、Gabor滤波器、Haar小波相关推荐

离散小波变换 python_用python中的“haar”小波对图像进行离散小波变换
我试图在python中的图像上应用haar小波.这是代码from pywt import dwt2, idwt2 img = cv2.imread('xyz.png') cA, (cH, cV, cD ...
【小波变换】小波变换入门----haar小波
小波可以认为是一个带通滤波器,只允许频率和小波基函数频率相近的信号通过.小波变换的基本思想是用一组小波函数和基函数表示一个函数或者信号. haar小波变换首先,以haar小波变换过程为例来理解小波变 ...
小波变换中的多贝西小波（DB小波函数）概述
内容均来源于维基百科对db小波函数的介绍多贝西小波(英语:Daubechies Wavelet),是以比利时女性物理暨数学家英格丽·多贝西(Ingrid Daubechies)的名字命名之一种小波函 ...
R语言从经济时间序列中用HP滤波器，小波滤波和经验模态分解等提取周期性成分分析
经济时间序列的分析通常需要提取其周期性成分.最近我们被客户要求撰写关于经济时间序列的研究报告,包括一些图形和统计输出.这篇文章介绍了一些方法,可用于将时间序列分解为它们的不同部分.它基于<宏观经 ...
matlab 小波滤波器,matlab小波滤波器使用
研究db小波,用分解和重构滤波器和上下采样函数实现多分辨分析,代码如下: %%% 小波分解与重构 clear;close all; load noissin;Sig=noissin; %% 滤波器分解 ...
【转】小波图像分解与重构
原文地址:[转]小波图像分解与重构作者:freedesert 小波图像融合综述(1) 分类: 小波变换 2007-11-14 11:15 4728人阅读评论(6) 收藏举报图像融合是将 ...
基于Gabor滤波器的人脸特征提取
Gabor小波变换的核函数,其最大的特点是其与人脑的皮层神经细胞的二维反射区具有相似的特征.基于Gabor小波变换的特征信息能够有效提取各种细节结构信息,因此Gabor小波变换在图像识别领域有着较为广 ...
二维小波变换_小波变换-哈尔小波（Haar）
最近比较忙,但是在网上无意中看到一篇哈尔小波的文章,觉得浅显易懂,就将其摘录下来,原文链接: https://blog.csdn.net/u014744118/article/details/5043 ...
数字图像处理与Python实现笔记之图像小波变换与多分辨率
数字图像处理与Python实现笔记摘要绪论 1 数字图像处理基础知识 2 彩色图像处理初步 3 空间滤波 4 频域滤波 5 图像特征提取 6 图像压缩 7 图像小波变换与多分辨率 7.1 从傅里叶 ...

图像算法五：【图像小波变换】多分辨率重构、Gabor滤波器、Haar小波

小波变换：

图像算法五：【图像小波变换】多分辨率重构、Gabor滤波器、Haar小波相关推荐

最新文章

热门文章