×

图像重建的全变分优化共轭梯度法。 (英语) Zbl 1430.94040号

反向问题。 34,第3号,文章ID 034001,26 p.(2018); 勘误表同上,36,第8号,文章编号089601,第2页(2020)。
摘要:共轭梯度(CG)方法通常用于最小二乘问题的相对快速求解。在图像重建中,该问题可能不存在,也可能受到噪声污染;因此,应该使用正则化等方法。总变差(TV)是一种有用的正则化惩罚,常用于图像重建,以生成具有尖锐边缘的图像。当选择非二次范数进行正则化时,如TV的情况,则不再可能使用CG。非线性CG是一种替代方法,但它不具备CG用最小二乘法显示的效率,对于TV范数问题,首选快速迭代收缩阈值算法(FISTA)等方法。包含先验信息的另一种方法是优势化。
本文证明了共轭梯度法的优越性。提出了五种不同的CG变体,包括预处理CG。提出了以总变分范数为优势的CG方法,并证明了它们在图像重建中的性能。结果表明,由于原始CG处理最小二乘问题的速度快,一些提出的高级CG方法变体可以生成比FISTA生成的重建质量更好、计算时间更少的重建。
在附录中,我们检查了一种高级CG方法(我们称之为S-CG)的行为;它的输入参数之一是正数\(\varepsilon\)。证明了,对于任何大于最小二乘解的半平方残差的给定(varepsilon),S-CG以有限步数终止,其输出的半平方余量小于或等于。重要的是,对于相同的半平方残差值(varepsilon),输出的TV值也会低于未经验证的CG提供的TV值。

MSC公司:

