CN105701775B - 一种基于改进自适应字典学习的图像去噪方法 - Google Patents
一种基于改进自适应字典学习的图像去噪方法 Download PDFInfo
- Publication number
- CN105701775B CN105701775B CN201610005896.1A CN201610005896A CN105701775B CN 105701775 B CN105701775 B CN 105701775B CN 201610005896 A CN201610005896 A CN 201610005896A CN 105701775 B CN105701775 B CN 105701775B
- Authority
- CN
- China
- Prior art keywords
- dictionary
- coefficient
- signal
- image
- rarefaction representation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 64
- 230000006872 improvement Effects 0.000 title description 4
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000008569 process Effects 0.000 claims description 7
- 239000000203 mixture Substances 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- HOWHQWFXSLOJEF-MGZLOUMQSA-N systemin Chemical compound NCCCC[C@H](N)C(=O)N[C@@H](CCSC)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H]([C@@H](C)O)C(=O)N[C@@H](CC(O)=O)C(=O)OC(=O)[C@@H]1CCCN1C(=O)[C@H]1N(C(=O)[C@H](CC(O)=O)NC(=O)[C@H](CCCN=C(N)N)NC(=O)[C@H](CCCCN)NC(=O)[C@H](CO)NC(=O)[C@H]2N(CCC2)C(=O)[C@H]2N(CCC2)C(=O)[C@H](CCCCN)NC(=O)[C@H](CO)NC(=O)[C@H](CCC(N)=O)NC(=O)[C@@H](NC(=O)[C@H](C)N)C(C)C)CCC1 HOWHQWFXSLOJEF-MGZLOUMQSA-N 0.000 claims description 3
- 108010050014 systemin Proteins 0.000 claims description 3
- 238000000205 computational method Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 8
- 230000013016 learning Effects 0.000 description 14
- 230000006870 function Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000001815 facial effect Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于改进自适应字典学习的图像去噪方法,包括:输入数据集合,初始化,设定字典D的初值;固定字典D,寻找信号yi在字典D第j列上的最大稀疏表示系数,以及字典D中相对应的列;更新信号yi;计算新的稀疏系数与已有稀疏系数均值的比值,得到信号yi的稀疏表示;更新字典D;得到适合的冗余字典D。本发明有益效果:相比于传统利用噪声方差作为阈值的方法,该方法不用估计噪声能量,而且在噪声变化的情况下,具有较好的鲁棒性。
Description
技术领域
本发明涉及图像处理领域的,特别涉及一种基于改进自适应字典学习的图像去噪方法。
背景技术
1.1图像去噪
在图像处理过程中,从传感器获得的图像信息常常被噪声污染,因此,为了提高图像的质量,需要将原始图像中掺杂的噪声去除掉。图像去噪的数学模型一般表示为:
y=y1+n
其中,y1为原始图像信息,n为噪声信息,y为获得的受污染的图像信息,它实际是图像信息和噪声信息的混合。图像去噪的目的就是从受污染的图像y中去除噪声n的信息,以提高图像的质量。
图像去噪有很多种方法,总的来说分为基于空域的方式和基于变换域的方式。基于字典学习的图像去噪方法属于变换域的方式,这种方法可以通过对被噪声污染的图像进行学习,得到能够反映图像特点的字典,然后利用得到的字典对图像进行稀疏表示,将噪声去除掉,该方法可以有效地抑制噪声,较好地保留图像边缘和纹理等重要信息,而且此方法能够根据不同的特点的图像,生成相应的字典,具有很好的适应性。
1.2噪声分类
高斯噪声
高斯噪声是指概率密度函数服从高斯分布(即正态分布)的一类噪声,宇宙噪声、热噪声和散粒噪声等噪声都可以表示为高斯噪声。图像噪声通常被假设为高斯噪声,其概率密度函数为
其中,μ为噪声均值,σ为方差。
非高斯噪声
虽然我们通常假设图像中的噪声为高斯噪声,但是在实际应用中,还有许多噪声不符合高斯分布,例如电磁噪声,大气噪声和人为噪声等。非高斯噪声主要有混合高斯分布,广义高斯分布和α稳定分布等几种类型。混合高斯分布是比较常用的一种类型,其概率密度函数为
其中,εi为第i个高斯分量的加权系数,且满足每个高斯分量的均值为μi,方差为σi,噪声总的方差为当C=2时,称为二元混合高斯分布。
1.3基于字典学习的图像去噪
根据稀疏表示理论,任意图像y1可以表示为
y1=Dx1
其中,为冗余字典,其列向量可以表示为为稀疏表示系数,可以通过下式求得:
其中,|| ||0表示向量的0范数,即向量中非零元素的个数。
根据稀疏表示理论,图像y1是稀疏的,则向量x1的0范数是一个有限值,即
||x1||0<L
对于含有噪声的图像y,在冗余字典下,可以表示为
y=y1+n=Dx1+Dx2=D(x1+x2)=Dx
其中,x2为噪声的稀疏表示系数,由于噪声是非稀疏信号,因此x2中各元素的值均较小,且分布均匀。因此我们可以用x中较大的L个元素,来精确逼近y1,即
||y1-DxL||p≤δ,xL=[x1,x2...xL],x1≥x2≥...xL≥xL+1≥...≥xK
其中,δ为噪声的方差,xL为x中较大的L个稀疏表示系数。
在稀疏表示过程中,相对于不同的字典,信号有不同的稀疏表示形式。因此字典的选择是信号稀疏表示的关键。字典学习的方法通过对一系列含噪声图像的训练,得到适应于特定图像的字典,从而提高稀疏表示的性能,达到去噪目的。
1.3传统字典学习方法
K-SVD是一种有效的传统字典学习方法,其利用数据集合来得到适合的冗余字典D,即
整个过程主要分两步:
第一步,固定字典D,寻找图像的稀疏表示
第二步,更新字典D
不断重复这两个过程,最终会得到一个较好的字典。
1.4传统K-SVD方法的缺点
传统基于字典学习的图像去噪方法中,在字典学习过程中和利用字典去噪的过程中,都涉及到稀疏表示问题,即用有限个稀疏系数,结合字典来表示信号。在表示的过程中都假定稀疏系数的个数为L个,但参数L实际上是未知的,因此在实际应用中,都是当稀疏表示的信号和原始信号的差值足够小时,就认为能够稀疏表示该信号了,即
||y1-DxL||p≤δ,xL=[x1,x2...xL],x1≥x2≥...xL≥xL+1≥...≥xK
其中δ为噪声的方差。
专利《基于稀疏自适应字典的图像去噪方法》(申请号:CN201310161279.7)在稀疏表示过程中,也是应用噪声方差作为判据。
这种利用噪声方差来作为稀疏表示是否完成的判据方式,存在以下问题:
(1)图像噪声实际上是未知的。需要用额外的方法来估计噪声,这不仅增加了去噪算法的复杂性,而且降低了稀疏表示的精度。
(2)图像噪声实际上是不稳定的。我们通常都假设图像噪声是稳定的高斯白噪声,有确定的方差δ,但是在实际环境中,噪声的来源千差万别,并不是都是稳定的高斯白噪声,因此,噪声的方差在整个图像范围内并不是均匀的,而且在某些情况下,是随时间变化的。在这种情况下,利用噪声的方差δ来作为稀疏表示是否完成的判据,存在很大的不确定性。
因此,本发明设计了一种与噪声方差无关的判定稀疏表示是否完成的方法,即,根据稀疏系数之间的比值来确定该系数是否属于较大系数,并将其用于字典学习中,大大提高了字典学习的性能和图像去噪的效果。
发明内容
本发明的目的就是为了解决上述问题,提供了一种基于改进自适应字典学习的图像去噪方法,该方法不依赖于噪声的方差,大大提高了字典学习的性能和图像去噪的效果。
为实现上述目的,本发明采用下述技术方案:
一种基于改进自适应字典学习的图像去噪方法,包括以下步骤:
(1)输入数据集合初始化,设定字典D的初值;
(2)固定字典D,设定寻找次数j=1,寻找其在字典D各个列上的最大稀疏表示系数x’j,以及字典D中相对应的列dj;
(3)更新信号yi,令j=j+1,计算更新后的信号yi,寻找其在字典D其余各个列上的新的最大稀疏表示系数x’j,以及字典D中相对应的列dj;
(4)计算新的稀疏系数与已有稀疏系数均值的比值,如果所述比值满足设定条件,则得到信号yi的一个新的稀疏表示系数,返回步骤(3),继续寻找下一个稀疏表示系数;否则,说明信号yi的稀疏表示系数已经寻找完毕,结束本过程,进行下一步;
(5)更新字典D;
(6)重复步骤(2)——步骤(5),得到适合的冗余字典D,使其满足:
其中,yi为输入的信号数据,xi为数据的稀疏表示系数,Y为信号yi的集合,X为稀疏表示系数xi的集合。
所述步骤(1)中字典D的初值设为余弦矩阵。
所述步骤(3)中更新信号yi的方法为:
yi=yi-<dj,xj>。
所述步骤(4)的具体方法为:
计算新的稀疏系数与已有稀疏系数均值的比值,即:
如果返回步骤(3);否则,说明寻找得到的新的稀疏系数过小,是噪声在字典上的表示,而不是图像在字典上的表示,说明信号yi的稀疏表示系数已经寻找完毕,结束本过程;
其中,xj为更新信号yi后计算得到的新的稀疏系数;xw为原有稀疏系数。
所述步骤(5)的具体方法为:
1)对于字典D中的第k列dk,稀疏表示系数的集合X中与之相对应的系数为计算其残差Ek;
2)数据集合中所有的稀疏表示系数组成的矩阵中,将与字典D中的第k列dk的相对应项的编号组成集合ωk;
3)定义Ωk为N×|ωk|的矩阵,所述矩阵Ωk中,第(ωk(i),i)个元素为1,其余元素为0,则根据矩阵Ωk与残差Ek的乘积得到矩阵
4)对矩阵进行奇异值分解,即
其中,Δ为的特征值的对角阵,U的列为的正交特征向量,V的列为的正交特征向量。
5)选择矩阵U的第一列作为dk,完成字典D中的第k列的更新。
所述步骤1)中,残差Ek的计算方法具体为:
其中,Y为信号yi的集合,dj为字典D中的第j列;为字典D中的第j列上的最大稀疏表示系数。
本发明的有益效果是:
本方法在字典学习的稀疏编码过程中,利用信号系数与噪声系数幅值差别较大的特性,通过计算新的稀疏系数与已有稀疏系数均值的比值,来判断新的稀疏系数是否属于信号。相比于传统利用噪声方差作为阈值的方法,该方法不用估计噪声能量,而且在噪声变化的情况下,具有较好的鲁棒性。
附图说明
图1为本发明方法流程图;
图2(a)为本发明具体实施例原始图像;
图2(b)为本发明具体实施例加入噪声后的图像;
图3(a)为本发明具体实施例采用K-SVD方法去噪的效果图;
图3(b)为本发明具体实施例采用本发明方法去噪的效果图。
下面结合附图和实施例对本发明进一步说明。
本发明对传统的K_SVD字典学习方法进行改进,在稀疏表示过程中,根据稀疏系数之间的比值来确定该系数是否属于较大系数。该方法不依赖于噪声的方差,大大提高了字典学习的性能和图像去噪的效果,我们称之为自适应字典学习方法(AK-SVD)。
一种基于改进自适应字典学习的图像去噪方法,如图1所示,包括以下步骤:
第一步,初始化,将字典D的初值设为余弦矩阵,
第二步,固定字典D,寻找信号yi的稀疏表示
具体方法为:
(1)对于每一个信号yi,设寻找次数j=1,寻找其在字典D各个列上的最大稀疏表示系数x’j,以及字典D中相对应的列dj。
(2)更新信号yi,即yi=yi-<dj,xj>,j=j+1。
(3)对于信号yi,寻找其在字典D各个列上的新的最大稀疏表示系数x’j,以及字典D中相对应的列dj。
(4)如果重复步骤(2)、(3)。如果不满足,说明所有系数都已找到,得到信号yi的稀疏表示:
第三步,更新字典D。
1)对于字典D中的第k列dk,X中与之相对应的系数为计算其残差Ek
2)数据集合中,用到字典D中的第k列dk的集合为ωk,即
3)定义Ωk为N×|ωk|的矩阵,其中第(ωk(i),i)个元素为1,其余元素为0,则
(4)对矩阵进行奇异值分解,即
(5)选择矩阵U的第一列作为dk,完成字典D中的第k列的更新。
重复上述第二步和第三步,优化字典D,重复次数在10到20次之间,根据图像的噪声来确定,噪声越重,需要的次数越多。
利用数据集合来得到适合的冗余字典D,即
本发明通过实际图像去噪实验对比了传统K-SVD字典学习方法和本发明提出的AK-SVD方法对于非高斯噪声的去噪效果。
如图2(a)和图2(b)所示,在原始图像中加入二元混合高斯噪声,参数为μ1=μ2=0,ε1=ε2=0.5,σ1=10,σ2=30则总的噪声方差加入噪声后,图像的峰值信噪比为24.2dB。
实验结果如图3(a)和图3(b)所示。试验结果表明,在非高斯噪声的环境下,K-SVD方法在已知总的噪声方差的条件下,去噪后图像的峰值信噪比为30.4dB,而利用本发明方法在未知噪声参数的条件下,去噪后图像的峰值信噪比为30.7dB。可见本发明方法较K-SVD方法理想情况具有略好的去噪效果,而且对比两幅去噪后的图像可以看出,在纹理比较丰富的部分,例如桌布,人物脸部和人物裤子部分,本发明方法的细节保留程度更好一些。
需要指出的是,虽然在本实验中本发明方法较K-SVD方法的去噪效果在峰值信噪比上仅改善了0.3dB,但是这是建立在K-SVD方法准确知道所加噪声方差的情况下的结果,是理想情况下的结果。而在实际应用中,图像噪声是未知的,需要通过其他算法来估计噪声的方差,其估计结果与真实值还会存在一定偏差,这会影响K-SVD方法的去噪效果,而本发明方法则不依赖于噪声的方差,具有计算简单和性能卓越的特点。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (6)
1.一种基于改进自适应字典学习的图像去噪方法,其特征是,包括以下步骤:
(1)输入数据集合初始化,设定字典D的初值;
(2)固定字典D,设定寻找次数j=1,寻找信号yi在字典D第j列上的最大稀疏表示系数x’j,以及字典D中相对应的列dj;
(3)更新信号yi,令j=j+1,计算更新后的信号yi在字典D第j列上的最大稀疏表示系数,以及字典D中相对应的列;
(4)计算新的稀疏系数与已有稀疏系数均值的比值,如果所述比值满足设定条件,则得到信号yi的一个新的稀疏表示系数,返回步骤(3),继续寻找下一个稀疏表示系数;否则,说明信号yi的稀疏表示系数已经寻找完毕,结束本过程,进行下一步;
(5)更新字典D;
(6)重复步骤(2)——步骤(5),得到适合的冗余字典D,使其满足:
其中,yi为输入的信号数据,xi为数据的稀疏表示系数,Y为信号yi的集合,X为稀疏表示系数xi的集合,L为稀疏系数的个数。
2.如权利要求1所述的一种基于改进自适应字典学习的图像去噪方法,其特征是,所述步骤(1)中字典D的初值设为余弦矩阵。
3.如权利要求1所述的一种基于改进自适应字典学习的图像去噪方法,其特征是,所述步骤(3)中更新信号yi的方法为:
yi=yi-<dj,x’j>。
4.如权利要求1所述的一种基于改进自适应字典学习的图像去噪方法,其特征是,所述步骤(4)的具体方法为:
计算新的稀疏系数与已有稀疏系数均值的比值,即:
如果返回步骤(3);否则,说明新的稀疏系数过小,是噪声在字典上的表示,而不是图像在字典上的表示,说明信号yi的稀疏表示系数已经寻找完毕,结束本过程;
其中,xj为更新信号yi后计算得到的新的稀疏系数;xw为原有稀疏系数。
5.如权利要求1所述的一种基于改进自适应字典学习的图像去噪方法,其特征是,所述步骤(5)的具体方法为:
1)对于字典D中的第k列dk,稀疏表示系数的集合X中与之相对应的系数为计算其残差Ek;
2)数据集合中所有的稀疏表示系数组成的矩阵中,将与字典D中的第k列dk的相对应项的编号组成集合ωk;
3)定义Ωk为N×|ωk|的矩阵,所述矩阵Ωk中,第(ωk(i),i)个元素为1,其余元素为0,则根据矩阵Ωk与残差Ek的乘积得到矩阵
4)对矩阵进行奇异值分解,即
其中,Δ为的特征值的对角阵,U的列为的正交特征向量,V的列为的正交特征向量;
5)选择矩阵U的第一列作为dk,完成字典D中的第k列的更新。
6.如权利要求5所述的一种基于改进自适应字典学习的图像去噪方法,其特征是,所述步骤1)中,残差Ek的计算方法具体为:
其中,Y为信号yi的集合,dj为字典D中的第j列;为字典D中的第j列上的最大稀疏表示系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610005896.1A CN105701775B (zh) | 2016-01-06 | 2016-01-06 | 一种基于改进自适应字典学习的图像去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610005896.1A CN105701775B (zh) | 2016-01-06 | 2016-01-06 | 一种基于改进自适应字典学习的图像去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105701775A CN105701775A (zh) | 2016-06-22 |
CN105701775B true CN105701775B (zh) | 2018-06-08 |
Family
ID=56226184
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610005896.1A Expired - Fee Related CN105701775B (zh) | 2016-01-06 | 2016-01-06 | 一种基于改进自适应字典学习的图像去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105701775B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106203414B (zh) * | 2016-07-01 | 2019-07-05 | 昆明理工大学 | 一种基于判别字典学习与稀疏表示的场景图片文字检测的方法 |
CN106228524A (zh) * | 2016-08-01 | 2016-12-14 | 广东工业大学 | 一种图像去噪方法及装置 |
CN107065006B (zh) * | 2017-01-23 | 2019-06-11 | 武汉大学 | 一种基于在线字典更新的地震信号编码方法 |
CN107451992B (zh) * | 2017-07-20 | 2020-08-11 | 广东工业大学 | 一种sar图像变化检测的方法与装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102088606B (zh) * | 2011-02-28 | 2012-12-05 | 西安电子科技大学 | 基于稀疏表示的去块效应方法 |
CN102427527B (zh) * | 2011-09-27 | 2013-11-06 | 西安电子科技大学 | 基于分布式视频压缩感知系统非关键帧的重构方法 |
CN103218791B (zh) * | 2013-05-05 | 2015-09-30 | 西安电子科技大学 | 基于稀疏自适应字典的图像去噪方法 |
-
2016
- 2016-01-06 CN CN201610005896.1A patent/CN105701775B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105701775A (zh) | 2016-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105701775B (zh) | 一种基于改进自适应字典学习的图像去噪方法 | |
TWI794157B (zh) | 自動多閾值特徵過濾方法及裝置 | |
CN110531313B (zh) | 一种基于深度神经网络回归模型的近场信号源定位方法 | |
Hyvarinen | Fast ICA for noisy data using Gaussian moments | |
Bhotto et al. | Robust shrinkage affine-projection sign adaptive-filtering algorithms for impulsive noise environments | |
CN107368799B (zh) | 基于多特征和自适应时延估计的泄漏检测与定位方法 | |
Li et al. | Signal denoising with random refined orthogonal matching pursuit | |
CN110084828A (zh) | 一种图像分割方法、装置及终端设备 | |
CN107292855B (zh) | 一种结合自适应非局部样本和低秩的图像去噪方法 | |
Li et al. | Direct-optimization-based DC dictionary learning with the MCP regularizer | |
CN111079691A (zh) | 一种基于双流网络的剪枝方法 | |
Li et al. | Blind separation and extraction of binary sources | |
CN110118958B (zh) | 基于变分编码-解码网络的宽带雷达复回波去噪方法 | |
Xi et al. | Recursive identification for Hammerstein systems with state-space model | |
CN109544593B (zh) | 一种基于自相似度量的通用图像噪声估计方法 | |
CN113452350B (zh) | 一种变步长块稀疏仿射投影自适应滤波器 | |
Novosadová et al. | Piecewise-polynomial signal segmentation using reweighted convex optimization | |
CN107085832A (zh) | 一种数字图像非局部平均去噪的快速实现方法 | |
Tekin et al. | Benefits of consistency in image denoising with steerable wavelets | |
Priorov et al. | Applications of image filtration based on principal component analysis and nonlocal image processing | |
Fang et al. | An efficient dictionary learning algorithm for sparse representation | |
Tang et al. | Image denoising via expected patch log likelihood with Gaussian model identification | |
Alhabsi | Improved SPSA optimization algorithm requiring a single measurement per iteration | |
Mowakeaa et al. | On the characterization, generation, and efficient estimation of the complex multivariate GGD | |
Li et al. | Image denoising algorithm based on independent component analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180608 Termination date: 20190106 |
|
CF01 | Termination of patent right due to non-payment of annual fee |