(19)中华人民共和国国家知识产权局
(12)发明专利说明书 | ||
(10)申请公布号 CN 102208100 A 正则化残差(43)申请公布日 2011.10.05 | ||
(21)申请号 CN201110144242.4
(22)申请日 2011.05.31
(71)申请人 重庆大学
地址 400030 重庆市沙坪坝区沙正街174号
(72)发明人 李伟红 李权利 龚卫国 唐述 李正浩 杜兴
(74)专利代理机构 重庆华科专利事务所
代理人 康海燕
(51)Int.CI
G06T5/00
权利要求说明书 说明书 幅图 |
(54)发明名称
基于Split Bregman 迭代的全变差正则化图像盲复原方法 | |
(57)摘要
本发明是一种基于SplitBregman迭代的全变差正则化图像盲复原方法，属于图像处理技术领域。其核心是通过图像先验信息及正则化的特点建立TV正则化盲复原代价函数，利用算子分裂技术将最小化问题转化为新的约束求解问题，通过加入惩罚项的方法将约束求解问题转化为分裂的代价函数，进而提出了扩展的SplitBregman迭代求解框架对分裂的代价函数进行求解。实验表明本发明方法能够有效、快速地复原图像，能针对多种模糊类型进行复原，克服传统TV正则化盲复原方法对复杂模糊类型或复杂图像复原效果不佳的缺点，且该发明方法对噪声具有很好的鲁棒性。 | |
法律状态
法律状态公告日 | 法律状态信息 | 法律状态 |
权 利 要 求 说 明 书
1.一种基于Split Bregman迭代的全变差正则化图像盲复原方法，该方法包括以下步骤：<br>（1）定义TV正则化盲复原最小化代价函数，采用算子替换的方法将盲去卷积模型中的梯度算子进行替换： <Image>，<Image>，将最小化问题转化为约束优化问题，其中，<Image>和<Image>分别为图像<Image>和点扩散函数<Image>的梯度，<Image>和<Image>为引入的两个辅助算子；<br>（2）引入惩罚项对<Image>和<Image>进行惩罚，将步聚（1）中的约束求解问题转化为分裂最小化代价函数；<br>（3）采用扩展的Split Bregman迭代最小化方法对步聚（2）中的分裂最小化代价函数进行求解，经过迭代最终复原出清晰图像<I>u</I>。<br><br>2.根据权利要求1所述的基于Split Bregman迭代的全变差正则化图像盲复原方法，其特征在于步骤（1）中：<br>所述TV
正则化盲复原最小化代价函数定义如下：<br><Image></Image><br>其中，<I>u</I>原清晰图像，<I>k</I>点扩散函数，<I>f</I>为已知的退化图像；<Image>和<Image>为两个大于0的正则化参数，控制图像和点扩散函数的正则化程度；<Image>和<Image>为TV范数，∑表示所有像素相加；<I>u</I><Sub><I>x</I></Sub>和<I>u</I><Sub><I>y</I></Sub>分别代表图像<I>u</I>在水平方向和垂直方向的一阶微分，在位置(<I>i</I>, <I>j</I>)处水平和垂直方向的一阶微分定义分别为：<br><Image></Image><br><I>k</I><Sub><I>x</I></Sub>和<I>k</I><Sub><I>y</I></Sub>分别代表点扩散函数<I>k</I>在水平方向和垂直方向的一阶微分，在位置(<I>i</I>, <I>j</I>)处的定义分别为：<br><Image></Image><br>引入两个辅助算子<I>b</I><Sub>1</Sub>和<I>b</I><Sub>2</Sub>，采用算子替换的方法将TV正则化盲去卷积模型中的梯度算子进行替换：<Image>，<Image>，将最小化问题转化为以下约束优化问题：<br><Image></Image>。<br><br><Claim>3.根据权利要求1或2所述的基于Split Bregman迭代的全变差正则化图像盲复原方法，其特征在于，所述步骤（2）中引入两个惩罚项对<Image>和<Image>进行惩罚，将步骤（1）约束求解问题转化为新的分裂的最小化代价函数如下：<br><Image></Image><br>其中后两项为惩罚项，<Image>和<Image>为大于0的惩罚参数，控制两个惩罚项的权重。<br><br>4.根据权利要求3所述的基于Split Br
egman迭代的全变差正则化图像盲复原方法，其特征在于，所述步骤（3）中为解决步骤（2）中分裂最小化代价函数，采用扩展的Split Bregman迭代求解框架如下：<br><Image></Image><br>其中，<I>t</I><Sub>1</Sub>和<I>t</I><Sub>2</Sub>为Bregman 迭代方法中强约束所引入的两个辅助变量。<br><br>5.根据权利要求4所述的基于Split Bregman迭代的全变差正则化图像盲复原方法，其特征在于，所述扩展的Split Bregman迭代求解框架分解为五个子问题进行交替最小化，在复原过程中，采用图像与PSF的非负约束条件，并假定PSF是归一化的，且其支持域的大小是已知的，所述方法采用复原过程中相邻的两次迭代所估计得的图像的相对差异<Image>，为较小的正数，作为迭代终止条件，或通过设定最大迭代次数来终止迭代；通过交替迭代最优化，最终复原出清晰的图像<I>u</I>。<br><br>
说 明 书
技术领域
本发明属于图像处理技术领域。
背景技术
图像是人们最主要的信息源之一，然而在图像的获取、传输等过程中，由于各种因素的干扰，会造成图像的退化降质。图像的退化会使大量的真实信息丢失，不仅会降低图像的科学价值，而且也会带来巨大的经济损失。因此，我们需要利用图像复原技术从退化的图像复原出原本面貌。目前，图像复原技术已应用到众多科学与技术领域，如天文观测、医学成像、多媒体、刑事侦察等。众多图像复原方法要求先验信息较多，或存在效果较差，算法复杂度高等缺点。至今，研究出有效、快速的图像复原方法仍是图像处理领域中最具有挑战性的难题之一。
图像复原技术的发展经历了约40年的历史，一些经典的复原方法假定造成图像模糊的点扩散函数(Point spread function, PSF)是已知的，如逆滤波、维纳滤波、R-L方法等等。但在实际中点扩散函数往往是未知的，因此图像的盲复原技术得到了大量研究。目前，盲复原方法有多种，按估计点扩散函数和图像的顺序可分为先验辨识法和联合辨识法。先验辨识法是首先估计出点扩散函数，然后根据所得到的PSF利用经典的非盲复原方法复原。先验辨识法只能针对一些较简单的模糊类型进行复原，且受噪声干扰较大。联合辨识法是采用迭代方式交替估计点扩散函数和图像，该方法又分为参数法和非参数法，参数法是把模糊过程看作具有一定参数的数学模型，通过估计模型中的参数得到复原图像。而非参数法不需要模糊过程的参
数化模型。
为了克服图像复原这一逆过程的病态性，许多正则化方法引入到了复原方法中。其中两个著名的正则化方法是Tikhonov正则化和全变差(Total variation, TV)正则化。 Tikhonov正则化具有各向同性的特点，缺点是会使图像过度平滑，对图像细节的保护不够，而全变差正则化方法因弥补了Tikhonov正则化方法的不足而受到了众多研究者的关注。全变差正则化盲复原方法利用全变差正则化的特点，将全变差范数作为正则项引入所构造的代价函数中，使算法克服病态性，且使求解过程变得稳定以得到较好的复原结果。
传统的全变差方法采用固定点法迭代求代价函数的最优解，针对一些背景简单的图像有较好的复原效果，但对复杂的图像或针对复杂的模糊类型效果不佳。本发明通过建立盲复原代价函数，采用Split Bregman（译为分裂布雷格曼，一种<I>L</I><Sub>1</Sub>正则化代价函数的优化方法）迭代方法进行最优化求解，很好地克服了全变差正则化盲复原方法的缺点。基于Split Bregman迭代的全变差正则化盲复原方法是我们提出的新的盲复原方法。
发明内容
本发明的目的在于提出一种针对多种模糊类型的退化图像进行复原的有效、快速图像复原方法，力求从已知的退化图像恢复出原本面貌，使模糊的图像变得清晰，旨在改善图像的质量。
本发明是基于Split Bregman迭代的全变差正则化盲复原方法，克服传统全变差盲复原方法对复杂模糊类型或复杂图像复原效果不佳的缺点，并使复原方法对噪声具有很好的鲁棒性。具体内容为将算子分裂技术引入全变差正则化盲复原中，采用算子分裂的方法对代价函数中的算子进行替换，进而提出新的约束求解问题。然后运用惩罚项的方法将所提出的约束求解问题变成新的无约束的分裂代价函数。进一步运用提出的扩展的Split Bregman迭代方法对分裂的代价函数进行交替求解获得原始图像的估计，从而对退化图像进行有效、快速的复原。
该盲复原方法的实现步骤如下：
（1）定义TV正则化盲复原最小化代价函数，采用算子替换的方法将盲去卷积模型中的梯度算子进行替换： <Image>，<Image>，将最小化问题转化为约束优化问题；
（2）引入惩罚项对<Image>和<Image>进行惩罚，将步聚（1）中的约束问题转化为分裂的最小化代价函数；
（3）采用扩展的Split Bregman迭代最优化方法对步聚（2）中的分裂最小化代价函数进行最优化求解，通过迭代运算最终复原出原清晰图像<I>u</I>。
版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。
发表评论