一种自适应参数更新的贝叶斯图像反卷积恢复方法
【专利摘要】本发明属于图像处理技术领域，公开了一种自适应参数更新的贝叶斯图像反卷积恢复方法，包括：观测模型建模；图像先验建模；求观测模型的下界：利用优化最小化方法获得观测模型的下界函数；求图像先验模型的下界：利用优化最小化方法获得图像先验模型的下界函数；利用最大化后验概率迭代更新观测模型超参数、图像先验超参数和图像；本发明提供的自适应参数更新的贝叶斯图像反卷积恢复方法自适应更新能力好，计算代价低，恢复效果显著，具有去除高斯噪声和椒盐噪声的能力，可应用于在高斯噪声或椒盐噪声污染情形下且受高斯模糊退化图像的恢复。
【专利说明】
-种自适应参数更新的贝叶斯图像反卷积恢复方法
技术领域
[0001] 本发明设及图像处理技术领域，尤其设及一种自适应参数更新的贝叶斯图像反卷 积恢复方法，用于当图像受到高斯噪声或者椒盐噪声污染的情形下，对图像进行贝叶斯反 卷积恢复。
【背景技术】
[0002] 图像反卷积是计算机视觉和图像处理中的一个基础和关键问题。在实际成像过程 中，由于受端流、运动、抖动等多种因素的影响，观测图像普遍包含端流、运动、抖动等模糊 退化，如相机在成像过程中由于手的抖动导致的运动模糊。
[0003] 此外，由于电磁干扰及通信系统的故障和缺陷等因素，图像数据常常受到噪声的 干扰，例如电视信号会由于大气层的干扰而产生脉冲噪声(椒盐噪声）、电子电路的热效应 和光电传感器的光子波动等会产生加性噪声（如高斯噪声），等等。很多成像应用系统的处 理结果很大程度上都依赖于反卷积后图像的质量，所W图像反卷积预处理是十分必须的、 也是很多图像预处理应用系统的一个重要组成部分。
[0004] 现有变分正则化图像反卷积方法往往能有效去除噪声和模糊，但方法中设及的关 键正则化参数需要手动选取，方法耗时且不实用。近年来，变分贝叶斯方法已经用于图像反 卷积，现有贝叶斯图像反卷积方法往往仅仅去除高斯噪声和模糊，性能受限，如D. Babacan 等人于2008年提出的变分贝叶斯总变分图像反卷积方法仅仅能处理高斯噪声和模糊退化， 对脉冲噪声污染退化无能为力。
【发明内容】
[0005] 针对上述问题，本发明的目的在于提供一种自适应参数更新的贝叶斯图像反卷积 恢复方法，既能处理高斯噪声又能处理脉冲噪声，解决统一模型下多种噪声和模糊的有效 去除问题。
[0006] 本发明实现的技术思路为：先将观测模型建模为广义高斯分布;再将图像先验模 型建模为广义拉普拉斯分布，通过优化最小化方法分别求他们的下界，最后通过最大化后 验概率迭代自适应更新超参数和求解图像，具有去除高斯噪声和椒盐噪声和自适应参数更 新能力。
[0007] 为达到上述目的，本发明的实施例采用如下技术方案予W实现。
[000引一种自适应参数更新的贝叶斯图像反卷积恢复方法，所述方法包括如下步骤：
[0009] 步骤1，获取观测图像y，所述观测图像包含高斯噪声和高斯模糊，或者所述观测图 像包含椒盐噪声和高斯模糊;所述观测图像y = Nim化X) +n，其中n表示高斯噪声，Nim表示椒 盐噪声，H表示模糊核函数矩阵，X是待恢复的清晰图像，y是观测图像；
[0010] 步骤2,采用零均值的广义高斯分布对观图像y中的噪声进行建模，从而得到观测 图像y的条件概率分布函数P(y|e，x)，p(〇| ?)表示已知.时O的条件概率分布，0为多变量 广义高斯观测模型的超参数，朗g从均匀分布；
[0011] 步骤3,对所述待恢复的清晰图像X进行先验建模，得到待恢复的清晰图像X的先验 概率分布函数p(x|a)，a为图像先验模型超参数，a服从均匀分布；
[0012] 步骤4,在分层贝叶斯模型下，构建联合全局概率分布p(a，e，x，y):
[0013] p(a，0，x，y)=p(a)p(0)p(y |0，x)p(x|a)
[0014]步骤5,计算观测图像y的条件概率分布函数p(y |0，x)的下界Mi(0，x，y) :p(y |0，x) >]/[1(0，义，7);计算待恢复的清晰图像义的先验概率分布函数9^|0)的下界12(0，义）：9^|口） >M2(a，x);
[0015] 步骤6,根据观测图像y的条件概率分布函数p(y|e，x)的下界、待恢复的清晰图像X 的先验概率分布函数p(x|a)的下界，得到所述联合全局概率分布p(a，e，x，y)的下界P(Q)P (0)Mi(0，x，y)M2(a，x):
[0016] p(a，0，x，y) =p(a)p(0)p(y |0，x)p(x |a)>p(a)p(0)Mi化，x，y)M2(a，x)
[0017] 步骤7,求解后验分布9(〇，0，^|7)=9(〇，0，^，7)/9(7)的下界，得到关于图像先验 模型超参数〇、多变量广义高斯观测模型的超参数e、待恢复的清晰图像X的泛函；
[0018] 步骤8,迭代求解关于图像先验模型超参数a、多变量广义高斯观测模型的超参数 e、待恢复的清晰图像X的泛函，直到满足迭代停止条件，并将最后一次迭代得到的待恢复的 清晰图像作为反卷积恢复的清晰图像。
[0019] 本发明技术方案的有益效果：
[0020] (1)本发明提供的贝叶斯图像反卷积方法，将观测模型和图像先验模型分别建模 为广义高斯分布和广义拉普拉斯分布，具有更大兼容性，既能处理高斯噪声又能处理椒盐 噪声；（2)本发明提供的贝叶斯图像反卷积方法，利用优化最小化方法分别求观测模型和图 像先验模型概率分布的下界，通过最大化后验概率，得到观测模型和图像先验模型超参数 的自适应更新公式，相比于现有文献中手动选取参数，本发明提供的方法在迭代过程中自 适应更新超参数，具有更强的实用价值;有利于贝叶斯图像反卷积方法走向实用。
【附图说明】
[0021] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可W 根据运些附图获得其他的附图。
[0022] 图1为本发明实施例提供的贝叶斯图像反卷积恢复方法的流程示意图；
[0023] 图2为本发明实施例中的退化观测图像和恢复图像一；
[0024] 其中，图2(a)和图2(c)是实施例中受高斯噪声污染的退化图像；图2(b)和图2(d) 是实施例中恢复的图像；
[0025] 图3为本发明实施例中的退化观测图像和恢复图像二；图3(a)和图3(c)是实施例 中受椒盐噪声污染的退化图像;图3(b)和图像3(d)是实施例中恢复的图像。
【具体实施方式】
[0026] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0027] -种自适应参数更新的贝叶斯图像反卷积恢复方法，参考如图1，所述方法包括如 下步骤：
[0028] 步骤1，获取观测图像y，所述观测图像包含高斯噪声和高斯模糊，或者所述观测图 像包含椒盐噪声和高斯模糊;所述观测图像y=Nimp化X) +n，其中n表示高斯噪声，Nimp表示椒 盐噪声，H表示模糊核函数矩阵，X是待恢复的清晰图像，y是观测图像。
[0029] 其中n大小:N X 1，H大小:N X N，X大小:N X 1，y是观测退化图y大小:N X 1，N是观测 图像的像素个数。
[0030] 当n声0，Nimp = l时，表示观测图像受到高斯噪声污染；
[0031] 当n = 0，Nimp声1时，表示观测图像受到椒盐噪声污染。
[0032] 步骤2,采用零均值的广义高斯分布对观图像y中的噪声进行建模，从而得到观测 图像y的条件概率分布函数P(y|e，x)，p(〇| ?)表示已知?时O的条件概率分布，0为多变量 广义高斯观测模型的超参数，朗g从均匀分布。
[0033] 采用零均值的广义高斯分布对观测图像y中的噪声进行建模，得到观测图像y的条 件概率分布函数p(y|e，x):
[0034]
[00对其中，Al是正的实数，0<P《2，N是观测图像y的像素个数，e邱（?）表示We为底的 指数函数，e是多变量广义高斯观测模型的超参数，I I ? I Ip表示P范数处理，If表示P范数的 P次方处理，P(0 I ?)表示已知?时O的条件概率分布，Ki为正的常数，H表示模糊核函数矩 阵，X是待恢复的清晰图像，y是观测图像。
[0036] 步骤3,对所述待恢复的清晰图像X进行先验建模，得到待恢复的清晰图像X的先验 概率分布函数p(x|a)，a为图像先验模型超参数，a服从均匀分布。
[0037] 对所述待恢复的清晰图像X进行先验建模，得到待恢复的清晰图像X的先验概率分 布函数P(x|a);
[003引
痒中，人2是正的实数，0<q《2，a是图像先 验模型超参数，exp( ?)表示We为底的指数函数，01、化、03、04、0日、1)6分别表示0°、45°、90°、 135°方向一阶差分算子[1-1]、[0 1;-1 0]、[1;-1]、[10;0-1]W及水平和竖直方向二阶差 分算子[1 -2 对应的差分矩阵，M ? Ilq表示q范数处理，K表示q范数的q次方 处理，K2为正的常数。
[0039] 需要补充的是，超参数巧€{斯判，使用扁平非正常超先验建模超参数分布使得 衣(巧)='const.;其中const.是constant的缩写，表不常数，扁平先验是指该先验在定义域上 是均匀分布。
[0040] 步骤4,在分层贝叶斯模型下，构建联合全局概率分布p(a，e，x，y):
[0041 ] p(a，0，x，y) =p(a)p(0)p(y |0，x)p(x |a)
[0042] 图像反卷积问题的贝叶斯推断是基于未知后验分布p(a，e，x|y)的估计p(a，e，x y)=p(a，e，x，y)/p(y)，由于观测模型和图像先验的非凸非平滑特点，直接最大化后验概率 p(a，0，x I y)比较困难，因此，本发明实施例使用优化最小化(majorization-minimization， MM)方法来解决上述问题，首先需要求非凸先验的下界。
[0043] 步骤5,计算观测图像y的条件概率分布函数p(y |0，x)的下界Mi(0，x，y) :p(y |0，x) >]/[1(0，义，7);计算待恢复的清晰图像义的先验概率分布函数9^|0)的下界12(0，义）：9^|口） >M2(a，x)。
[0044] 已知不等；!^
，其等价表达式为i
，z >0,0<p《2;
[0045] i
[0046] 其中，化x-y)i表示向量化-y的第i个元素，Zi表示N维变分列向量Z的第i个元素，且 设有函数：
[0047]
[004引
[0049] N是观测图像y或待恢复的清晰图像X的像素个数，Ui表示N维变分列向量U的第i个 元素，M是正的实数，0<p《2，e邱（?）表示We为底的指数函数，0是多变量广义高斯观测 模型的超参数，H表示模糊核函数矩阵，A2是正的实数，0<q《2，a是图像先验模型超参数，Dj 表示差分矩阵；
[00加]观細Il巧傻V的签化舰莖麻巧狱n f V I B. 一的下巧责Mi化.X，y):
[0化1 ]
[0化2] 3M2(a，x):
[0化3]
[0054]步骤6,根据观测图像y的条件概率分布函数p(y|e，x)的下界、待恢复的清晰图像X 的先验概率分布函数p(x|a)的下界，得到所述联合全局概率分布p(a，e，x，y)的下界P(Q)P (0)Mi(0，x，y)M2(a，x):
[0化5] p(a，0，x，y) =p(a)p(0)p(y |0，x)p(x |a)>p(a)p(0)Mi化，x，y)M2(a，x)
[0056] 步骤7，求解后验分布P(a，e，X I y) = P(a，e，X，y)/p(y)的下界，得到关于图像先验 模型超参数〇、多变量广义高斯观测模型的超参数e、待恢复的清晰图像X的泛函。
[0057] 最大化后验分布p(a，0，x |y) =p(a，P，x，y)/p(y)的下界等价于：
[0化引
[0059] 其中0 ={a，0，x} ,log表示对数运算符，? = ui'gmin表示求泛函最小时对应的变量 贷 0的值?，N是观测图像y或待恢复的清晰图像X的像素个数，Ui表示N维变分列向量U的第i个 元素，M是正的实数，0<p《2，e邱（?）表示We为底的指数函数，0是多变量广义高斯观测 模型的超参数，H表示模糊核函数矩阵，A2是正的实数，0<q《2，a是图像先验模型超参数，Dj 表不差分矩阵。
[0060] 由上述最小化泛函可看出，利用优化最小化准则将原始观测模型和先验模型中含 有的非凸先验用一系列凸泛函代替，并引入变分向量Z和U。于是通过交替最小化方法可迭 代地求解上述泛函中的未知量a、0、X、Z和U。
[0061] 步骤8,迭代求解关于图像先验模型超参数a、多变量广义高斯观测模型的超参数 e、待恢复的清晰图像X的泛函，直到满足迭代停止条件，并将最后一次迭代得到的待恢复的 清晰图像作为反卷积恢复的清晰图像。
[0062] 步骤8具体包括如下子步骤：
[0063] (8a)初始化参数p，q，Ai，A2，0〇，a〇，
'i = l，...，N， 迭代次数k初值为0;
[0064] (8b)对于第k次迭代的结果zk，uk，利用下式求解关于待恢复的清晰图像X的最小化 问题：
[00 化]
[0066] 可得关于X的更新式为：
[0067] V, 间
J
[006引其中k表示迭代次数：
，diag (?)表示对角矩阵处理;在迭代的过程中，为了防止Z和U取值过小而使得Wf和Wr取到无穷大 的不合理估，化化下处理：
[0069]
[0070] 其中
[0071]
[0072] 运里 Ei = O.001。
[0073] (8c)利用下式求解关于变分向量Z的最小化问题：
[0074]
[00巧]可得关于Z的更新式夫
[0076] (8d)利用下式求解关于变分向量U的最小化问题：
[0077]
[0078] 可得关于U的更新式为
[0079] (Se)利用下式更新多变量广义高斯观测模型的超参数0:
[0080]
[0081]
[0082]
[0083] (8f)利用下式更新图像先验模型超参数a:
[0084]
[0085]
[0086]
[0087] (Sg)令迭代次数k加1，并依次重复子步骤(8b)至(8f)，直到||xk-xk"| |f/| |刮If 《e ;或迭代次数k达到预设的最大迭代次数；
[0088] 示例性的，迭代终止条件E = 1〇-6。
[0089] (她)获取最后一次迭代得到的待恢复的清晰图像，并将其作为反卷积恢复的清晰 图像。
[0090] 在本实施例中，参数Al和A2需要去逼近与先验分布p(y|e，x)和p(x|a)相应的划分 函数。但是，与先验分布P(y|e，x)和p(x|a)相应的划分函数在解析上难W处理，因此，有必 要近似处理先验分布P (y I e，X)和P (X I a)的划分函数。近似地确定先验分布P (y 10，X)和P (X a)的划分函数后，就可W近似地确定、和人2的值。且、和的值一旦选定，就保持不变。
[0091 ]示例性的，
[0092] 当n声0 ,Nimp = 1时，表示图像受高斯噪声污染。此时，可W设置P = 2, q = 1，入1 = 1.2,入 2 = 0.4。
[0093] 当n = 0 ,Nimp声1时，表示图像受椒盐脉冲噪声污染。此时，可W设置P = 1，Q= 1，入1 =1.2，入2 = 0.4。
[0094] 在本实施例中，部分参数的初始化值如下:迭代开始时，首先迭代初始化图像/为 输入的退化图像y。初始佔
，i = l,...，N;护=1〇-3，日0=1〇 ^;最大迭代次数maxiter = 500;上述参数初始化值对所有测试图像均适用。
[00%]为了验证本发明提供的贝叶斯图像反卷积方法的有效性，图2中的图2(a)和图2 (C)表示受高斯模糊和高斯噪声污染的退化图像，提出方法应用到图2(a)和图2(c)获得相 应的恢复图像为图2(b)和图2(d)，其中dB是decibel的缩写，表示分贝，BSNR是blur signa]_-t〇-noise ratio的缩写，表示模糊信噪比，PSNR是peak si即al-t〇-noise ratio的 缩写，表示峰值信噪比，BSNR和PSNR的单位均是地。图帥的图3(a)和图3(c)表示受高斯模 糊和椒盐噪声污染的退化图像，提出方法应用到图3(a)和图3(c)获得相应的恢复图像为图 3(b)和图 3(d)。
[0096] 从图2和图3的运2组退化图像可看出，图像被模糊和噪声退化后，图像本身细苄基 本被淹没，从视觉上难W分辨，运为图像后续处理带来很大的困难。从图2(b)和图2(d) W及 图3 (b)和图3 (d)可看出，本发明提供的贝叶斯反卷积方法有效地去除图像模糊，较好地保 留了图像本身的细节结构，也最大程度地抑制图像强噪声。运从主观视觉验证了本发明提 供的方法的有效性。
[0097] 表1列出了本发明给出的方法应用到四种不同噪声级别污染的退化图像恢复的 PSNR结果对比。从运个表可看出，提出的方法恢复结果的PSNR相比退化图像的PSNR获得了 明显的提高，此外，本发明提出的方法的PSNR也高于2008年Babacan等人提出的方法，且 Babacan等人提出的方法很难去除椒盐噪声。运再一次从客观评价指标上验证了本发明方 法的有效性。
[009引
[0099] 综上所述，本发明公开了一种分层贝叶斯图像反卷积恢复方法，首先对观测模型 和图像进行先验建模，然后再对超参数进行建模;其次为了克服观测模型和图像先验模型 中非凸泛函数值上难W求解的难题，本发明利用优化最小化方法将上述问题转化为易于求 解的非凸先验的下界，然后再最大化后验分布结合交替最小化的思想W迭代的方式得到想 要的解。上述实验结果表明本发明提出的自适应参数更新的方法不需要任何人工手动调整 参数，而且可W有效地抑制高斯噪声和椒盐噪声同时去除模糊。
[0100] W上所述，仅为本发明的【具体实施方式】，但本发明的保护范围并不局限于此，任何 熟悉本技术领域的技术人员在本发明掲露的技术范围内，可轻易想到变化或替换，都应涵 盖在本发明的保护范围之内。因此，本发明的保护范围应W所述权利要求的保护范围为准。
【主权项】
1. 一种自适应参数更新的贝叶斯图像反卷积恢复方法，其特征在于，所述方法包括如 下步骤： 步骤1，获取观测图像y，所述观测图像包含高斯噪声和高斯模糊，或者所述观测图像包 含椒盐噪声和高斯模糊;所述观测图像y=Nimp(Hx)+n，其中η表示高斯噪声，Nimp表示椒盐噪 声，H表示模糊核函数矩阵，X是待恢复的清晰图像，y是观测图像； 步骤2,采用零均值的广义高斯分布对观测图像y中的噪声进行建模，从而得到观测图 像y的条件概率分布函数P(y I β，χ)，β为多变量广义高斯观测模型的超参数，β服从均匀分 布； 步骤3,对所述待恢复的清晰图像X进行先验建模，得到待恢复的清晰图像X的先验概率 分布函数Ρ(Χ|α)，α为图像先验模型超参数，α服从均匀分布； 步骤4,在分层贝叶斯模型下，构建联合全局概率分布p(a，i3，X，y): ρ(α,β,χ,γ)=ρ(α)ρ(β)ρ(γ |β,χ)ρ(χ|α) 步骤5，计算观测图像y的条件概率分布函数ρ (y I β, χ)的下界Mi (0，x，y):p(y|0，x)彡Mi (P，x，y);计算待恢复的清晰图像x的先验概率分布函数ρ(χ |α)的下界Μ2(α，χ) :ρ(χ |α)彡M2 (α,χ)； 步骤6,根据观测图像y的条件概率分布函数p(y|i3，x)的下界、待恢复的清晰图像x的先 验概率分布函数p(x I α)的下界，得到所述联合全局概率分布ρ(α，β，χ，γ)的下界P(Ct)p(P)M 1 (β,χ,γ)Μ2(α,χ)： ρ(α,β,χ,γ)=ρ(α)ρ(β)ρ(γ |β,χ)ρ(χ|α)^：ρ(α)ρ(β)Μι(β,χ,γ)Μ2(α,χ) 步骤7,求解后验分布？(〇，0，^7)=?(〇，00，7)/?( 7)的下界，得到关于图像先验模型 超参数α、多变量广义高斯观测模型的超参数β、待恢复的清晰图像X的泛函； 步骤8,迭代求解关于图像先验模型超参数α、多变量广义高斯观测模型的超参数β、待 恢复的清晰图像X的泛函，直到满足迭代停止条件，并将最后一次迭代得到的待恢复的清晰 图像作为反卷积恢复的清晰图像。2. 根据权利要求1所述的一种自适应参数更新的贝叶斯图像反卷积恢复方法，其特征 在于，步骤1中所述观测图像y=Nimp (Hx) +η: 当η辛0，NimP=l时，表示观测图像包含高斯噪声； 当n = O，Nimp辛1时，表示观测图像包含椒盐噪声。3. 根据权利要求1所述的一种自适应参数更新的贝叶斯图像反卷积恢复方法，其特征 在于，步骤2具体为： 采用零均值的广义高斯分布对观测图像y中的噪声进行建模，得到观测图像y的条件概 率分布函数P(y |β，χ):其中，A1是正的实数，0〈ρ$2，Ν是观测图像y的像素个数，exp( ·)表示以e为底的指数函 数，β是多变量广义高斯观测模型的超参数，I I · I Ip表示P范数处理，It表示P范数的P次方 处理，K1为正的常数，H表示模糊核函数矩阵，X是待恢复的清晰图像。4. 根据权利要求1所述的一种自适应参数更新的贝叶斯图像反卷积恢复方法，其特征 在于，步骤3具体为： 对所述待恢复的清晰图像X进行先验建模，得到待恢复的清晰图像X的先验概率分布函 数P(x|a):其中，λ2是正的实数，0〈q<2，a是图像先验模型超参数，N是观测图像y的像素个数， 示差分矩阵，I I · I U表示q范数处理，If表示q范数的q次方处理，K2为正的常数。 5 .根据权利要求1所述的一种自适应参数更新的贝叶斯图像反卷积恢复方法，其特征 在于，步骤5具体为：Ρ^Ξ2；其中，（版^^表示向量Hx-y的第i个元素，Zl表示N维变分列向量ζ的第i个元素，且设有 函数：N是观测图像y或待恢复的清晰图像X的像素个数，U1表示N维变分列向量u的第i个元素， 入:是正的实数，0〈p<2，eXp( ·)表示以e为底的指数函数，β是多变量广义高斯观测模型的 超参数，H表示模糊核函数矩阵，λ2是正的实数，〇〈q< 2, α是图像先验模型超参数，Dj表示差 分矩阵； 观测图像y的条件概率分布函数P (y I β，X)的下界为Mi (β，X，y):表示P范数处理，|'|$表示 P范数的P次方处理； 待恢复的清晰图像X的条件概率分布函数P (X I α)的下界为M2 (α，X):I I · I Iq表示q范数处理，表示q范数的q次方处理。 6 .根据权利要求1所述的一种自适应参数更新的贝叶斯图像反卷积恢复方法，其特征 在于，步骤7具体为： 最大化后验分布口(〇，04|7)=口(€[，04,7)/^(7)的下界等价于：其中0 ={α，β，χ}，l〇g表示对数运算符，0 = arg._min表示求泛函最小时对应的变量@的 Θ 值忌，N是观测图像y或待恢复的清晰图像X的像素个数，U1表示N维变分列向量u的第i个元 素，A1是正的实数，0〈p<2，eXp( ·)表示以e为底的指数函数，β是多变量广义高斯观测模型 的超参数，H表示模糊核函数矩阵，λ2是正的实数，〇〈q<2，a是图像先验模型超参数，Dj表示 差分矩阵。7.根据权利要求6所述的一种自适应参数更新的贝叶斯图像反卷积恢复方法，其特征 在于，步骤8具体包括如下子步骤： (8a)初始化参彳迭代次 数k初值为0; (8b)对于第k次迭代的结果zk，uk，利用下式求解关于待恢复的清晰图像X的最小化问 题：对角矩阵处理； (8c)利用下式求解关于变分向量z的最小化问题：(8d)利用下式求解关于变分向量u的最小化问题： 厂认双K刀μι，开TO认里反丁灰3來、〇1^王、〇1;，且到I |x:-xk+1| |f/| |xk| |f彡ε;或 迭代次数k达到预设的最大迭代次数； (8h)获取最后一次迭代得到的待恢复的清晰图像，并将其作为反卷积恢复的清晰图 像。
【文档编号】G06T5/00GK106023106SQ201610330261
【公开日】2016年10月12日
【申请日】2016年5月18日
【发明人】方厚章, 刘宏伟, 潘东辉, 时愈, 许述文, 刘军
【申请人】西安电子科技大学