CN108596994B - 一种基于深度学习和数据自洽的磁共振弥散加权成像方法 - Google Patents
一种基于深度学习和数据自洽的磁共振弥散加权成像方法 Download PDFInfo
- Publication number
- CN108596994B CN108596994B CN201810376066.9A CN201810376066A CN108596994B CN 108596994 B CN108596994 B CN 108596994B CN 201810376066 A CN201810376066 A CN 201810376066A CN 108596994 B CN108596994 B CN 108596994B
- Authority
- CN
- China
- Prior art keywords
- data
- network
- image
- consistent
- phase
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000013135 deep learning Methods 0.000 title claims abstract description 23
- 238000002597 diffusion-weighted imaging Methods 0.000 title claims abstract description 21
- 238000003384 imaging method Methods 0.000 claims abstract description 54
- 238000012549 training Methods 0.000 claims abstract description 30
- 238000013507 mapping Methods 0.000 claims abstract description 18
- 238000002592 echocardiography Methods 0.000 claims abstract description 6
- 238000009792 diffusion process Methods 0.000 claims description 30
- 230000005284 excitation Effects 0.000 claims description 25
- 238000013527 convolutional neural network Methods 0.000 claims description 24
- 230000008569 process Effects 0.000 claims description 21
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000010587 phase diagram Methods 0.000 claims description 10
- 230000035945 sensitivity Effects 0.000 claims description 10
- 230000006870 function Effects 0.000 claims description 8
- 238000013528 artificial neural network Methods 0.000 claims description 7
- 239000006185 dispersion Substances 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 5
- 238000012935 Averaging Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 238000012512 characterization method Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 4
- 238000011084 recovery Methods 0.000 claims description 4
- 238000003786 synthesis reaction Methods 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 3
- 239000003550 marker Substances 0.000 claims description 3
- 230000004913 activation Effects 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 description 16
- 230000033001 locomotion Effects 0.000 description 10
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 6
- 238000012360 testing method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 229920002521 macromolecule Polymers 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于深度学习和数据自洽的磁共振弥散加权成像方法,涉及磁共振弥散加权成像领域;其包括:1:通过由数据自洽层、CNN网络和相位约束层组成的网络模块重复叠加构建网络后,执行包含预扫描和导航回波的序列确定网络的输入数据和训练标记数据;2:将训练标记数据作为目标,将输入数据对应的图像输入已构建的网络通过反向传播训练得到网络的输入输出映射关系;3:执行包含预扫描和无导航回波的序列获得成像信号和数据自洽卷积核,将其输入网络通过前向传播映射获得输出图像完成重建;解决了现有基于传统的弥散加权成像受限各因素导致成像分辨率低的问题,达到了提高网络学习能力,提高成像分辨率的效果。
Description
技术领域
本发明涉及磁共振弥散加权成像领域,尤其是一种基于深度学习和数据自洽的磁共振弥散加权成像方法。
背景技术
磁共振弥散成像技术,是一种依赖体内水分子随机运动而提供图像对比度的崭新技术。组织中水分子的扩散符合随机的热运动模型,扩散的幅度和方向受到生物膜和组织中生物大分子的影响。当梯度磁场存在时,水分子的扩散运动会引起磁化矢量的失相位,导致磁共振信号的降低。磁共振信号降低的程度依赖于组织类型、结构、物理和生理的状态及微环境。上述过程中,专门用来影响水分子热运动的梯度磁场称为弥散敏感梯度。弥散敏感梯度能够显著提高各种序列对水分子随机热运动的敏感性,进而帮助观察水分子的扩散特性,但是这种梯度磁场也对其他类型的运动如头部运动十分敏感。单次激发弥散成像技术,在一次信号激发后,采集用于成像的所有数据;这种方式能有效的缩短扫描时间,有利于避免引入更多宏观运动对图像造成影响。但是,单次激发扫描方式采用的回波链较长,容易引起磁敏感伪影和几何形变;单次采集所能够得到的数据有限,因此图像的分辨率较低,不利于细微结构的诊断;为了减少图像伪影和几何形变,以及有效的提高图像空间分辨率,采用的策略是使用多次激发磁共振弥散成像技术;该技术面临的主要挑战是有效的处理不同激发所采集数据之间存在的由于宏观运动导致的相位误差。根据数据采集方式的不同,可以将宏观运动矫正分为两类:第一类需要在正常数据采集之前,采集导航回波信号,这个信号会被用来矫正接下来每个激发所采集的成像数据;第二类不需要采集导航回波信号,而是通过每次激发数据之间的关系来矫正相互的相位。同采集导航回波的方式相比,不需要导航回波的采样方式有更高的数据采集效率,也能避免导航回波与实际成像数据之间失配的问题。
2002年,M.A.Griswold等人提出了GRAPPA技术。GRAPPA技术同样利用多通道相控阵线圈,但不需要额外计算线圈的空间灵敏度分布,而是利用K空间数据的相关性恢复欠采样的数据,GRAPPA技术在对K空间进行欠采样的同时,会保留K空间中心的数据为全采样,这部分数据叫做ACS信号;然后,根据ACS信号计算出代表K空间数据之间相关性的卷积核;最后,使用估计出的卷积核以及已采集数据,重建欠采样的数据,进而得到完整的K空间数据,对应的就是没有卷褶的图像。GRAPPA技术避免了精确的估计线圈灵敏度的空间分布,不过其重建质量也直接受到接收线圈以及卷积核的影响。
2010年,Michael Lustig等人提出了基于信号域数据自洽的SPIRiT技术。该技术将欠采样K空间的重建问题定义为一个优化问题,优化的目标是保证ACS信号与欠采样信号有最好的自洽性,SPIRiT技术在K空间欠采样的同时,也需要保留K空间中心的ACS数据为全采样,不过SPIRiT可以兼容任意的欠采样模式;然后,根据ACS信号计算出代表K空间局部数据自洽性卷积核;最后,使用估计出的卷积核以及已经采集的数据,通过迭代计算的方式重建欠采样的数据,进而得到完整的K空间数据。该技术对于任意形式的K空间填充都是适用的,而且能够在求解过程中方便的融合先验知识,能够有效的提高重建的精确性以及获得更高的扫描加速因子。依赖于并行成像的无导航回波多次激发弥散成像,其临床不足在于受到并行成像的性能的限制,其中激发次数越多,则对于单次激发数据而言欠采样率就越高,这样对于并行成像的要求就越高,并行成像只能处理一定程度欠采样的数据,因此整体上就限制了激发次数的提升,另一方面并行成像数据恢复能力本身收到接收通道的数目、几何形状基于算法本身的限制,因此现有并行成像激发次数不能太多从而导致最终图像的分辨率难以继续提升。
近些年,以卷积神经网络为代表的深度学习在计算机视觉、语言理解等领域获得了重大进展。近两年,为了获得更高的加速比以及更好的重建效果,将深度学习的思想用于磁共振稀疏重建进而更有效的恢复欠采样数据的技术不断涌现。2016年,Wang等人提出将深度学习用于磁共振稀疏重建的技术。该技术首先构造一个卷积神经网络,然后使用端到端的训练模式使得该网络能够将输入的欠采样数据转化为全采样数据,最后,将网络的输出结果作为压缩感知重建的初始值或者作为压缩感知重建的方程中新引入的正则化项,优点在于首次将基于神经网络的重建优势应用到实际中,缺点在于,只是将神经网络重建作为辅助的手段,后期仍然要基于数学模型的重建,这就从根本上限制了算法的性能提升;网络结构较浅、且简单,难以获得更高的性能;2017年,Jo Schlemper等人提出一种串联的深度学习网络模式用于磁共振稀疏重建。该技术将磁共振稀疏重建定义为图像域消除伪影的学习过程,构造的卷积神经网络能够通过训练过程学会如何消除由于欠采样而导致的图像伪影。该技术首先构造出较浅的卷积神经网络,然后在该网络后面添加专门设置的数据一致性层,用于提供已采样数据的信息;最后,将上述两个基本结构重复的串联起来,构造出更深的卷积网络,但是技术尚不能支持多通道信号
基于深度学习的磁共振稀疏重建技术显示出了明显优于传统并行成像、稀疏低秩或者压缩感知技术的性能。但是,目前尚没有相关研究将深度学习的强大数据恢复能力与多次激发磁共振弥散技术相结合,以更好的处理多次激发之间的相位误差,获得更高的分辨率。
发明内容
本发明的目的在于:本发明公开了一种基于深度学习和数据自洽的磁共振弥散加权成像方法,解决了现有基于传统的弥散加权成像受限并行成像性能因素导致成像分辨率低的问题。
本发明采用的技术方案如下:
综上所述,由于采用了上述技术方案,本发明的有益效果是:
1.本发明使用数据自洽层使得该网络结构能够处理来自多通道的磁共振数据且更有效地利用磁共振多通道数据本身的冗余特性,帮助基于深度学习的卷积神经网络建立更稳定、精准的端到端映射关系,从根本上提高多次激发磁共振重建的质量,使用相位约束层处理来自多次激发的磁共振弥散加权数据,进一步利用多次激发数据之间的冗余信息,提升网络重建的性能;将数据自洽层、相位约束层以串联的形式融入序贯模型的卷积神经网络结构能够改善网络结构的泛化能力、提升网络训练及测试的稳定性;解决了现有基于传统的弥散加权成像受限各因素导致成像分辨率低的问题,达到了提高网络学习能力,提高成像分辨率的效果;
2.本发明采用卷积神经网络处理单次激发的多通道信号相比传统的基于模型的方法采用卷积网络的方法能够更好的恢复欠采样数据,能够获得更好的图像质量;多个数据自洽层共享基于磁共振先验知识计算出的自洽卷积核和已采集K空间数据数据,自洽层保证了整个网络的多通道信号处理能力,更有效的利用多个通道之间的信号冗余,生成更为精准的图像;
3.本发明能够支持无导航回波磁共振弥散成像提高激发次数,并且有效的处理不同激发数据之间的相位误差,获得更高的图像分辨率、更高的信噪比;
4.本发明使用同样的磁共振弥散信号,相对于传统的方法能够生成更为精准的图像;
5.本发明提出的卷积神经网络中的所有超参数都不是人为设定,而是经过大量数据训练得到的,对于临床复杂的结构,包括扫描部位、图像信噪比等,都能够提供更为稳定的图像重建质量。
6.数据自洽层、相位约束层和已采集K空间数据按照序贯模型排列,一方面通过增加网络结构的深度提高网路的学习能力,另一方面,重复的使用已知信息监督整个学习过程,能够保证学习的稳定性,提高网络的泛化能力。
附图说明
本发明将通过例子并参照附图的方式说明,其中:
图1是本发明的网络示意图;
图2是本发明步骤1的方法流程图;
图3是本发明序列示意图;
图4是本发明网络输入输出示意图;
图5是本发明采集序列生成网络标记数据的流程图。
具体实施方式
本说明书中公开的所有特征,或公开的所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以以任何方式组合。
下面结合图1-5对本发明作详细说明。
步骤1:通过由数据自洽层、CNN网络和相位约束层组成的网络模块重复叠加构建网络后,执行包含预扫描和导航回波的多次激发弥散加权序列获得多次激发成像信号和数据自洽卷积核并确定网络的输入数据和训练标记数据;
步骤2:将训练标记数据作为目标,将输入数据对应的图像、多次激发成像信号和数据自洽卷积核输入已构建的网络通过反向传播训练得到网络的输入输出映射关系;
步骤3:执行包含预扫描和无导航回波的多次激发弥散加权序列获得成像信号和数据自洽卷积核,将成像信号对应的图像输入已完成训练的网络并融入数据自洽卷积核通过前向传播映射获得输出图像完成重建。
每个CBR单元至少包含一个卷积层、一个规范化层和一个非线性激活层;CNN网络前半部分的CBR单元还加入汇聚层,后半部分加入反汇聚层。反向传播过程:如图4所示,一个神经网络结构包含三个层次,输入层,包含三个输入i1,i2,b1;中间层,包括h1,h2,b2;输出层,包括o1,o2;层与层之间连接的权重为标记的w1,w2,反向传播:在神经网络训练阶段,给定已知的输入和输出,以最小化损失函数为目标,通过链式法则递归求解权重w参数和b参数;因此,反向传播即为通过链式法则计算参数w和b相对于损失函数的梯度确定最优的参数。前向传播过程:已知的测试集数据为多次激发的磁共振弥散加权图像,以及网络结构各个层次的共享数据;已知权重w参数和b参数即前期训练阶段求得的,给定输入i1和i2即测试集数据,通过网络的逐层计算,得到精确的输出结果即未知映射数据完成磁共振高分辨率弥散加权成像的重建,因此前向传播即为已经完成训练的神经网络将给定的输入映射为输出的过程。
本发明将数据自洽层、相位约束层以串联的形式融入序贯模型的卷积神经网络结构能够改善网络结构的泛化能力、提升网络训练及测试的稳定性;解决了现有基于传统的弥散加权成像受限各因素导致成像分辨率低的问题,达到了提高网络学习能力,提高成像分辨率的效果。
实施例1
步骤1包括如下步骤:
步骤1.1:将由数据自洽层、CNN网络和相位约束层组成的网络模块依次叠加初步完成网络构建,每个网络模块中数据自洽层与CNN网络多对多连接,CNN网络与相位约束层多对一连接;
步骤1.2:采集多次激发弥散加权序列中的预扫描数据,基于与扫描数据生成矫正矩阵和数据自洽方程,如图2中的101-102;
步骤1.3:基于矫正矩阵计算所有接收线圈的灵敏度分布如图2中的104,基于数据自洽方程计算卷积核如图2中标记的103;
步骤1.4:采集多次激发弥散加权序列中的导航回波信号和成像信号,成像信号在第一个180度脉冲后采集,导航回波在第二个180度脉冲后采集,虚线表示弥散敏感梯度,如图3所示,将成像信号作为已构建网络的输入数据,通过导航回波信号矫正多次激发成像信号之间的相位差获得磁共振弥散模值图像并将其作为已构建网络的训练标记数据,如图5所示,矫正相位误差如公式1所示;
其中,g为多次激发成像信号对应的多组图像,E为由线圈灵敏度和导航回波对应的相位图计算出来的编码矩阵,Ψ为所有接收通道的相关矩阵,H为厄米特矩阵运算符,为磁共振弥散加权图像。矩阵E经过厄米特算符后,相当于对矩阵E进行了共轭转置操作,属于数学领域的公知,在此不赘述。
步骤2包括如下步骤:
步骤2.1:将输入数据转换为图像输入已构建的网络如图1中标记的401,将已采集的多次激发成像信号和数据自洽卷积核融入已构建网络如图1中的402和403,通过数据自洽层进行数据自洽将每次激发的多通道图像映射为卷褶伪影减弱的多通道图像如图1中标记的404;
步骤2.1中的数据自洽过程包括如下步骤:
步骤2.1.1:第一层数据自洽层输入为输入图像,其他数据自洽层输入为相位约束层输出即每次激发的多通道图像Ipcon,将其通过傅里叶变换获得K空间信号Spcon,计算如公式4所示:
Spcon=FIpcon;
步骤2.1.2:将K空间信号Spcon进行数据投影,若K空间位置上的数据已被采集,则填充为已采集数据Sacq,若未被采集则计算Spcon后构造新的完整K空间数据Sdp,计算如公式5所示:
步骤2.1.3:将已采集数据和融入的数据自洽卷积核进行卷积运算获得新K空间数据Ssc,计算如公式6所示:
Ssc=fsc(Sdp,G)=G*Sdp;
步骤2.1.4:将新的K空间数据Ssc通过傅里叶逆变换转换为图像域的多通道图像Isc完成数据自洽,计算如公式7所示:
Isc=F-1Ssc。
步骤2.2:将卷褶伪影减弱的多通道图像通过CNN网络映射获得无卷褶或卷褶减弱的图像Icnn如图1中标记的405;
步骤2.3:将无卷褶或卷褶减弱的图像Icnn输入相位约束层经过激发平均和相位恢复完成相位约束如图1中标记的406,将训练标记数据作为目标代入最小化损失函数进行反向传播训练获得输入输出的映射关系,最小化损失函数如公式2所示:
其中,Xi为成像信号对应图像,Yi为已知的训练标记图像,θ为网络表征参数。
步骤2.3中的相位约束包括如下步骤:
其中,Icnn,i为每次激发对应图像,Ns为激发次数;
步骤2.3.2:将不同激发图像Icnn,i减去对应的具有平滑性质的相位信息后,采用取平均的方式进行激发合,计算如公式9所示:
步骤2.3.3:将每次激发对应的具有平滑性质的相位图重新附加到激发合成后的图像中,得到相位约束后的各次激发的图像Ipcon,i,计算如公式10所示:;
步骤3包括如下步骤:
步骤3.1:根据步骤2得到的映射关系确定网络的结构和运算特征,扫描包含无导航回波的多次激发弥散加权序列得到成像信号和数据自洽卷积核;
步骤3.2:将成像信号转换为图像后输入已完成训练的网络进行前向传播根据输入输出映射关系即公式3得到磁共振弥散模值图像完成重建,如图1中标记的407,公式3计算如下:
f(θ,X):X→Y
其中,X为成像信号对应的图像,Y为磁共振弥散模值图像,θ为网络表征参数。
本发明将深度学习网络应用于多次激发的磁共振弥散加权成像领域中,需要克服的问题:一、多次激发信号的合成问题,通过训练得到超参数设定克服了将相位约束算法转换为网络结构的形式的技术难点,通过设置梯度传递和反向传播克服了如何将相位约束层融合至卷积网络中的技术难点,通过设置相关约束层实现有效地合成多次激发信号并产生了意外的效果:超参数通过训练得到并非人为设定,利于提高网络的泛化能力使其具有稳定的性能;二、通过将数据自洽层、CNN网络和相关约束层层串联组成序贯模型,实现多通道的磁共振弥散加权成像,同时实现快速高倍欠采样,克服了处理图像域和K空间信号域的转换、前向传播和反向传播过程要在两个域中交替进行的问题,通过数据自洽层实现多个通道的处理的问题,解决了现有基于深度学习的成像方法只能处理单通道数据的问题。另一方面在训练阶段采用包含无导航回波的序列,测试阶段采用包含导航回波的序列,通过线下调试保证导航回波获得好的重建结果,然后把这个结果作为训练的目标;用无导航回波信号作为网络输入,通过重建神经网络前向传播,使得无导航回波网络重建方案能够得到又好有稳定的重建效果,同时具有较短的扫描时间,克服了采用无导航回波无法达到采用导航回波的效果的偏见,保证了基于深度学习网络的得到的重建效果的稳定性,同时缩短了扫描时间。
Claims (8)
1.一种基于深度学习和数据自洽的磁共振弥散加权成像方法,其特征在于:包括如下步骤:
步骤1:通过由数据自洽层、CNN网络和相位约束层组成的网络模块重复叠加构建网络后,执行包含预扫描和导航回波的多次激发弥散加权序列,获得多次激发成像信号和数据自洽卷积核,并确定网络的输入数据和训练标记数据;
步骤2:将训练标记数据作为目标,将输入数据对应的图像、多次激发成像信号和数据自洽卷积核输入已构建的网络,通过反向传播训练得到网络的输入输出映射关系;
步骤3:执行包含预扫描和无导航回波的多次激发弥散加权序列,获得成像信号和数据自洽卷积核,将成像信号对应的图像输入已完成训练的网络并融入数据自洽卷积核,通过前向传播映射获得输出图像,完成重建。
2.根据权利要求1所述的一种基于深度学习和数据自洽的磁共振弥散加权成像方法,其特征在于:所述步骤1包括如下步骤:
步骤1.1:将由数据自洽层、CNN网络和相位约束层组成的网络模块依次叠加初步完成网络构建,每个网络模块中数据自洽层与CNN网络多对多连接,CNN网络与相位约束层多对一连接;
步骤1.2:采集多次激发弥散加权序列中的预扫描数据,基于预扫描数据生成矫正矩阵和数据自洽方程;
步骤1.3:基于矫正矩阵计算所有接收线圈的灵敏度分布,基于数据自洽方程计算卷积核;
步骤1.4:采集多次激发弥散加权序列中的导航回波信号和成像信号,将成像信号作为已构建网络的输入数据,通过导航回波信号矫正多次激发成像信号之间的相位差获得磁共振弥散模值图像并将其作为已构建网络的训练标记数据,矫正相位差如公式1所示;
3.根据权利要求1所述的一种基于深度学习和数据自洽的磁共振弥散加权成像方法,其特征在于:所述步骤2包括如下步骤:
步骤2.1:将输入数据转换为图像输入已构建的网络,将已采集的多次激发成像信号和数据自洽卷积核融入已构建网络,通过数据自洽层进行数据自洽,将每次激发的多通道图像映射为卷褶伪影减弱的多通道图像;
步骤2.2:将卷褶伪影减弱的多通道图像通过CNN网络映射获得无卷褶或卷褶减弱的图像Icnn;
步骤2.3:将无卷褶或卷褶减弱的图像Icnn输入相位约束层,经过激发平均和相位恢复完成相位约束,将训练标记数据作为目标代入最小化损失函数,进行反向传播训练获得输入输出的映射关系,最小化损失函数如公式2所示:
其中,Xi为成像信号对应的图像,Yi为已知的训练标记图像,θ为网络表征参数,n为单次训练过程中,输入的图像的组,即表示此次训练过程使用了n组图像,f(θ,Xi)表示以Xi为输入,θ为神经网络过程的输出结果。
4.根据权利要求1所述的一种基于深度学习和数据自洽的磁共振弥散加权成像方法,其特征在于:所述步骤3包括如下步骤:
步骤3.1:根据步骤2得到的映射关系确定网络的结构和运算特征,扫描包含无导航回波的多次激发弥散加权序列得到成像信号和数据自洽卷积核;
步骤3.2:将成像信号转换为图像后输入已完成训练的网络进行前向传播,根据输入输出映射关系即公式3得到磁共振弥散模值图像,完成重建,公式3计算如下:
f(θ,X):X→Y
其中,X为成像信号对应的图像,Y为磁共振弥散模值图像,θ为网络表征参数,f(θ,X)为以X为输入,θ为神经网络过程的输出。
5.根据权利要求3所述的一种基于深度学习和数据自洽的磁共振弥散加权成像方法,其特征在于:所述步骤2.1中的数据自洽过程包括如下步骤:
步骤2.1.1:第一层数据自洽层的输入为输入图像,其他数据自洽层的输入为相位约束层输出,即每次激发的多通道图像Ipcon,将其通过傅里叶变换获得K空间信号Spcon,计算如公式4所示:
Spcon=FIpcon
其中,F表示傅里叶变换;
步骤2.1.2:将K空间信号Spcon进行数据投影,若K空间位置上的数据已被采集,则填充为已采集数据Sacq,若未被采集则计算Spcon后构造第一K空间数据Sdp,计算如公式5所示:
其中,k表示K空间中的某个具体位置,Ω表示K空间中多个位置构成的集合,fdp(Spcon,Sacq,k)表示以Spcon、Sacq、k为输入的数据投影函数,输出为Sdp;
步骤2.1.3:将已采集数据和融入的数据自洽卷积核进行卷积运算获得第二K空间数据Ssc,计算如公式6所示:
Ssc=fsc(Sdp,G)=G*Sdp
其中,G表示数据自洽卷积核,fsc(Sdp,G)表示以K空间数据Sdp和数据自洽卷积核G为输入的函数,其输出为第二K空间数据Ssc;
步骤2.1.4:将第二K空间数据Ssc通过傅里叶逆变换转换为图像域的多通道图像Isc完成数据自洽,计算如公式7所示:
Isc=F-1Ssc
其中,F-1表示傅里叶逆变换。
6.根据权利要求5所述的一种基于深度学习和数据自洽的磁共振弥散加权成像方法,其特征在于:所述步骤2.3中的相位约束包括如下步骤:
其中,Icnn,i为每次激发对应图像,Ns为激发次数,TV(Icnn,i)表示对图像Icnn,i进行平滑处理后的图像;
步骤2.3.2:将不同激发图像Icnn,i减去对应的具有平滑性质的相位信息后,采用取平均的方式进行合并,计算如公式9所示:
步骤2.3.3:将每次激发对应的具有平滑性质的相位图重新附加到激发合成后的图像中,得到相位约束后的各次激发的图像Ipcon,i,计算如公式10所示:
7.根据权利要求3所述的一种基于深度学习和数据自洽的磁共振弥散加权成像方法,其特征在于:所述CNN网络包括多个CBR单元、成对汇聚层和反汇聚层以及残差链接,残差链接的一端连接CNN网络的输入,另一端与最后一个CBR单元连接构成最终的输出图像。
8.根据权利要求7所述的一种基于深度学习和数据自洽的磁共振弥散加权成像方法,其特征在于:每个所述CBR单元至少包含一个卷积层、一个规范化层和一个非线性激活层;所述CNN网络前半部分的CBR单元还加入汇聚层,后半部分加入反汇聚层。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810376066.9A CN108596994B (zh) | 2018-04-24 | 2018-04-24 | 一种基于深度学习和数据自洽的磁共振弥散加权成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810376066.9A CN108596994B (zh) | 2018-04-24 | 2018-04-24 | 一种基于深度学习和数据自洽的磁共振弥散加权成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108596994A CN108596994A (zh) | 2018-09-28 |
CN108596994B true CN108596994B (zh) | 2022-05-03 |
Family
ID=63609326
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810376066.9A Active CN108596994B (zh) | 2018-04-24 | 2018-04-24 | 一种基于深度学习和数据自洽的磁共振弥散加权成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108596994B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111856364B (zh) * | 2019-04-24 | 2023-03-28 | 深圳先进技术研究院 | 一种磁共振成像方法、装置、系统及存储介质 |
CN111856365B (zh) * | 2019-04-24 | 2023-03-14 | 深圳先进技术研究院 | 磁共振成像方法、装置、系统及存储介质 |
CN110346743B (zh) * | 2019-07-22 | 2021-09-14 | 上海东软医疗科技有限公司 | 一种磁共振弥散加权成像方法和装置 |
US11586982B2 (en) * | 2019-09-18 | 2023-02-21 | Samsung Electronics Co., Ltd. | Electronic and atomic structure computation utilizing machine learning |
US11537898B2 (en) | 2019-10-02 | 2022-12-27 | Samsung Electronics Co., Ltd. | Generative structure-property inverse computational co-design of materials |
CN112763958B (zh) * | 2020-12-10 | 2022-06-21 | 复旦大学 | 一种基于神经网络的多次激发平面回波磁共振成像方法 |
CN113854995B (zh) * | 2021-10-19 | 2023-11-24 | 复旦大学 | 一种基于单次激发的弥散加权成像扫描重建方法及系统 |
CN113936073B (zh) * | 2021-11-02 | 2024-05-14 | 哈尔滨理工大学 | 一种基于注意力机制的AttISTANet压缩感知磁共振重建方法 |
CN114114117A (zh) * | 2021-11-05 | 2022-03-01 | 中国科学院精密测量科学与技术创新研究院 | 一种基于深度学习的磁共振波谱相位校正方法 |
CN116152367B (zh) * | 2022-11-02 | 2023-10-20 | 佛山瑞加图医疗科技有限公司 | 一种基于弥散加权的自旋回波磁共振成像方法及装置 |
CN118229821A (zh) * | 2024-05-22 | 2024-06-21 | 山东奥新医疗科技有限公司 | 磁敏感伪影校正方法、装置、设备及介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103646410A (zh) * | 2013-11-27 | 2014-03-19 | 中国科学院深圳先进技术研究院 | 磁共振快速参数成像方法和系统 |
CN104020430A (zh) * | 2014-04-15 | 2014-09-03 | 清华大学 | 磁共振成像运动伪影的校正方法及系统 |
KR101659578B1 (ko) * | 2015-09-01 | 2016-09-23 | 삼성전자주식회사 | 자기 공명 영상 처리 방법 및 장치 |
CN106096616A (zh) * | 2016-06-08 | 2016-11-09 | 四川大学华西医院 | 一种基于深度学习的磁共振影像特征提取及分类方法 |
WO2017113205A1 (zh) * | 2015-12-30 | 2017-07-06 | 中国科学院深圳先进技术研究院 | 一种基于深度卷积神经网络的快速磁共振成像方法及装置 |
CN106970343A (zh) * | 2017-04-11 | 2017-07-21 | 深圳先进技术研究院 | 一种磁共振成像方法及装置 |
CN107369189A (zh) * | 2017-07-21 | 2017-11-21 | 成都信息工程大学 | 基于特征损失的医学图像超分辨率重建方法 |
CN107610194A (zh) * | 2017-08-14 | 2018-01-19 | 成都大学 | 基于多尺度融合cnn的磁共振图像超分辨率重建方法 |
CN107633486A (zh) * | 2017-08-14 | 2018-01-26 | 成都大学 | 基于三维全卷积神经网络的结构磁共振图像去噪方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9279873B2 (en) * | 2011-07-28 | 2016-03-08 | Siemens Corporation | Self consistent parallel imaging with temporal sensitivity estimation in real-time magnetic resonance imaging |
DE102016204145A1 (de) * | 2016-03-14 | 2017-09-14 | Siemens Healthcare Gmbh | Verfahren zur Aufnahme eines Magnetresonanzdatensatzes, Magnetresonanzeinrichtung, Computerprogramm und elektronisch lesbarer Datenträger |
-
2018
- 2018-04-24 CN CN201810376066.9A patent/CN108596994B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103646410A (zh) * | 2013-11-27 | 2014-03-19 | 中国科学院深圳先进技术研究院 | 磁共振快速参数成像方法和系统 |
CN104020430A (zh) * | 2014-04-15 | 2014-09-03 | 清华大学 | 磁共振成像运动伪影的校正方法及系统 |
KR101659578B1 (ko) * | 2015-09-01 | 2016-09-23 | 삼성전자주식회사 | 자기 공명 영상 처리 방법 및 장치 |
WO2017113205A1 (zh) * | 2015-12-30 | 2017-07-06 | 中国科学院深圳先进技术研究院 | 一种基于深度卷积神经网络的快速磁共振成像方法及装置 |
CN107182216A (zh) * | 2015-12-30 | 2017-09-19 | 中国科学院深圳先进技术研究院 | 一种基于深度卷积神经网络的快速磁共振成像方法及装置 |
CN106096616A (zh) * | 2016-06-08 | 2016-11-09 | 四川大学华西医院 | 一种基于深度学习的磁共振影像特征提取及分类方法 |
CN106970343A (zh) * | 2017-04-11 | 2017-07-21 | 深圳先进技术研究院 | 一种磁共振成像方法及装置 |
CN107369189A (zh) * | 2017-07-21 | 2017-11-21 | 成都信息工程大学 | 基于特征损失的医学图像超分辨率重建方法 |
CN107610194A (zh) * | 2017-08-14 | 2018-01-19 | 成都大学 | 基于多尺度融合cnn的磁共振图像超分辨率重建方法 |
CN107633486A (zh) * | 2017-08-14 | 2018-01-26 | 成都大学 | 基于三维全卷积神经网络的结构磁共振图像去噪方法 |
Non-Patent Citations (2)
Title |
---|
Synthesizing CT from Ultrashort Echo-Time MR Images via Convolutional Neural Networks;Roy et al.;《Simulation and Synthesis in Medical Imaging Second International Workshop》;20170910;第24-32页 * |
磁共振图像处理中部分傅里叶重建算法的比较;何汶静 等;《重庆医学》;20160720;第45卷(第20期);第2804-2806页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108596994A (zh) | 2018-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108596994B (zh) | 一种基于深度学习和数据自洽的磁共振弥散加权成像方法 | |
CN108594146B (zh) | 一种基于深度学习和凸集投影的磁共振弥散加权成像方法 | |
CN108535675B (zh) | 一种基于深度学习和数据自洽的磁共振多通道重建方法 | |
CN108335339B (zh) | 一种基于深度学习和凸集投影的磁共振重建方法 | |
CN110378980B (zh) | 一种基于深度学习的多通道磁共振图像重建方法 | |
CN108828481B (zh) | 一种基于深度学习和数据一致性的磁共振重建方法 | |
US10527699B1 (en) | Unsupervised deep learning for multi-channel MRI model estimation | |
US9229082B2 (en) | Method for diffusion magnetic resonance imaging | |
CN109615675B (zh) | 一种多通道磁共振成像的图像重建方法 | |
US11327137B2 (en) | One-dimensional partial Fourier parallel magnetic resonance imaging method based on deep convolutional network | |
CN108459289B (zh) | 一种基于数据一致性的多次激发弥散加权磁共振成像方法 | |
US20100244825A1 (en) | System And Method Of Parallel Imaging With Calibration To A Virtual Coil | |
CN106772167A (zh) | 核磁共振成像方法及装置 | |
CN113971706B (zh) | 一种快速磁共振智能成像方法 | |
US11948676B2 (en) | Qualitative and quantitative MRI using deep learning | |
CN110728624B (zh) | 一种高分辨率扩散加权图像重建方法 | |
CN113298902B (zh) | 卷褶视野磁共振图像的重建方法、计算机设备及存储介质 | |
CN110095742A (zh) | 一种基于神经网络的平面回波成像方法和装置 | |
CN110782452B (zh) | 一种t2定量图像成像方法及系统 | |
CN103091656B (zh) | 基于正则化约束多项式拟合磁共振线圈灵敏度的计算方法 | |
CN113129401B (zh) | 一种参数化磁共振成像的图像重建方法 | |
CN114255291A (zh) | 用于磁共振参数定量成像的重建方法、系统 | |
CN117094897B (zh) | 一种相衬光学相干弹性成像的超分辨相位梯度估计方法 | |
Dar et al. | Learning deep mri reconstruction models from scratch in low-data regimes | |
Gan et al. | SS-JIRCS: Self-supervised joint image reconstruction and coil sensitivity calibration in parallel MRI without ground truth |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |