(19)中华人民共和国国家知识产权局
(12)发明专利说明书
(10)申请公布号 CN 108629741 A
(43)申请公布日 2018.10.09
(21)申请号 CN201810252207.6
(22)申请日 2018.03.26
(71)申请人 中南大学
    地址 410083 湖南省长沙市岳麓区麓山南路932号
(72)发明人 谢永芳 张骞 桂卫华 徐德刚 蒋朝辉 唐朝晖
(74)专利代理机构 北京路浩知识产权代理有限公司
    代理人 王莹
(51)Int.CI
     
                                                                  权利要求说明书 说明书 幅图
(54)发明名称
      一种基于L0和L1正则项的模糊核估计方法
(57)摘要
l1正则化的作用
      本发明提供了一种基于L0和L1正则项的模糊核估计方法,利用L0正则项作为中间清晰图像求解时的约束条件,有效抑制振铃效应;利用L1正则项作为估计模糊核求解时的正则化约束条件,使得估计模糊核足够稀疏,根据建立的模糊核估计模型,将基于L0和L1正则项的模糊核估计方法,应用于图像金字塔理论中,由粗到细地从图像金字塔的顶层至底层求解估计模糊核,将所求得的模糊核用于图像非盲复原算法,复原出最终清晰图像。本发明能够较为准确地估计出模糊图像的模糊核,从而用估计出的模糊图像的模糊核,将初始模糊图像复原出最终清晰图像。
法律状态
法律状态公告日
法律状态信息
法律状态
权 利 要 求 说 明 书
1.一种基于L0和L1正则项的模糊核估计方法,其特征在于,包括:
S1、确立模糊核估计模型,输入初始模糊图像,对所述初始模糊图像进行下采样至图像金字塔的顶层,得到初始模糊核;
S2、根据所述模糊核估计模型、所述初始模糊图像和所述初始模糊核,得到中间清晰图像;
S3、根据所述模糊核估计模型、所述初始模糊图像和所述中间清晰图像,得到估计模糊核;
S4、判断是否达到所述图像金字塔的底层;若是,则得到最终模糊核;若否,则重复所述步骤S2和S3,对所述中间清晰图像和所述估计模糊核进行上采样至所述图像金字塔的下一层,直至达
到所述图像金字塔的底层并得到所述最终模糊核;
S5、将所述最终模糊核用于图像非盲复原算法,复原出最终清晰图像;
其中,所述步骤S2中,选用L0正则项作为所述中间清晰图像的正则化约束条件;所述步骤S3中,选用L1正则项作为所述估计模糊核的正则化约束条件;所述步骤S4中,所述图像金字塔的底层的层数为指定迭代次数,所述指定迭代次数为所述步骤S1中输入所述初始模糊图像时设置的预设值。
2.根据权利要求1所述的一种基于L0和L1正则项的模糊核估计方法,其特征在于,所述步骤S1中,所述模糊核估计模型为:
<Image>
其中,y为初始模糊图像,k为模糊核,*为卷积算子,x为中间清晰图像,p(x)为中间清晰图像的先验知识,p(k)为模糊核的先验知识,||k*x-y||为数据拟合项。
3.根据权利要求2所述的一种基于L0和L1正则项的模糊核估计方法,其特征在于,所述模糊核估
计模型根据图像模糊模型和贝叶斯理论得到,所述图像模糊模型为y=k*x+n,所述贝叶斯理论为p(x,k|y)∝p(y|x,k)p(x)p(k);
其中,y为初始模糊图像,k为模糊核,*为卷积算子,x为中间清晰图像,p(x)为中间清晰图像的先验知识,p(k)为模糊核的先验知识。
4.根据权利要求1所述的一种基于L0和L1正则项的模糊核估计方法,其特征在于,所述步骤S5中,所述图像金字塔的底层值为:
N=floor(2log(5/k<Sub>s</Sub>)/log0.5)+1
其中,N为图像金字塔的底层值,k<Sub>s</Sub>为模糊核的尺寸。
5.根据权利要求1所述的一种基于L0和L1正则项的模糊核估计方法,其特征在于,中间清晰图像的求解方程为:
<Image>
其中,y为初始模糊图像,k为模糊核,*为卷积算子,x为中间清晰图像,<Image>为L<Sub>0</Sub
>正则项,λ为L<Sub>0</Sub>正则项参数,||k*x-y||为数据拟合项。
6.根据权利要求5所述的一种基于L0和L1正则项的模糊核估计方法,其特征在于,中间清晰图像的方程基于二次分裂的方法求解,引入辅助变量g=(gh,gv)<Sup>T</Sup>,使得中间清晰图像的求解方程为:
<Image>
其中,y为初始模糊图像,k为模糊核,*为卷积算子,x为中间清晰图像,<Image>为x的L<Sub>0</Sub>正则项,λ为x的L<Sub>0</Sub>正则项参数,||k*x-y||为数据拟合项,g为辅助变量,g=(g<Sub>h</Sub>,g<Sub>v</Sub>)<Sup>T</Sup>。
当g初始化为0时,中间清晰图像x的求解方程为:
<Image>
其中,F<Sup>-1</Sup>(·)、F(·)分别代表快速傅里叶反变换和快速傅里叶变换,<Image>为F(·)的复共轭算子,<Image><Image><Image>分别为水平微分算子和垂直微分算子。
当x确定时,通过逼近估计得到的g的解为:
<Image>
其中,λ为x的L<Sub>0</Sub>正则化参数。
7.根据权利要求1所述的一种基于L0和L1正则项的模糊核估计方法,其特征在于,模糊核的求解方程为:
其中,y为初始模糊图像,k为模糊核,*为卷积算子,x为中间清晰图像,||k||<Sub>1</Sub>为k的L<Sub>1</Sub>正则项,β为k的L<Sub>1</Sub>正则项参数,||k*x-y||为数据拟合项。
8.根据权利要求1所述的一种基于L0和L1正则项的模糊核估计方法,其特征在于,得到所述估计模糊核之后,对模糊核进行优化处理,所述优化处理包括:
对模糊核进行非负处理:将模糊核中负值元素设为0;
对模糊核进行连续性处理:保留模糊核的主要结构,剔除分支,使得模糊核的非0区域为一个连通域;
对模糊核进行归一化处理:对模糊核进行归一化处理,使得∑<Sub>i</Sub>k<Sub>i</Sub>=1。

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。