94A08型 信息和通信理论中的图像处理(压缩、重建等)
65J20型 抽象空间中不适定问题的数值解;正则化
68单位10 图像处理的计算方法
90C20个 二次规划
90摄氏52度 减少梯度类型的方法
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Björck,A.,最小二乘问题的数值方法,(1996),费城,宾夕法尼亚州:SIAM,费城·Zbl 0847.65023号
[2] Barrett,H.H。;Myers,K.J.,《图像科学基础》(2004),新泽西州霍博肯:新泽西州威利
[3] Kim,K。;Shevlyakov,G.,《为什么是Gaussianity?》?,IEEE信号处理。Mag.,25,102-113,(2008)·doi:10.1109/MSP.2007.913700
[4] Vogel,C.R.,《反问题的计算方法》,(2002),费城,宾夕法尼亚州:SIAM,费城·Zbl 1008.65103号
[5] Bovik,A.C.,《图像和视频处理手册》,(2000),加利福尼亚州圣地亚哥:学术出版社,加利福尼亚州圣迭戈·Zbl 0967.68155号
[6] Romberg,J.K.,《压缩采样成像》,IEEE信号处理。Mag.,25,14-20,(2008年)·doi:10.1109/MSP.2007.914729
[7] 戴,Y。;Han,J。;刘,G。;Sun,D。;尹,H。;Yuan,Y.,非线性共轭梯度法的收敛性,SIAM J.Optim。,10, 345-358, (2000) ·Zbl 0957.65062号 ·doi:10.1137/S10526223494268443
[8] Hager,W。;Zhang,H.,非线性共轭梯度法综述,Pac。J.Optim。,2, 35-58, (2006) ·兹比尔1117.90048
[9] Rockafellar,R.T。;Wets,R.J B.,变分分析,(1998),柏林:施普林格出版社,柏林·Zbl 0888.49001号
[10] 赵,F。;诺尔,哥伦比亚特区。;尼尔森,J-F;Fessler,J.A.,《通过压缩传感分离幅度和相位正则化》,IEEE Trans。医学成像,311713-1723,(2012)·doi:10.1109/TMI.2012.2196707
[11] Zibetti,M.V W。;Pipa,D.R。;De Pierro,A.R.,作为迭代重建算法加速器的快速准确一维L2-L1优化,Digit。信号处理。,48, 178-187, (2016) ·doi:10.1016/j.dsp.2015.09.009
[12] Combettes公司。;Pesquet,J.C。;Burachik,H.H.,《信号处理中的近距离分裂方法》,《科学与工程中反问题的定点算法》,185-212,(2011),纽约:Springer,纽约·Zbl 1242.90160号
[13] 齐布列夫斯基,M。;Elad,M.,《信号和图像处理中的L1-L2优化》,IEEE信号处理。Mag.,27,76-88,(2010)·doi:10.1109/MSP.2010.936023
[14] Bioucas-Dias,J。;Figueiredo,M.,《新的TwIST:图像恢复的两步迭代收缩/阈值算法》,IEEE Trans。图像处理。,16, 2992-3004, (2007) ·doi:10.1109/TIP.2007.909319
[15] 贝克,A。;Teboulle,M.,线性反问题的快速迭代收缩阈值算法,SIAM J.成像科学。,2183-202年,(2009年)·Zbl 1175.94009号 ·doi:10.1137/080716542
[16] 北卡罗来纳州帕里赫。;Boyd,S.,《近似算法》,Found。最佳趋势。,1, 127-239, (2014) ·数字对象标识代码:10.1561/24000003
[17] Chambolle,A.,《总变异最小化算法及其应用》,J.Math。成像视觉。,20, 89-97, (2004) ·Zbl 1366.94048号 ·doi:10.1023/B:JMIV.00000011320.81911.38
[18] Chambolle,A。;Pock,T.,凸问题的一阶原对偶算法及其在成像中的应用,J.Math。成像视觉。,40, 120-145, (2011) ·Zbl 1255.68217号 ·doi:10.1007/s10851-010-0251-1
[19] 贝克,A。;Teboulle,M.,约束全变差图像去噪和去模糊问题的快速梯度算法,IEEE Trans。图像处理。,182419-244,(2009年)·Zbl 1371.94049号 ·doi:10.1109/TIP.2009.2028250
[20] Censor,Y。;Davidi,R。;Herman,G.T.,扰动弹性和迭代算法的优越性,反问题,26,(2010)·Zbl 1193.65019号 ·doi:10.1088/0266-5611/26/6/065008
[21] Herman,G.T。;加杜尼奥,E。;Davidi,R。;Censor,Y.,《高级化:医学物理的优化启发式》,医学物理。,39, 5532-5546, (2012) ·数字对象标识代码:10.1118/1.4745566
[22] Hansen,P.C.,秩亏和离散不适定问题:线性反演的数值方面,(1998),费城,宾夕法尼亚州:SIAM,费城
[23] Censor,Y。;Davidi,R。;Herman,G.T。;舒尔特,R.W。;Tetraushvili,L.,投影次梯度最小化与优化,J.Optim。理论应用。,160, 730-747, (2014) ·Zbl 1298.90104号 ·doi:10.1007/s10957-013-0408-3
[24] 加杜尼奥,E。;Herman,G.T.,ML-EM算法的优化,IEEE Trans。编号。科学。,61, 162-172, (2014) ·doi:10.1109/TNS.2013.2283529
[25] 施拉普,M.J。;Herman,G.T.,《使用优势方法的x射线计算机断层扫描数据融合》,《科学评论》。仪器。,85, (2014) ·数字对象标识代码:10.1063/1.4872378
[26] Helou,E.S。;Zibetti,M.V W。;Miqueles,E.X.,统计层析图像重建增量优化算法的优化,反问题,33,(2017)·Zbl 1367.65029号 ·doi:10.1088/1361-6420/3/34/044010
[27] Björck,A。;Elfving,T。;Strakos,Z.,线性最小二乘问题的共轭梯度稳定性和Lanczos方法,SIAM J.矩阵分析。申请。,19, 720-736, (1998) ·Zbl 0924.65035号 ·doi:10.1137/S089547989631202X
[28] Herman,G.T.,《计算机层析成像基础:投影图像重建》(2009),伦敦:施普林格出版社,伦敦·Zbl 1280.92002年
[29] 戴,Y-H;Kou,C-X,具有最优性质的非线性共轭梯度算法和改进的Wolfe线搜索,SIAM J.Optim。,23, 296-320, (2013) ·Zbl 1266.49065号 ·数字对象标识代码:10.1137/100813026
[30] Luenberger,D.G。;Ye,Y.,线性和非线性规划,(2008),纽约:Springer,纽约·Zbl 1207.90003号
[31] 弗莱塞,V。;Giraud,L.,实数和复数算术的共轭梯度例程集,技术报告TR/PA/00/47,(2000)
[32] Bolz,J。;农民,I。;格林斯彭,E。;Schröder,P.,《GPU上的稀疏矩阵解算器:共轭梯度和多重网格》,ACM Trans。图表。,22, 917-924, (2003) ·doi:10.1145/882262.882364
[33] Shewchuk,J.R.,无痛苦共轭梯度法简介,技术报告,(1994)
[34] 汉弗莱斯,T。;Winn,J。;Faridani,A.,从稀疏视图和有限角度多能量数据重建CT图像的优化算法,Phys。医学生物学。,62, 6762-6783, (2017) ·doi:10.1088/1361-6560/aa7c2d
[35] 董,Q。;Gibali,A。;江,D。;唐寅,超梯度型方法的有界扰动弹性及其应用,J.不等式应用。,280, (2017) ·Zbl 1378.49026号 ·doi:10.1186/s13660-017-1555-0
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。