非盲图像复原综述

非盲图像复原综述

摘要：作为目前图像处理领域的研究重点，图像复原可移除图像中的模糊与噪声，具有重要的理论价值和广阔的应用前景。为使图像复原的研究被人们所了解，该文首先对图像复原做了简单的描述，接着介绍了近年来出现的一些非盲图像复原算法，包括基于总变分模型的算法、基于bregman迭代的算法和基于稀疏表示的算法等，最后基于对现有算法的了解与分析，总结了图像复原研究的难点与趋势。

关键词：图像复原；总变分模型； bregman迭代；稀疏表示；优化问题

中图分类号：tp391 文献标识码：a 文章编号：1009-3044（2013）07-1642-03

由于噪声、模糊等不利因素的影响，数字图像的质量通常难以令人满意，无法对其进一步进行研究和利用。观测到的退化与理想的原始图像之间的关系可表示为：

式中，y表示观测图像；a表示模糊算子；x表示原始图像；n表示加性噪声。因此，图像复原的目的是在模型（1）的框架下，估计最优的原始图像x，这是一种典型的线性逆问题。图像复原已成为图像处理领域乃至计算机领域的研究热点，其研究可追溯到上个世纪六十年代，经过逾半个世纪的发展，不断有新的算法和技术涌现。目前，图像复原算法主要分为非盲（non-blind）图像复原算法和图像盲复原算法两大类。当模型（1）中的模糊算子a是已知

的，图像复原被称为非盲的图像复原。非盲图像复原算法主要有：基于总变分（tv，total variation）模型的算法、基于bregman

方法的算法和基于稀疏表示的算法等。

1 基于tv模型的图像复原算法

tv模型亦称rof模型，可表示为：

由l. i. rudin、s. osher和e. fatemi提出[1]，其最突出的优点是抑噪的同时可保留图像边缘等重要信息，该模型已成为图像复原、图像去噪[2]和图像修复[3]等领域使用最广泛的先验模型之一。由于tv模型是不可微分的，基于tv模型的算法需重点考虑图像复原的数值解问题。为求解经典的tv/l2问题，y. l. wang等人[4]建立了一种半二次化（quadratic minimization）模型，采用二次最小化方法计算该模型以提高图像复原的效率。j. f. wang等人[5]扩展了y. l. wang等人的算法，将其用于彩色图像复原问题的处理。a. beck等人[6]提出了一种基于nesterov梯度法的tv图像复原算法，该算法具有全局收敛速率。基于变量分离和交替方向法（adm，alternating direction method），m. k. ng等人[7]将tv/l2模型分解为更简单的子问题进行求解，adm方法的优点是通过有限步迭代即可得到问题的精确解。结合增广拉格朗日法，s. h. chan等人[8]提出了与m. k. ng等人类似的算法，该算法主要用于复原视频图像序列。尽管tv模型可保留图像细节，但复原图像的平滑区域产生阶梯效应（staircase effect），因此一些研究提出在tv模型中加入高阶项以克服其缺点。例如，f. li等人[9]将tv

滤波器和四阶pde滤波器结合复原图像，避免了阶梯效应的出现。原始-对偶法（primal-dual method）是一类求解经典tv图像复原比较有效的方法。它最早由t. f. chan等人提出[10]，为了消除模型（2）目标函数的奇异性和不可微性，t. f. chan等人用新变量替换了其欧拉-拉格朗日方程中的?x/|?x|，得到如下的模型：

及其对偶模型，并通过半光滑牛顿法（semismooth newton method）对模型（4）及其对偶模型进行求解。m. hintermuller等人的算法具有超线性的收敛速度，相比t. f. chan等人的算法，该算法对对偶变量没有过多的约束。e. esser等人[12]修改了原始对偶混合梯度（pdhg，primal-dual hybrid gradient）算法，将tv最小化问题转化为等价的鞍点问题（saddle point problem），并利用不精确uzawa算法求解该鞍点问题。由最终的实验结果，e. esser等人的算法略优于pdhg算法。a. chambolle等人[13]将tv 优化问题转化为更通用的鞍点问题，采用了与pdhg类似的算法处理该鞍点问题。针对不同的具体问题，a. chambolle等人的算法有着不同的收敛速率。

2 基于bregman迭代的图像复原算法

近年来，在图像复原领域，bregman迭代因其简单、稳定、速度快和效率高而受到越来越多的关注。bregman迭代是一系列以bregman距离为基础的方法的总称，包括了经典的bregman方法、线性（linerized）bregman方法和分裂（split）bregman方法等。bregman迭代方法的基本思想是将优化问题分解为等价的非约束优

化子问题，其中某些子问题的目标函数由bregman距离定义。该方法最早s. osher等人[14]引入图像复原领域，用以改进传统方法对tv模型的处理。为提高经典bregman方法的性能，将该方法中的二次项0.5×||ax-y||2 2替换为，并加入误差项0.5×μ

||x-xk||2 2，w. t. yin等人提出了线性bregman方法[15]，并将其应用于压缩传感问题[16]。为适应模糊矩阵的变化，在对线性bregman方法修改的基础上，结合稀疏表示技术，j. f. cai等人[17]提出了一种紧框架域的图像复原算法，该算法可以快速找到图像复原问题的稀疏解。为进一步提高经典bregman方法的性能，t. goldstein等人[18]提出了适用更通用的l1正则化问题的分裂bregman方法。它可看作bregman方法和算子分裂（operator splitting）相结合的一种方法，且与douglas-rachford分裂方法和向前-向后（forward-backward）分裂方法本质上属于同一种方法[19]。j. f. cai等人[20]将分裂bregman方法用于基于分析的图像复原问题，分析算子为紧框架小波。相比传统的罚函数法和连续法（continuation method），该算法收敛速度更快，更稳定，参数在迭代过程保持不变。s. setzer等人[21]为泊松模糊图像复原建立新的能量泛函，它由i-散度（divergence）和tv正则项构成，分裂bregman方法用以求解该能量泛函。该算法无需内部迭代，相比同类算法更加高效。针对基于有界总变分的图像复原问题，x. w. liu等人[22]引入扩展的bregman迭代方法以快速获得最优解。除此之外，x. q. zhang等人[23]还提出了一种bregman化的算法分