基于卷积盲降噪的混合式核磁共振成像.pdf

基于卷积盲降噪的混合式核磁共振成像宗春梅1,张月琴2,郝耀军11(忻州师范学院计算机系,忻州034000)2(太原理工大学计算机科学与技术学院,太原030024)通信作者:宗春梅,E-mail:摘要:为了解决图像压缩感知重建研究领域中通过有效的图像先验信息重构与原图相似性高且保留细节消除伪影的高质量图像的问题,针对不足采样的 K 空间数据,在经典的 CNN 算法 CBDNet 算法的基础上,通过融合深度学习先验信息及传统图像恢复各自优势的方法,研究了基于深度神经网络去噪先验和 BM3D 块压缩感知算法的混合式重构算法.该算法采用交互式方法训练多尺度残差网络抑制噪声水平,借优化选择的方式将深度学习与传统块匹配多尺度结合以提取图像不同尺度的特征数据从而实现抑制伪影、快速重建高质量 MRI.实结果表明深度学习结合 BM3D 在 MR 图像重构领域能够有效降低伪影保留细节信息,加强重构效果.与此同时,通过采用 GPU 的加速运算,算法的计算复杂度较使用单一算法并未增加很多.可见基于卷积盲降噪的混合式核磁共振成像效果更佳.关键词:压缩感知;卷积盲降噪;图像重建;深度学习;非局部相似性引用格式:宗春梅,张月琴,郝耀军.基于卷积盲降噪的混合式核磁共振成像.计算机系统应用,2023,32(12):1220.http:/www.c-s- Magnetic Resonance Imaging Based on Convolutional Blind DenoisingZONGChun-Mei1,ZHANGYue-Qin2,HAOYao-Jun11(DepartmentofComputerScience,XinzhouNormalUniversity,Xinzhou034000,China)2(CollegeofComputerScienceandTechnology,TaiyuanUniversityofTechnology,Taiyuan030024,China)Abstract:Inthefieldofimagecompressionperceptualreconstruction,high-qualityimagesarereconstructedwithhighsimilaritytotheoriginalimage,anddetailsareretainedtoeliminateartifactsthrougheffectiveimagepriorinformationreconstruction.Thus,aimingattheK-spacedatawithinsufficientsampling,basedontheclassicCNNalgorithmCBDNetalgorithm,thisstudyadoptsthemethodtocombinetheadvantagesoffusingdeeplearningpriorinformationandtraditionalimagerestoration.Meanwhile,ahybridreconstructionalgorithmbasedonpriordenoisingofdeepneuralnetworkandcompressedsensingalgorithmofBM3Dblockisstudied.Thealgorithmemploysaninteractivemethodtotrainamulti-scaleresidualnetworktosuppressnoiselevelsandcombinesdeeplearningwiththemulti-scalematchingoftraditionalblockstoextractimagefeaturedataatdifferentscalesthroughoptimalselection,thussuppressingartifactsandquicklyreconstructinghigh-qualityMRI.TheexperimentalresultsshowthatdeeplearningcombinedwithBM3DcanreduceartifactsandretaindetailsinMRimagereconstruction,enhancingthereconstructioneffect.Additionally,thecomputationalcomplexityofthealgorithmisnotmuchmorethanthatofthesinglealgorithmbytheGPUacceleratedoperation.ItcanbeseenthatthehybridMRIbasedonconvolutionblinddenoisinghasabettereffect.Key words:compressedsensing;convolutionalblinddenoising;imagereconstruction;deeplearning;non-localsimilarity计算机系统应用ISSN1003-3254,CODENCSAOBNE-mail:ComputerSystems&Applications,2023,32(12):1220doi:10.15888/ki.csa.009316http:/www.c-s-中国科学院软件研究所版权所有.Tel:[phone]基金项目:国家自然科学基金(61876124);教育部人文社科青年基金(20YJC630034)收稿时间:2023-05-28;修改时间:2023-06-28;采用时间:2023-07-03;csa 在线出版时间:2023-09-21CNKI 网络首发时间:2023-09-2212专论综述SpecialIssue20 世纪 80 年代初引入医学影像诊断的核磁共振成像(MRI)技术是近年来一种新型的具有高度软组织分辨能力的高科技影像学检查方法,磁共振图像的获取是通过进动的质子切割磁感线,产生电信号,信号频率的空间信息的存储位置存放在 K 空间的值上,再通过反傅立叶变换,就得到了我们看到的磁共振影像,但是各因素间的相互制约导致核磁共振成像时间长、信噪比低、分辨率差等,寻找一种采集时间短重构效果好的成像方式是临床核磁共振成像所考虑的最重要因素之一1.压缩感知技术是在满足测量矩阵与稀疏基的不相关性的前提下,通过对 K 空间信号进行远低于奈奎斯特率值的高度欠采样,边采样边压缩,进而通过获得的测量矩阵与稀疏化后的信号所构成的非线性重构优化算法来恢复信号,显著减少传感元数量、减少了采集时间、降低了存储空间等,并获得很好的重建效果.因此针对压缩感知的采样和重构学者们提出了众多算法.完全填充的 K 空间获得最高质量图像的同时也会导致较长扫描时间,在压缩感知里,通过采集较少数据可以加快采集速度,但同时会导致图像信噪比降低,因此信号稀疏采集是关键,可以忽略其余信号只捕获最重要的 MR 信号.Yuan 等人2提出了使用非下采样碎片变换对核磁共振图像实现了一个灵活的分解,在每个级别上都有任意偶数的方向子带,提供了足够的稀疏表示和任何方向的分解,此算法的高方向灵敏度及最优近似特性提高了 CS-MRI 重建应用.Huang 等人3提出在多个小波帧上具有平均稀疏度的基础上提高平均信号稀疏性,利用稀疏度平均先验的快速迭代收缩阈值算法进行 CS-MRI 重建.Kojima 等人4提出随机获取数据需要一个最优的随机欠采样模式,在 K 空间中心区域估计了不同加速度因子和全采样范围的图像恢复精度,指出可以通过增加高频区域的相位编码数来优化高斯分布参数从而提高空间分辨率.Kazama等人5研究了非随机信号欠采样 CS-MRI以稳定重建图像质量和避免采样点选择的任意性的可行性,提出了一种利用在迭代重建过程中使用曲线变换作为稀疏变换函数的方法增加了傅里叶基和稀疏变换函数之间的非相干性.即便采样技术相对成熟,采样后数据集进行傅里变换时得到的图像仍然是有噪声的,需要通过滤波等技术去除图像中的噪声以提高图像质量.随着稀疏采集技术的逐步成熟,采集后的重建算法也在不断创新.Yuan 等人6研究了两种不同的分析L1/2 正则化的迭代半阈值算法求解非凸优化问题,利用图像在紧帧下的稀疏 K 空间数据重建高分辨率的图像,从而有效地缩短成像时间.但是为获得最优图像非线性最优化算法需要一个迭代最优化过程,这样就造成重建时间过长,并且迭代可能陷入局部收敛.Mathew等人7提出了采用“使用参数值处理稀疏正则化”与“估计后续迭代中使用的参数值”两个步骤交替最小化方法得到解.此类基于最优化模型的重建方法一般需要设置各种优化参数,而不恰当的参数设置会导致过约束,从而造成重建图像看起来不自然,如过平滑或残留欠采样伪影.He 等人8提出了一种能够在幅度和相位图像中捕获更多细节的基于双树复小波稀疏性的新方法,该方法消除了由于相位跳跃而造成的幅度伪影.但是目前广泛应用的基于整体变分(totalvariation)、小波等稀疏变换对于处理具有复杂结构的生物组织图像仍显得过于简单,会引入阶梯状或块状伪影.复杂的正则化与过程迭代会加大耗时和成本,参数无法保证调到最优,通过模糊的先验信息转化成的正则化不能确保符合通用需求.近年来,深度学习在图像修复领域的成功为 MRI 高速重建提升了可能性,广泛的关注促使产生了多种将深度学习应用于 CS-MRI的算法,并且已经证明了深度学习能够在扫描时间差距不大的情况下同时提升图像重建质量.Hashimoto 等人9提出了利用卷积神经网络(CNN)作为先验信息,实现了用原始 K 空间数据替换未填充区域的 K 空间修正,使用迭代的方式实现了 CS-MRI 重建,减少了混叠伪影.Ran 等人10提出了一种同时操作 K 空间和空间域数据的包含两个并行和交互式分支的 MRI 双域重建网络,进一步提高了重构性能.Zhao 等人11提出了一种以高分辨率参考 MR 图像作为卷积神经网络的输入,然后加入 K 空间数据校正步骤,从欠采样的K 空间数据中获得了更准确的重建结果.Njeh 等人12提出了基于自动编码器的体系结构的 CS-MRI 重建深度卷积编码器-解码器架构,提供了更好地重建并保留了纹理图像的细节,同时缩短了运行时间.Mardani 等人13提出了使用生成对抗网络(GAN)来建模高质量MR 图像,与传统的基于小波和基于字典学习的 CS 方案以及使用像素级训练的基于深度学习的方案相比,此算法以更精细的纹理细节检索到更高质量的图像.在本文中,为了提高重建效果,我们提出了一种新2023年第32卷第12期http:/www.c-s-计 算 机 系 统 应 用SpecialIssue专论综述13的深度学习结合 BM3D 优势的核磁共振成像算法.与目前的方法不同,盲降噪深度学习具有更好的去除自然噪声的能力,BM3D 在图像细节保留有突出的优势,所提出的算法将两者结合,可以将深度网络的学习能力强但可解释性差和 BM3D 的去噪效果具有可解释性相结合以实现互相弥补.1基于稀疏基和先验信息的压缩感知核磁共振成像图像在空间上的离散化称为采样,也就是用空间上部分点的灰度值代表图像,这些点称为采样点.二维信号一般是采用采样矩阵对图像进行亚采样,压缩采样在图像的稀疏域进行,被称为采样即压缩.许多真实信号本身并不是稀疏的,而是在一些变换域中有一个稀疏的表示.理论上 MR 压缩感知欠采样原理14如下:y=Fux+g(1)x CNFuCMN(M N)yg N(0,2)其中,表示具有 N 个像素的 MR 原图(磁共振最终图像),M 是测量值的长度,表示部分欠采样算子如傅里叶编码矩阵,依采样函数而定.观测在 K 空间进行,表示测量数据,使用高斯噪声向量模拟真实场景的自然数据噪声.x=sxxsx磁共振信号采集是逐行填充的,例如一幅 256256的磁共振图像,则至少要扫描 256 次(相位编码),每一次采集 256 个采样点(频率编码),形成了 K 空间.视野及成像的区域对固定大小的成像样品视野越大,使样品在图像中越小,大多数情况下使用正方形视野,也就是数字化采样一般是按照正方形点阵取样的.原始的MR 图像存在但未被表示出来,医疗检查过程中通过硬件设备获得了它的部分图像并对这部分图像进行了采样产生测量值.本文主要研究实 MR 图像.表示信号 在变换域 中可以被稀疏表示成 s.观测矩阵 与稀疏基 的乘积称为恢复矩阵.核磁共振成像 MRI信号重构的关键是在满足有限等距定理 RIP 的前提下“找到非稀疏信号 在中的稀疏表示”.通过观测矩阵求解稀疏系数 s 进而重构原始图像 是压缩感知进行图像恢复的关键,重建算法的限定性结构通常构建成以下优化问题:minssps.t.y=x,x=s(2)其中,p 是范数,当 p=1 时式(2)是一个凸优化问题,可以用梯度下降法、变量分离法等许多算法来求解.欠采样傅立叶编码矩阵算子有很多种,填充 K 空间的数据直接来源于MR信号,但它可以以任意的顺序进行采集填充.在过去的几十年里,K 空间的填充方式主要是逐条线性填充的笛卡尔方式,现在随机和辐射状K 空间填充方式越来越流行.高斯随机采样更接近于真实噪声,非常适合用于基于压缩传感(CS)的快速成像,有利于实现细节重建、图像边缘留存及抑制成像伪影.径向采样的 K 空间可以提供改善视场时间分辨率,减少总扫描时间,对于从径向采样的 K 空间重建图像,插值起着重要的作用,它可以在图像域或在频域(如重构方法)中实现.基于全局稀疏性的方法能够快速进行图像重建,但由于全局变换是固定的,其重建质量有待提升.近年来出现的有着出众重建效果的 BM3D 算法15利用其去噪器隐含的先验知识进行非局部相似性图像重建,基于寻找相似块进行协同滤波最后把同一个像素的不同滤波结果整合的思想16,提出构建以下 CSMRI 优化框架,对于第 i 次迭代,重构图像的优化模型如下:x(i)=argminxFuxy22+?xx(i1)?22(3)x(i)=argminx?x x(i)?22+x2(4)x(i1)x(i)其中,表示当前的估计值,、是正则化参数,表示 BM3D 去噪中使用的解析稀疏算子,表示中间图像.式(3)可通过求导近似求其闭式解.求解式(4)时,直接利用 BM3D 去噪器近似求解.但 BM3D 算法去噪先验有其自己的局限性,如果待处理的图像具有非自相似性,则使用 BM3D 去噪先验会带来过于光滑和不规则的结构等.因此,采用混合先验可能是一个更好的选择.2基于深度学习的稀疏表示模型 2.1 判别式学习先验去噪与目前已有的多种手动设计的正则先验相比,使用判别式学习方法自动学习出来的隐含在图像中的先验信息是一种更好的选择17.基于深度学习的卷积神经网络(CNN)采用数据驱动的策略,通过对数据的学习探究图像中隐含的先验,具有良好的特征提取能力,其性能优于许多手动设计的去噪算子如 BM3D.另外,由于 CNN 网络的运行往往结合 GPU(graphicsprocessingunit)的并行运算能力使其具有高效的运行速度,因此计 算 机 系 统 应 用http:/www.c-s-2023年第32卷第12期14专论综述SpecialIssue可加速其判别式学习能力.由此,从图像重构的性能和速度的角度考虑,可以尝试采用基于判别式策略学习的深度 CNN 去噪器18作为一种深度先验加入到 MRI图像重构框架.这种方法能够高效恢复图像依赖于迭代求解,因此负面结果是极大的消耗时间.2.2 网络结构模式虽然仍处于早期阶段,但深度神经网络的最新进展已经显示出了它们在计算机视觉和图像处理方面的潜力.Quan 等人19提出了一种新的基于深度学习的生成逆向模型用于快速和准确的 CS-MRI 重建.该模型是全残余卷积自编码器和生成式对抗网络(GANs)的一种变体,是单向部署在前馈网络,重建过程非常迅速,低至几十毫秒重建 256256 图像,即使采样率低至10%左右图像质量也是优越的.基于深度学习的 CS-MRI 重建过程将 CNN 与 CS相结合20构造出如下优化模型:minx12Fuxy22+R(x)+?x fCNN(xu?)?22(5)fCNNxu=FHuy其中,是通过参数 的 CNN 正向传播,是另一个正则化参数,CNN 网络被训练用以优化重建效果,表示训练后的 CNN 优化参数.表明对 y 在 K 空间随机 0 填充欠采样重建后的图像,H 表示 Hermitian转置.MRI 数据通常以复数形式编码了幅值和相位信息,在深度学习中,通常有两种表示方式,一种是将虚部置零,只保留实部信息,另一种是将实部和虚部当作两个通道共同作为输入.该论文提到为了减少计算量,使用第 1 种方式.深度神经网络的压缩感知 MRI 重建是一个由欠采样的 K 空间数据生成全重建图像的过程,在这样的网络结构模式中,这个能量最小化过程被神经网络的训练过程所取代.具体来说,压缩感知 MRI 的深度学习网络结构模式的基本思想通过网络训练使经验损失最小化21,优化模型如下.min1TTt=1?vvi?2=1TTt=1?viB(w,)(zi)?2(6)vB(w,)wv其中,z 为网络的输入数据,T 是训练数据个数,是输出数据即目标向量,是关于权值和可训练网络参数的神经网络,为标签图像.式(6)能够以损失最小化为约束通过使用非线性权值自适应求取网络的输出数据.为得到质量较好的图像,一般可采用“细量化粗采样图像缓变部分以避免假轮廓,粗量化细采样图像细节部分以避免模糊”的准则进行重构.将深度学习应用于图像重构需要考虑选择合适的网络架构实现网络复杂性和性能之间的权衡,可以在 K 空间域和图像域中进行端到端训练实现这一平衡.3融合深度学习和传统方法的压缩感知核磁共振成像深度技术通过学习对最终重建提供快速恢复,然而由于数据驱动的重构行为不能保证根据领域知识的原则约束进行重建,其重建过程的可靠性有待论证.在本文中,我们应用了一个组装 CS-MRI 原则模块的深度框架,该框架将学习策略与传统重建能量的迭代求解器相结合.该框架嵌入了一个最优的条件检查机制,促进了高效和可靠的重建.由多种先验知识构建的去噪器所表达的稀疏表示模型能够将多种去噪器蕴含的先验知识融合到图像重建中,从而提升低采样率下 MRI的重建质量.MR 过程中会产生有多个来源的高度复杂且不可预测的噪声,因此消除现实世界的盲噪声即进行盲去噪是图像恢复中一个重要的问题.CBDNet 深度学习算法22可以从有噪声的样本学习图像中的噪声分布并实现盲去噪.CBDNet 由 A 部分噪声估计子网络(CNNE)和 B 部分非盲去噪子网络(CNND)两部分组成,非对称学习的噪声估计子网络可以有效抑制噪声,利用这种非对称损失来提高对噪声的泛化能力.CNND 结构包括 5 层全卷积(Conv+ReLU),32 个卷积核,尺寸33.CNNE 使用了一个没有 pooling(池化)和 BN(batchnormalization)操作的 5 层全卷积网络;CNND则使用了一个 16 层的 U-Net,学习残差映射来预测结果.经实验验证,结合使用 CNNE 和 CNND 残差学习对高斯去噪操作不仅能保证网络训练快速稳定,而且还能获得效果较好的去噪结果.CBDNet 在训练过程中同时结合了合成噪声和真实噪声图像进行学习,实现了交互式去噪,提高了网络的去噪能力,非常适合应用于医学核磁共振成像领域.深度去噪器和传统去噪器各有优势,因此本文试图利用基于两种去噪器联合解决现有 CSMRI 算法重建质量低的问题.本文所研究的融合盲去噪 CNN23算法 CBDNet 和传统降噪器 BM3D 算法的总体构造如2023年第32卷第12期http:/www.c-s-计 算 机 系 统 应 用SpecialIssue专论综述15图 1 所示,AB 两部分构成了 CBDNet 降噪器,BM3D降噪器由基础估计和最终估计两部分构成,在联合降噪的过程中通过某种优选机制以及滤波处理等流程获得最终的估计图像.测量数据基础估计估计图像优选机制滤波处理BM3D 降噪图最终估计ABCBDNet 降噪图图 1融合深度学习和传统方法的压缩感知核磁共振成像原理流程本文所研究的联合降噪具体构建以下对应的优化模型:minx,sFuxy22+Wxs2+s p(7)xyWWTW=E I(x)Ws式(7)中第 1 项为保证重建图像 与测量数据 相匹配的数据保真项;第 2 项、第 3 项为构建的基于深度去噪器和传统去噪器的稀疏表示模型.为正则化参数用以平衡两项的比重,表示对应的解析稀疏算子(),是正则化参数,是经去噪后的图像在下的稀疏系数,与 是相似的.因此本文研究的混合降噪过程中稀疏表示模型如下:minsWIs22+s p(8)G=max?Ic(x(i1)G(i1)?2,?Ib(x(i1)G(i1)?2(9)=TWI(x(i1),(10)p IIc(x(i1)IcG(i1)G(0)ITWI(x),s式(8)中 值本研究中取为 0,稀疏系数 的解需要依托迭代硬阈值法获取.依据式(9)求取,其中表示用深度降噪器,用其获得降噪图像,表示本轮图像恢复前的恢复图像,表示初始恢复图像.表示本轮降噪过程中由式(9)获得的恢复图像.式(10)中根据实验效果选用硬阈值算子进行求解.稀疏表示系数 确定后,采用式(11)求解估计图像:minxFuxy22+Wxs2(11)对最小化代价函数式(11)进行求导推演出式(12),表示对去噪器去噪后的图像进行某稀疏变换,代表对图像进行某稀疏变换,求出估计图像的系数变换,在进行反变换求出最终的估计图像.Fx(i)=f(I(x(i1)+f(G0)f(I(x(i1)(+1)I(12)本文算法如算法 1 所示.算法 1.CBDNet-BM3DyG(0)=FHuy 输入:测量数据,初始化估计图像,参数,xi输出:由迭代产生的重构图像G(i1)Ic(x(i1)IcIb(x(i1)Ib1)对估计图像利用深度去噪器进行滤波获得降噪图像,同时利用传统去噪器进行滤波获得降噪图像;IcIbG(i1)I2)由式(9)对降噪图像与依其求其与的相似性,二者 max择优获取综合降噪恢复图像;3)由式(10)更新稀疏系数;G(i1)G(0)Fx(i1)4)在其他变量视为已知量的情况下,在 K 空间由式(12)通过和求取估计图像的傅立叶变换得到:;Fx(i1)xiG(i)5)对进行 IFFT 变换计算新的估计图像,也即;6)返回步骤(1),直到迭代次数 i 符合设计要求.4实验结果实验将对本文提出的算法(后续用 CBDNet-BM3D算法指代)与 5 种不同 MR 图像重构算法的性能进行程序运行结果对比评价,实验将对不同类型的 MR 图像在不同采样方式的不同欠采样率进行图像恢复,以论证该方法的高效性、鲁棒性及泛化性等.仿真实验软件环境:Windows10 系统,MatlabR2019b,Matconvnet 深度学习框架.仿真实验硬件环境:Intel(R)Core(TM)i9-9900CPU,主频 3.10GHz,内存 16.0GB,显卡为 NVIDIAGeForceRTX2060.图像的重构质量通过采用信号重建质量测量指数峰值信噪比(peaksignaltonoiseratio,PSNR)和衡量幅图像相似度的指标结构相似性(structuralsimilarity,SSIM)作为客观指标视觉效果作为主观指标来进行评价.本文将计 算 机 系 统 应 用http:/www.c-s-2023年第32卷第12期16专论综述SpecialIssue所提算法 CBDNet-BM3D 与现有的 4 种算法进行对比.这 4 种算法分别为经典的传统算法(BM3D-MRI24)和 3 种基于深度学习的算法(FFDNet18、CSMRI-FFDNet25、CBDNet23).本文实验所使用的图像块大小为,算法中的所有测量值分别采用多种采样频率的高斯随机采样矩阵、径向采样矩阵、笛卡尔采样矩阵产生,再使用 ifftshift去中心化将低频和高频的位置还原(将白色的、亮的移动到图像四周,将黑色的、暗的移动到图像中心),实践证明变换后所产生的采样矩阵恢复效果比未变换前要好,采样矩阵如图 2 所示,图 2(a)分别为 3 类采样矩阵,图 2(b)分别为 3 类采样矩阵的去中心化结果,采样率均为 20%.图 3 是实验中所用到的 6 幅 256256的核磁共振图像.(a1)Random sampling(b1)图(a1)去中心化(a2)Ratial sampling(b2)图(a2)去中心化(a3)Cartesian sampling(b3)图(a3)去中心化图 2高斯随机采样、径向采样和笛卡尔采样(sampling20%)(a)brain2d(b)shouder(c)head(d)brain3d(e)bust(f)bone图 3实验数据集表 1 是图像数据集中 brain 在不同采样率和采样算法下的 PSNR 和 SSIM 值,实验结果表明高斯随机采样的图像恢复效果是最好的.表 1 显示了在采样率分别为 10%、20%和 30%的高斯随机采样矩阵、径向采样矩阵、笛卡尔采样矩阵时使用 5 种不同重构算法得到的 PSNR 和 SSIM 值的均值和标准差.表 1 数据显示,当采样算子使用高斯随机矩阵时,CBDNet-BM3D 算法的图像重构结果相对于其他 4 个算法呈现较优的状态.在 10%的采样率下,CBDNet-BM3D 算法分别高于 BM3D-MRI、FFDNet、CSMRI-FFDNet 和CBDNet 算法 1.9383dB、3.1657dB、3.0823dB 和0.8092dB;在 20%的采样率下,CBDNet-BM3D 算法分别高于 BM3D-MRI、FFDNet、CSMRI-FFDNet 和CBDNet 算法 1.1389dB、0.9457dB、0.8401dB 和0.0888dB;在 30%的采样率下,CBDNet-BM3D 算法分别高于 BM3D-MRI、FFDNet、CSMRI-FFDNet 和CBDNet 算法 1.818dB、1.3148dB、1.1024dB 和0.1622dB.当计算混合降噪结果时,先分别求基于深度去噪的滤波图像与基于 BM3D 去噪的滤波图像与上一轮循环中估计图像间的距离,再使用 max 函数取两者中差异大的视为修复度更高.实验中外循环取 66 次,基于深度去噪修复度高的次数与基于 BM3D 去噪修复度高的次数比约为 47:19,数据表明深度去噪修复效果相对优于 BM3D 去噪修复效果,并且二者结合应用修复效果更好.为了评估本文提出算法对测量噪声的鲁棒性,我们在具有加性高斯噪声的 brain3d 图上进行了实验,如图 4 所示.CT 讲的是密度,细节清晰非常重要,brain3d图是 T1 加权像,骨皮质是黑色暗的低信号,皮下脂肪是白色亮的高信号,脑室内的液体(脑脊液)是黑色的.对本文所提出 CBDNet-BM3D 算法及其他 4 种算法的重构图像进行视觉对比,其中 8%的原始 K 空间样本被保留.图 4(b)图 4(f)显示了不同算法对同一副图像进行采样后恢复的结果图,黄色框及蓝色框是对应的感兴趣区域的放大状态.从图 4 中可看出,基于传统的BM3D-MRI 算法以及深度学习 CBDNet 算法重构效果具有一定的优势,重构图像的边缘以及纹理保留了更多的细节信息,且低伪影使得图像更为清晰,实验结果进一步支持了我们提出的方法在存在较高测量噪声的情况下重建效果有很好的竞争能力.2023年第32卷第12期http:/www.c-s-计 算 机 系 统 应 用SpecialIssue专论综述17表 1不同测量矩阵不同采样率下 5 种算法的 PSNR 和 SSIMSamplingApproachessampling:10%sampling:20%sampling:30%PSNR(dB)SSIMPSNR(dB)SSIMPSNR(dB)SSIMRandomBM3D-MRI21.24690.61370.66880.045433.31941.92410.96120.015937.87401.17380.98550.0043FFDNet20.67490.04170.62490.001533.02242.41430.95720.020837.99071.56020.98600.0052CSMRI-FFDNet21.85321.13660.66410.037733.09162.45080.95760.020938.17541.5880.98650.0051CBDNet23.05430.06460.70460.002933.83532.45830.96340.018539.02021.68330.98860.0045CBDNet-BM3D23.37350.42540.66910.03633.92492.45760.96420.018139.17421.69150.98890.0044RadialBM3D-MRI24.28550.61370.72250.045431.93301.92410.96130.015934.66391.17380.98410.0043FFDNet25.72420.04170.73430.001531.45472.41430.95220.020834.49351.56020.98160.0052CSMRI-FFDNet25.77951.13660.73670.037731.48252.45080.95320.020934.57851.5880.98250.0051CBDNet25.55740.06460.74350.002931.86832.45830.95850.018534.95631.68330.98420.0045CBDNet-BM3D25.43740.42540.74270.03631.92262.45760.96020.018135.01601.69150.98500.0044CartesianBM3D-MRI26.19655.62120.81180.131826.90115.08690.85160.104929.29663.19860.93000.0488FFDNet25.65105.31570.77410.108326.51245.04880.82690.099328.85853.4590.91500.053CSMRI-FFDNet25.69225.37670.77650.112526.51975.08120.82800.102428.86903.45540.91540.0534CBDNet25.83185.49640.79630.127726.74415.03470.84340.102729.11693.20240.92290.0475CBDNet-BM3D25.97525.54460.80450.128726.78925.05190.84770.102929.18403.19560.92590.0483(a)原图(b)BM3D-MRI 算法(c)FFDNet 算法(d)CSMRI-FFDNet 算法(e)CBDNet 算法(f)CBDNet-BM3D 算法302520151050302520151050302520151050302520151050302520151050302520151050图 45 种算法 brain3d 重建图、重建结果局部放大图及重建误差图(Randomsampling5%)当采样率为 30%高斯随机采样、迭代次数为660 时运行时长如表 2 所示.虽然数据显示各算法中运行时间 BM3D 算法最少,但如图 5 所示,取外循环为2166 次所对应的文中所提 6 种算法的 PSNR 值形成曲线图,BM3D 算法峰值信噪比值程下抛物线状,说明此算法图像恢复效果不稳定,FFDNet 及 CSMRI-FFDNet算法 PSNR 值上升后成稳定状态,但其图像恢复效果不够理想.数据统计可以看出 CBDNet-BM3D 算法优于其他 4 种算法且非常稳定.重构算法获得的估计图像PSNR值差距大都在3,3范围内,数值非常近似,为更好地显示此类型的数据,图 6 使用了堆积折线图展示随着采样率的提升各算法PSNR 值的上升率.在高斯随机采样率分别为 1%、5%、10%、15%、20%、25%、30%、35%、50%的取值下,从折线图可以看出,最上面代表本文所提出算法的曲线斜率是最大的,说明采样率越高本文所提出算法的图像恢复效果会更加显示其优势.表 2采样率 30%迭代次数 660 时的 runtime 对比(s)Methods BM3D FFDNet CSMRI-FFDNet CBDNet CBDNet-BM3DRuntime 39.3822 78.1435235.8293363.6485419.4434计 算 机 系 统 应 用http:/www.c-s-2023年第32卷第12期18专论综述SpecialIssue2138.038.539.039.540.040.541.041.524 27 30 33 36 39 42外循环次数PSNR(dB)45 48 51 54 57 60 63 66BM3DFFDNETCSMRI-FFDNETCBDNETPaper图 530%高斯随机采样、外循环为 2166 次所对应的文中所提 5 种算法的 PSNR 值曲线图箱线图显示 CBDNet-BM3D 算法与其他 4 种对比算法的 PSNR 值,方框内的线表示中值,方框代表第 25 到第 75 个百分位数.方框外线代表第 10 个和第 90 个百分位数,上下水平线分别代表最小值和最大值.实验数据表明,深度学习算法与传统重构算法结合的算法效果均要优于单纯使用深度学习算法的重构结果.表 1 和图 7(b)显示 BM3D 算法的图像重构效果从数据上看非常有优势,本文提出的方法是深度学习与传统重构算法的混合式重构算法,图像恢复效果均优于其他算法,这些结果证明了所提出的算法的有效性.105010015020025051015PSNR(dB)20高斯随机采样率(%)25303550BM3D-MRIFFDNetCSMRI-FFDNetCBDNetPAPER图 6实验数据集 6 幅图在 5 种高斯随机采样下的 PSNRBM3D-MRIFFDNetCSMRI-FFDNetCBDNetPAPERBM3D-MRIFFDNetCSMRI-FFDNetCBDNetPAPER26283032343638404244460.9600.9650.9700.975(a)PSNR 箱线图(b)SSIM 箱线图0.9800.9850.9900.9951.000PSNR(dB)SSIM图 7实验数据集 6 幅图高斯随机采样 8%的 PSNR 和 SSIM5结论本研究提出了利用 CNN 结合传统方式进行欠采图像的重构,针对不足采样的 K 空间数据,在经典的CNN 算法 CBDNet 算法的基础上,提出了基于深度神经网络去噪先验和 BM3D 算法的块压缩感知重构法,得到以下结论.(1)该算法实现了抑制伪影、快速重建高质量 MRI.(2)实验结果表明,该理论有效降低了伪影并保留了细节信息,加强了重构效果,印证了我们提出的方法比单纯的深度学习方法或传统的块匹配方法更接近目标图像.(3)与此同时,通过采用 GPU 的加速运算,算法的计算复杂度较使用单一算法并未增加很多.综上所述,证明混合式 CNN 算法在 MR 图像领域非常具有现实应用价值.参考文献QiuJX,LiuJ,BiZX,et al.Aninvestigationof2Dspinemagneticresonanceimaging(MRI)withcompressedsensing(CS).SkeletalRadiology,2022,51(6):12731283.doi:10.1007/s00256-021-03954-x1YuanM,Ya