CN107767351A - 一种盲去模糊的图像恢复方法 - Google Patents
一种盲去模糊的图像恢复方法 Download PDFInfo
- Publication number
- CN107767351A CN107767351A CN201711069974.5A CN201711069974A CN107767351A CN 107767351 A CN107767351 A CN 107767351A CN 201711069974 A CN201711069974 A CN 201711069974A CN 107767351 A CN107767351 A CN 107767351A
- Authority
- CN
- China
- Prior art keywords
- image
- model
- norm
- solving
- fuzzy
- 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 31
- 238000011084 recovery Methods 0.000 title claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 235000009508 confectionery Nutrition 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 abstract description 3
- 230000000750 progressive effect Effects 0.000 abstract 1
- 230000015556 catabolic process Effects 0.000 description 3
- 238000006731 degradation reaction Methods 0.000 description 3
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明公开了一种盲去模糊的图像恢复方法,属于图像处理领域。利用lq/l2范数作为正则化先验项,采用多尺度方法由粗尺度到细尺度逐步迭代求解,用***法求解模型时,用l1范数保真项更新估计图像的高频信息,在清晰图像恢复时,采用封闭阈值公式以解析解的形式给出,提高了算法速度,同时,更新模糊核时,提出线性递增权重参数,对模糊核按多尺度方法由粗到细逐步进行估计,进一步改善了模糊核的收敛性,提高了去模糊图像质量。
Description
技术领域
本发明属于图像处理技术领域,具体涉及一种盲去模糊的图像恢复方法。
背景技术
图像恢复是图像处理中的一大领域,有着广泛的应用,正成为当前研究的热点。图像恢复的主要目的是使退化图像经过一定的加工处理,去掉退化因素,以最大的保真度恢复成原来的图像。传统的图像恢复假设图像的降质模型是己知的。
而在实际应用中,图像退化***的点扩展函数一般是未知的,只能凭退化图像的观测数据,再附加很少的关于***与原图像先验知识来估计原图像,称之为盲图像复原。
现有技术在去模糊图像的质量上还不够高,有待于进一步提高,同时在算法上运算速度也有待于进一步提高,以及模糊核的收敛性还不够好。
发明内容
针对上述现有技术中描述的不足,本发明提供一种盲去模糊的图像恢复方法,该方法能进一步改善模糊核的收敛性,提高去模糊图像质量。
为解决上述技术问题,本发明所采用的技术方案如下:
一种盲去模糊的图像恢复方法,步骤如下:
步骤1,建立图像模糊模型为:
其中,f表示模糊图像,k表示模糊核矩阵,u表示清晰图像,n表示成像过程中的噪声;
步骤2,利用离散滤波器 产生高频信息构造凹凸范数比值正则化模型为:
其中,k≥0,∑iki=1,ki表示模糊核矩阵k中的元素,λ和λ1为权重参数,ψ为常数值,y为模糊图像的高频信息,x为清晰图像的高频信息,q为方程次方数且0<q<1;
步骤3,求解步骤2的凹凸范数比值正则化模型,得到清晰图像;
步骤3-1,更新清晰图像的高频信息x;
步骤3-1-1,保真项选择l1范数,求解模型为:
步骤3-1-2,将步骤3-1-1中的||x||2视为常数,则步骤3-1-1的求解模型转换为关于非凸||x||q范数的正则化模型:
步骤3-1-3,引入辅助变量v和权重参数θ,并将步骤3-1-2中的正则化模型转换为:
其中,
θi=ciθ0 (7);
式中,θ0表示权重参数θ的初始值,ci与尺度i成线性关系,ci=2i;
步骤3-1-4,采用β代替常数项λ||x||2,并对步骤3-1-3的正则化模型转换中的v和x进行分别求解,求解公式为:
步骤3-1-5,对步骤3-1-4的求解公式进行分别求导:
式中,xk+1表示第k+1步的x值,Δt为迭代步长;
步骤3-2,更新模糊核k;
步骤3-2-1,保真项选择l2范数的平方项,求解模型为:
步骤3-2-2,在模糊核k的更新过程中,利用IRLS方法计算模糊核k的权重:
其中,λ1和ψ为常数值,k0为初始模糊核,wk为更新的模糊核k的权重;
步骤3-2-3,根据步骤3-2-1和步骤3-2-2,在最细尺度求出模糊核函数k;
步骤3-3,得到恢复图像;
步骤3-3-1,在模糊核已知的情况下,图像模型求解变为非盲图像去卷积,求解模型为:
步骤3-3-2,令d=Du,其中D表示运算,将变量u与D分离,同时增加l2范数的平方项,引入相应的正则化因子β,步骤3-3-1的模型转换为:
步骤3-3-3,对步骤3-3-2中的模型进行求导,并用二维快速傅里叶法求得最优解,最优解对应的就是恢复后的清晰图像:
求解d时,阈值公式为:
其中,
本发明利用lq/l2范数作为正则化先验项,采用多尺度方法由粗尺度到细尺度逐步迭代求解,用***法求解模型时,用l1范数保真项更新估计图像的高频信息,在清晰图像恢复时,采用封闭阈值公式以解析解的形式给出,提高了算法速度,同时,更新模糊核时,提出线性递增权重参数,对模糊核按多尺度方法由粗到细逐步进行估计,进一步改善了模糊核的收敛性,提高了去模糊图像质量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明可选实施例的一种盲去模糊的图像恢复方法的流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,一种盲去模糊的图像恢复方法,步骤如下:
步骤1,建立图像模糊模型为:
其中,f表示模糊图像,k表示模糊核矩阵,u表示清晰图像,n表示成像过程中的噪声;
步骤2,利用离散滤波器 产生高频信息构造凹凸范数比值正则化模型为:
其中,k≥0,∑iki=1,ki表示模糊核矩阵k中的元素,λ和λ1为权重参数,ψ为常数值,y为模糊图像的高频信息,x为清晰图像的高频信息,q为方程次方数且0<q<1;
步骤3,求解步骤2的凹凸范数比值正则化模型,得到清晰图像;
步骤3-1,更新清晰图像的高频信息x;
步骤3-1-1,保真项选择l1范数,求解模型为:
步骤3-1-2,将步骤3-1-1中的||x||2视为常数,则步骤3-1-1的求解模型转换为关于非凸||x||q范数的正则化模型:
步骤3-1-3,引入辅助变量v和权重参数θ,并将步骤3-1-2中的正则化模型转换为:
其中,
θi=ciθ0 (7);
式中,θ0表示权重参数θ的初始值,ci与尺度i成线性关系,ci=2i;
步骤3-1-4,采用β代替常数项λ||x||2,并对步骤3-1-3的正则化模型转换中的v和x进行分别求解,求解公式为:
步骤3-1-5,对步骤3-1-4的求解公式进行分别求导:
式中,xk+1表示第k+1步的x值,Δt为迭代步长;
步骤3-2,更新模糊核k;
步骤3-2-1,保真项选择l2范数的平方项,求解模型为:
步骤3-2-2,在模糊核k的更新过程中,利用迭代最小二乘法IRLS方法计算模糊核k的权重:
其中,λ1和ψ为常数值,k0为初始模糊核,wk为更新的模糊核k的权重;
步骤3-2-3,根据步骤3-2-1和步骤3-2-2,在最细尺度求出模糊核函数k;
步骤3-3,得到恢复图像;
步骤3-3-1,在模糊核已知的情况下,图像模型求解变为非盲图像去卷积,求解模型为:
步骤3-3-2,令d=Du,其中D表示运算,将变量u与D分离,同时增加l2范数的平方项,引入相应的正则化因子β,步骤3-3-1的模型转换为:
步骤3-3-3,对步骤3-3-2中的模型进行求导,并用二维快速傅里叶法求得最优解,最优解对应的就是恢复后的清晰图像:
求解d时,阈值公式为:
其中,
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种盲去模糊的图像恢复方法,其特征在于,步骤如下:
步骤1,建立图像模糊模型为:
<mrow>
<mi>f</mi>
<mo>=</mo>
<mi>u</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>+</mo>
<mi>n</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中,f表示模糊图像,k表示模糊核矩阵,u表示清晰图像,n表示成像过程中的噪声;
步骤2,利用离散滤波器产生高频信息构造凹凸范数比值正则化模型为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<munder>
<mi>min</mi>
<mi>x</mi>
</munder>
<mi>&lambda;</mi>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>y</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>|</mo>
<msup>
<mo>|</mo>
<mi>q</mi>
</msup>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munder>
<mi>min</mi>
<mi>k</mi>
</munder>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>y</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mi>&psi;</mi>
<mo>|</mo>
<mo>|</mo>
<mi>k</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中,k≥0,∑iki=1,ki表示模糊核矩阵k中的元素,λ和λ1为权重参数,ψ为常数值,y为模糊图像的高频信息,x为清晰图像的高频信息,q为方程次方数且0<q<1;
步骤3,求解步骤2的凹凸范数比值正则化模型,得到清晰图像;
步骤3-1,更新清晰图像的高频信息x;
步骤3-2,更新模糊核k;
步骤3-3,得到恢复图像。
2.根据权利要求1所述的盲去模糊的图像恢复方法,其特征在于,在步骤3-1中,具体步骤为:步骤3-1-1,保真项选择l1范数,求解模型为:
<mrow>
<munder>
<mi>min</mi>
<mi>x</mi>
</munder>
<mi>&lambda;</mi>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>y</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>|</mo>
<msup>
<mo>|</mo>
<mi>q</mi>
</msup>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
步骤3-1-2,将步骤3-1-1中的||x||2视为常数,则步骤3-1-1的求解模型转换为关于非凸||x||q范数的正则化模型:
<mrow>
<munder>
<mi>min</mi>
<mi>x</mi>
</munder>
<mi>&lambda;</mi>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>y</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
<mo>+</mo>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>|</mo>
<msup>
<mo>|</mo>
<mi>q</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
步骤3-1-3,引入辅助变量v和权重参数θ,并将步骤3-1-2中的正则化模型转换为:
<mrow>
<munder>
<mi>min</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>v</mi>
</mrow>
</munder>
<mi>&lambda;</mi>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>v</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
<mo>+</mo>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>|</mo>
<msup>
<mo>|</mo>
<mi>q</mi>
</msup>
<mo>+</mo>
<mfrac>
<mi>&theta;</mi>
<mn>2</mn>
</mfrac>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>y</mi>
<mo>-</mo>
<mi>v</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中,
<mrow>
<mi>v</mi>
<mo>=</mo>
<mi>x</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>y</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
θi=ciθ0 (7);
式中,θ0表示权重参数θ的初始值,ci与尺度i成线性关系,ci=2i;
步骤3-1-4,采用β代替常数项λ||x||2,并对步骤3-1-3的正则化模型转换中的v和x进行分别求解,求解公式为:
<mrow>
<munder>
<mi>min</mi>
<mi>v</mi>
</munder>
<mi>&beta;</mi>
<mo>|</mo>
<mo>|</mo>
<mi>v</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
<mo>+</mo>
<mfrac>
<mi>&theta;</mi>
<mn>2</mn>
</mfrac>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>y</mi>
<mo>-</mo>
<mi>v</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<munder>
<mi>min</mi>
<mi>v</mi>
</munder>
<mo>|</mo>
<mo>|</mo>
<mi>v</mi>
<mo>|</mo>
<msup>
<mo>|</mo>
<mi>q</mi>
</msup>
<mo>+</mo>
<mfrac>
<mi>&theta;</mi>
<mn>2</mn>
</mfrac>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>y</mi>
<mo>-</mo>
<mi>v</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
步骤3-1-5,对步骤3-1-4的求解公式进行分别求导:
<mrow>
<mi>v</mi>
<mo>=</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mi>x</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>y</mi>
<mo>|</mo>
<mo>-</mo>
<mfrac>
<mi>&beta;</mi>
<mi>&theta;</mi>
</mfrac>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
<mrow>
<msup>
<mi>x</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>=</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mo>|</mo>
<msup>
<mi>x</mi>
<mi>k</mi>
</msup>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mi>&theta;</mi>
<mo>(</mo>
<mi>x</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>y</mi>
<mo>-</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>k</mi>
<mi>T</mi>
</msup>
<mo>|</mo>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mi>q</mi>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>|</mo>
<msup>
<mo>|</mo>
<mrow>
<mi>q</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
<mo>&CenterDot;</mo>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mo>&lsqb;</mo>
<msup>
<mi>x</mi>
<mi>k</mi>
</msup>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mi>&theta;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>y</mi>
<mo>-</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mi>k</mi>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
式中,xk+1表示第k+1步的x值,Δt为迭代步长。
3.根据权利要求1所述的盲去模糊的图像恢复方法,其特征在于,在步骤3-2中,具体步骤为:
步骤3-2-1,保真项选择l2范数的平方项,求解模型为:
<mrow>
<munder>
<mi>min</mi>
<mi>k</mi>
</munder>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>y</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mi>&psi;</mi>
<mo>|</mo>
<mo>|</mo>
<mi>k</mi>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
步骤3-2-2,在模糊核k的更新过程中,利用IRLS方法计算模糊核k的权重:
其中,λ1和ψ为常数值,k0为初始模糊核,wk为更新的模糊核k的权重;
步骤3-2-3,根据步骤3-2-1和步骤3-2-2,在最细尺度求出模糊核函数k。
4.根据权利要求1所述的盲去模糊的图像恢复方法,其特征在于,在步骤3-2中,具体步骤为:
步骤3-3-1,在模糊核已知的情况下,图像模型求解变为非盲图像去卷积,求解模型为:
<mrow>
<munder>
<mi>min</mi>
<mi>u</mi>
</munder>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mn>2</mn>
</mfrac>
<mo>|</mo>
<mo>|</mo>
<mi>u</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>f</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mo>|</mo>
<mo>|</mo>
<mo>&dtri;</mo>
<mi>u</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mi>q</mi>
<mi>q</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
步骤3-3-2,令d=Du,其中D表示运算,将变量u与D分离,同时增加l2范数的平方项,引入相应的正则化因子β,步骤3-3-1的模型转换为:
<mrow>
<munder>
<mi>min</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>d</mi>
</mrow>
</munder>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mn>2</mn>
</mfrac>
<mo>|</mo>
<mo>|</mo>
<mi>u</mi>
<mo>&CircleTimes;</mo>
<mi>k</mi>
<mo>-</mo>
<mi>f</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mfrac>
<mi>&beta;</mi>
<mn>2</mn>
</mfrac>
<mo>|</mo>
<mo>|</mo>
<mi>D</mi>
<mi>u</mi>
<mo>-</mo>
<mi>d</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mo>|</mo>
<mo>|</mo>
<mi>d</mi>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mi>q</mi>
<mi>q</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
步骤3-3-3,对步骤3-3-2中的模型进行求导,并用二维快速傅里叶法求得最优解,最优解对应的就是恢复后的清晰图像:
求解d时,阈值公式为:
<mrow>
<mi>d</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mi>A</mi>
<mo>|</mo>
<mo>+</mo>
<msqrt>
<mrow>
<mfrac>
<mrow>
<mn>2</mn>
<mo>|</mo>
<mi>a</mi>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<mi>A</mi>
<mo>|</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>|</mo>
<mi>A</mi>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>a</mi>
<mo>></mo>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<mi>a</mi>
<mo>|</mo>
<mo>></mo>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mi>A</mi>
<mo>|</mo>
<mo>+</mo>
<msqrt>
<mrow>
<mfrac>
<mrow>
<mn>2</mn>
<mo>|</mo>
<mi>a</mi>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<mi>A</mi>
<mo>|</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>|</mo>
<mi>A</mi>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>a</mi>
<mo>></mo>
<mo>-</mo>
<mi>p</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中,
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2017107967298 | 2017-09-06 | ||
CN201710796729 | 2017-09-06 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107767351A true CN107767351A (zh) | 2018-03-06 |
Family
ID=61273146
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711069974.5A Pending CN107767351A (zh) | 2017-09-06 | 2017-11-03 | 一种盲去模糊的图像恢复方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107767351A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113191413A (zh) * | 2021-04-25 | 2021-07-30 | 华中科技大学 | 基于中央凹残差网络的***多模mr图像分类方法和*** |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103761710A (zh) * | 2014-01-08 | 2014-04-30 | 西安电子科技大学 | 基于边缘自适应的高效图像盲去模糊方法 |
CN104112261A (zh) * | 2014-07-17 | 2014-10-22 | 五邑大学 | 基于范数比值正则化的快速图像盲去模糊方法 |
CN104599242A (zh) * | 2014-12-09 | 2015-05-06 | 西安电子科技大学 | 使用多尺度非局部正则的模糊核估计方法 |
CN106097256A (zh) * | 2016-05-31 | 2016-11-09 | 南京邮电大学 | 一种基于图像盲去模糊的视频图像模糊度检测方法 |
CN106204472A (zh) * | 2016-06-30 | 2016-12-07 | 北京大学 | 基于稀疏特性的视频图像去模糊方法 |
CN106875349A (zh) * | 2016-12-30 | 2017-06-20 | 无锡高新兴智能交通技术有限公司 | 盲图像复原方法中模糊核的计算方法及盲图像复原方法 |
-
2017
- 2017-11-03 CN CN201711069974.5A patent/CN107767351A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103761710A (zh) * | 2014-01-08 | 2014-04-30 | 西安电子科技大学 | 基于边缘自适应的高效图像盲去模糊方法 |
CN104112261A (zh) * | 2014-07-17 | 2014-10-22 | 五邑大学 | 基于范数比值正则化的快速图像盲去模糊方法 |
CN104599242A (zh) * | 2014-12-09 | 2015-05-06 | 西安电子科技大学 | 使用多尺度非局部正则的模糊核估计方法 |
CN106097256A (zh) * | 2016-05-31 | 2016-11-09 | 南京邮电大学 | 一种基于图像盲去模糊的视频图像模糊度检测方法 |
CN106204472A (zh) * | 2016-06-30 | 2016-12-07 | 北京大学 | 基于稀疏特性的视频图像去模糊方法 |
CN106875349A (zh) * | 2016-12-30 | 2017-06-20 | 无锡高新兴智能交通技术有限公司 | 盲图像复原方法中模糊核的计算方法及盲图像复原方法 |
Non-Patent Citations (1)
Title |
---|
余义斌 等: "凹凸范数比值正则化的快速图像盲去模糊", 《电子学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113191413A (zh) * | 2021-04-25 | 2021-07-30 | 华中科技大学 | 基于中央凹残差网络的***多模mr图像分类方法和*** |
CN113191413B (zh) * | 2021-04-25 | 2022-06-21 | 华中科技大学 | 基于中央凹残差网络的***多模mr图像分类方法和*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103761710B (zh) | 基于边缘自适应的高效图像盲去模糊方法 | |
CN108921786B (zh) | 基于残差卷积神经网络的图像超分辨率重构方法 | |
WO2020015167A1 (zh) | 一种基于融合网络的图像超分辨率及去非均匀模糊方法 | |
Pan et al. | Fast $\ell^{0} $-Regularized Kernel Estimation for Robust Motion Deblurring | |
CN108154487B (zh) | 基于多通道分解的屏摄图像摩尔纹消除方法 | |
CN109035146B (zh) | 一种基于深度学习的低质量图像超分方法 | |
CN104103052B (zh) | 一种基于稀疏表示的图像超分辨率重建方法 | |
CN104112261A (zh) | 基于范数比值正则化的快速图像盲去模糊方法 | |
CN108564544A (zh) | 基于边缘感知的图像盲去模糊组合稀疏优化方法 | |
CN112164011A (zh) | 基于自适应残差与递归交叉注意力的运动图像去模糊方法 | |
CN109360157B (zh) | 基于tv和小波正则化的空间变化模糊图像复原方法 | |
Yoshikawa et al. | Super resolution image reconstruction using total variation regularization and learning-based method | |
CN107292842A (zh) | 基于先验约束和离群值抑制的图像去模糊方法 | |
CN104021529A (zh) | 一种模糊图像非盲复原方法 | |
CN111986122B (zh) | 基于混合型总变分正则化的模糊图像非盲复原方法 | |
CN105590296B (zh) | 一种基于双字典学习的单帧图像超分辨率复原方法 | |
CN107194889A (zh) | 分块双边全变分正则化图像噪声消除方法 | |
CN114092416A (zh) | 一种dr模糊图像盲反卷积复原方法及*** | |
CN107146202B (zh) | 基于l0正则化和模糊核后处理的图像盲去模糊的方法 | |
CN107767351A (zh) | 一种盲去模糊的图像恢复方法 | |
CN111986121B (zh) | 基于Framelet l0范数约束的模糊图像非盲复原方法 | |
CN111626943A (zh) | 基于一阶前向后向算法的全变分图像去噪方法 | |
CN107590781B (zh) | 基于原始对偶算法的自适应加权tgv图像去模糊方法 | |
CN107292836B (zh) | 基于外部图像块先验信息与稀疏表示的图像盲去模糊方法 | |
CN101567079B (zh) | 基于Hopfield神经网络的运动模糊图像恢复方法 |
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: 20180306 |
|
RJ01 | Rejection of invention patent application after publication |