CN109727200A - 基于贝叶斯张量分解的相似块堆聚图像消噪方法及系统 - Google Patents
基于贝叶斯张量分解的相似块堆聚图像消噪方法及系统 Download PDFInfo
- Publication number
- CN109727200A CN109727200A CN201711030903.4A CN201711030903A CN109727200A CN 109727200 A CN109727200 A CN 109727200A CN 201711030903 A CN201711030903 A CN 201711030903A CN 109727200 A CN109727200 A CN 109727200A
- Authority
- CN
- China
- Prior art keywords
- tensor
- noising
- image
- preliminary
- parameter
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000011159 matrix material Substances 0.000 claims abstract description 23
- 238000011156 evaluation Methods 0.000 claims abstract description 4
- 239000000203 mixture Substances 0.000 claims abstract description 4
- 239000013598 vector Substances 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 230000008030 elimination Effects 0.000 claims description 4
- 238000003379 elimination reaction Methods 0.000 claims description 4
- 239000000284 extract Substances 0.000 claims description 4
- 238000011084 recovery Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Landscapes
- Image Processing (AREA)
Abstract
一种基于贝叶斯张量分解的相似块堆聚图像消噪方法及系统,通过从原始图像中划分出若干图像块,将其中相似的图像块堆聚组成更高阶的张量,对高阶张量进行张量分解得到张量矩阵和噪声精度的联合概率,通过交替迭代法近似估计基于联合概率的极大对数似然模型参数和超参数并得到消噪高阶张量,最后将消噪高阶张量中每个图像块复原后得到消噪后的低阶图像。
Description
技术领域
本发明涉及的是一种图像处理领域的技术,具体是一种基于贝叶斯张量并行因子分解的相似块堆聚图像消噪方法及系统。
背景技术
现代社会数字图像的应用已经相当普及,人们对高质量图像要求越来越高。但在图像获取过程中,由于环境的影响和设备的精度问题,获取的图像数据会包含各种噪声,影响了图像的视觉效果。减少图像噪音的影响,尽可能地恢复图像原来的信息具有广泛的应用前景。对应各种不同的噪声,传统的消噪方法有:高斯低通滤波、均值滤波、中值滤波、小波变换消噪等。近年来基于图像特征提取的方法越来越多地应用在这上面,譬如基于稀疏编码的消噪算法、基于图像非局部相似性的算法、基于深层神经网路的算法等。然而,这些方法都需要把噪音具体统计指标作为输入变量,这个在实际图像中一般是很难直接得到的。这些方法中很多需要复杂的参数调节或者需要使用大量数据进行训练才能较好地完成某些特定场景下的消噪任务。
发明内容
本发明针对现有技术存在的上述不足,提出一种基于贝叶斯张量分解的相似块堆聚图像消噪方法及系统,只需要输入含有噪音的图像,在做了一个初步的预处理后之后,程序会把输入图像中的相似的小块组合成高阶张量,实现“堆聚”,然后利用的贝叶斯张量并行因子分解模型对张量的秩进行自动估计并进行张量分解,算出消噪结果并估计噪音方差,最后根据张量分解的结果对图像进行重构可以有效估计无噪音图像。
本发明是通过以下技术方案实现的:
本发明通过从原始图像中划分出若干图像块,将其中相似的图像块堆聚组成更高阶的张量,对高阶张量进行张量分解得到张量矩阵和噪声精度的联合概率,通过交替迭代法近似估计基于联合概率的极大对数似然模型参数和超参数并得到消噪高阶张量,最后将消噪高阶张量中每个图像块复原后得到消噪后的低阶图像。
所述的原始图像,具有三阶张量,即空间水平、垂直方向和颜色维度。
所述的初步消噪图像,是指在原图像上进行高斯滤波等简单操作后初步去除噪声的图像。
所述的图像块,通过用滑动窗口在图像中取出的若干小图像块。
所述的堆聚是指:首先用相对大一点的滑动距离生成初步消噪图像的参考图像块。对于每个参考图像块Yi,在图像中的一个更小的区域,以较小滑动距离的滑动窗口生成新的图像块,并从其中选择与之最相似的图像块以组成初步消噪图像块组原始图像相同位置的图像块则组成原始图像块组当原来的图像块是三阶张量,现在加上相似的非局部结合,则得到四阶张量;使用双滑动窗口是为了选出与Yi更相似的图像块。
所述的张量分解,即对得到的多个四阶张量原始图像块组按张量并行因子模式进行分解。张量分解算法采用贝叶斯模型,模型中参数和超参数具有条件概率关系,算法中对模型中每个参数和超参数采用先固定其他参数计算其极大对数似然的最优值的交替迭代法计算出因子,其中每个原始图像块组的噪音精度参数τ根据原始图像块组与初步消噪图像块组之差设定初始值。
所述的张量重构,当张量分解迭代收敛时,根据得到的因子通过重构得到消噪高阶张量
所述的复原是指:将消噪高阶张量中每一个图像块提取出来并置于原始位置得到消噪后的低阶图像。优选地,当存在多个图像块重叠,则取平均值作为消噪后的图像。
本发明涉及一种实现上述方法的系统,包括:初步消噪模块、堆聚模块、张量分解模块、张量重构模块以及复原模块,其中:初步消噪模块与堆聚相连并传输初步消噪后的图像,堆聚模块与张量模块相连并传输相似图像块堆聚组,张量分解模块与张量重构模块相连并传输因子分解结果,张量重构模块复原模块相连并传输张量重构的低秩图像堆聚块组。
技术效果
与现有技术相比,本发明张量分解技术能把图像数据在各个维度上的结构信息很好地提取出来,有助于更好地恢复被噪声污染的图像信息。本发明利用图像的非局部特性,使相似图像块数据堆聚,增强了此类方法在消噪应用中的效果。本发明可以提取带噪声图像在各个维度上的结构特征,对于图像或视频处理的后续步骤,具有很大的作用。
在消噪过程中,噪音参数τ等参数都是自动估计,不需要事先输入正确的数字,这是其他的消噪音算法所不能做到的,而实际应用中图像的真实噪声参数是很难得到的。本方法中有许多可以并行处理的过程,如果通过适当的并行计算优化,可以很大程度上减低计算时间。这些特性使得本发明在智能设备已经普及的现在具有较大实际应用价值。
附图说明
图1为本发明图像超立方体组成示意图;
图2为本发明方法示意图;
图3为实施例效果示意图;
图中:a为原图;b为噪声标准差为75的噪声图;c为消噪结果示意图,PSNR=25.0476。
具体实施方式
如图2所示,本实施例包括初步消噪、堆聚、张量分解、张量重构以及复原几个步骤。:
其中,张量分解包括以下步骤:
步骤1)设定图像块为三阶张量,即:三个维度分别是空间水平、垂直方向和颜色维度;然后将其中相似的图像块沿新的维度组合成一个四阶张量具体步骤包括:输入的带噪音图像立方体可分解成其中:张量ε中的元素为为0均值高斯白噪音,τ为噪声精度(precision),即方差的倒数。
步骤2)将四阶张量按张量并行因子分解为:
其中:o表示向量的外积,是一维向量,张量矩阵是Kruskal算子表示中间式子的运算,R是张量可以分解成的秩一向量外积项的个数为即张量的并行因子秩,N是张量的阶数。
对每一个张量矩阵A(n)的概率满足:其中:是A(n)的第in行,Λ=diag(λ)是精度矩阵,即协方差矩阵的逆矩阵。
对所有的张量矩阵A(n),其精度矩阵Λ均相同,因为该精度矩阵为对角矩阵,对角线上的元素为λ=[λ1,...,λR],其概率满足:其中:表示γ分布;对应噪声精度的概率也满足p(τ)=Ga(τ|a0,b0),则有张量矩阵和噪声精度的联合概率其中:参数Θ={A(1),...,A(N),λ,τ}。
步骤3)为降低后验概率的计算复杂度,采用可解耦的分布q(Θ)近似估计联合概率即:对于每一个参数Θj∈Θ有:极大对数似然模型参数其中:表示在q分布下除去Θj外所有变量的均值,Θ为超参数。
采用交替迭代法求解极大对数似然模型参数和超参数,即:固定除Θj外的其他参数值,极大化对数似然lnq(Θj)得到更新后的参数Θj,循环此过程直到迭代收敛,从而求得消噪后的张量。
所述的迭代是指:
①对于A(n),可以推导得其中的参数通过以下方式更新:其中:⊙表示张量Khatri-Rao连乘。
②关于A的超参数λ,可以推得其中的参数通过求解以下得到:可以得到
③对于噪声的参数τ,同理可得:
④对于秩R,初始化时选取可能的最大值;当迭代进行时,可能会出现某个λr会得到比较大的值,导致对应外积项趋于0,于是这个外积项就自动去除,估计的秩减一,迭代到最后就可以自动计算出秩R。
本方法基于贝叶斯张量分解模型将带噪音的图像分解成清晰图像张量部分和噪音部分,清晰图像通过利用秩一张量的线性分解表示,其秩通过自适应算法估计。在自动估计张量秩之后可以有效地自动估计噪音图的噪音统计数据,从而使得算法不需要外部输入噪音的统计量也能够自动完成消噪过程。由于张量的多维结构,本方法也可支持多颜色等通道的图像数据。
本实施例使用经典的Lena彩色图像作为消噪实例的输入:选取图像大小为512×512像素,在图像上加上标准差为75的高斯白噪声作为本发明的输入图像,然后选取大小为32×32像素的图像块,首次采样时滑动窗口滑动距离为15,小区域采样滑动窗口滑动距离为3,每个图像块组有20个图像块,得到的结果如图3所示;可以看到本发明可以从高噪声中得出很好的结果。
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式对其进行局部调整,本发明的保护范围以权利要求书为准且不由上述具体实施所限,在其范围内的各个实现方案均受本发明之约束。
Claims (9)
1.一种基于贝叶斯张量分解的相似块堆聚图像消噪方法,其特征在于,通过从原始图像中划分出若干图像块,将其中相似的图像块堆聚组成更高阶的张量,对高阶张量进行张量分解得到张量矩阵和噪声精度的联合概率,通过交替迭代法近似估计基于联合概率的极大对数似然模型参数和超参数并得到消噪高阶张量,最后将消噪高阶张量中每个图像块复原后得到消噪后的低阶图像。
2.根据权利要求1所述的方法,其特征是,所述的初步消噪图像块,通过对原始图像经过初步消噪操作,得到初步消噪图,然后在初步消噪图上用滑动窗口生成若干图像块。
3.根据权利要求1所述的方法,其特征是,所述的堆聚是指:先用较大滑动距离生成参考初步消噪图像块,对于每个参考初步消噪图像块Yi,在初步消噪图像中的一个更小的区域,以不同的滑动距离选取初步消噪图像块,并从其中选择与之最相似的图像块以组成初步消噪图像块组原始图像相同位置的图像块则组成原始图像块组当原来的图像块是三阶张量,现在加上相似的非局部结合,则得到四阶张量。
4.根据权利要求1所述的方法,其特征是,所述的张量分解,即对得到的多个四阶张量原始图像块组按张量并行因子模式进行分解;张量分解算法采用贝叶斯模型,模型中参数和超参数具有条件概率关系,算法中对模型中每个参数和超参数采用先固定其他参数计算其极大对数似然的最优值的交替迭代法计算出因子,其中每个原始图像块组的噪音精度参数τ根据原始图像块组与初步消噪图像块组之差设定初始值。
5.根据权利要求1所述的方法,其特征是,所述的复原是指:当张量分解迭代收敛时,根据得到的因子通过重构得到消噪高阶张量将消噪高阶张量中每一个图像块提取出来并置于原始位置得到消噪后的低阶图像。
6.根据权利要求5所述的方法,其特征是,当存在多个图像块重叠,则取平均值作为消噪后的图像。
7.根据权利要求3或4所述的方法,其特征是,所述的四阶张量为:
其中:o表示向量的外积,是一维向量,张量矩阵是Kruskal算子表示中间式子的运算,R是张量分解成的秩一向量外积项的个数为即张量的并行因子秩,N是张量的阶数;
对每一个张量矩阵A(n)的概率满足:其中:是A(n)的第in行,Λ=diag(λ)是精度矩阵,即协方差矩阵的逆矩阵;对所有的张量矩阵A(n),其精度矩阵Λ均相同,因为该精度矩阵为对角矩阵,对角线上的元素为λ=[λ1,...,λR],其概率满足:其中:表示γ分布;对应噪声精度的概率也满足p(τ)=Ga(τ|a0,b0),则有张量矩阵和噪声精度的联合概率其中:参数Θ={A(1),...,A(N),λ,τ}。
8.根据权利要求4所述的方法,其特征是,所述的交替迭代法是指:
①对于A(n),可以推导得其中的参数通过以下方式更新:其中:⊙表示张量Khatri-Rao连乘;
②关于A的超参数λ,可以推得其中的参数通过求解以下得到:得到
③对于噪声的参数τ,同理可得: τ=aM/bM;
④对于秩R,初始化时选取可能的最大值;当迭代出现某个λr会得到比较大的值,导致对应外积项趋于0,于是这个外积项就自动去除,估计的秩减一,迭代到最后就可以自动计算出秩R。
9.一种实现上述任一权利要求所述方法的系统,其特征在于,包括:初步消噪模块、堆聚模块、张量分解模块、张量重构模块以及复原模块,其中:初步消噪模块与堆聚相连并传输初步消噪后的图像,堆聚模块与张量模块相连并传输相似图像块堆聚组,张量分解模块与张量重构模块相连并传输因子分解结果,张量重构模块复原模块相连并传输张量重构的低秩图像堆聚块组。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711030903.4A CN109727200A (zh) | 2017-10-30 | 2017-10-30 | 基于贝叶斯张量分解的相似块堆聚图像消噪方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711030903.4A CN109727200A (zh) | 2017-10-30 | 2017-10-30 | 基于贝叶斯张量分解的相似块堆聚图像消噪方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109727200A true CN109727200A (zh) | 2019-05-07 |
Family
ID=66291198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711030903.4A Pending CN109727200A (zh) | 2017-10-30 | 2017-10-30 | 基于贝叶斯张量分解的相似块堆聚图像消噪方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109727200A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111276183A (zh) * | 2020-02-25 | 2020-06-12 | 云南大学 | 一种基于参数估计的张量分解处理海量基因序列的方法 |
CN116233697A (zh) * | 2022-04-19 | 2023-06-06 | 中国科学院声学研究所 | 一种基于深度学习的声反馈的抑制方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040008904A1 (en) * | 2003-07-10 | 2004-01-15 | Samsung Electronics Co., Ltd. | Method and apparatus for noise reduction using discrete wavelet transform |
US20080166064A1 (en) * | 2007-01-05 | 2008-07-10 | Guoyi Fu | Method And Apparatus For Reducing Noise In An Image Using Wavelet Decomposition |
CN103049892A (zh) * | 2013-01-27 | 2013-04-17 | 西安电子科技大学 | 基于相似块矩阵秩最小化的非局部图像去噪方法 |
CN105957026A (zh) * | 2016-04-22 | 2016-09-21 | 温州大学 | 基于非局部相似图像块内部和块间隐性低秩结构的去噪方法 |
CN106846268A (zh) * | 2017-01-04 | 2017-06-13 | 温州大学 | 一种高斯‑脉冲混合图像噪声去除方法 |
CN107292337A (zh) * | 2017-06-13 | 2017-10-24 | 西北工业大学 | 超低秩张量数据填充方法 |
-
2017
- 2017-10-30 CN CN201711030903.4A patent/CN109727200A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040008904A1 (en) * | 2003-07-10 | 2004-01-15 | Samsung Electronics Co., Ltd. | Method and apparatus for noise reduction using discrete wavelet transform |
US20080166064A1 (en) * | 2007-01-05 | 2008-07-10 | Guoyi Fu | Method And Apparatus For Reducing Noise In An Image Using Wavelet Decomposition |
CN103049892A (zh) * | 2013-01-27 | 2013-04-17 | 西安电子科技大学 | 基于相似块矩阵秩最小化的非局部图像去噪方法 |
CN105957026A (zh) * | 2016-04-22 | 2016-09-21 | 温州大学 | 基于非局部相似图像块内部和块间隐性低秩结构的去噪方法 |
CN106846268A (zh) * | 2017-01-04 | 2017-06-13 | 温州大学 | 一种高斯‑脉冲混合图像噪声去除方法 |
CN107292337A (zh) * | 2017-06-13 | 2017-10-24 | 西北工业大学 | 超低秩张量数据填充方法 |
Non-Patent Citations (3)
Title |
---|
QIBIN ZHAO ET AL: "Bayesian CP Factorization of Incomplete Tensors with Automatic Rank Determination", 《IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE》, vol. 37, no. 9, 15 January 2015 (2015-01-15), pages 1752 - 1757 * |
QIBIN ZHAO ET AL: "Bayesian Robust Tensor Factorization for Incomplete Multiway Data", 《 IEEE TRANSACTIONS ON NEURAL NETWORKS AND LEARNING SYSTEMS》, vol. 27, no. 4, 30 April 2016 (2016-04-30), pages 736 - 748, XP011603274, DOI: 10.1109/TNNLS.2015.2423694 * |
王忠美: "张量组稀疏表示的高光谱图像去噪算法", 《测绘学报》, vol. 46, no. 5, 31 May 2017 (2017-05-31), pages 614 - 622 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111276183A (zh) * | 2020-02-25 | 2020-06-12 | 云南大学 | 一种基于参数估计的张量分解处理海量基因序列的方法 |
CN111276183B (zh) * | 2020-02-25 | 2023-03-21 | 云南大学 | 一种基于参数估计的张量分解处理海量基因序列的方法 |
CN116233697A (zh) * | 2022-04-19 | 2023-06-06 | 中国科学院声学研究所 | 一种基于深度学习的声反馈的抑制方法及系统 |
CN116233697B (zh) * | 2022-04-19 | 2023-09-05 | 中国科学院声学研究所 | 一种基于深度学习的声反馈的抑制方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shi et al. | Deep networks for compressed image sensing | |
Liu et al. | Single image dehazing with depth-aware non-local total variation regularization | |
CN106204467B (zh) | 一种基于级联残差神经网络的图像去噪方法 | |
CN109859147A (zh) | 一种基于生成对抗网络噪声建模的真实图像去噪方法 | |
CN103873743B (zh) | 一种基于结构张量和卡尔曼滤波的视频去噪方法 | |
CN105913393A (zh) | 一种自适应小波阈值图像去噪算法及装置 | |
Xu et al. | Patch ordering-based SAR image despeckling via transform-domain filtering | |
CN109345474A (zh) | 基于梯度域和深度学习的图像运动模糊盲去除方法 | |
CN105427264B (zh) | 一种基于群稀疏系数估计的图像重构方法 | |
US11880903B2 (en) | Bayesian image denoising method based on distribution constraint of noisy images | |
CN103020918B (zh) | 基于形状自适应邻域均值的非局部均值去噪方法 | |
CN104680510A (zh) | Radar视差图优化方法、立体匹配视差图优化方法及系统 | |
CN101944230B (zh) | 基于多尺度的自然图像非局部均值去噪方法 | |
CN104156918B (zh) | 基于联合稀疏表示与残差融合的sar图像噪声抑制方法 | |
CN103455991A (zh) | 一种多聚焦图像融合方法 | |
CN106846253A (zh) | 一种基于反向传播神经网络的图像超分辨率重建方法 | |
CN107123094A (zh) | 一种混合泊松、高斯和脉冲噪声的视频去噪方法 | |
Cozzolino et al. | Nonlocal SAR image despeckling by convolutional neural networks | |
CN107301631B (zh) | 一种基于非凸加权稀疏约束的sar图像降斑方法 | |
Zha et al. | The power of triply complementary priors for image compressive sensing | |
CN109727200A (zh) | 基于贝叶斯张量分解的相似块堆聚图像消噪方法及系统 | |
Liu et al. | Generic model-agnostic convolutional neural network for single image dehazing | |
CN107292855A (zh) | 一种结合自适应非局部样本和低秩的图像去噪方法 | |
CN110533608A (zh) | 图像条带噪声抑制方法及其装置、电子设备、存储介质 | |
CN109003247A (zh) | 彩色图像混合噪声的去除方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190507 |