杰西·亚当斯;马蒂亚斯·莫斯菲尔德;凯文·乔伊斯;玛丽丽莎·霍华德;亚伦·路特曼 一种用于大规模图像去模糊中的降维Gibbs采样的分块方案。 (英语) Zbl 07480113号 反向探测。科学。工程师。 29,第12期,1789-1810(2021). 摘要:使用马尔可夫链蒙特卡罗(MCMC)方法从贝叶斯反问题的后验分布中进行抽样时,最重要的挑战之一是抽样变得难以计算的速度,这是估计参数数量的函数。在图像去模糊中,文献中有许多MCMC算法,但很少尝试重建大于(512次512)像素(10^5参数)的图像。在用于诊断动态材料实验的定量X射线照相术中,图像可能会大得多,导致数百万参数的问题。我们解决了这个问题,并通过一种分块方案构造了一个Gibbs采样器,该分块方案产生一个稀疏且高度结构化的后验精度矩阵。Gibbs取样器在采样过程中自然利用了特殊的矩阵结构,使其成为“维数-分辨率”,因此其混合特性几乎与图像大小无关,生成一个样本在计算上是可行的。维度丰满性能够在适度的计算平台上表征大规模图像去模糊问题的后验概率。我们通过对美国能源部内华达州国家安全局Cygnus双束X射线照相设备拍摄的尺寸为(4096times4096)像素((10^7)参数)的射线照片进行去模糊处理,证明了该方法的适用性。 MSC公司: 15年29日 线性代数中的反问题 65层22 数值线性代数中的病态性和正则化问题 65立方厘米05 蒙特卡罗方法 65C60个 统计中的计算问题(MSC2010) 94A08级 信息与通信理论中的图像处理(压缩、重建等) 关键词:反褶积;马尔科夫蒙特卡洛;吉布斯采样 软件:贝叶斯DA;GMRF库 PDF格式BibTeX公司 XML格式引用 \textit{J.Adams}等人,《反问题》。科学。工程29,编号12,1789--1810(2021;Zbl 07480113) 全文: 内政部 链接 参考文献: [1] 霍尔,BR;罗伊,S。;Gord,JR,使用宽带闪光X射线照相术对喷雾中单次喷射液体分布进行定量成像,国际J多相流,87,241-249(2016) [2] Hanson,KM;GS.坎宁安,《贝叶斯推理引擎》。In:最大熵和贝叶斯方法,125-134(1996),荷兰多德雷赫特:荷兰多德莱赫特Kluwer学院·Zbl 0900.62147号 [3] 霍华德,M。;Fowler先生。;Luttman,A.,定量X射线照相术中的Bayesian Abel反演,SIAM科学计算杂志,38,B396-B413(2016)·Zbl 1339.92036号 [4] 梅尔,E。;PJ.威瑟斯。,定量X射线断层成像,国际材料评论,59,1-43(2014) [5] 福勒,M。;霍华德,M。;Luttman,A.,《高能X射线成像系统用不确定性估计量化模糊的随机方法》,《逆向概率科学工程》,24,353-371(2016) [6] SVS纳格什;拉纳,R。;Russ,M。;爱奥尼亚,中国;Bednarek博士;Rudin,S.,高分辨率直接转换x射线探测器的焦点去模糊,Proc SPIE Int Soc Opt Eng,9783,97833R(2016) [7] von Wittenau,AES;洛根,CM;Aufderheide,MB,《使用平板成像系统的兆伏射线照相术中的模糊伪影:蒙特卡罗模拟与测量的比较》,《医学物理学》,29,2559-2570(2002) [8] 汉森,PC;Nagy,JG;奥利里,DP。,去模糊图像:矩阵、光谱和滤波(2006),费城,宾夕法尼亚州:SIAM,费城·Zbl 1112.68127号 [9] Bardsley,JM.,反问题的计算不确定性量化(2018),费城(PA):费城SIAM(PA)·Zbl 1435.60001号 [10] 巴兹利,JM;Luttman,A.,《处理基于MCMC的反褶积中的边界伪影》,《线性代数应用》,47339-358(2015)·Zbl 1312.94008号 [11] 巴兹利,JM;Luttman,A.,《泊松似然和高斯先验线性反问题的Metropolis-Hastings方法》,《国际不确定量子化杂志》,6,35-55(2016)·Zbl 1498.94009号 [12] 福克斯,C。;拉脱维亚诺顿。,线性高斯反问题中的快速采样,SIAM/ASA J Uncertain Quantif,41191-1218(2016)·Zbl 1398.94081号 [13] 福克斯,C。;Parker,A.,使用矩阵分裂和多项式对正态分布进行加速吉布斯采样,Bernoulli,233711-3743(2017)·Zbl 1457.62085号 [14] 霍华德,M。;福勒,M。;Luttman,A.,X射线照片反褶积中基于采样的不确定性量化,计算机应用数学杂志,270,43-51(2014)·Zbl 1321.65031号 [15] 乔伊斯,KT;吉米·巴德斯利;Luttman,A.,部分坍塌吉布斯采样X射线成像中的点扩散函数估计,SIAM J Sci-Comput,40,B766-B787(2018)·Zbl 1387.65005号 [16] 王,Z。;巴兹利,J。;Solonen,A.,《具有(####)先验的贝叶斯反问题:随机优化方法》,SIAM科学计算杂志,39,S140-S166(2017)·Zbl 1422.65028号 [17] A.帕克。;皮特斯,B。;Lorenz,L.,《生物膜成像用大高斯模型的多项式加速解:理论和有限精度》,美国统计学会杂志,1131431-1442(2018)·Zbl 1409.62226号 [18] 陈,J。;Anitescu,M。;Saad,Y.,《通过最小二乘多项式近似进行计算》,SIAM科学计算杂志,33,195-222(2011)·Zbl 1234.65027号 [19] 周,E。;Saad,Y.,采样多元高斯分布的预处理krylov子空间方法,SIAM科学计算杂志,36,A588-A608(2014)·Zbl 1296.60087号 [20] A.帕克。;Fox,C.,带共轭梯度的krylov空间中的抽样高斯分布,SIAM科学计算杂志,34,B312-B334(2012)·Zbl 1246.65057号 [21] 吉尔克斯,W。;理查森,S。;Spiegelhalter,D.,《马尔可夫链蒙特卡罗实践》(1996),佛罗里达州博卡拉顿:斯普林格,博卡拉顿,佛罗里达州·Zbl 0832.00018号 [22] 鲁宾斯坦,RY;Kroese博士。,模拟和蒙特卡罗方法(2017),新泽西州霍博肯:新泽西州威利 [23] Morzfeld,M。;汤,X。;马祖,YM。,MCMC的定位:用局部结构抽样高维后验分布,《计算物理杂志》,380,1-28(2018)·Zbl 1451.65011号 [24] 罗伯茨,GO;Sahu,S.,吉布斯采样器的更新方案、相关结构、阻塞和参数化,J R Stat Soc Ser B,59,291-317(1997)·Zbl 0886.62083号 [25] Chen,V,Dunlop,MM,Papaspiliopoulos,O,et al.贝叶斯反问题中的尺寸鲁棒MCMC。提交;2018 [26] Beskos,A。;罗伯茨,G。;Stuart,A.,《高维非产品目标上本地Metropolis-Hastings链的最优缩放》,Ann Appl Probab,19863-898(2009)·Zbl 1172.60328号 [27] 罗伯茨,G。;Gelman,A。;Gilks,W.,《随机行走Metropolis算法的弱收敛性和最优缩放》,Ann Appl Probab,7110-120(1997)·Zbl 0876.60015号 [28] 罗伯茨,G。;Rosenthal,J.,《Langevin扩散离散近似的最佳缩放》,J R Stat Soc,Ser B(Stat Methodol),60,255-268(1998)·Zbl 0913.60060号 [29] 布奇尼,A。;多纳特利,M。;Reichel,L.,带一般惩罚项的迭代Tikhonov正则化,数值线性代数应用,24,1-19(2017)·Zbl 1424.65045号 [30] 陈,H。;王,C。;Song,Y.,用于图像去模糊的分裂Bregmaized各向异性全变分模型,J Vis Commun image Representation,31282-293(2015) [31] 焦,Y。;金,Q。;Lu,X.,线性反问题乘数的交替方向法,SIAM J Numer Ana,54,2114-2137(2016)·Zbl 1342.65137号 [32] 马,L。;徐,L。;Zeng,T.,图像去模糊的低秩先验和全变差正则化,科学计算杂志,701336-1357(2017)·Zbl 1366.65041号 [33] 陶,S。;Dong,W。;Xu,Z.,泊松噪声污染模糊图像的快速全变差反褶积,J Vis Commun image Representation,38,582-594(2016) [34] 徐,J。;Chang,HB;Qin,J.,图像去模糊的区域分解方法,J Comput Appl Math,271,401-414(2014)·Zbl 1334.68287号 [35] 贝特罗,M。;Boccacci,P.,《Richardson-Lucy图像反褶积方法中减少边界效应的简单方法》,Astron Astrophys,437369-374(2005) [36] 维奥·R。;巴兹利,JM;Donatelli,M.,处理最小二乘图像反褶积问题中的边缘效应,Astron Astrophys,442397-403(2005) [37] Bardsley,JM.,逆问题的高斯马尔可夫随机场先验,逆概率成像,7397-416(2013)·Zbl 1268.65011号 [38] H街。;Held,L.,《高斯马尔可夫随机场:理论与应用》(2005),纽约:查普曼和霍尔/CRC,纽约·邮编1093.60003 [39] Gelman,A。;卡林,JB;Stern,HS,贝叶斯数据分析(2013),佛罗里达州博卡拉顿:查普曼和霍尔CRC出版社,佛罗里达州波卡拉顿 [40] Geman,S。;Geman,D.,《随机松弛、吉布斯分布和图像的贝叶斯恢复》,IEEE Trans-Pattern Ana Mach Intell,PAMI-6721-741(1984)·Zbl 0573.62030号 [41] Adams,J.X射线照相术中图像去模糊的可缩放块-Gibbs采样[博士论文]。图森(亚利桑那州):亚利桑纳大学;2019 [42] Asaki,TJ;Chartrand,R。;Vixie,KR,使用全变分正则化的Abel反演,反演问题,211895-1903(2005)·Zbl 1094.65129号 [43] Asaki,TJ;公关部坎贝尔;Chartrand,R.,《使用总变差正则化的阿贝尔反演:应用》,《反演问题科学与工程》,第14873-885页(2006年)·Zbl 1200.65115号 [44] 戴维斯,G。;Jain,N。;Elliott,J.,《光束硬化修正的建模方法》,Proc SPIE Int Soc Opt Eng,7078,70781E(2008) [45] 关天俊;伯宁格,M。;Snell,C.,使用射线照相链模型模拟天鹅座棒箍缩二极管,IEEE Trans Plasma Sci,37530-537(2009) [46] 西曼,HE;Roth,B.,《射线照相用新型阶梯楔》,《放射学报》,第53期,第215-226页(1960年) [47] 英国,H。;汉克,M。;Neubauer,A.,反问题的正则化(1996),多德雷赫特,荷兰:Kluwer Academic,荷兰多德雷希特·Zbl 0859.65054号 [48] 戈卢布,GH;希思,M。;Wahba,G.,《广义交叉验证作为选择良好岭参数的方法》,《技术计量学》,21,215-223(1979)·Zbl 0461.62059号 [49] 雷诺特,RA;舵手,AW;Vatankhah,S.,tikhonov正则化参数的无偏预测风险估计:奇异值分解的递增秩近似收敛,BIT数值数学,59,1031-1061(2019)·Zbl 1434.65039号 [50] Vogel,CR,反问题的计算方法(2002),费城(PA):SIAM,费城·Zbl 1008.65103号 [51] Sokal,A。;德维特·莫雷特,C。;卡地亚,P。;Folacci,A.,《功能整合》。北约ASI系列(B系列:物理)。,361,《统计力学中的蒙特卡罗方法:基础和新算法》(1998),波士顿(马萨诸塞州):斯普林格,波士顿(麻省) [52] Wolff,U.,Monte Carlo errors with less errors,计算物理通讯,156,145-153(2004)·Zbl 1196.65138号 此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。