CN107292838A - 基于模糊区域分割的图像去模糊方法 - Google Patents
基于模糊区域分割的图像去模糊方法 Download PDFInfo
- Publication number
- CN107292838A CN107292838A CN201710420889.2A CN201710420889A CN107292838A CN 107292838 A CN107292838 A CN 107292838A CN 201710420889 A CN201710420889 A CN 201710420889A CN 107292838 A CN107292838 A CN 107292838A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mover
- fuzzy region
- 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 20
- 238000009826 distribution Methods 0.000 claims description 15
- 230000009466 transformation Effects 0.000 claims description 10
- 230000021615 conjugation Effects 0.000 claims description 8
- 238000013507 mapping Methods 0.000 claims description 3
- 238000004904 shortening Methods 0.000 abstract description 2
- 230000015556 catabolic process Effects 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000011218 segmentation Effects 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/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- 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/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明涉及一种基于模糊区域分割的图像去模糊方法，包括以下步骤：S1为表征特征差异，采用尖峰措施区分模糊区域B和非模糊区域U；S2利用Graph‑cut算法将模糊区域B和非模糊区域U分割；S3对模糊区域B进行原始图像区域L的估计；S4对模糊区域B进行模糊核k的估计；S5对模糊区域B进行图像反卷积，得到清晰区域；S6将清晰区域与非模糊区域U重新融合，得到最后的去模糊结果。本发明先将一幅模糊图像的模糊区域与非模糊区域进行分割，对于模糊区域进行原始图像区域和模糊核的估计，然后再进行图像反卷积，从而将模糊区域去模糊得到清晰区域，再将清晰区域与之前的非模糊区域相结合，得到模糊图像去模糊的最终结果，有效提高效率和缩短时间。
Description
技术领域
本发明涉及图像处理领域，具体涉及一种基于模糊区域分割的图像去模糊方法。
背景技术
模糊图像是图像退化导致细节丢失的其中一种类型。造成图像模糊的原因是在成像过程中，由于外部环境条件的限制和成像设备自身的物理局限等众多因素的影响，不可避免地会导致拍摄得到的图像出现模糊。运动模糊图像复原的关键是找到图像的退化模型，并采取逆过程求解原始图像，这个过程就称为模糊图像的去模糊过程。去模糊方法可以分为非盲源图像去卷积方法和图像盲反卷积方法。非盲源图像去卷积方法是假设在模糊核已知的情况下，对模糊图像进行复原；图像盲反卷积方法是假设模糊核未知的情况下，对模糊图像进行复原。
在模糊图像去模糊方法的研究中，人们更经常是直接解决图像非盲源去卷积和盲反卷积问题，即人们直接明确地估计原始图像和模糊核，从而复原图像。模糊图像往往仅仅存在局部模糊，将完整一幅图像进行图像去模糊处理，必然会花费更多的时间，不能做到时间和效率并行。
发明内容
本发明的目的在于提供一种基于模糊区域分割的图像去模糊方法，以解决现有技术对局部模糊的图像进行去模糊处理需要花费很多时间，效率低下的问题。
为了实现上述的目的，采用如下的技术方案。一种基于模糊区域分割的图像去模糊方法，包括以下步骤：
S1为表征特征差异，采用尖峰措施区分模糊区域B和非模糊区域U；
S2利用Graph-cut算法将模糊区域B和非模糊区域U分割；
S3对模糊区域B进行原始图像区域L的估计；
S4对模糊区域B进行模糊核k的估计；
S5对模糊区域B进行图像反卷积，得到清晰区域；
S6将清晰区域与非模糊区域U重新融合，得到最后的去模糊结果。
步骤S1中，尖峰措施为采用kurtosis衡量尖峰分布：
式中，K(·)表示kurtosis分布，用来表示尖峰分布情况；a表示输入数据向量；E(·)表示输入数据a的期望值；
对于模糊图像I，定义第一个特征为：
f1＝min(ln(K(Ix)+3),ln(K(Iy)+3))；
式中，f1即为求得的第一个特征值；ln(·)将特征映射到一个合适的范围中；min(·)运算符在x轴和y轴的值之间选择一个更小的值；K(·)即为kurtosis尖峰分布求得的值；
对于kurtosis，模糊区域B对应更小的值，非模糊区域U对应更大的值，以此将模糊图像I分割成模糊区域B和非模糊区域U。
步骤S3中，模糊区域B表示为模糊核k卷积原始图像区域L后再加噪声η，其数学模型为：
原始图像区域L的能量函数为：
式中，表示使得后面对应公式成立的L的最小值；表示不同方向和顺序的偏导数，ω*∈{ω0,ω1,ω2}是对每一类偏导数的权重，λ是正则项的权重；
式中的第二项是一个正则项，正则项中的和可以通过下面的公式计算得到：
式中，根据经验，能够使得公式达到最优值；
为求解上面的公式，引入新的变量lxi和lyi：
式中，是对相应项求偏导；
则最优化能量函数后，根据Parseval定理，由傅里叶变换可以得到原始图像区域L的能量函数的最优解：
式中，F(·)表示快速傅里叶变换；而F-1(·)和分别表示快速傅里叶变换的逆变换和快速傅里叶变换的共轭取值；λ是正则项的权重；d由以下公式求得：
式中，F(·)和分别表示快速傅里叶变换和快速傅里叶变换的共轭取值；ω0，ω1，ω2分别表示每一项对应的权重。步骤S3中，模糊区域B表示为模糊核k卷积原始图像区域L后再加噪声η，其数学模型为：
原始图像区域L的能量函数为：
式中，表示使得后面对应公式成立的L的最小值；表示不同方向和顺序的偏导数，ω*∈{ω0,ω1,ω2}是对每一类偏导数的权重，λ是正则项的权重；
式中的第二项是一个正则项，正则项中的和可以通过下面的公式计算得到：
式中，根据经验，能够使得公式达到最优值；
为求解上面的公式，引入新的变量lxi和lyi：
式中，是对相应项求偏导；
则最优化能量函数后，根据Parseval定理，由傅里叶变换可以得到原始图像区域L的能量函数的最优解：
式中，F(·)表示快速傅里叶变换；而F-1(·)和分别表示快速傅里叶变换的逆变换和快速傅里叶变换的共轭取值；λ是正则项的权重；d由以下公式求得：
式中，F(·)和分别表示快速傅里叶变换和快速傅里叶变换的共轭取值；ω0，ω1，ω2分别表示每一项对应的权重。
步骤S4中，由原始图像区域L得到模糊核k的能量函数：
式中，表示使得后面对应公式成立的k的最小值；表示x,y方向的偏导数算子；λ是正则项的权重；
由傅里叶变换得到模糊核k的能量函数的最优解为：
式中，表示使得后面对应公式成立的k的最小值；表示x,y方向的偏导数算子；λ是正则项的权重。
步骤S5中，由于自然图像的梯度基本服从“重尾”分布，所以在反卷积过程中，采用L0.5为正则项的超拉普拉斯先验算法进行反卷积。
与现有技术相比，本发明先将一幅模糊图像的模糊区域与非模糊区域进行分割，对于模糊区域进行原始图像区域和模糊核的估计，然后再进行图像反卷积，从而将模糊区域去模糊得到清晰区域，再将清晰区域与之前的非模糊区域相结合，得到模糊图像去模糊的最终结果，有效提高效率和缩短时间。
附图说明
图1为本发明的流程图；
图2为实施例选取的模糊图像；
图3为实施例利用尖峰措施得到的模糊区域和非模糊区域；
图4为实施例利用Graph-cut算法得到的分割图像；
图5为实施例模糊区域的模糊核；
图6为实施例模糊图像的去模糊结果。
具体实施方式
下面结合附图和实施例对本发明作进一步的描述。
本发明的流程如图1所示，包括以下步骤：
S1为表征特征差异，采用尖峰措施区分模糊区域B和非模糊区域U；
S2利用Graph-cut算法将模糊区域B和非模糊区域U分割；
S3对模糊区域B进行原始图像区域L的估计；
S4对模糊区域B进行模糊核k的估计；
S5对模糊区域B进行图像反卷积，得到清晰区域；
S6将清晰区域与非模糊区域U重新融合，得到最后的去模糊结果。
实施例采用本发明对图2进行去模糊处理。
(1)自然图像服从重尾分布。一般而言，模糊区域一般不包含尖锐边缘，则其分布一般包含较小的值。因此，为表征特征差异，采用尖峰措施来区分模糊区域B和非模糊区域U。采用kurtosis来衡量尖峰分布：
式中，K(·)表示kurtosis分布，用来表示尖峰分布情况；a表示输入数据向量；E(·)表示输入数据a的期望值；
对于一幅模糊图像I，假设它可以被分为模糊区域和非模糊区域，则定义第一个特征为：
f1＝min(ln(K(Ix)+3),ln(K(Iy)+3))；
式中，f1即为求得的第一个特征值；ln(·)将特征映射到一个合适的范围中；min(·)运算符在x轴和y轴的值之间选择一个更小的值；K(·)即为kurtosis尖峰分布求得的值；
在上面的公式中，对于kurtosis而言，模糊区域对应于更小的值，而非模糊区域对应于更大的值。根据此来将模糊图像I分割成模糊区域和非模糊区域。如图3所示，图3中的黑色部分表示非模糊区域，白色区域表示模糊区域。
(2)利用Graph-cut算法将模糊区域B和非模糊区域U分割。利用Graph-cut算法分割后的非模糊区域U如图4所示。将非模糊区域U保留，将模糊区域B转入到步骤(3)继续处理。
(3)对模糊区域B进行原始图像区域L的估计。模糊区域B可以表示为模糊核k卷积原始图像区域L后再加噪声η，其数学模型为：
为了对模糊区域B去模糊，首先需要先求原始图像区域L。原始图像L的能量函数为：
式中，表示使得后面对应公式成立的L的最小值；表示不同方向和顺序的偏导数，ω*∈{ω0,ω1,ω2}是对每一类偏导数的权重，λ是正则项的权重；
式中的第二项是一个正则项，正则项中的和可以通过下面的公式计算得到：
式中，根据经验，能够使得公式达到最优值；
为了更好的求解上面的公式，引入新的变量lxi和lyi：
式中，是对相应项求偏导；
则最优化能量函数后，根据Parseval定理，由傅里叶变换可以得到L的最优解：
式中，F(·)表示快速傅里叶变换；而F-1(·)和分别表示快速傅里叶变换的逆变换和快速傅里叶变换的共轭取值；λ是正则项的权重；d由以下公式求得：
式中，F(·)和分别表示快速傅里叶变换和快速傅里叶变换的共轭取值；ω0，ω1，ω2分别表示每一项对应的权重。
(4)对模糊区域B进行模糊核k的估计。得到L的最优解后，对模糊核k进行估计，其能量函数为：
式中，表示使得后面对应公式成立的k的最小值；表示x,y方向的偏导数算子；λ是正则项的权重；
同样，通过傅里叶变换后，最优化上面的公式，可得
式中，表示使得后面对应公式成立的k的最小值；表示x,y方向的偏导数算子；λ是正则项的权重。图5展示的是两幅不同模糊图像的模糊区域分别得到的模糊核。
(5)对模糊区域B进行图像反卷积。由于自然图像的梯度基本服从“重尾”分布，所以在反卷积过程中，将采用L0.5为正则项的超拉普拉斯先验算法进行反卷积。
(6)将模糊区域B处理后得到的清晰区域与非模糊区域U重新融合，得到最后的去模糊结果。图6展示的是最后的去模糊结果。
Claims (5)
1.一种基于模糊区域分割的图像去模糊方法，其特征在于，包括以下步骤：
S1为表征特征差异，采用尖峰措施区分模糊区域B和非模糊区域U；
S2利用Graph-cut算法将模糊区域B和非模糊区域U分割；
S3对模糊区域B进行原始图像区域L的估计；
S4对模糊区域B进行模糊核k的估计；
S5对模糊区域B进行图像反卷积，得到清晰区域；
S6将清晰区域与非模糊区域U重新融合，得到最后的去模糊结果。
2.根据权利要求1所述的基于模糊区域分割的图像去模糊方法，其特征在于，步骤S1中，尖峰措施为采用kurtosis衡量尖峰分布：
<mrow>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>a</mi>
<mn>4</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mn>3</mn>
<mo>;</mo>
</mrow>
式中，K(·)表示kurtosis分布，用来表示尖峰分布情况；a表示输入数据向量；E(·)表示输入数据a的期望值；
对于模糊图像I，定义第一个特征为：
f1＝min(ln(K(Ix)+3),ln(K(Iy)+3))；
式中，f1即为求得的第一个特征值；ln(·)将特征映射到一个合适的范围中；min(·)运算符在x轴和y轴的值之间选择一个更小的值；K(·)即为kurtosis尖峰分布求得的值；
对于kurtosis，模糊区域B对应更小的值，非模糊区域U对应更大的值，以此将模糊图像I分割成模糊区域B和非模糊区域U。
3.根据权利要求1所述的基于模糊区域分割的图像去模糊方法，其特征在于，步骤S3中，模糊区域B表示为模糊核k卷积原始图像区域L后再加噪声η，其数学模型为：
<mrow>
<mi>B</mi>
<mo>=</mo>
<mi>L</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>+</mo>
<mi>&eta;</mi>
</mrow>
原始图像区域L的能量函数为：
<mrow>
<munder>
<mrow>
<mi>arg</mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>L</mi>
</munder>
<munder>
<mo>&Sigma;</mo>
<msub>
<mo>&part;</mo>
<mo>*</mo>
</msub>
</munder>
<msub>
<mi>&omega;</mi>
<mo>*</mo>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>k</mi>
<mo>&CircleTimes;</mo>
<msub>
<mo>&part;</mo>
<mo>*</mo>
</msub>
<mi>L</mi>
<mo>-</mo>
<msub>
<mo>&part;</mo>
<mo>*</mo>
</msub>
<mi>B</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>&phi;</mi>
<mo>(</mo>
<mrow>
<msub>
<mo>&part;</mo>
<mi>x</mi>
</msub>
<mi>L</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mi>&phi;</mi>
<mo>(</mo>
<mrow>
<msub>
<mo>&part;</mo>
<mi>y</mi>
</msub>
<mi>L</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
式中，表示使得后面对应公式成立的L的最小值；表示不同方向和顺序的偏导数，ω*∈{ω0,ω1,ω2}是对每一类偏导数的权重，λ是正则项的权重；
式中的第二项是一个正则项，正则项中的和可以通过下面的公式计算得到：
<mrow>
<msub>
<mi>&phi;</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<msup>
<mi>&epsiv;</mi>
<mn>2</mn>
</msup>
</mfrac>
<mo>|</mo>
<mi>L</mi>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>L</mi>
<mo><</mo>
<mi>&epsiv;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mi>o</mi>
<mi>t</mi>
<mi>h</mi>
<mi>e</mi>
<mi>r</mi>
<mi>w</mi>
<mi>i</mi>
<mi>s</mi>
<mi>e</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
式中，根据经验，能够使得公式达到最优值；
为求解上面的公式，引入新的变量lxi和lyi：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>l</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<msub>
<mo>&part;</mo>
<mi>x</mi>
</msub>
<msub>
<mi>L</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>&le;</mo>
<mi>&epsiv;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mo>&part;</mo>
<mi>x</mi>
</msub>
<msub>
<mi>L</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>o</mi>
<mi>t</mi>
<mi>h</mi>
<mi>e</mi>
<mi>r</mi>
<mi>w</mi>
<mi>i</mi>
<mi>s</mi>
<mi>e</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>l</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<msub>
<mo>&part;</mo>
<mi>y</mi>
</msub>
<msub>
<mi>L</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>&le;</mo>
<mi>&epsiv;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mo>&part;</mo>
<mi>y</mi>
</msub>
<msub>
<mi>L</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>o</mi>
<mi>t</mi>
<mi>h</mi>
<mi>e</mi>
<mi>r</mi>
<mi>w</mi>
<mi>i</mi>
<mi>s</mi>
<mi>e</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
式中，是对相应项求偏导；
则最优化能量函数后，根据Parseval定理，由傅里叶变换可以得到原始图像区域L的能量函数的最优解：
<mrow>
<mi>L</mi>
<mo>=</mo>
<msup>
<mi>F</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>B</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>d</mi>
<mo>+</mo>
<mfrac>
<mi>&lambda;</mi>
<msup>
<mi>&epsiv;</mi>
<mn>2</mn>
</msup>
</mfrac>
<mrow>
<mo>(</mo>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>x</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mo>(</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>+</mo>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>y</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mo>(</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>d</mi>
<mo>+</mo>
<mfrac>
<mi>&lambda;</mi>
<msup>
<mi>&epsiv;</mi>
<mn>2</mn>
</msup>
</mfrac>
<mrow>
<mo>(</mo>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>x</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>x</mi>
</msub>
<mo>)</mo>
<mo>+</mo>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>y</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>y</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
式中，F(·)表示快速傅里叶变换；而F-1(·)和分别表示快速傅里叶变换的逆变换和快速傅里叶变换的共轭取值；λ是正则项的权重；d由以下公式求得：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>d</mi>
<mo>=</mo>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<mo>&CenterDot;</mo>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&omega;</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>x</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>x</mi>
</msub>
<mo>)</mo>
<mo>+</mo>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>y</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>y</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>+</mo>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>+</mo>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中，F(·)和分别表示快速傅里叶变换和快速傅里叶变换的共轭取值；ω0，ω1，ω2分别表示每一项对应的权重。
4.根据权利要求1所述的基于模糊区域分割的图像去模糊方法，其特征在于，步骤S4中，由原始图像区域L得到模糊核k的能量函数：
<mrow>
<munder>
<mrow>
<mi>arg</mi>
<mi>min</mi>
</mrow>
<mi>k</mi>
</munder>
<mo>|</mo>
<mo>|</mo>
<mi>k</mi>
<mo>&CircleTimes;</mo>
<mo>&dtri;</mo>
<mi>L</mi>
<mo>-</mo>
<mo>&dtri;</mo>
<mi>B</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mi>&lambda;</mi>
<mo>|</mo>
<mo>|</mo>
<mi>k</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>;</mo>
</mrow>
式中，表示使得后面对应公式成立的k的最小值；表示x,y方向的偏导数算子；λ是正则项的权重；
由傅里叶变换得到模糊核k的能量函数的最优解为：
<mrow>
<mi>k</mi>
<mo>=</mo>
<msup>
<mi>F</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>x</mi>
</msub>
<mi>L</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>x</mi>
</msub>
<mi>B</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>y</mi>
</msub>
<mi>L</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>y</mi>
</msub>
<mi>B</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>x</mi>
</msub>
<mi>L</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>x</mi>
</msub>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mover>
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>y</mi>
</msub>
<mi>L</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
<mi>F</mi>
<mrow>
<mo>(</mo>
<msub>
<mo>&part;</mo>
<mi>y</mi>
</msub>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&lambda;</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
式中，表示使得后面对应公式成立的k的最小值；表示x,y方向的偏导数算子；λ是正则项的权重。
5.根据权利要求1所述的基于模糊区域分割的图像去模糊方法，其特征在于，步骤S5中，采用L0.5为正则项的超拉普拉斯先验算法进行反卷积。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710420889.2A CN107292838A (zh) | 2017-06-07 | 2017-06-07 | 基于模糊区域分割的图像去模糊方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710420889.2A CN107292838A (zh) | 2017-06-07 | 2017-06-07 | 基于模糊区域分割的图像去模糊方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107292838A true CN107292838A (zh) | 2017-10-24 |
Family
ID=60096080
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710420889.2A Pending CN107292838A (zh) | 2017-06-07 | 2017-06-07 | 基于模糊区域分割的图像去模糊方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107292838A (zh) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108230276A (zh) * | 2018-02-06 | 2018-06-29 | 江西理工大学 | 一种融合无参图像质量评价的自然场景图像去模糊方法 |
| CN108564552A (zh) * | 2018-04-28 | 2018-09-21 | 深圳市商汤科技有限公司 | 图像去模糊的方法及装置 |
| CN108898563A (zh) * | 2018-07-02 | 2018-11-27 | 京东方科技集团股份有限公司 | 一种显示面板光学检测图像的处理方法及计算机可读介质 |
| CN109089040A (zh) * | 2018-08-20 | 2018-12-25 | Oppo广东移动通信有限公司 | 图像处理方法、图像处理装置及终端设备 |
| CN109636733A (zh) * | 2018-10-26 | 2019-04-16 | 华中科技大学 | 基于深度神经网络的荧光图像解卷积方法及系统 |
| CN110111282A (zh) * | 2019-05-09 | 2019-08-09 | 杭州电子科技大学上虞科学与工程研究院有限公司 | 一种基于运动矢量和cnn的视频去模糊方法 |
| CN111340722A (zh) * | 2020-02-20 | 2020-06-26 | Oppo广东移动通信有限公司 | 图像处理方法、处理装置、终端设备及可读存储介质 |
| WO2020182230A3 (zh) * | 2019-03-11 | 2020-10-22 | 影石创新科技股份有限公司 | 一种图像融合方法及便携式终端 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103186894A (zh) * | 2013-03-22 | 2013-07-03 | 南京信息工程大学 | 一种自适应分块的多聚焦图像融合方法 |
| US20150086127A1 (en) * | 2013-09-20 | 2015-03-26 | Samsung Electronics Co., Ltd | Method and image capturing device for generating artificially defocused blurred image |
| CN104867111A (zh) * | 2015-03-27 | 2015-08-26 | 北京理工大学 | 一种基于分块模糊核集的非均匀视频盲去模糊方法 |
| CN106157268A (zh) * | 2016-07-28 | 2016-11-23 | 浙江工业大学 | 一种基于l0凸近似的退化图像复原方法 |
-
2017
- 2017-06-07 CN CN201710420889.2A patent/CN107292838A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103186894A (zh) * | 2013-03-22 | 2013-07-03 | 南京信息工程大学 | 一种自适应分块的多聚焦图像融合方法 |
| US20150086127A1 (en) * | 2013-09-20 | 2015-03-26 | Samsung Electronics Co., Ltd | Method and image capturing device for generating artificially defocused blurred image |
| CN104867111A (zh) * | 2015-03-27 | 2015-08-26 | 北京理工大学 | 一种基于分块模糊核集的非均匀视频盲去模糊方法 |
| CN106157268A (zh) * | 2016-07-28 | 2016-11-23 | 浙江工业大学 | 一种基于l0凸近似的退化图像复原方法 |
Non-Patent Citations (2)
| Title |
|---|
| JIANPING SHI等: "Discriminative Blur Detection Features", 《2014 IEEE CONFERENCE ON COMPUTER VISION AND PATTERN RECOGNITION》 * |
| 闫敬文等: "基于L0范数正则项的运动图像去模糊", 《激光与光电子学进展》 * |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108230276B (zh) * | 2018-02-06 | 2020-06-02 | 江西理工大学 | 一种融合无参图像质量评价的自然场景图像去模糊方法 |
| CN108230276A (zh) * | 2018-02-06 | 2018-06-29 | 江西理工大学 | 一种融合无参图像质量评价的自然场景图像去模糊方法 |
| CN108564552A (zh) * | 2018-04-28 | 2018-09-21 | 深圳市商汤科技有限公司 | 图像去模糊的方法及装置 |
| CN108898563A (zh) * | 2018-07-02 | 2018-11-27 | 京东方科技集团股份有限公司 | 一种显示面板光学检测图像的处理方法及计算机可读介质 |
| CN109089040A (zh) * | 2018-08-20 | 2018-12-25 | Oppo广东移动通信有限公司 | 图像处理方法、图像处理装置及终端设备 |
| CN109089040B (zh) * | 2018-08-20 | 2021-05-14 | Oppo广东移动通信有限公司 | 图像处理方法、图像处理装置及终端设备 |
| CN109636733A (zh) * | 2018-10-26 | 2019-04-16 | 华中科技大学 | 基于深度神经网络的荧光图像解卷积方法及系统 |
| CN109636733B (zh) * | 2018-10-26 | 2020-07-24 | 华中科技大学 | 基于深度神经网络的荧光图像解卷积方法及系统 |
| WO2020182230A3 (zh) * | 2019-03-11 | 2020-10-22 | 影石创新科技股份有限公司 | 一种图像融合方法及便携式终端 |
| JP2022524806A (ja) * | 2019-03-11 | 2022-05-10 | 影石創新科技股▲ふん▼有限公司 | 画像融合方法及び携帯端末 |
| JP7253635B2 (ja) | 2019-03-11 | 2023-04-06 | 影石創新科技股▲ふん▼有限公司 | 画像融合方法及び携帯端末 |
| US11967051B2 (en) | 2019-03-11 | 2024-04-23 | Arashi Vision Inc. | Image fusion method and portable terminal |
| CN110111282A (zh) * | 2019-05-09 | 2019-08-09 | 杭州电子科技大学上虞科学与工程研究院有限公司 | 一种基于运动矢量和cnn的视频去模糊方法 |
| CN110111282B (zh) * | 2019-05-09 | 2021-05-11 | 杭州电子科技大学上虞科学与工程研究院有限公司 | 一种基于运动矢量和cnn的视频去模糊方法 |
| CN111340722A (zh) * | 2020-02-20 | 2020-06-26 | Oppo广东移动通信有限公司 | 图像处理方法、处理装置、终端设备及可读存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107292838A (zh) | 基于模糊区域分割的图像去模糊方法 | |
| Teodoro et al. | Image restoration and reconstruction using variable splitting and class-adapted image priors | |
| Chierchia et al. | Epigraphical projection and proximal tools for solving constrained convex optimization problems | |
| Sun et al. | Edge-based blur kernel estimation using patch priors | |
| JP6409872B2 (ja) | ぶれ画像をぶれ修正する方法およびシステム | |
| Zhang et al. | Bregmanized nonlocal regularization for deconvolution and sparse reconstruction | |
| Yuan et al. | l0tv: A new method for image restoration in the presence of impulse noise | |
| US11263728B2 (en) | Priori constraint and outlier suppression based image deblurring method | |
| Barbu | PDE-based restoration model using nonlinear second and fourth order diffusions | |
| CN108257108A (zh) | 一种超分辨率图像重建方法及系统 | |
| Liu et al. | Learning collaborative generation correction modules for blind image deblurring and beyond | |
| Kim et al. | An optimized first-order method for image restoration | |
| Zhou et al. | Towards the automation of deep image prior | |
| CN108898557B (zh) | 图像恢复方法及装置、电子设备、计算机程序及存储介质 | |
| Yang et al. | Guided filter based edge-preserving image non-blind deconvolution | |
| Tian et al. | A modeling method for face image deblurring | |
| Innerhofer et al. | A convex approach for image hallucination | |
| JP5784472B2 (ja) | 画像処理装置、方法及びプログラム | |
| Sahragard et al. | Image restoration by variable splitting based on total variant regularizer | |
| Hong et al. | Single image motion deblurring using anisotropic regularization | |
| Ono et al. | Total variation-wavelet-curvelet regularized optimization for image restoration | |
| Choudhury et al. | A denoising inspired deblurring framework for regularized image restoration | |
| El Abbadi et al. | Estimating Blur Parameters to Reconstruct the Motion Blurred Images | |
| He et al. | Wavelet frame truncated ℓ 1-regularized image deblurring | |
| Zhi et al. | An image inpainting algorithm based on energy minimization |
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 |
Application publication date: 20171024 |
|
| RJ01 | Rejection of invention patent application after publication |