CN107292855A - 一种结合自适应非局部样本和低秩的图像去噪方法 - Google Patents
一种结合自适应非局部样本和低秩的图像去噪方法 Download PDFInfo
- Publication number
- CN107292855A CN107292855A CN201710651882.1A CN201710651882A CN107292855A CN 107292855 A CN107292855 A CN 107292855A CN 201710651882 A CN201710651882 A CN 201710651882A CN 107292855 A CN107292855 A CN 107292855A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- image
- msubsup
- msub
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 13
- 230000009466 transformation Effects 0.000 claims abstract description 9
- 238000000354 decomposition reaction Methods 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 6
- 235000013399 edible fruits Nutrition 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000012360 testing method Methods 0.000 abstract description 6
- 230000000007 visual effect Effects 0.000 abstract description 4
- 239000000654 additive Substances 0.000 abstract description 2
- 230000000996 additive effect Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000007812 deficiency Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000003909 pattern recognition 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/74—Image or video pattern matching; Proximity measures in feature spaces
- G06V10/75—Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
- G06V10/751—Comparing pixel values or logical combinations thereof, or feature values having positional relevance, e.g. template matching
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20004—Adaptive image processing
- G06T2207/20012—Locally adaptive
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Databases & Information Systems (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Artificial Intelligence (AREA)
- Health & Medical Sciences (AREA)
- Multimedia (AREA)
- Image Processing (AREA)
Abstract
本发明公开一种结合自适应非局部样本和低秩的图像去噪方法,首先用对数变换将图像转换到对数域,将乘性噪声模型转换成加性噪声模型;将图像分块并且按照相似度分组,得到具有相似块的图像组;然后对图像组做低秩逼近处理,得到初始的估计值;再对初始的估计值用自适应非局部样本模型处理,得到对数域恢复结果;最后用指数变换将对数域图像还原到实数域并且进行修正,得到最终去噪图像。实验结果表明,本发明对乘性噪声有较好的鲁棒性,针对含有乘性噪声的图像不仅能得到很好的峰值信噪比和结构相似度,还较好地改善图像的视觉质量。
Description
技术领域
本发明涉及数字图像处理技术领域,具体涉及一种结合自适应非局部样本和低秩的图像去噪方法。
背景技术
图像去噪技术的发展历史最早从50年代开始,由于图像采集和传递中的技术缺陷和存在外界干扰,所以需要寻求新的技术来弥补这些弊端,提高图像的质量。随着科学技术的飞速发展,图像处理技术在人民生活、天文学、生物医学、VR技术、人工智能、公安司法、文化艺术等应用领域受到广泛重视并取得了开拓性成就,这些技术的预处理阶段的好坏都会直接关系到图像的后续处理。近十几年来,许多研究者都在寻求各种各样的思路对图像进行处理,探究如何利用图像自身的结构特点和数据规律对目标图像进行分析。
基于图像自身信息的图像恢复算法包括非局部均值算法“Buades A,Coll B,Morel J M.A non-local algorithm for image denoising.Computer Vision andPattern Recognition,2005:60-65”、块匹配3D滤波算法“Dabov K,Foi A,Katkovnik V,etal.Image denoising by sparse 3-D transform domain collaborative filtering[J].IEEE Transactions on Image Processing,2007,16(8):2080-2095”和低秩逼近方法“Cai J F,Candès E J,Shen Z.A Singular Value Thresholding algorithm for matrixcompletion.Siam Journal on Optimization,2010,20(4):1956-1982.”等。上述方法虽然在图像去噪上取得较好的效果,但是在保持边缘及轮廓等方面存在一定的不足。
由于自然图像中具有一定的冗余信息和相似性结构,因此,图像包含相似的块矩阵,而理论上是低秩的,但是在噪声的干扰下,却破坏了图像的低秩性。
发明内容
本发明所要解决的是现有图像去噪方法在保持边缘及轮廓等方面存在的不足,提供一种结合自适应非局部样本和低秩的图像去噪方法,其能提高去噪图像的质量,达到理想的实用效果。
为解决上述问题,本发明是通过以下技术方案实现的:
一种结合自适应非局部样本和低秩的图像去噪方法,包括步骤如下:
步骤1,读取实数域中的含噪图像y;
步骤2,将实数域中的含噪图像转化到对数域中,得到对数域中的含噪图像Y;
步骤3,将对数域中的含噪图像Y进行分块,并对每个图像块进行非局部相似匹配,将与之最相似的m个相似图像块放到一组,得到每个图像块所对应的图像组YI;
步骤4,初始化迭代次数t,X(0),Y(0),令t=0和X(0)=Y(0)=Y;其中为待恢复图像X(0)的第I个图像组,为含噪图像Y(0)的第I个图像组,Y为输入的含噪图像;
步骤5,对每个含噪图像组求均值并且训练它的PCA变换矩阵
步骤6,对每个待恢复图像组进行奇异值分解,并据此计算该图像组所对应的辅助变量即
步骤7,根据下式更新每个待恢复图像组的第k个稀疏编码即
步骤8,根据下式更新每个待恢复图像组的第k个稀疏编码即
步骤9,根据下式更新每个待恢复图像组的第k个图像相似块即
步骤10,将每个待恢复图像组的m个图像相似块集合在一起,得到待恢复图像组所有的待恢复图像组进行组合,得到目标图像
步骤11,根据下式更新噪声图,即
步骤12,判断迭代次数t是否达到预先设定的迭代次数T;如果达到,则转至步骤13;如果未达到,则返回步骤5进行继续迭代;
步骤13,利用指数变换和误差修正将步骤10所得到的目标图像还原到实数域中,再进行误差修正,得到最终实数域的去噪图像
上述各式中,上标t表示第t次迭代,上标t+1表示第t+1次迭代,上标t+1 2表示第t次和t+1次两次迭代的中间迭代;I=1,2,…,N,N为图像组的总数;k=1,2,…,m,m为图像组中相似图像块的个数;上标T代表转置;Sλ为第一软阈值算子,λ为设定的第一软阈值参数;为辅助变量,A(·)为奇异值分解所得的左酉矩阵,Σ(·)为奇异值分解所得的对角矩阵,B(·)为奇异值分解所得的右酉矩阵;Y为含噪图像,λ2为设定的正参数,(αI k)(·)为第k个稀疏编码,为PCA变换矩阵,E为单位向量,σI为标准差向量;Sτ为第二软阈值算子,τ为设定的第二阈值,为第二软阈值参数,c为一个常数,为噪声的方差;为均值;为第k个图像相似块;ρ为一个大于0的常数。
上述步骤13中,通过步骤11得到对数域中的去噪图像再利用指数变换和误差修正,得到最终实数域的去噪图像为
其中,L表示噪声的视数,exp(·)表示指数函数。
与现有技术相比,本发明的有益效果为:
1、引入了低秩逼近算法,克服了传统算法对图像边缘和轮廓的保持不足,使得图像具有清晰的结构,提高视觉效果。
2、使用了自适应非局部样本模型,充分使用了图像的自身信息和自身结构,可以很好的保持图像的细节和纹理,克服了传统稀疏表示的自适应性不足,使得去噪图像内容丰富,视觉质量高。
附图说明
图1是本发明结合自适应非局部样本和低秩的图像去噪方法的流程图。
图2是几种图像乘性噪声移除方法的仿真效果图。其中(a)为原图,(b)为CC算法,(c)为HMNZ算法,(d)为本发明算法。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实例,并参照附图,对本发明进一步详细说明。
参照图1,一种结合自适应非局部样本和低秩的图像去噪方法,具体实施步骤如下:
步骤1、读取实数域中的含噪图像y。
步骤2、将实数域中的含噪图y进行对数变换,得到对数域中的含噪图Y。
步骤3、利用块分配技术将对数域中的含噪图像Y分成7×7的小块,并利用非局部相似匹配算法寻找每个图像块的相似块,再通过欧式距离计算相似度,将最相似的m个图像块组成一个图像组其中是含噪图像块的yI第k个相似块,I指第I组含噪图像。
先利用块分配技术将大小为256×256的含噪图像Y分成有重叠的小块,每块的大小为7×7,定义为n为图像块的总数。然后,对于每个图像块yi,利用非局部相似匹配算法在搜索窗内寻找每个图像块的相似块,再通过欧式距离计算相似度,将最相似的m个图像块组成一个图像组其中是含噪图像块yI的第k个相似块,I指第I组含噪图像,I=1,2,…,N,N为图像组的总数。
步骤4、初始化迭代次数t,X(0),Y(0),
令t=0,X(0)=Y(0)=Y。为待恢复图像X(0)的第I个图像组,为含噪图像Y(0)的第I个图像组,Y为输入的含噪图像。
步骤5、对每个含噪图像组求均值并且训练它的PCA变换矩阵
对每个含噪图像组求该图像组的均值得到将每个图像组减去其均值得到中心化向量则中心化矩阵对的协方差矩阵Ω进行奇异值分解:Ω=ΦΛΦT,则PCA变换矩阵为Φ,其中Φ是酉矩阵,Λ是由Ω的特征值为对角线组成的对角矩阵,ΦT是Φ的转置。
步骤6、对每个待恢复图像组进行奇异值分解,即得到左,右酉矩阵A(t)和(B(t))T,其中是对角矩阵,是矩阵的正奇异值且σ1 (k-1)>σ2 (k-1)>…>σr (k-1),r是矩阵的秩,利用第一软阈值算子Sw(Σ(t))=max(Σ(t)-w,0),求出Sw(∑(t)),最后利用更新辅助变量其中为软阈值参数,d是参数,本实施例中取值为d=0.54,n是相似块的个数,本实施例中取值为n=8,ε是一个大于0的常数,防止分母为0,本实施例中取值为ε=0.00001,表示的奇异值,在计算时是未知的,用估计。
步骤7、利用最小二乘法迭代
更新每个待恢复图像组的第k个稀疏编码其中t+1/2是第t次和第t+1次两次迭代的中间值,计算出所有的得到相应图像组的稀疏编码其中λ2是正的参数,本实施例中取值为λ2=0.8,σI是标准差向量,E是单位向量。
步骤8、利用第二软阈值算子迭代更新每个待恢复图像组的第k个稀疏编码计算出所有的(αI k)(t+1),得到相应图像组的稀疏编码其中 是噪声的方差,Sτ(·)是第二软阈值算子,是第二软阈值参数。
步骤9、利用更新每个待恢复图像组的第k个图像相似块
步骤10、将每个待恢复图像组的m个图像相似块集合在一起,得到待恢复图像组所有待恢复图像组进行组合,得到目标图像
步骤11、利用更新噪声图,其中ρ是一个大于0的常数,在本实施例中取值为ρ=0.65。
步骤12、判断迭代次数t是否达到预先设定的迭代次数T,本实施例中噪声视数L=9,迭代次数T=4,如果未达到,返回步骤5进行继续迭代;如果达到,则转至步骤13。
步骤13、利用指数变换将还原到实数域中,并进行误差修正,得到最终实数域的去噪图像
通过步骤10得到的待恢复的图像再利用指数变换将还原到实数域中,并进行误差修正,得到最终的去噪图像,其中,L表示噪声的视数,exp(·)表示指数函数。
本发明的效果由以下仿真进一步说明。
1.仿真条件:
(1)选取大小为256×256的标准测试图像库中自然图像进行测试;
(2)将一幅图像分为7×7的图像块,每组相似块的个数为15,测试的图像噪声视数L=9。
2.仿真内容与结果:
仿真内容:利用256×256的House图、Barbara图和Pirate图,用本发明与已有技术的CC算法,HMNZ算法分别对其进行乘性噪声的移除。
3.实验结果:
本发明比CC算法和HMNZ算法在峰值信噪比PSNR和SSIM指标上取得更好的噪声移除效果。实验结果如表1和表2所示。
表1实验结果峰值信噪比(PSNR/dB)
表2实验结果结构相似度(SSIM)
从表中可以看出,本发明与CC算法和HMNZ算法相比,在不同的视数和测试图像上,均取得了较高的PSNR值和SSIM值。
图2是噪声视数为9时,对测试图Pirate采用不同的去噪算法所得到的实验效果图,其中(a)为原图,(b)为CC算法,(c)为HMNZ算法,(d)为本发明算法。从放大区域可以明显看出本发明恢复出的衣袖边缘清晰可见,而对比算法明显失真。由此可见本发明的视觉效果高于对比算法。
本发明首先用对数变换将图像转换到对数域,将乘性噪声模型转换成加性噪声模型;将图像分块并且按照相似度分组,得到具有相似块的图像组;然后对图像组做低秩逼近处理,得到初始的估计值;再对初始的估计值用自适应非局部样本模型处理,得到对数域恢复结果;最后用指数变换将对数域图像还原到实数域并且进行修正,得到最终去噪图像。实验结果表明,本发明对乘性噪声有较好的鲁棒性,针对含有乘性噪声的图像不仅能得到很好的峰值信噪比和结构相似度,还较好地改善图像的视觉质量。
需要说明的是,尽管以上本发明所述的实施例是说明性的,但这并非是对本发明的限制,因此本发明并不局限于上述具体实施方式中。在不脱离本发明原理的情况下,凡是本领域技术人员在本发明的启示下获得的其它实施方式,均视为在本发明的保护之内。
Claims (2)
1.一种结合自适应非局部样本和低秩的图像去噪方法,其特征是,包括步骤如下:
步骤1,读取实数域中的含噪图像y;
步骤2,将实数域中的含噪图像转化到对数域中,得到对数域中的含噪图像Y;
步骤3,将对数域中的含噪图像Y进行分块,并对每个图像块进行非局部相似匹配,将与之最相似的m个相似图像块放到一组,得到每个图像块所对应的图像组YI;
步骤4,初始化迭代次数t,X(0),XI (0),Y(0),令t=0和X(0)=Y(0)=Y;其中为待恢复图像X(0)的第I个图像组,为含噪图像Y(0)的第I个图像组,Y为输入的含噪图像;
步骤5,对每个含噪图像组求均值并且训练它的PCA变换矩阵
步骤6,对每个待恢复图像组进行奇异值分解,并据此计算该图像组所对应的辅助变量即
<mrow>
<msubsup>
<mi>U</mi>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>=</mo>
<msup>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msup>
<msub>
<mi>S</mi>
<mi>&lambda;</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>&Sigma;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>B</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>;</mo>
</mrow>
步骤7,根据下式更新每个待恢复图像组的第k个稀疏编码(αI k)(t+1/2),即
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>&alpha;</mi>
<mi>I</mi>
</msub>
<mi>k</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<mn>2</mn>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&Phi;</mi>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<msub>
<mi>&Phi;</mi>
<mi>I</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>&Phi;</mi>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>(</mo>
<msubsup>
<mi>&Phi;</mi>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<msubsup>
<mi>&Phi;</mi>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mi>Y</mi>
<mo>+</mo>
<mn>2</mn>
<msubsup>
<mi>&Phi;</mi>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msubsup>
<msubsup>
<mi>U</mi>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msubsup>
<mo>-</mo>
<mfrac>
<mi>E</mi>
<msub>
<mi>&sigma;</mi>
<mi>I</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
步骤8,根据下式更新每个待恢复图像组的第k个稀疏编码即
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>&alpha;</mi>
<mi>I</mi>
</msub>
<mi>k</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msub>
<mi>S</mi>
<mi>&tau;</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<msub>
<mi>&alpha;</mi>
<mi>I</mi>
</msub>
<mi>k</mi>
</msup>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>-</mo>
<msubsup>
<mi>&mu;</mi>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>&mu;</mi>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</msubsup>
<mo>;</mo>
</mrow>
步骤9,根据下式更新每个待恢复图像组的第k个图像相似块即
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mi>I</mi>
</msub>
<mi>k</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>&Phi;</mi>
<mi>I</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>I</mi>
<mi>k</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>;</mo>
</mrow>
步骤10,将每个待恢复图像组的m个图像相似块集合在一起,得到待恢复图像组所有的待恢复图像组进行组合,得到目标图像
步骤11,根据下式更新噪声图,即
<mrow>
<msup>
<mi>Y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>=</mo>
<msup>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>+</mo>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<mi>Y</mi>
<mo>-</mo>
<msup>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
步骤12,判断迭代次数t是否达到预先设定的迭代次数T;如果达到,则转至步骤13;如果未达到,则返回步骤5进行继续迭代;
步骤13,利用指数变换和误差修正将步骤10所得到的目标图像还原到实数域中,再进行误差修正,得到最终实数域的去噪图像
上述各式中,上标t表示第t次迭代,上标t+1表示第t+1次迭代,上标t+1/2表示第t次和t+1次两次迭代的中间迭代;I=1,2,…,N,N为图像组的总数;k=1,2,…,m,m为图像组中相似图像块的个数;上标T代表转置;Sλ为第一软阈值算子,λ为设定的第一软阈值参数;为辅助变量,A(·)为奇异值分解所得的左酉矩阵,Σ(·)为奇异值分解所得的对角矩阵,B(·)为奇异值分解所得的右酉矩阵;Y为含噪图像,λ2为设定的正参数,(αI k)(·)为第k个稀疏编码,为PCA变换矩阵,E为单位向量,σI为标准差向量;Sτ为第二软阈值算子,τ为设定的第二阈值,为第二软阈值参数,c为一个常数,为噪声的方差;为均值;为第k个图像相似块;ρ为一个大于0的常数。
2.根据权利要求1所述的一种结合自适应非局部样本和低秩的图像去噪方法,其特征是,步骤13中,通过步骤11得到对数域中的去噪图像再利用指数变换和误差修正,得到最终实数域的去噪图像为
<mrow>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<msup>
<mover>
<mi>X</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<mrow>
<mfrac>
<msup>
<mi>d</mi>
<mn>2</mn>
</msup>
<mrow>
<msup>
<mi>dL</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mo>(</mo>
<mrow>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mo>+</mo>
<mi>&infin;</mi>
</mrow>
</msubsup>
<mi>exp</mi>
<mrow>
<mo>(</mo>
<mrow>
<mo>-</mo>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
</mrow>
<msup>
<mi>&delta;</mi>
<mrow>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>d</mi>
<mi>&delta;</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mn>2</mn>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
其中,L表示噪声的视数,exp(·)表示指数函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710651882.1A CN107292855B (zh) | 2017-08-02 | 2017-08-02 | 一种结合自适应非局部样本和低秩的图像去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710651882.1A CN107292855B (zh) | 2017-08-02 | 2017-08-02 | 一种结合自适应非局部样本和低秩的图像去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107292855A true CN107292855A (zh) | 2017-10-24 |
CN107292855B CN107292855B (zh) | 2020-08-07 |
Family
ID=60104271
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710651882.1A Active CN107292855B (zh) | 2017-08-02 | 2017-08-02 | 一种结合自适应非局部样本和低秩的图像去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107292855B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110827223A (zh) * | 2019-11-05 | 2020-02-21 | 郑州轻工业学院 | 结合分数阶全变差的cs高噪声天文图像去噪重建方法 |
CN111405206A (zh) * | 2020-02-25 | 2020-07-10 | 南京邮电大学 | 基于帧平均及权值调整二阶盲识别的泊松-高斯联合噪声图像序列分离降噪方法 |
CN113436093A (zh) * | 2021-06-16 | 2021-09-24 | 中国电子科技集团公司第五十四研究所 | 一种基于能量域变换的遥感影像去噪方法 |
CN116188333A (zh) * | 2023-05-04 | 2023-05-30 | 北京师范大学 | 基于结构化整体最小二乘的乘性加性混合噪声去除方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103049892A (zh) * | 2013-01-27 | 2013-04-17 | 西安电子科技大学 | 基于相似块矩阵秩最小化的非局部图像去噪方法 |
CN103473740A (zh) * | 2013-08-31 | 2013-12-25 | 西安电子科技大学 | 基于稀疏表示和低秩双重约束的非局部去噪方法 |
CN104657951A (zh) * | 2015-03-02 | 2015-05-27 | 桂林电子科技大学 | 图像乘性噪声移除方法 |
US9111359B2 (en) * | 2012-05-07 | 2015-08-18 | Samsung Techwin Co., Ltd. | Method and system for detecting matrix-based motion including frequency transform and filtering |
CN105513026A (zh) * | 2015-12-17 | 2016-04-20 | 浙江大学 | 一种基于图像非局部相似的压缩感知重构方法 |
CN106204483A (zh) * | 2016-07-08 | 2016-12-07 | 桂林电子科技大学 | 基于非局部自适应字典的乘性噪声去除方法 |
-
2017
- 2017-08-02 CN CN201710651882.1A patent/CN107292855B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9111359B2 (en) * | 2012-05-07 | 2015-08-18 | Samsung Techwin Co., Ltd. | Method and system for detecting matrix-based motion including frequency transform and filtering |
CN103049892A (zh) * | 2013-01-27 | 2013-04-17 | 西安电子科技大学 | 基于相似块矩阵秩最小化的非局部图像去噪方法 |
CN103473740A (zh) * | 2013-08-31 | 2013-12-25 | 西安电子科技大学 | 基于稀疏表示和低秩双重约束的非局部去噪方法 |
CN104657951A (zh) * | 2015-03-02 | 2015-05-27 | 桂林电子科技大学 | 图像乘性噪声移除方法 |
CN105513026A (zh) * | 2015-12-17 | 2016-04-20 | 浙江大学 | 一种基于图像非局部相似的压缩感知重构方法 |
CN106204483A (zh) * | 2016-07-08 | 2016-12-07 | 桂林电子科技大学 | 基于非局部自适应字典的乘性噪声去除方法 |
Non-Patent Citations (3)
Title |
---|
RUI ZHU: "Spectral Nonlocal Restoration of Hyperspectral Images With Low-Rank Property", 《IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING 》 * |
唐中和: "低秩逼近理论及其在自然图像去噪中的应用", 《中国优秀硕士学位论文全文数据库》 * |
姚波旭: "基于低秩矩阵分解的非局部稀疏模型图像去噪方法研究", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110827223A (zh) * | 2019-11-05 | 2020-02-21 | 郑州轻工业学院 | 结合分数阶全变差的cs高噪声天文图像去噪重建方法 |
CN111405206A (zh) * | 2020-02-25 | 2020-07-10 | 南京邮电大学 | 基于帧平均及权值调整二阶盲识别的泊松-高斯联合噪声图像序列分离降噪方法 |
CN111405206B (zh) * | 2020-02-25 | 2022-07-26 | 南京邮电大学 | 一种泊松-高斯联合噪声图像序列分离降噪方法 |
CN113436093A (zh) * | 2021-06-16 | 2021-09-24 | 中国电子科技集团公司第五十四研究所 | 一种基于能量域变换的遥感影像去噪方法 |
CN116188333A (zh) * | 2023-05-04 | 2023-05-30 | 北京师范大学 | 基于结构化整体最小二乘的乘性加性混合噪声去除方法 |
CN116188333B (zh) * | 2023-05-04 | 2023-07-11 | 北京师范大学 | 基于结构化整体最小二乘的乘性加性混合噪声去除方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107292855B (zh) | 2020-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yeh et al. | Multi-scale deep residual learning-based single image haze removal via image decomposition | |
CN109859147B (zh) | 一种基于生成对抗网络噪声建模的真实图像去噪方法 | |
Lin et al. | Hyperspectral image denoising via matrix factorization and deep prior regularization | |
CN112233026A (zh) | 一种基于多尺度残差注意力网络的sar图像去噪方法 | |
CN107292855A (zh) | 一种结合自适应非局部样本和低秩的图像去噪方法 | |
CN106204482B (zh) | 基于加权稀疏的混合噪声去除方法 | |
Wang et al. | Deep convolutional architecture for natural image denoising | |
CN109636722B (zh) | 一种基于稀疏表示的在线字典学习超分辨率重建的方法 | |
Singh et al. | ResDNN: deep residual learning for natural image denoising | |
CN104657951A (zh) | 图像乘性噪声移除方法 | |
Wang et al. | New image restoration method associated with tetrolets shrinkage and weighted anisotropic total variation | |
CN110349112A (zh) | 一种基于自适应奇异值阈值的两阶段图像去噪方法 | |
CN112967210B (zh) | 一种基于全卷积孪生网络的无人机图像去噪方法 | |
CN106934398B (zh) | 基于超像素聚类和稀疏表示的图像去噪方法 | |
CN113378620B (zh) | 监控视频噪声环境下跨摄像头行人重识别方法 | |
Hu et al. | Patch-based low-rank minimization for image denoising | |
CN114005046A (zh) | 基于Gabor滤波器和协方差池化的遥感场景分类方法 | |
CN106169180A (zh) | 一种基于组的非局部稀疏表示加性噪声去除方法 | |
CN111340741B (zh) | 基于四元数与l1范数的粒子群优化灰度图像增强方法 | |
CN116051444A (zh) | 一种有效的红外与可见光图像自适应融合方法 | |
CN113781333A (zh) | 一种基于引导滤波的gan网络对水下图像进行处理的方法 | |
CN108090907B (zh) | 基于非高斯hmt模型的图像分割方法 | |
Xie et al. | Model‐guided boosting for image denoising | |
CN114723614B (zh) | 基于小波域admm深度网络的图像复原方法 | |
Jia et al. | Detachable image decomposition and illumination mapping search for low-light image enhancement |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |