CN104068857B - 一种用于磁共振弹性成像的回波采集方法及装置 - Google Patents
一种用于磁共振弹性成像的回波采集方法及装置 Download PDFInfo
- Publication number
- CN104068857B CN104068857B CN201310106824.2A CN201310106824A CN104068857B CN 104068857 B CN104068857 B CN 104068857B CN 201310106824 A CN201310106824 A CN 201310106824A CN 104068857 B CN104068857 B CN 104068857B
- Authority
- CN
- China
- Prior art keywords
- echo
- data acquisition
- acquisition time
- magnetic resonance
- 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
- 238000002592 echocardiography Methods 0.000 title claims abstract description 89
- 238000002091 elastography Methods 0.000 title claims abstract description 33
- 238000010587 phase diagram Methods 0.000 claims description 32
- 238000010586 diagram Methods 0.000 claims description 24
- 230000000875 corresponding Effects 0.000 claims description 9
- 230000005284 excitation Effects 0.000 description 5
- 238000000034 method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 210000000056 organs Anatomy 0.000 description 2
- 230000000737 periodic Effects 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 210000001519 tissues Anatomy 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000012447 hatching Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000035812 respiration Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000000630 rising Effects 0.000 description 1
- 231100000817 safety factor Toxicity 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000001131 transforming Effects 0.000 description 1
Abstract
本发明公开的一种用于磁共振弹性成像的回波采集方法及装置,其在数据采集梯度轴上设置多个数据采集时间段。本发明提出的一种用于磁共振弹性成像的回波采集方法及装置,能够充分的利用现有技术没有利用的空闲时间,减少扫描时间,提高回波信号的信噪比。
Description
技术领域
本发明涉及医学设备技术领域,特别涉及磁共振弹性成像技术领域。
背景技术
磁共振弹性成像(MRE)的基本原理是利用磁共振技术检测体内的组织或器官在外力作用下产生的质点位移。通过在磁共振成像(MRI)设备中附加一套产生机械运动的激发装置,该激发装置于体表产生低频率剪切波,只要在梯度场的x、y或z轴上施加运动敏感梯度脉冲序列,剪切波导致的介质内的周期性移位就会使得接收到的信号产生周期性相位偏移,以此为基础得出组织或器官内部各点的弹性系数的分布图(即弹性图),由此作为医学诊断的依据。图1显示了现有的MRE回波数据采集的步骤,包括:先通过读出梯度轴采集回波数据得到K空间图;再通过K空间图得到磁共振的相位图;最后由相位图得到弹性图。
在现有的基于快速梯度回波序列(FGRE)的磁共振弹性成像(MRE)中,脉冲序列的重复时间(TR)一般设置成上述的激发装置产生的低频波周期的整数倍,比如:当激发装置产生的低频波的频率是60Hz时,TR为16.6毫秒。此外,再考虑到其他的一些安全因素,在实际配置MRE产品参数的时候,TR最少需要配置的更长一些,比如:25毫秒。这就意味着最小的TR值应当设置成上述低频波周期的两倍,即:33.3毫秒。这就导致每一个TR时间内有大约8毫秒的时间没有得到有效利用。图2显示了现有的MRE回波数据采集时的上述时长关系,其中,x轴是读出梯度方向,y轴是相位编码方向,z轴是选层梯度方向,正弦波202表示激发装置发出的低频波,从图1可以看出,现有技术在x轴上只有一个回波获取时间段201。
此外,为了保证足够的信噪比(SNR)水平,通常会将TR设置成上述低频波周期的三倍,即:50毫秒。这就导致每个TR内有更多没有得到利用的空闲时间。在这种情况下,如果没有采用阵列空间敏感性编码技术(ASSET),总的扫描时间为26秒;如果采用了阵列空间敏感性编码技术(ASSET),总的扫描时间为13秒。在不采用ASSET的情况下,对于病人而言,要在一个20多秒的扫描周期屏住呼吸是有困难的,而采用ASSET会导致相当的信噪比损失。上述扫描时间过长的问题,使得现有的MRE产品实现时倾向于采用平面回波成像(EPI)而非基于快速梯度回波序列(FGRE)。然而,EPI存在相当的图像变形且实现复杂。
发明内容
本发明所要解决的技术问题是提供一种新的磁共振弹性成像的回波采集方法,能够提升信噪比,降低MRE的扫描时间,以解决现有技术中MRE扫描时间过长,信噪比不高的问题。
为了解决这个问题,本发明提供了一种磁共振弹性成像的回波采集方法,包括如下步骤:在数据采集梯度轴上设置多个数据采集时间段。
为了解决这个问题,本发明还提供一种用于磁共振弹性成像的回波采集装置,包括用于在数据采集梯度轴上设置多个数据采集时间段的装置。
因此,与现有技术相比,本发明提供的一种磁共振弹性成像的回波采集方法和装置,其有益的技术效果是:
1、由于在一个TR内设置了多个用于获取回波的数据采集时间段,就能够充分的利用现有技术没有利用的空闲时间。
2、来自多个回波的相位图通过求平均的方式合成在一起,使得MRE扫描的信噪比提高,从而缩短了扫描时间。
3、由于缩短了扫描时间,使得可以不采用ASSET就可以完成MRE扫描,就避免了ASSET在进行校准扫描时由于病人呼吸运动带来的麻烦。
4、由于扫描时间的减少,使得基于FGRE的MRE和基于EPI的MRE扫描具有可比性。
总的来说,本发明提出的技术方案在没有额外的成本开销的情况下,提升了信噪比。
附图说明
图1是现有技术用于磁共振弹性成像的回波采集方法的步骤示意图;
图2是现有技术磁共振弹性成像的脉冲序列和低频波的关系的示意图;
图3是本发明提出的用于磁共振弹性成像的回波采集方法的步骤示意图;
图4是本发明提出的磁共振弹性成像的脉冲序列和低频波的关系的示意图;
图5是本发明提出的用于磁共振弹性成像的回波采集装置的示意图;
图6是对比现有技术和本发明采集到的数据的相位图及其噪声的示意图;
图7是在选定区域定量对比现有技术和本发明采集到的数据的相位差的示意图。
具体实施方式
下面结合附图对本发明作进一步说明。
为简单起见,以下描述中省略了本领域技术人员公知的某些技术特征。
图3示出了用于磁共振弹性成像回波数据采集方法可以包括的步骤。
其中的步骤301为在数据采集梯度轴上设置多个数据采集时间段。
如图4所示,假定x、y、z是磁共振弹性成像中的梯度轴,x轴是读出梯度方向,即:x轴是数据采集梯度轴,y轴是相位编码方向,z轴是选层梯度方向。低频剪切波用图中的正弦波202来表示,脉冲序列的重复时间TR被设置成正弦波202周期的两倍。在现有的数据采集时间段201之后,增加设置另一个回波数据采集时间段401。根据本发明的一个实施例,新增加的回波数据采集时间段401可以与原有的回波数据采集时间段201在时间上是相互邻接,即:在原有的回波数据采集时间段201结束以后,紧接着就会开始新增加的回波数据采集时间段401。两者之间没有时间间隔。新增加的回波数据采集时间段401的极性既可以和原有的数据采集时间段201的极性相同,也可以相反。图4给出的是极性相反的例子。新增加的回波数据采集时间段401的时长不小于该数据采集时间段401内需要采集的数据点个数除以接收带宽得到的值,比如,当接收带宽为31.25KHz时,如果一个数据采集时间段内需要采集96个点,则该时间段401的时长不小于96/31.25毫秒,即:3.072毫秒。在实际应用中,考虑到回波数据采集时间段401是一个梯形,而在我们只能在梯形的顶部所在的时间内进行数据采集,因此,实际的时长可以设置为大约4毫秒。
根据本发明的一个实施例,由于新增加了回波数据采集时间段401,y轴上的回聚梯度脉冲402和z轴上的反向梯度脉冲403的位置应当挪动到x轴上的最后一个数据采集时间段401结束以后。
需要说明的是,本实施例在此仅给出了增加一个回波数据采集时间段401的示例,实际上还可以增加更多的回波数据采集时间段。这些回波数据采集时间段既可以相互邻接,也可以有时间间隔;它们的极性既可以设置成全部相同,比如:全部为正或者全部为负,也可以设置成相邻的两个回波数据采集时间段的极性相反(一正一负)。每个数据采集时间段的时长都不小于该数据采集时间段内需要采集的数据点个数除以接收带宽得到的值。
其中的步骤302为获取每个数据采集时间段内采集到的数据的K空间图。根据本发明的一个实施例,可以分别获取两个数据采集时间段201和401内采集到的数据的K空间图。
其中的步骤303为计算每个K空间图对应的单回波相位图。根据本发明的一个实施例,分别将两个回波数据采集时间段内的K空间图做傅里叶变换,然后进行相位取模、相位相减和相位解缠绕(phase unwrapping),就可以得到每个回波数据采集时间段内的单回波相位图。相位相减的目的是为了消除B0场的不一致性、本地磁感性(localsusceptibility)和涡流(eddy currents)导致的相位积累误差。根据本发明的一个实施例,相位相减可以对于每一个相位偏移(phase offset),先做两次数据采集,这两次数据采集的运动编码梯度的极性相反;然后将两次采集得到的相位相减。
其中的步骤304为将单回波相位图进行加权叠加得到多回波相位图。根据本发明的一个实施例,可以将单回波相位图进行加权叠加,得到多回波相位图。
由于通过上述处理以后的单回波相位图只包含有运动带来的相位变化,和回波时间(TE)无关,这些相位图就可以直接加权叠加到一起来得到多回波对应的相位图。
这里所述的加权叠加,既可以是将所有单回波相位图简单的求算术平均,也可以对每个单回波相位图赋予不同的权重,进行加权叠加。比如:对某些质量较好的单回波相位图赋予较大的权重,对某些质量较差的单回波相位图赋予较小的权重。
其中的步骤305为用多回波相位图计算得到组织体的弹性分布图。根据本发明的一个实施例,可以利用加权叠加得到的多回波相位图就可以计算得到组织体的弹性分布图,即:刚度图。
图5给出了本发明提出的用于磁共振弹性成像的回波采集装置501的示意图,该装置可以包括:用于在数据采集梯度轴上设置多个数据采集时间段的装置502;用于获取每个数据采集时间段内采集到的数据的K空间图的装置503;用于计算每个K空间图对应的单回波相位图的装置504;用于将单回波相位图进行加权叠加得到多回波相位图的装置505;用于通过多回波相位图计算弹性分布图的装置506。
根据本发明的一个实施例,多个数据采集时间段既可以是相互邻接的,也可以有时间间隔。
根据本发明的一个实施例,每个数据采集时间段的时长不小于该数据采集时间段内需要采集的数据点个数除以接收带宽得到的值。
根据本发明的一个实施例,多个数据采集时间段的极性既可以完全相同,即:全部为正或者全部为负,也可以是相邻两个数据采集时间段的极性相反。
根据本发明的一个实施例,用于计算每个K空间图对应的单回波相位图的装置504可以进一步包括用于通过进行两次数据采集来计算单回波相位图的装置,其中,这两次数据采集的运动编码梯度的极性相反。
根据本发明的一个实施例,用于计算每个K空间图对应的单回波相位图的装置504可以进一步包括用于通过将两次采集得到的数据的相位相减来计算所述的单回波相位图的装置。
下面给出本发明技术方案在真实的MRE设备上测试的结果。
首先,这台基于FGRE的MRE回波序列配置成图4所示,使得其可以采集多回波数据。然后采用了如下扫描参数进行两个回波的数据采集:FOV=30cm,矩阵大小matrix=256*64,RBW=31.25kHz,TR/TE=33.3/22.6ms,相位偏移有4个,低频波的频率为60Hz,扫描时间是17秒。由此得到的其中一个回波的相位图及其噪声如图6b所示,相位差如图7b所示;由此得到的两个回波求平均叠加以后的多回波相位图及其噪声如图6c所示,相位差如图7c所示。作为对比,还采用现有技术进行了回波数据获取,扫描参数与上面的相同,只是TR改为50毫秒,总的扫描时间变成了26秒,由此得到的相位图及其噪声如图6a所示,相位差如图7a所示。上述三种不同类型的扫描都重复的进行了两次采集,这两次数据采集的运动编码梯度的极性相反。这两次采集之间的相位差就被认作是噪声,这个噪声通过平均相位进行归一化并被用作一个基准尺度。
可以看出,图6b相对于图6a表现出了更高的噪声波动,因为相对于图6a而言,其TR减少了;图6c中的多回波平均相位的噪声水平相对于图6b有所改善。相位图上面还标识了剖面线来说明上述观点。
分别在图7a,图7b和图7c上选择两个区域A和B,定量计算这两个选定区域归一化以后的噪声,结果如下:
图7a中的区域A的噪声为0.15,图7a中的区域B的噪声为0.13;
图7b中的区域A的噪声为0.20,图7b中的区域B的噪声为0.25;
图7c中的区域A的噪声为0.14,图7c中的区域B的噪声为0.11;
上述定量分析的结果与图6的结果一致,即:减少TR会导致噪声上升;而采用本发明提出的多回波获取不仅能够有效地补偿这个信噪比损失,还能在扫描时间减少了30%的情况下,比现有技术的信噪比有所上升。
应当注意,上面所描述的实施例仅仅是示例性而非限制性的,并且本领域技术人员在不偏离所附权利要求书的范围的情况下可以设计出很多备选实施例。所使用的动词“包括”并不排除权利要求书或说明书中所记载的元件和步骤之外的元件和步骤。在元件之前使用的词语“一个”并不排除存在多个这种元件。
Claims (14)
1.一种用于磁共振弹性成像的回波采集方法,其特征是包括如下步骤:
在数据采集梯度轴的一个重复时间内设置多个数据采集时间段;
获取每个数据采集时间段内采集到的数据的K空间图;
基于相位相减计算每个K空间图对应的单回波相位图;
将单回波相位图进行加权叠加得到多回波相位图;以及
通过多回波相位图计算弹性分布图。
2.根据权利要求1所述的一种用于磁共振弹性成像的回波采集方法,其特征是:所述的多个数据采集时间段是相互邻接的。
3.根据权利要求1所述的一种用于磁共振弹性成像的回波采集方法,其特征是:所述多个数据采集时间段中的每个数据采集时间段的时长不小于该数据采集时间段内需要采集的数据点个数除以接收带宽得到的值。
4.根据权利要求1所述的一种用于磁共振弹性成像的回波采集方法,其特征是:所述的多个数据采集时间段的极性相同。
5.根据权利要求1所述的一种用于磁共振弹性成像的回波采集方法,其特征是:所述的多个数据采集时间段中,相邻两个数据采集时间段的极性相反。
6.根据权利要求1所述的一种用于磁共振弹性成像的回波采集方法,其特征是:通过进行两次数据采集来计算所述单回波相位图,这两次数据采集的运动编码梯度的极性相反。
7.根据权利要求6所述的一种用于磁共振弹性成像的回波采集方法,其特征是:通过将两次采集得到的数据的相位相减来计算所述的单回波相位图。
8.一种用于磁共振弹性成像的回波采集装置,其特征是包括:
用于在数据采集梯度轴的一个重复时间内设置多个数据采集时间段的装置;
用于获取每个数据采集时间段内采集到的数据的K空间图的装置;
用于基于相位相减计算每个K空间图对应的单回波相位图的装置;
用于将单回波相位图进行加权叠加得到多回波相位图的装置;以及
用于通过多回波相位图计算弹性分布图的装置。
9.根据权利要求8所述的一种用于磁共振弹性成像的回波采集装置,其特征是:所述的多个数据采集时间段是相互邻接的。
10.根据权利要求8所述的一种用于磁共振弹性成像的回波采集装置,其特征是:所述多个数据采集时间段中的每个数据采集时间段的时长不小于该数据采集时间段内需要采集的数据点个数除以接收带宽得到的值。
11.根据权利要求8所述的一种用于磁共振弹性成像的回波采集装置,其特征是:所述的多个数据采集时间段的极性相同。
12.根据权利要求8所述的一种用于磁共振弹性成像的回波采集装置,其特征是:所述的多个数据采集时间段中,相邻两个数据采集时间段的极性相反。
13.根据权利要求8所述的一种用于磁共振弹性成像的回波采集装置,其特征是:所述用于计算每个K空间图对应的单回波相位图的装置进一步包括用于通过进行两次数据采集来计算单回波相位图的装置,其中,这两次数据采集的运动编码梯度的极性相反。
14.根据权利要求13所述的一种用于磁共振弹性成像的回波采集装置,其特征是:所述的用于计算每个K空间图对应的单回波相位图的装置进一步包括用于通过将两次采集得到的数据的相位相减来计算所述的单回波相位图的装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310106824.2A CN104068857B (zh) | 2013-03-29 | 一种用于磁共振弹性成像的回波采集方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310106824.2A CN104068857B (zh) | 2013-03-29 | 一种用于磁共振弹性成像的回波采集方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104068857A CN104068857A (zh) | 2014-10-01 |
CN104068857B true CN104068857B (zh) | 2019-07-16 |
Family
ID=
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5899858A (en) * | 1997-04-10 | 1999-05-04 | Mayo Foundation For Medical Education And Research | MR imaging with enhanced sensitivity of specific spin motion |
CN101708123A (zh) * | 2009-10-28 | 2010-05-19 | 上海理工大学 | 肝纤维化分级研究的磁共振弹性成像检测系统及其方法 |
CN102782518A (zh) * | 2009-12-21 | 2012-11-14 | 皇家飞利浦电子股份有限公司 | 磁共振弹性成像 |
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5899858A (en) * | 1997-04-10 | 1999-05-04 | Mayo Foundation For Medical Education And Research | MR imaging with enhanced sensitivity of specific spin motion |
CN101708123A (zh) * | 2009-10-28 | 2010-05-19 | 上海理工大学 | 肝纤维化分级研究的磁共振弹性成像检测系统及其方法 |
CN102782518A (zh) * | 2009-12-21 | 2012-11-14 | 皇家飞利浦电子股份有限公司 | 磁共振弹性成像 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103608693B (zh) | 具有用于脂肪抑制的化学位移编码的对比度增强磁共振血管造影 | |
JP6998218B2 (ja) | 動き検出を用いるmr撮像 | |
US7511494B2 (en) | Method for measuring the microarchitecture of complex tissue with MRI | |
JP5193447B2 (ja) | セグメント分割マルチショット放射方向ファンビーム・エンコード順序でmrデータを収集するための装置 | |
CN105074491B (zh) | 具有使用压缩感知的图像重建的动态mri | |
CN104204838B (zh) | 具有Dixon类型水/脂肪分离和关于主磁场的不均匀性的先验知识的MRI | |
JP5719968B2 (ja) | Mrデータを収集する方法及び装置 | |
CN105934683B (zh) | 具有对主磁场变化的估计的dixon型水/脂肪分离的mri | |
US9746539B2 (en) | MR imaging with suppresion of flow artifacts | |
CN102413762B (zh) | 具有运动补偿的介入式mr成像 | |
CN110031786B (zh) | 磁共振图像重建方法、磁共振成像方法、设备及介质 | |
CN105283125B (zh) | 具有增强的磁化率对比度的mr成像 | |
EP2623028A1 (en) | Magnetic resonance imaging device and magnetic resonance imaging method | |
CN103099618A (zh) | 用于生成磁共振血管造影图像的方法和相应的磁共振设备 | |
CN103957785A (zh) | Mr电学性质断层摄影 | |
US20120281901A1 (en) | Magnetic resonance imaging apparatus and blood vessel image capturing method | |
CN106796274B (zh) | 具有伪迹抑制的propeller-mr成像 | |
CN107076818B (zh) | 零回波时间mr成像 | |
CN102772207A (zh) | 磁共振成像方法和装置 | |
JP5735916B2 (ja) | 磁気共鳴イメージング装置及び同期計測方法 | |
JP2018033691A (ja) | 磁気共鳴測定装置および画像処理方法 | |
Ibrahim et al. | Real‐time MR imaging of myocardial regional function using strain‐encoding (SENC) with tissue through‐plane motion tracking | |
CN109716155A (zh) | 具有迪克逊型水/脂肪分离的mr成像 | |
JP6433425B2 (ja) | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 | |
CN104068857B (zh) | 一种用于磁共振弹性成像的回波采集方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |