CN113970716A - 梯度回波序列压脂方法及磁共振成像方法及设备 - Google Patents
梯度回波序列压脂方法及磁共振成像方法及设备 Download PDFInfo
- Publication number
- CN113970716A CN113970716A CN202111212023.5A CN202111212023A CN113970716A CN 113970716 A CN113970716 A CN 113970716A CN 202111212023 A CN202111212023 A CN 202111212023A CN 113970716 A CN113970716 A CN 113970716A
- Authority
- CN
- China
- Prior art keywords
- gradient
- fat
- pulse
- dephasing
- imaging
- 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.)
- Granted
Links
- 238000003825 pressing Methods 0.000 title claims abstract description 89
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000002595 magnetic resonance imaging Methods 0.000 title claims abstract description 28
- 238000003384 imaging method Methods 0.000 claims abstract description 76
- 230000005284 excitation Effects 0.000 claims abstract description 9
- 238000007443 liposuction Methods 0.000 claims description 28
- 238000004590 computer program Methods 0.000 claims description 15
- 230000007246 mechanism Effects 0.000 claims description 11
- 230000009467 reduction Effects 0.000 claims description 10
- 238000001514 detection method Methods 0.000 claims description 9
- 239000004519 grease Substances 0.000 claims description 9
- 230000003094 perturbing effect Effects 0.000 claims description 2
- 239000010410 layer Substances 0.000 abstract description 23
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 18
- 239000002356 single layer Substances 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 23
- 210000001519 tissue Anatomy 0.000 description 12
- 230000006870 function Effects 0.000 description 9
- 238000004088 simulation Methods 0.000 description 7
- 241000446313 Lamella Species 0.000 description 6
- 230000006399 behavior Effects 0.000 description 4
- 239000006185 dispersion Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000005415 magnetization Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 208000035126 Facies Diseases 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 210000000577 adipose tissue Anatomy 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 150000002632 lipids Chemical class 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 210000004761 scalp Anatomy 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 241000321520 Leptomitales Species 0.000 description 1
- 240000004371 Panax ginseng Species 0.000 description 1
- 235000002789 Panax ginseng Nutrition 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本申请提供了一种梯度回波序列压脂方法、磁共振成像方法、磁共振成像设备、压脂装置及计算机设备。在压脂脉冲激发被测组织的脂肪信号后,向不同的片层组施加不同的梯度散相组合,在压脂的同时使得不同片层组之间脂肪信号不串扰。每个梯度散相组合中的梯度散相的大小相同,保证了水和脂肪稳态信号按稳定的路径进行演化,稳态信号的回聚只发生在射频脉冲激发之前,消除了磁共振成像的伪影和斜纹。本申请提供的梯度回波序列压脂方法解决了单层、多层扫描场景下常规压脂和快速压脂带来的成像问题。
Description
技术领域
本申请涉及医学成像技术领域,特别涉及一种梯度回波序列压脂方法、磁共振成像方法、磁共振成像设备、压脂装置及计算机设备。
背景技术
常规临床协议大部分需要通过压脂来消除脂肪高亮信号对诊断的影响。梯度回波(GRE)序列为快速成像的常用序列。梯度回波序列采集的是自由感应衰减信号(FID)。为了获得T1对比,梯度回波序列必须同时具有梯度扰相和射频扰相。一种传统的常规压脂的实现方式是一个压脂脉冲之后跟一个成像脉冲,重复多次采满整个K空间。另一种传统的快速压脂的实现方式是一个压脂脉冲之后跟多个(两个及其以上)成像脉冲。然而,以上两种压脂方式形成的自旋回波或刺激回波都容易混入到成像回波中,多层扫描还会把不同层的脂肪信号混入到其他层当中。
为了解决上述问题,现有技术通常采用的技术方案为:在压脂脉冲后施加散相梯度,且散相梯度的大小按照设定的规律改变。然而,不断改变的梯度一方面会影响水质子的稳态信号;另一方面,变化的梯度会在演化路径上某个时刻因为梯度矩相加为0而形成回波,如果被信号采集窗接收了该回波信号,就会形成伪影。基于此,有必要对现有压脂成像方法进行改进。
发明内容
基于此,为了解决传统技术中单层、多层扫描场景下常规压脂和快速压脂成像时的水的稳态信号较差问题,本发明提出一种梯度回波序列压脂方法、磁共振成像设备、压脂装置及计算机设备。
一种梯度回波序列压脂方法,包括:
发出压脂脉冲,激发被测组织的脂肪信号;
向不同的片层或片层组施加不同的梯度散相组合,以压脂所述脂肪信号,其中,每个所述梯度散相组合中的梯度散相的大小相同。
在一个实施例中,所述梯度回波序列压脂方法还包括:设置用于提高被压脂的所述脂肪信号的稳定性的所述压脂脉冲的射频扰相。
在一个实施例中,所述向不同的片层组施加不同的梯度散相组合,以压脂所述脂肪信号,其中,每个所述梯度散相组合中的梯度散相的大小相同中,所述梯度散相组合包含梯度轴和极性,通过使所述梯度散相施加于不同的所述梯度轴和/或具有不同的所述极性,和/或使每个所述梯度散相组合施加所述梯度散相的所述梯度轴的数量不同,以使不同的所述片层组施加的所述梯度散相组合不同。
在一个实施例中,所述设置用于提高被压脂的所述脂肪信号的稳定性的所述压脂脉冲的射频扰相包括:
设置每个所述片层组的所述压脂脉冲额外添加的相位;
所述设置用于提高被压脂的所述脂肪信号的稳定性的所述压脂脉冲的射频扰相之后包括:发出成像脉冲,激发成像信号。
在一个实施例中,所述射频扰相满足以下三个条件:所述成像脉冲之间的射频扰相满足第一设定的射频扰相机制,所述成像脉冲和上一个紧挨的所述压脂脉冲之间满足第二设定的射频扰相机制,所述压脂脉冲和上一个紧挨的所述成像脉冲之间满足第三设定的射频扰相机制。
一种磁共振成像方法,所述磁共振成像方法包括:
获取第一扫描序列,所述第一扫描序列包括第一压脂脉冲、所述第一压脂脉冲之后施加的第一梯度散相脉冲和第一成像序列;
将所述第一扫描序列施加在检测对象,以获取所述检测对象的第一片层组的磁共振信号;
获取第二扫描序列,所述第二扫描序列包括第二压脂脉冲、所述第二压脂脉冲之后施加的第二梯度散相脉冲和第二成像序列,所述第二梯度散相脉冲与所述第一梯度散相脉冲的位置和/或极性不同;
将所述第二扫描序列施加在所述检测对象,以获取所述检测对象的第二片层组的磁共振信号;
重建所述第一片层组的磁共振信号、所述第二片层组的磁共振信号,以获取所述检测对象的磁共振图像。在一个实施例中,所述第一压脂脉冲的相位与所述第二压脂脉冲的相位非相干,和/或所述第一梯度散相脉冲与所述第二梯度散相脉冲的大小相同。
一种磁共振成像设备,所述磁共振成像设备包括:
发射线圈,用于发出压脂脉冲,激发被测组织的脂肪信号;
梯度线圈,用于向不同的片层或片层组施加不同的梯度散相组合,每个所述梯度散相组合中的梯度散相的大小相同。
一种压脂装置,所述装置包括:激发模块,用于发出压脂脉冲,激发被测组织的脂肪信号;梯度散相组合模块,用于向不同的片层或片层组施加不同的梯度散相组合,每个所述梯度散相组合中的梯度散相的大小相同。
一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序。所述处理器执行所述计算机程序时实现上述实施例中任一项所述方法的步骤。
本申请提供的梯度回波序列压脂方法,在压脂脉冲激发被测组织的脂肪信号后,向不同的片层或片层组施加不同的梯度散相组合,在压脂的同时使得不同片层组之间脂肪信号不串扰。每个梯度散相组合中的梯度散相的大小相同,保证了水和脂肪稳态信号按稳定的路径进行演化,稳态信号的回聚只发生在射频脉冲激发之前,消除了磁共振成像的伪影和斜纹。本申请提供的梯度回波序列压脂方法解决了单层、多层扫描场景下常规压脂和快速压脂带来的成像问题。
附图说明
图1为传统技术中梯度回波序列的梯度散相示意图。
图2为传统技术磁共振成像中头皮脂肪伪影示意图。
图3为传统技术磁共振成像中条状伪影示意图。
图4为本申请提供的一实施例中梯度回波序列压脂方法的流程示意图。
图5为本申请提供的一实施例中油膜和水膜与压脂脉冲和成像脉冲之间影响关系的示意图。
图6为本申请提供的一实施例中梯度回波序列的梯度散相示意图。
图7为本申请提供的一实施例中射频扰相方法的EPG仿真图。
图8为传统技术与本申请提供的梯度扰相实施例产生的其他层的脂肪伪影对比示意图。
图9为传统技术与本申请提供的梯度扰相实施例产生的条状伪影对比示意图。
图10为传统技术与本申请提供的射频扰相实施例产生的油膜与水膜对比示意图。
图11为本申请提供的另一实施例中射频扰相方法的EPG仿真图。
图12为本申请提供的另一实施例中射频扰相方法的EPG仿真图。
图13为本申请一实施所使用的扫描序列时序示意图。
图14为本申请另一实施所使用的扫描序列时序示意图。
具体实施方式
为使本申请的上述目的、特征和优点能够更加明显易懂,下面结合附图对本申请的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本申请。但是本申请能够以很多不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本申请内涵的情况下做类似改进,因此本申请不受下面公开的具体实施的限制。
本文中为部件所编序号本身,例如“第一”、“第二”等,仅用于分所描述的对象,不具有任何顺序或技术含义。而本申请所说“连接”、“联接”,如无特别说明,均包括直接和间接连接(联接)。在本申请的描述中,需要理解的是,术语“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”、“顺时针”、“逆时针”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请的限制。
在本申请中,除非另有明确的规定和限定,第一特征在第二特征“上”或“下”可以是第一和第二特征直接接触,或第一和第二特征通过中间媒介间接接触。而且,第一特征在第二特征“之上”、“上方”和“上面”可是第一特征在第二特征正上方或斜上方,或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”可以是第一特征在第二特征正下方或斜下方,或仅仅表示第一特征水平高度小于第二特征。
图1为传统技术中梯度回波序列的梯度散相示意图。申请人/发明人研究发现,传统压脂方法改变梯度散相的大小会影响水的稳态信号。改变梯度散相的大小会在信号的演化路径上某个时刻因为梯度零阶矩相加为0形成回波,即散相和回聚刚好相互抵消。如果模数转换(ADC)看到了回波信号,就会在图像中形成伪影。如果在成像回波中心之外看到其他回波信号,就能在图像上看到斜纹。如图2所示为临床场景所看到的其他层的头皮脂肪伪影。如图3所示,当前层采集到其他层的脂肪信号,可以在图像中见明显的条状伪影。
参见图4,本申请实施例提供一种梯度回波序列压脂方法。所述梯度回波序列压脂方法包括:
S10,发出压脂脉冲,激发被测组织的脂肪信号;
S20,向不同的片层或片层组施加不同的梯度散相组合,以压脂所述脂肪信号。其中,每个所述梯度散相组合中的梯度散相的大小相同。
所述S10中,所述压脂脉冲可以为一种射频脉冲。所述压脂脉冲可以由发射线圈发射。所述发射线圈可以是体发射线圈或局部发射线圈。所述被测组织是磁共振成像过程中患者需要检查病变的人体部位组织。
所述S20中,片层是磁共振成像过程中选定的层。所述片层组是多个片层的集合。所述片层组可以是一个,也可以是多个。一个所述片层组代表常规压脂,多个所述片层组代表快速压脂。压脂所述脂肪信号可以是抑制所述脂肪信号。所述梯度散相组合包含梯度轴、极性和大小相同的梯度散相。所述梯度散相的大小即为梯度强度对时间的积分。因此本实施例提供的所述梯度回波序列压脂方法可以用于解决单层、多层扫描场景下常规压脂和快速压脂带来的成像问题。
固定大小的所述梯度散相可以保证水和脂肪的信号按稳定的路径进行演化,信号回聚也都只发生在成像脉冲激发之前。不同所述片层组不同的所述梯度散相组合可以使得不同所述片层组之间脂肪信号不串扰。本实施例提供的所述梯度回波序列压脂方法通过在不同所述片层组设置独立不同的但大小相等的所述梯度散相组合,可以在消除脂肪带来的伪影同时不对水的稳态信号产生影响。
在一个实施例中,所述梯度回波序列压脂方法还包括:S30,设置用于提高被压脂的所述脂肪信号的稳定性的所述压脂脉冲的射频扰相。
所述S30中,射频扰相即为所述压脂脉冲额外添加的相位。图5为本申请提供的一实施例中油膜和水膜与压脂脉冲和成像脉冲之间影响关系的示意图。如图5所示,单独的油模和水模中,水和残留的脂肪都可以受到所述成像脉冲的作用,但无法受到所述压脂脉冲的作用。被压脂的脂肪可以同时感受到所述压脂脉冲和所述成像脉冲,图中的对号(√)代表能感受到对应的脉冲,叉号(X)代表不能感受到对应的脉冲。因为所述压脂脉冲和所述成像脉冲的翻转角不一样,所述梯度散相的大小也不一样,容易导致稳态信号被破坏,仍会造成被压脂的脂肪信号的扰动。本实施例设置所述压脂脉冲的射频扰相可以保证被压脂信号的稳定性,同时不影响水信号和残留脂肪的稳定性。
在一个实施例中,所述S20中,所述梯度散相组合包含梯度轴和极性。所述极性可以为梯度极性。通过使所述梯度散相施加于不同的所述梯度轴和/或具有不同的所述极性,和/或使每个所述梯度散相组合施加所述梯度散相的所述梯度轴的数量不同,以使不同的所述片层组施加的所述梯度散相组合不同。
所述梯度轴包括层面选择梯度轴(GSS)、相位编码梯度轴(GPE)和频率编码梯度轴(GRO)。所述极性可以为正梯度或负梯度。
图6为本申请提供的一实施例中梯度回波序列的梯度散相示意图。在一个实施例中,第1片层组可以只施加层面选择梯度轴方向、正梯度的梯度散相。第2片层组可以只施加相位编码梯度轴方向、正梯度的梯度散相。第3片层组可以只施加频率编码梯度轴方向、正梯度的梯度散相。第4片层组可以同时施加层面选择梯度轴方向和相位编码梯度轴方向、正梯度的梯度散相。第5片层组可以同时施加层面选择梯度轴方向和频率编码梯度轴方向、正梯度的梯度散相。第N-1片层组可以施加层面选择梯度轴方向、正梯度和相位编码梯度轴方向、负梯度的梯度散相,N为大于6的正整数。第N片层组可以施加层面选择梯度轴方向、正梯度和相位编码梯度轴方向、负梯度和频率编码梯度轴方向、正梯度的梯度散相,N为大于6的正整数。不同所述片层组使用不同形式的固定大小的所述梯度散相组合。所述梯度散相组合的组合形式包含不同轴、不同极性。在磁共振成像过程中,脂肪组织T2*不长,所述梯度散相组合的所有组合形式用完,脂肪组织T2*已经可以完成。因此,若所述梯度散相组合的所有组合形式用完,接下来就重复使用已排列好的组合形式,不会影响下一个重复使用所述梯度散相组合的所述片层组信号。
在一个实施例中,所述S30包括:设置每个所述片层组的所述压脂脉冲额外添加的相位。所述S30之后包括:发出成像脉冲,激发成像信号。所述成像脉冲可以为一种射频脉冲。所述成像脉冲可以由发射线圈发射。所述发射线圈可以是体发射线圈或局部发射线圈。一个所述压脂脉冲之后跟一个所述成像脉冲,重复多次采满整个K空间为常规压脂。一个所述压脂脉冲之后跟两个或两个以上所述成像脉冲为快速压脂。本实施例中所述的梯度回波序列压脂方法,不同所述片层组独立不同的固定大小所述梯度散相组合策略和特殊设计的所述射频扰相,可以完全消除脂肪带来的伪影,同时不对水信号产生影响。本实施例中所述的梯度回波序列压脂方法可以用于解决单层、多层扫描场景下常规压脂和快速压脂带来的成像问题。
在一个实施例中,所述射频扰相满足以下三个条件:所述成像脉冲之间的射频扰相满足第一设定的射频扰相机制,所述成像脉冲和上一个紧挨的所述压脂脉冲之间满足第二设定的射频扰相机制,所述压脂脉冲和上一个紧挨的所述成像脉冲之间满足第三设定的射频扰相机制。第一设定的射频扰相机制可以满足公式θ(n)-θ(n-2)=(n-1)/2×θA。第二设定的射频扰相机制可以满足公式θ(n)-θ(n-1)=n×θB+θC。第三设定的射频扰相机制可以满足公式θ(n+1)-θ(n)=(n+1)×θB+θC。
其中,θ(n)是所述片层组的所述压脂脉冲和成像脉冲额外添加的相位,θA、θB、θC为任意指定值。n的计数包含压脂脉冲和成像脉冲,如图5所示,从0开始计数,θ(0)为压脂脉冲,θ(1)为成像脉冲,θ(2)为压脂脉冲,θ(3)为成像脉冲,一直循环往后。上述公式是所述梯度回波序列的射频扰相要求,可以消除横向磁矢。
在一个实施例中,θA为117°或50°或123°。为了保证磁共振成像过程中射频扰相之后拥有正常的T1对比度,θA可以是常用的117°、50°、123°等。当θA为117°或50°或123°时,可以通过θ(n)-θ(n-2)=(n-1)/2×θA,θ(n)-θ(n-1)=n×θB+θC,θ(n+1)-θ(n)=(n+1)×θB+θC三个公式推导出:θ(n)=(n+1)×(n-1)/8×θA。所以,θB=1/4×θA,θC=-1/8×θA。图7为本申请提供的一实施例中射频扰相方法的EPG仿真图。参见图7,是在θA为117°或50°或123°,θ(n)=(n+1)×(n-1)/8×θA,θB=1/4×θA,θC=-1/8×θA的条件下,扩展相位图EPG(ExtendedPhase Graph)仿真的一串稳态信号幅值、相位、点扩散函数PSF(point spread function)分析。观察点扩散函数PSF图像,当出现峰值时,所述射频扰相效果最好。上述所述射频扰相方案可以完美地保证被压脂信号的稳定性,同时不影响水信号和残留脂肪的稳定性。
参见图8,左侧图像是传统压脂方法的图像,右侧图像是本申请实施例所述的梯度回波序列压脂方法的图像。可见所述梯度回波序列压脂方法可以很好解决其他层串入的脂肪伪影。参见图9,左侧图像是传统压脂方法的图像,右侧图像是本申请实施例所述的梯度回波序列压脂方法的图像。可见所述梯度回波序列压脂方法可以很好消除条状伪影。参见图10,左侧图像是传统压脂方法产生的油膜与水膜对比示意图,右侧图像是本申请实施例所述的梯度回波序列压脂方法的油膜与水膜对比示意图。可见所述梯度回波序列压脂方法不影响水信号和残留脂肪的稳定性。
在一个实施例中,所述压脂脉冲也可以采用随机的射频扰相角。所述射频扰相角即为所述压脂脉冲的激发角度。图11为本实施例中射频扰相方法的EPG仿真图。参见图11,是所述压脂脉冲采用随机的射频扰相角条件下,扩展相位图EPG(Extended Phase Graph)仿真的一串稳态信号幅值、相位、点扩散函数PSF(point spread function)分析。观察点扩散函数PSF图像,当出现峰值时,所述射频扰相效果最好。
在一个实施例中,所述压脂脉冲射频扰相满足θ(n)-θ(n-2)=n×θB+θC的任何θB和θC的组合,n为偶数,n为偶数代表的就是压脂脉冲。图12为本实施例中射频扰相方法的EPG仿真图。参见图12,是所述射频扰相满足θ(n)-θ(n-2)=n×θB+θC(当n偶数时)的任何θB和θC的组合条件下,扩展相位图EPG(Extended Phase Graph)仿真的一串稳态信号幅值、相位、点扩散函数PSF(point spread function)分析。观察点扩散函数PSF图像,当出现峰值时,所述射频扰相效果最好。上述两个实施例中的射频扰相方案在磁共振成像中未发现比较明显的伪影,也可以满足临床诊断对图像的要求。
应该理解的是,虽然图4的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图4中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行.
本申请实施例还提供一种磁共振成像方法,其包括:
获取第一扫描序列,所述第一扫描序列包括第一压脂脉冲(第一饱和脉冲)、所述第一压脂脉冲之后施加的第一梯度散相脉冲和第一成像序列;
将所述第一扫描序列施加在检测对象,以获取所述检测对象的第一片层组的磁共振信号;
获取第二扫描序列,所述第二扫描序列包括第二压脂脉冲(第二饱和脉冲)、所述第二压脂脉冲之后施加的第二梯度散相脉冲和第二成像序列,所述第二梯度散相脉冲与所述第一梯度散相脉冲的位置和/或极性不同;
将所述第二扫描序列施加在所述检测对象,以获取所述检测对象的第二片层组的磁共振信号;
重建所述第一片层组的磁共振信号、所述第二片层组的磁共振信号,以获取所述检测对象的磁共振图像。
进一步的,所述第一压脂脉冲的相位与所述第二压脂脉冲非相干。即第一饱和脉冲的相位与第二饱和脉冲的相位的大小不相同。
进一步的,所述第二梯度散相脉冲与所述第一梯度散相脉冲的大小相同。
进一步的,所述第一成像序列、所述第二成像序列为GRE(梯度回波)成像序列。
如图13所示,为本申请所使用的扫描序列时序示意图,该扫描序列可以是第一扫描序列或第二扫描序列,且每一压脂脉冲之后采集一片层的磁共振信号。其中,第一行为射频脉冲(RF);第二行为层面选择梯度(GSS);第三行为相位编码梯度(GPE);第四行为频率编码梯度(GRO);第五行为回波信号(Echo)。在此实施例中,首先施加压脂脉冲(β)1201;随后施加梯度散相脉冲1202,用以在压脂脉冲之后,成像射频脉冲施加之前,消除压脂脉冲1201之后残留的脂肪横向磁化矢量;成像序列可包括梯度散相脉冲1202后施加的成像射频脉冲(α)1203、与成像射频脉冲1203同时施加在GSS轴的层面选择梯度脉冲1204、施加在GPE轴的相位编码梯度1205、施加在GRO轴的离相位梯度1206与聚相位梯度1207,以采集一个片层的回波信号。进一步的,在回波信号采集完之后,该扫描序列还包括施加在GRO轴的扰相梯度1208。在成像射频脉冲1203的回波采集之后,下一次成像射频脉冲来临之前,通过施加扰相梯度1208,能够人为造成磁场的不均匀,加快质子失相位,以消除成像射频脉冲1203的回波采集之后残留的横向磁化矢量。
如图14所示,为本申请所使用的扫描序列时序示意图,该扫描序列同样可以是第一扫描序列或第二扫描序列,且每一压脂脉冲之后采集多片层的磁共振信号。其中,第一行为射频脉冲(RF);第二行为层面选择梯度(GSS);第三行为相位编码梯度(GPE);第四行为频率编码梯度(GRO);第五行为回波信号(Echo)。在此实施例中,首先施加压脂脉冲1301;随后施加梯度散相脉冲1311,用以在压脂脉冲之后,成像射频脉冲施加之前,消除压脂脉冲1301之后残留的脂肪横向磁化矢量;成像序列可包括梯度散相脉冲1311后施加的多个成像射频脉冲1302、1303和1304,且成像射频脉冲1302、1303和1304分别对应第一片层组(片层组1)的片层1-片层n。与1302、1303和1304分别对应同时施加在GSS轴的层面选择梯度脉冲1312、1313和1314。施加在GPE轴的相位编码梯度1321、1322和1323,施加在GRO轴的离相位梯度与聚相位梯度形成的频率编码梯度1331、1332和1333,以采集片层1-片层3的回波信号。进一步的,在每个回波信号采集完之后,该扫描序列还包括施加在GRO轴的扰相梯度1334、1335和1336。在成像射频脉冲的回波采集之后,下一次成像射频脉冲来临之前,通过施加扰相梯度,能够人为造成磁场的不均匀,加快质子失相位,以消除成像射频脉冲的回波采集之后残留的横向磁化矢量。本申请实施例中,通过在每一压脂脉冲之后采集多片层的磁共振信号,能够提高成像速度。本申请实施例还提供一种磁共振成像设备。所述磁共振成像设备包括发射线圈和梯度线圈。所述发射线圈用于发出压脂脉冲,激发被测组织的脂肪信号。所述梯度线圈用于向不同的片层或片层组施加不同的梯度散相组合。每个所述梯度散相组合中的梯度散相的大小相同。所述发射线圈可以是体发射线圈或局部发射线圈。所述发射线圈还可以用于发射成像脉冲。
所述磁共振成像设备还可以包括控制装置。通过所述控制装置控制所述发射线圈发出压脂脉冲,激发被测组织的脂肪信号。通过所述控制装置控制所述梯度线圈向不同的片层或片层组施加不同的梯度散相组合,压脂所述脂肪信号。每个所述梯度散相组合中的梯度散相的大小相同。通过所述控制装置还可以控制所述发射线圈发出成像脉冲,激发成像信号。
本申请实施例还提供一种压脂装置。所述压脂装置包括激发模块和梯度散相组合模块。所述激发模块用于发出压脂脉冲,激发被测组织的脂肪信号。所述梯度散相组合模块用于向不同的片层或片层组施加不同的梯度散相组合。每个所述梯度散相组合中的梯度散相的大小相同。
关于所述压脂装置的具体限定可以参见上文中对于所述梯度回波序列压脂方法的限定。所述压脂装置在所述压脂脉冲激发被测组织的脂肪信号后,向不同的片层或片层组施加不同的梯度散相组合,在压脂的同时使得不同片层组之间脂肪信号不串扰。每个梯度散相组合中的梯度散相的大小相同,保证了水和脂肪稳态信号按稳定的路径进行演化,稳态信号的回聚只发生在射频脉冲激发之前,消除了磁共振成像的伪影和斜纹。所述压脂装置实现的所述梯度回波序列压脂方法解决了单层、多层扫描场景下常规压脂和快速压脂带来的成像问题。所述压脂装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。在一个实施例中,提供了一种计算机设备。所述计算机设备包括存储器和处理器。所述存储器中存储有计算机程序。所述处理器执行计算机程序时实现以下步骤:
控制发射线圈发出压脂脉冲,激发被测组织的脂肪信号;
控制梯度线圈向不同的片层或片层组施加不同的梯度散相组合,以压脂所述脂肪信号。每个所述梯度散相组合中的梯度散相的大小相同。
所述计算机设备可以是终端。所述计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中,所述计算机设备的处理器用于提供计算和控制能力。所述计算机设备的存储器包括非易失性存储介质、内存储器。所述非易失性存储介质存储有操作系统和计算机程序。所述内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。所述计算机设备的网络接口用于与外部的终端通过网络连接通信。所述计算机程序被处理器执行时以实现一种梯度回波序列压脂方法。所述计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏,所述计算机设备的输入装置可以是显示屏上覆盖的触摸层,也可以是计算机设备外壳上设置的按键、轨迹球或触控板,还可以是外接的键盘、触控板或鼠标等。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:设置用于提高被压脂的所述脂肪信号的稳定性的所述压脂脉冲的射频扰相。具体地,可以设置每个所述片层组的所述压脂脉冲额外添加的相位。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:控制梯度线圈将所述梯度散相施加于不同的梯度轴和/或具有不同的极性,和/或使每个所述梯度散相组合施加所述梯度散相的所述梯度轴的数量不同,以使不同的所述片层组施加的所述梯度散相组合不同。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:控制发射线圈发出成像脉冲,激发成像信号。
在一个实施例中,提供了一种计算机可读存储介质。所述计算机可读存储介质存储有计算机程序。所述计算机程序被处理器执行时实现上述实施例中任一所述梯度回波序列方法的步骤。
具体地,所述处理器可以包括中央处理器(CPU),或者特定集成电路(ApplicationSpecific Integrated Circuit,简称为ASIC),或者可以被配置成实施本申请实施例的一个或多个集成电路。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,但并不能因此而理解为对本申请专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种梯度回波序列压脂方法,其特征在于,包括:
发出压脂脉冲,激发被测组织的脂肪信号;
向不同的片层或片层组施加不同的梯度散相组合,以压脂所述脂肪信号,其中,每个所述梯度散相组合中的梯度散相的大小相同。
2.如权利要求1所述的梯度回波序列压脂方法,其特征在于,还包括:
设置用于提高被压脂的所述脂肪信号的稳定性的所述压脂脉冲的射频扰相。
3.如权利要求2所述的梯度回波序列压脂方法,其特征在于,所述梯度散相组合包含梯度轴和极性,通过使所述梯度散相施加于不同的所述梯度轴和/或具有不同的所述极性,和/或使每个所述梯度散相组合施加所述梯度散相的所述梯度轴的数量不同,以使不同的所述片层组施加的所述梯度散相组合不同。
4.如权利要求2所述的梯度回波序列压脂方法,其特征在于,所述设置用于提高被压脂的所述脂肪信号的稳定性的所述压脂脉冲的射频扰相包括:
设置每个所述片层组的所述压脂脉冲额外添加的相位;
所述设置用于提高被压脂的所述脂肪信号的稳定性的所述压脂脉冲的射频扰相之后包括:
发出成像脉冲,激发成像信号。
5.如权利要求4所述的梯度回波序列压脂方法,其特征在于,所述射频扰相满足以下三个条件:所述成像脉冲之间的射频扰相满足第一设定的射频扰相机制,所述成像脉冲和上一个紧挨的所述压脂脉冲之间满足第二设定的射频扰相机制,所述压脂脉冲和上一个紧挨的所述成像脉冲之间满足第三设定的射频扰相机制。
6.一种磁共振成像方法,其特征在于,包括:
获取第一扫描序列,所述第一扫描序列包括第一压脂脉冲、所述第一压脂脉冲之后施加的第一梯度散相脉冲和第一成像序列;
将所述第一扫描序列施加在检测对象,以获取所述检测对象的第一片层组的磁共振信号;
获取第二扫描序列,所述第二扫描序列包括第二压脂脉冲、所述第二压脂脉冲之后施加的第二梯度散相脉冲和第二成像序列,所述第二梯度散相脉冲与所述第一梯度散相脉冲的位置和/或极性不同;
将所述第二扫描序列施加在所述检测对象,以获取所述检测对象的第二片层组的磁共振信号;
重建所述第一片层组的磁共振信号、所述第二片层组的磁共振信号,以获取所述检测对象的磁共振图像。
7.如权利要求6所述的磁共振成像方法,其特征在于,所述第一压脂脉冲的相位与所述第二压脂脉冲的相位非相干,和/或所述第一梯度散相脉冲与所述第二梯度散相脉冲的大小相同。
8.一种磁共振成像设备,其特征在于,包括:
发射线圈,用于发出压脂脉冲,激发被测组织的脂肪信号;
梯度线圈,用于向不同的片层或片层组施加不同的梯度散相组合,每个所述梯度散相组合中的梯度散相的大小相同。
9.一种压脂装置,其特征在于,所述装置包括:
激发模块,用于发出压脂脉冲,激发被测组织的脂肪信号;
梯度散相组合模块,用于向不同的片层或片层组施加不同的梯度散相组合,每个所述梯度散相组合中的梯度散相的大小相同。
10.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至7中任一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111212023.5A CN113970716B (zh) | 2021-10-18 | 2021-10-18 | 梯度回波序列压脂方法及磁共振成像方法及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111212023.5A CN113970716B (zh) | 2021-10-18 | 2021-10-18 | 梯度回波序列压脂方法及磁共振成像方法及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113970716A true CN113970716A (zh) | 2022-01-25 |
CN113970716B CN113970716B (zh) | 2024-04-23 |
Family
ID=79587594
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111212023.5A Active CN113970716B (zh) | 2021-10-18 | 2021-10-18 | 梯度回波序列压脂方法及磁共振成像方法及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113970716B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115616461A (zh) * | 2022-10-17 | 2023-01-17 | 浙江大学 | 一种基于可调回波时间编码的磁共振磁场测量方法及装置 |
CN116942134A (zh) * | 2023-09-19 | 2023-10-27 | 山东奥新医疗科技有限公司 | 一种磁共振成像方法、装置、设备及存储介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5757188A (en) * | 1995-01-19 | 1998-05-26 | Kabushiki Kaisha Toshiba | Magnetic resonance imaging apparatus with fat signal suppression |
CN102188245A (zh) * | 2010-03-17 | 2011-09-21 | 美国西门子医疗解决公司 | 用于在mr图像获取中进行脂肪抑制的系统 |
US20140303482A1 (en) * | 2013-04-03 | 2014-10-09 | Francesco Santini | Magnetic resonance imaging method for imaging components with short transverse relaxation times (t2) in a human or an animal heart |
CN105929350A (zh) * | 2016-05-05 | 2016-09-07 | 大连锐谱科技有限责任公司 | 一种单次激发水脂分离成像误差校正系统及方法 |
CN107997763A (zh) * | 2017-12-26 | 2018-05-08 | 中国科学院生物物理研究所 | 一种水脂分离磁共振成像方法 |
CN108957375A (zh) * | 2018-07-23 | 2018-12-07 | 深圳先进技术研究院 | 磁共振成像方法、装置、设备及存储介质 |
CN109620228A (zh) * | 2018-12-20 | 2019-04-16 | 上海联影医疗科技有限公司 | 磁共振成像中脂肪零点偏差校正方法及磁共振成像方法 |
CN111157932A (zh) * | 2020-01-02 | 2020-05-15 | 华东师范大学 | 一种快速自旋回波脉冲序列中射频脉冲的优化方法 |
WO2021196865A1 (zh) * | 2020-03-31 | 2021-10-07 | 浙江大学 | 一种振荡梯度准备的3d梯度自旋回波成像方法及设备 |
-
2021
- 2021-10-18 CN CN202111212023.5A patent/CN113970716B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5757188A (en) * | 1995-01-19 | 1998-05-26 | Kabushiki Kaisha Toshiba | Magnetic resonance imaging apparatus with fat signal suppression |
CN102188245A (zh) * | 2010-03-17 | 2011-09-21 | 美国西门子医疗解决公司 | 用于在mr图像获取中进行脂肪抑制的系统 |
US20140303482A1 (en) * | 2013-04-03 | 2014-10-09 | Francesco Santini | Magnetic resonance imaging method for imaging components with short transverse relaxation times (t2) in a human or an animal heart |
CN105929350A (zh) * | 2016-05-05 | 2016-09-07 | 大连锐谱科技有限责任公司 | 一种单次激发水脂分离成像误差校正系统及方法 |
CN107997763A (zh) * | 2017-12-26 | 2018-05-08 | 中国科学院生物物理研究所 | 一种水脂分离磁共振成像方法 |
CN108957375A (zh) * | 2018-07-23 | 2018-12-07 | 深圳先进技术研究院 | 磁共振成像方法、装置、设备及存储介质 |
CN109620228A (zh) * | 2018-12-20 | 2019-04-16 | 上海联影医疗科技有限公司 | 磁共振成像中脂肪零点偏差校正方法及磁共振成像方法 |
CN111157932A (zh) * | 2020-01-02 | 2020-05-15 | 华东师范大学 | 一种快速自旋回波脉冲序列中射频脉冲的优化方法 |
WO2021196865A1 (zh) * | 2020-03-31 | 2021-10-07 | 浙江大学 | 一种振荡梯度准备的3d梯度自旋回波成像方法及设备 |
Non-Patent Citations (1)
Title |
---|
丁庆国: "快速扰相梯度回波同相或反相位成像技术在含脂病变检查中的应用", 放射学实践, no. 5, pages 507 - 510 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115616461A (zh) * | 2022-10-17 | 2023-01-17 | 浙江大学 | 一种基于可调回波时间编码的磁共振磁场测量方法及装置 |
CN116942134A (zh) * | 2023-09-19 | 2023-10-27 | 山东奥新医疗科技有限公司 | 一种磁共振成像方法、装置、设备及存储介质 |
CN116942134B (zh) * | 2023-09-19 | 2023-11-28 | 山东奥新医疗科技有限公司 | 一种磁共振成像方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113970716B (zh) | 2024-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11880962B2 (en) | System and method for synthesizing magnetic resonance images | |
Jaubert et al. | Water–fat Dixon cardiac magnetic resonance fingerprinting | |
US9339239B2 (en) | Methods and devices for optimization of magnetic resonance imaging protocols | |
CN106997034B (zh) | 基于以高斯模型为实例整合重建的磁共振扩散成像方法 | |
CN111505553B (zh) | 磁共振成像系统和方法 | |
EP2870912A1 (en) | Magnetic resonance imaging device | |
US11181598B2 (en) | Multi-contrast MRI image reconstruction using machine learning | |
Dai et al. | Simultaneous multislice accelerated interleaved EPI DWI using generalized blipped‐CAIPI acquisition and 3D K‐space reconstruction | |
CN113970716A (zh) | 梯度回波序列压脂方法及磁共振成像方法及设备 | |
Haris et al. | Free‐breathing fetal cardiac MRI with doppler ultrasound gating, compressed sensing, and motion compensation | |
KR20150131124A (ko) | 신호 불균일 정정 및 성능 평가를 위한 방법 및 장치 | |
US8587305B2 (en) | System and method for fat suppression in chemical shift based water-fat separated MR imaging | |
Koolstra et al. | Cartesian MR fingerprinting in the eye at 7T using compressed sensing and matrix completion‐based reconstructions | |
US9971007B2 (en) | Method and apparatus for accelerated magnetic resonance imaging | |
JP6762284B2 (ja) | 磁気共鳴イメージング装置およびノイズ除去方法 | |
Ma et al. | Motion‐tolerant diffusion mapping based on single‐shot overlapping‐echo detachment (OLED) planar imaging | |
US10324150B2 (en) | Method and apparatus for generating a magnetic resonance image data set | |
CN115639510A (zh) | 磁共振成像方法、波谱成像方法、装置、设备和存储介质 | |
CN114076911B (zh) | 磁共振成像的方法及系统 | |
Dietrich et al. | Motion‐compensated fat‐water imaging for 3D cardiac MRI at ultra‐high fields | |
US20210166384A1 (en) | System, method and computer-accessible medium for tissue fingerprinting | |
CN111012349A (zh) | 磁共振成像中信号表征确定的系统与方法 | |
US20120194188A1 (en) | System and Method for Controlling Apparent Timing Dependencies for T2-Weighted MRI Imaging | |
CN116299108A (zh) | 基于3d-buda采集和联合低秩约束重建的mr成像方法及系统 | |
Ma et al. | Mitigating transmit‐B1 artifacts by predicting parallel transmission images with deep learning: A feasibility study using high‐resolution whole‐brain diffusion at 7 Tesla |
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 |