CN114137462A - 一种用于低场磁共振的多对比度成像方法和设备 - Google Patents
一种用于低场磁共振的多对比度成像方法和设备 Download PDFInfo
- Publication number
- CN114137462A CN114137462A CN202111376303.XA CN202111376303A CN114137462A CN 114137462 A CN114137462 A CN 114137462A CN 202111376303 A CN202111376303 A CN 202111376303A CN 114137462 A CN114137462 A CN 114137462A
- Authority
- CN
- China
- Prior art keywords
- echo
- radio frequency
- pulse sequence
- frequency pulse
- echo signals
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 38
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims abstract description 38
- 238000002592 echocardiography Methods 0.000 claims abstract description 15
- 238000003860 storage Methods 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 6
- 230000035945 sensitivity Effects 0.000 claims description 6
- 230000010363 phase shift Effects 0.000 claims description 5
- 230000009977 dual effect Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 238000000264 spin echo pulse sequence Methods 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 12
- 230000006870 function Effects 0.000 description 9
- 238000002595 magnetic resonance imaging Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 7
- 230000008901 benefit Effects 0.000 description 6
- 230000005284 excitation Effects 0.000 description 4
- 210000004885 white matter Anatomy 0.000 description 4
- 210000004556 brain Anatomy 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 210000004884 grey matter Anatomy 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004065 semiconductor Substances 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]
- G01R33/48—NMR imaging systems
- G01R33/50—NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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
技术领域
本发明涉及磁共振成像技术领域，更具体地，涉及一种用于低场磁共振的多对比度成像方法。
背景技术
在传统的磁共振成像(MRI)系统中，一次扫描仅仅可以重建出一种对比度加权图像。如果是定量图像则需要设置不同的重复时间(TR)和回波时间(TE)进行多次扫描，拟合出定量图像。这些不同的定性加权图和定量的图像通常需消耗大量的扫描时间，并且有可能得到错误的配准图像。此外，具有不同翻转角(FA)的GRE序列的T1加权图(T1W)图像产生有效的灰质(GM)和白质(WM)对比度，并且有来自射频发射(B1)和接收线圈信号响应变化的固有图像不均匀性。GM/WM对比度不足和RF不均匀性增加了信号的T1W自动脑结构分割的复杂性。此外，与没有任何磁化准备的常规GRE(梯度重聚回波成像)序列相比，T1W的生成过程常用梯度回波快速采集(T1MPRAGE)，但是T1MPRAGE中使用较小的FA，比使用较大的FA的GRE序列具有更多的B1场的不均匀性，最后长回波时间TE的QSM(定量磁敏感成像)可能会出现相位的混叠。由于长回波时间的数据对敏感性的微小变化更为敏感，因此需要更短的TE来更好的量化敏感性。
由于易于收集数据，可变翻转角(VFA)方法在评估T1定量图(T1 map)和PD(质子密度)方面得到了普及，该技术基于采集具有不同翻转角的GRE序列编码后的信号，然后对T1和PD进行适当量化。然而VFA方法需要准确的知道射频发射场或者翻转角。对于质子密度和T1 map的准确获取，前者取决于射频接收场，后者取决于射频发射场，可以通过使用三种不同的FA，在没有任何约束的情况下逐像素计算射频接收场，尽管这种方法需要较大的翻转角，但是还是可以计算出PD和T1，否则会成为不适定问题。最近的研究表明，偏置场矫正算法可用于射频发射场映射，但这在临床中难以实现。在主磁场强度从0.2T～7T上的MRI设备中，T1(纵向弛豫时间)的驰豫时间随着磁场的增加而增加，而T2/T2*驰豫时间随着磁场的增加而减少，相对于高场，这对于同时实现更多参数的同时成像提供了机会。然而，现有技术中还没有比较完善的用于低场磁共振的多对比度成像方法。
发明内容
本发明的目的是克服上述现有技术的缺陷，提供一种用于低场磁共振的多对比度成像方法和设备。
根据本发明的第一方面，提供一种用于低场磁共振的多对比度成像方法。该方法包括以下步骤：
依次采用设定的第一射频脉冲序列和第二射频脉冲序列应用于目标成像区域，其中第一射频脉冲序列采用第一翻转角，第二射频脉冲序列采用第二翻转角；
针对第一射频脉冲序列采集第一组回波信号，并针对第二射频脉冲序列采集第二组回波信号，其中第一组回波信号和第二组回波信号中均包含多个回波；
对所获得的回波信号进行调制和欠采样，进而重建得到目标成像区域的多种对比度的磁共振图像。
根据本发明的第二方面，提供一种用于低场磁共振的多对比度成像设备。该设备包括：
扫描单元：用于依次采用设定的第一射频脉冲序列和第二射频脉冲序列应用于目标成像区域，其中第一射频脉冲序列采用第一翻转角，第二射频脉冲序列采用第二翻转角；
信号采集单元：用于针对第一射频脉冲序列采集第一组回波信号，并针对第二射频脉冲序列采集第二组回波信号，其中第一组回波信号和第二组回波信号中均包含多个回波；
图像重建单元：用于对所获得的回波信号进行调制和欠采样，进而重建得到目标成像区域的多种对比度的磁共振图像。
与现有技术相比，本发明的优点在于，相对于高场的多对比图像生成方法，本发明能够有效利用低场磁共振中更长的T2/T2*(横向弛豫时间)，在实现可变翻转角对T1mapping(测量组织的T1值)的同时，有效利用两个翻转角之间的空隙进行T2图像生成。此外，还提出了基于波浪梯度场的编码方法对信号进行调制和欠采样，以利用基于多测量压缩感知的图像重建算法对高倍欠采样的数据进行高质量的重建，从而降低了系统的条件数并增加了重建图像信噪比。本发明提供的成像序列可以单次扫描实现多种对比度图像的生成和重建。
通过以下参照附图对本发明的示例性实施例的详细描述，本发明的其它特征及其优点将会变得清楚。
附图说明
被结合在说明书中并构成说明书的一部分的附图示出了本发明的实施例，并且连同其说明一起用于解释本发明的原理。
图1是根据本发明一个实施例的用于低场磁共振的多对比度成像方法的流程图；
图2是根据本发明一个实施例的快速多对比度磁共振成像序列的示意图。
具体实施方式
现在将参照附图来详细描述本发明的各种示例性实施例。应注意到：除非另外具体说明，否则在这些实施例中阐述的部件和步骤的相对布置、数字表达式和数值不限制本发明的范围。
以下对至少一个示例性实施例的描述实际上仅仅是说明性的，决不作为对本发明及其应用或使用的任何限制。
对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论，但在适当情况下，所述技术、方法和设备应当被视为说明书的一部分。
在这里示出和讨论的所有例子中，任何具体值应被解释为仅仅是示例性的，而不是作为限制。因此，示例性实施例的其它例子可以具有不同的值。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步讨论。
相对于高场以及超高场，低场磁共振下的组织T2/T2*值更长，更长的T2/T2*值使得组织信号衰减变慢，为信号的采集提供了更长的时间窗。因此，为了更加充分利用低场磁共振下信号衰减速度慢这一特点，本发明结合多回波GRE序列和FSE(快速自旋回波脉冲序列)序列的原理，设计了适用于低场磁共振的快速多对比度MRI序列和相应的多对比度成像方法。
参见图1所示，所提供的用于低场磁共振的多对比度成像方法包括以下步骤。
步骤S110，对于目标成像区域进行两次射频脉冲扫描，并利用低场磁共振中横向弛豫时间增加的优势，设计快速多对比度磁共振成像序列。
以第一次扫描的射频脉冲的翻转角设置为α1，第二次扫描的射频脉冲翻转角设置为α2为例，所设计的多对比度磁共振成像序列如图2所示。在翻转角为α1的非对称射频脉冲对目标成像区域组织进行激发后，首先利用双回波GRE对，读出2组梯度回波信号(标记为回波#1和#2)。随后，采用回波链长为2的FSE对2组自旋回波信号进行读出(回波#3和#4)。然后，使用翻转角为α2的非对称射频脉冲再次对组织进行激发，并使用双回波GRE采集2组梯度回波信号(回波#5和#6)。为了缩短回波时间TE以获得高信号的回波，图2实施例使用了非对称射频脉冲对组织进行激发。
在采集到上述多组回波之后，对于磁共振图像的生成，在一个实施例中，回波1和回波2通过采用非线性拟合的方式用于T2*map的生成；回波3和回波4利用非线性拟合的方式生成T2map；回波5和回波6用于生成T1map。需要说明的是，在α1和α2射频激发中间在高场中是无法使用的，然而对于低场磁共振，T2相对于高场磁共振会相对大的很多，因此在低场磁共振中存在固有的T2变长的优势，可以进一步生成T2map图像。并且，在回波1和回波3通过傅里叶变换的方式可以生成T1加权图像(T1W)和T2加权图像(T2W)，回波1和回波2可以生成T2*加权图像(T2*W)，此时一共生成了6种图像。应理解的是，α1和α2相位可以相同，但是不仅局限于相同相位。180°硬脉冲的相位相差π，并且第一个硬脉冲和α1具有π相位差。
具体地，通过求解Bloch方程可以得到：
其中，S表示快速自旋回波(FSE)总体信号强度，ρ0为质子密度，TR表示重复时间，T1表示纵向驰豫时间，T2表示横向弛豫时间，TE表示回波时间。设置长TR和长TE可以得到T2W图像，在一个实施例中，结合FSE的基本原理，可以调整回波3和回波4的时间，进行T2W图像的生成。当仅留有T2W的图像之后，公式(1)变成：
S′＝ρ0*e-TE/T2 (2)
通过设置不同的TE值，得到不同的信号响度(强度)，然后通过线性拟合的方式得到T2map图像。在本发明中，回波3和回波4的不同时间代表不同的TE值，因此可以直接使用回波3和回波4得到的T2W图像进行拟合T2map。
在梯度回波中，通过调整长TR和长TE值，可以得到类似于公式(2)的显示的信号模型，表示为：
S″＝ρ0*e-TE/T2 (3)
对于不同的TE值，可以得到不同的T2*加权图，通过不同的T2*加权图像，进行线性拟合后可以得到对应的T2*map图像。
进一步地，通过两个翻转角α1和α2，可以生成质子密度加权图(PD)和真实质子密度加权图(truePD，tPD)，表示为：
其中ρ0表示PD空间分布图像，ρi表示PDW图像，PD map可以从回波1和回波5进行线性拟合方式得到，T2*表示有效横向驰豫时间，i表示不同的梯度回波。
在射频场不均匀的情况下，射频衰减梯度回波数据采集的信号可以作为可变翻转角θ的函数，例如由下式得到：
其中，TR表示重复时间，bias表示结合线圈灵敏度的偏置，PD表示质子密度，TE表示回波时间，k表示发射场的缩放因子，为了使T1map的计算更加直接，将公式(5)简化为：
其中，E1＝e-TR/T1，
通过收集不同翻转角的激发数据，可以将转换后的数据拟合到线性，其中斜率为E1，截距为PDeff*(1-E1)，为了求解出T1map和PD map，至少需要两个翻转角。
对于理想的稳态RF损毁梯度GRE数据采集中，翻转角θ的恩斯特(Ernst)方程为：
其中，ρ0为质子密度，E1＝e-TR/T1，
TR为重复时间，TE表示回波时间，T1和T2*分别表示纵向和横向驰豫时间。结合公式(7)可知，由于射频发射场和接收场的变化，给定组织可能在不同的位置显示不同的信号强度，从公式(7)中可以得到小于白质的Ernst角的FA比大于Ernst的FA承受更多的射频发射场的变化。在一个实施例中，通过较大翻转角的回波信号减去较小翻转角的信号可以得到增强T1加权图(aT1W)，表示为：
aT1＝s(θ2，TEn)-λs(θ1，TEn) (8)
其中，θ2表示本发明中的α2，例如θ2＝24°提供aT1W图像，θ1＝6°提供PDW图像，θ1表示α1，aT1表示测量到的增强T1加权信号。应理解的是，两个RF的翻转角度不局限于6°和24°，
k为提取的射频发射场的变化，TEn为第n个回波，由于两次扫描都是双回波扫描，最终的aT1由第一次和第二次(此处的第一次和第二次指的是α2和α1两次激发的得到的回波)回波根据等式(8)计算得出的两个T1W图像之间的平均值给出。在GRE序列中，由Bloch方程和Ernst方程可知，大于Ernst角度FA给出T1W图像。
对于不同翻转角的射频脉冲(α1和α2)，后面紧跟的两个梯度回波，一共四个梯度回波(回波1、回波2、回波5和回波6)可以选择作为计算磁化率敏感图(QSM)加权图。例如，使用3D相位展开算法进行求解QSM加权图像。在本发明实施例中，使用回波1和回波2，回波5和回波6，TE时间相差如1.25ms(不局限于1.25ms)能够有效的防止相位混叠，在回波2处能够有效的解混叠，例如通过下式：
其中，
表示展开相位，
为每个回波处的原始相位，通过增加在回波1和回波2的展开相位到回波5处，可以得到在回波5(17.5ms)前面1.25ms处的相位，也就是16.25ms处的相位。最终，回波6处的相位，可以使用回波5处的相位结合公式(9)展开的相位，展开回波6处的相位，通过使用大脑对应的mask(掩膜)对生成的脑部相位展开图进行掩膜。随后使用截断k-space分割算法为每个回波生成QSM图像。
步骤S120，针对所提出的快速多对比度磁共振成像序列，采用波浪梯度磁场和相位偏移梯度磁场对回波信号进行调制和欠采样，以进行高质量的图像重建。
在一个实施例中，通过在相位和选层两个方向中添加Wave-CAIPI编码梯度场和相位偏移策略，通过在读出方向造成混叠的方式以及在相位两个方向造成相位偏移，降低重建系统条件数，增大重建图像的信噪比(SNR)，添加wave后的信号模型为：
其中，S表示在并行重建中使用的线圈敏感度，m[x，y，z]为待求解图像，Psf是点扩散函数，M表示采样模板，kx表示Psf沿读出方向的傅里叶变换，Fx表示傅里叶变换，重建可以使用传统的SENSE重建算法，在本发明中，通过使用wave-CAIPI的方式，可以增加重建图像的(SNR)，使得前面的各种定性和定量图像的求解更加准确。具体的图像重建方法可现有技术，在此不再赘述。
在此步骤中，针对所提出的快速多对比度MRI序列采用了波浪梯度磁场和相位偏移梯度磁场(wave-CAIPI)技术对回波信号进行调制和欠采样，以便利用基于多测量压缩感知的图像重建算法对高倍欠采样的数据进行高质量的重建。
相应地，本发明还提供一种用于低场磁共振的多对比度成像设备，用于实现上述方法的一个方面或多个方面。例如，该设备包括：扫描单元，其用于依次采用设定的第一射频脉冲序列和第二射频脉冲序列应用于目标成像区域，其中第一射频脉冲序列采用第一翻转角，第二射频脉冲序列采用第二翻转角；信号采集单元，其用于针对第一射频脉冲序列采集第一组回波信号，并针对第二射频脉冲序列采集第二组回波信号，其中第一组回波信号和第二组回波信号中均包含多个回波；图像重建单元，其用于对所获得的回波信号进行调制和欠采样，进而重建得到目标成像区域的多种对比度的磁共振图像。
综上，本发明设计的磁共振成像序列图，可以利用低场磁共振中T2/T2*相对于高场中变大的固有优势，在充分利用序列空闲时间(如回波2和回波5之间)生成T2W图像，并且进行线性拟合出T2map图像，然后通过可变翻转角射频激发脉冲实现T1图像和aT1图像，然后通过线性拟合的方式得到T2、T2*、PD参数图像。
综上所述，本发明提供的用于低场磁共振的单序列多对比度图像生成方法，利用了低场磁共振中T2/T2*变长的特点，在实现增强T1加权图像的同时，中间插入自旋回波梯度，进行T2加权图和定量图的生成；通过使用不同的翻转角生成T1加权和T1map图像的方法；通过使用Wave-CAIPI的编码方式，使得重建图像的信噪比更低。此外，本发明设计的成像序列除了上述方案外，还包括任何在低场磁共振中使用T2/T2*加长的优势的单序列多对比度图像生成方法。利用本发明提供的多对比度成像序列和对应的图像生成方法，通过单次扫描结合计算可以生成至少10种对比度图像，包括但不限于T1、T2、T2*和PD参数图像，以及T1、T1增强、T2*、T2、PD和QSM加权图像等，并显著缩短了传统的定性和定量图像生成的时间。此外，还通过引入wave-CAIPI技术进一步提高了定量和定性图像的信噪比，增强了重建图像的清晰度。
本发明可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质，其上载有用于使处理器实现本发明的各个方面的计算机可读程序指令。
计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是但不限于电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身，诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波(例如，通过光纤电缆的光脉冲)、或者通过电线传输的电信号。
这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备，或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令，并转发该计算机可读程序指令，以供存储在各个计算/处理设备中的计算机可读存储介质中。
用于执行本发明操作的计算机程序指令可以是汇编指令、指令集架构(ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码，所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++、Python等，以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中，通过利用计算机可读程序指令的状态信息来个性化定制电子电路，例如可编程逻辑电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA)，该电子电路可以执行计算机可读程序指令，从而实现本发明的各个方面。
这里参照根据本发明实施例的方法、装置(系统)和计算机程序产品的流程图和/或框图描述了本发明的各个方面。应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机可读程序指令实现。
这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器，从而生产出一种机器，使得这些指令在通过计算机或其它可编程数据处理装置的处理器执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令存储在计算机可读存储介质中，这些指令使得计算机、可编程数据处理装置和/或其他设备以特定方式工作，从而，存储有指令的计算机可读介质则包括一个制造品，其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的各个方面的指令。
也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、或其它设备上，使得在计算机、其它可编程数据处理装置或其它设备上执行一系列操作步骤，以产生计算机实现的过程，从而使得在计算机、其它可编程数据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方框中规定的功能/动作。
附图中的流程图和框图显示了根据本发明的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或指令的一部分，所述模块、程序段或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。对于本领域技术人员来说公知的是，通过硬件方式实现、通过软件方式实现以及通过软件和硬件结合的方式实现都是等价的。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。本发明的范围由所附权利要求来限定。
Claims (10)
1.一种用于低场磁共振的多对比度成像方法，包括以下步骤：
步骤S1：依次采用设定的第一射频脉冲序列和第二射频脉冲序列应用于目标成像区域，其中第一射频脉冲序列采用第一翻转角，第二射频脉冲序列采用第二翻转角；
步骤S2：针对第一射频脉冲序列采集第一组回波信号，并针对第二射频脉冲序列采集第二组回波信号，其中第一组回波信号和第二组回波信号中均包含多个回波；
步骤S3：对所获得的回波信号进行调制和欠采样，进而重建得到目标成像区域的多种对比度的磁共振图像。
2.根据权利要求1所述的方法，其特征在于，针对第一射频脉冲序列采集第一组回波信号包括：利用双回波梯度重聚回波成像序列对读出两组梯度回波信号，标记为第一回波和第二回波；采用回波链长为2的快速自旋回波脉冲序列读出两组自旋回波信号，标记为第三回波和第四回波；
针对第二射频脉冲序列采集第二组回波信号包括：使用双回波梯度重聚回波成像序列采集两组梯度回波信号，标记为第五回波和第六回波。
3.根据权利要求2所述的方法，其特征在于，对于第一回波和第二回波采用非线性拟合方式生成T2*图像；对于第三回波和第四回波采用非线性拟合方式生成T2图像；第五回波和第六回波用于生成T1图像；对于第一回波和第三回波，通过傅里叶变换方式生成T1加权图像和T2加权图像；对于第一回波和第二回波通过傅里叶变换方式生成T2*加权图像。
4.根据权利要求3所述的方法，其特征在于，还包括：利用所述第一翻转角和所述第二翻转角，通过从第一回波和第五回波生成质子密度加权图和真实质子密度加权图。
5.根据权利要求3所述的方法，其特征在于，还包括：利用第一回波和第五回波获得增强T1加权图，表示为：
aT1＝s(θ2,TEn)-λs(θ1,TEn)
其中，θ2第二射频脉冲序列对应的翻转角，θ1表示第一射频脉冲序列对应的翻转角，
k为提取的射频发射场的变化，TEn为第n个回波。
6.根据权利要求3所述的方法，其特征在于，在所述第一翻转角和所述第二翻转角设置为不同值的情况下，利用第一回波、第二回波、第五回波和第六回波使用3D相位展开算法计算磁化率敏感加权图QSM，其中第二回波处的解混频公式表示为：
其中，
表示展开相位，
为每个回波处的原始相位，TEn为第n个回波。
7.根据权利要求1所述的方法，其特征在于，在步骤S3中，通过在相位和选层两个方向中添加相位偏移梯度磁场编码梯度场和相位偏移，信号模型表示为：
其中，S表示在并行重建中使用的线圈敏感度，m[x,y,z]为待求解图像，Psf是点扩散函数，M表示采样模板，kx表示Psf沿读出方向的傅里叶变换，Fx表示傅里叶变换。
8.根据权利要求1所述的方法，其特征在于，所述第一翻转角和所述第二翻转角设置为相同或不同值。
9.一种用于低场磁共振的多对比度成像设备，包括：
扫描单元：用于依次采用设定的第一射频脉冲序列和第二射频脉冲序列应用于目标成像区域，其中第一射频脉冲序列采用第一翻转角，第二射频脉冲序列采用第二翻转角；
信号采集单元：用于针对第一射频脉冲序列采集第一组回波信号，并针对第二射频脉冲序列采集第二组回波信号，其中第一组回波信号和第二组回波信号中均包含多个回波；
图像重建单元：用于对所获得的回波信号进行调制和欠采样，进而重建得到目标成像区域的多种对比度的磁共振图像。
10.一种计算机可读存储介质，其上存储有计算机程序，其中，该程序被处理器执行时实现根据权利要求1至8中任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111376303.XA CN114137462B (zh) | 2021-11-19 | 2021-11-19 | 一种用于低场磁共振的多对比度成像方法和设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111376303.XA CN114137462B (zh) | 2021-11-19 | 2021-11-19 | 一种用于低场磁共振的多对比度成像方法和设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114137462A true CN114137462A (zh) | 2022-03-04 |
| CN114137462B CN114137462B (zh) | 2023-10-20 |
Family
ID=80390320
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111376303.XA Active CN114137462B (zh) | 2021-11-19 | 2021-11-19 | 一种用于低场磁共振的多对比度成像方法和设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114137462B (zh) |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080319301A1 (en) * | 2007-06-25 | 2008-12-25 | General Electric Company | Method and apparatus for generating a flip angle schedule for a spin echo train pulse sequence |
| US20100045288A1 (en) * | 2008-08-20 | 2010-02-25 | Wilhelm Horger | Magnetic resonance method and apparatus for generating different weighted images from the same magnetic resonance echo signal evolution |
| CN105249967A (zh) * | 2015-11-18 | 2016-01-20 | 深圳先进技术研究院 | 一种正对比磁共振成像方法及装置 |
| DE102015222833A1 (de) * | 2015-11-19 | 2017-05-24 | Siemens Healthcare Gmbh | Verfahren und Vorrichtung zur Ansteuerung eines Magnetresonanzbildgebungssystems zur simultanen Aufnahme von mehreren Teilvolumina |
| CN107271937A (zh) * | 2017-07-04 | 2017-10-20 | 大连锐谱科技有限责任公司 | 一种三维多参数加权磁共振成像的同步采集和校准方法 |
| CN108872897A (zh) * | 2018-04-19 | 2018-11-23 | 上海市东方医院 | 核磁共振t2图像成像方法 |
| WO2019038192A1 (en) * | 2017-08-24 | 2019-02-28 | Koninklijke Philips N.V. | MAGNETIC RESONANCE IMAGING (MRI) WITH WATER SEPARATION / DIXON TYPE GREASE |
| CN109549645A (zh) * | 2017-09-27 | 2019-04-02 | 三星电子株式会社 | 磁共振成像设备以及产生磁共振图像的方法 |
| CN110174632A (zh) * | 2019-06-10 | 2019-08-27 | 上海东软医疗科技有限公司 | 磁共振成像方法、装置、成像设备及磁共振成像系统 |
| CN111505553A (zh) * | 2019-07-12 | 2020-08-07 | 上海联影医疗科技有限公司 | 磁共振成像系统和方法 |
| CN113030817A (zh) * | 2021-03-02 | 2021-06-25 | 深圳市儿童医院 | 一种磁共振成像方法、设备及存储介质 |
| CN113017596A (zh) * | 2021-03-09 | 2021-06-25 | 深圳高性能医疗器械国家研究院有限公司 | 一种磁共振多参数定量方法及其应用 |
| CN113366330A (zh) * | 2018-07-31 | 2021-09-07 | 海珀菲纳股份有限公司 | 低场扩散加权磁共振成像 |
| CN113534032A (zh) * | 2020-04-14 | 2021-10-22 | 上海联影医疗科技股份有限公司 | 磁共振成像方法和系统 |
-
2021
- 2021-11-19 CN CN202111376303.XA patent/CN114137462B/zh active Active
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080319301A1 (en) * | 2007-06-25 | 2008-12-25 | General Electric Company | Method and apparatus for generating a flip angle schedule for a spin echo train pulse sequence |
| US20100045288A1 (en) * | 2008-08-20 | 2010-02-25 | Wilhelm Horger | Magnetic resonance method and apparatus for generating different weighted images from the same magnetic resonance echo signal evolution |
| CN105249967A (zh) * | 2015-11-18 | 2016-01-20 | 深圳先进技术研究院 | 一种正对比磁共振成像方法及装置 |
| DE102015222833A1 (de) * | 2015-11-19 | 2017-05-24 | Siemens Healthcare Gmbh | Verfahren und Vorrichtung zur Ansteuerung eines Magnetresonanzbildgebungssystems zur simultanen Aufnahme von mehreren Teilvolumina |
| CN107271937A (zh) * | 2017-07-04 | 2017-10-20 | 大连锐谱科技有限责任公司 | 一种三维多参数加权磁共振成像的同步采集和校准方法 |
| WO2019038192A1 (en) * | 2017-08-24 | 2019-02-28 | Koninklijke Philips N.V. | MAGNETIC RESONANCE IMAGING (MRI) WITH WATER SEPARATION / DIXON TYPE GREASE |
| CN111133327A (zh) * | 2017-08-24 | 2020-05-08 | 皇家飞利浦有限公司 | Dixon型水/脂肪分离MR成像 |
| CN109549645A (zh) * | 2017-09-27 | 2019-04-02 | 三星电子株式会社 | 磁共振成像设备以及产生磁共振图像的方法 |
| CN108872897A (zh) * | 2018-04-19 | 2018-11-23 | 上海市东方医院 | 核磁共振t2图像成像方法 |
| CN113366330A (zh) * | 2018-07-31 | 2021-09-07 | 海珀菲纳股份有限公司 | 低场扩散加权磁共振成像 |
| CN110174632A (zh) * | 2019-06-10 | 2019-08-27 | 上海东软医疗科技有限公司 | 磁共振成像方法、装置、成像设备及磁共振成像系统 |
| CN111505553A (zh) * | 2019-07-12 | 2020-08-07 | 上海联影医疗科技有限公司 | 磁共振成像系统和方法 |
| CN113534032A (zh) * | 2020-04-14 | 2021-10-22 | 上海联影医疗科技股份有限公司 | 磁共振成像方法和系统 |
| CN113030817A (zh) * | 2021-03-02 | 2021-06-25 | 深圳市儿童医院 | 一种磁共振成像方法、设备及存储介质 |
| CN113017596A (zh) * | 2021-03-09 | 2021-06-25 | 深圳高性能医疗器械国家研究院有限公司 | 一种磁共振多参数定量方法及其应用 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114137462B (zh) | 2023-10-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107750338B (zh) | 使用星形叠层采集的mr成像方法和设备 | |
| Knoll et al. | Second order total generalized variation (TGV) for MRI | |
| Hargreaves et al. | Variable‐rate selective excitation for rapid MRI sequences | |
| Zaitsev et al. | Point spread function mapping with parallel imaging techniques and high acceleration factors: Fast, robust, and flexible method for echo‐planar imaging distortion correction | |
| Lee et al. | On the signal‐to‐noise ratio benefit of spiral acquisition in diffusion MRI | |
| Chen et al. | Partial Fourier transform reconstruction for single‐shot MRI with linear frequency‐swept excitation | |
| Zahneisen et al. | Extended hybrid-space SENSE for EPI: Off-resonance and eddy current corrected joint interleaved blip-up/down reconstruction | |
| Dai et al. | Simultaneous multislice accelerated interleaved EPI DWI using generalized blipped‐CAIPI acquisition and 3D K‐space reconstruction | |
| US7683618B2 (en) | Slice-selective tunable-flip adiabatic low peak power excitation | |
| US7741842B2 (en) | Calibration maps for parallel imaging free of chemical shift artifact | |
| CN102713658B (zh) | 磁化率梯度绘图 | |
| CN104939828A (zh) | 磁共振成像方法 | |
| Chen et al. | An aliasing artifacts reducing approach with random undersampling for spatiotemporally encoded single-shot MRI | |
| Balaban et al. | Basic principles of cardiovascular magnetic resonance | |
| Chatnuntawech et al. | Vectorial total generalized variation for accelerated multi-channel multi-contrast MRI | |
| Dong et al. | Model‐based reconstruction for simultaneous multislice and parallel imaging accelerated multishot diffusion tensor imaging | |
| Liao et al. | Highly efficient MRI through multi-shot echo planar imaging | |
| US9841481B2 (en) | Multislice acquisition with incoherent aliasing | |
| Hilbert et al. | True constructive interference in the steady state (trueCISS) | |
| Cao et al. | Efficient T2 mapping with blip‐up/down EPI and gSlider‐SMS (T2‐BUDA‐gSlider) | |
| Zhang et al. | HF-SENSE: an improved partially parallel imaging using a high-pass filter | |
| Benjamin et al. | Multi-shot Echo Planar Imaging for accelerated Cartesian MR Fingerprinting: an alternative to conventional spiral MR Fingerprinting | |
| Chen et al. | Sequential combination of parallel imaging and dynamic artificial sparsity framework for rapid free‐breathing golden‐angle radial dynamic MRI: K‐T ARTS‐GROWL | |
| CN114137462B (zh) | 一种用于低场磁共振的多对比度成像方法和设备 | |
| Zhou et al. | Neural network enhanced 3D turbo spin echo for MR intracranial vessel wall imaging |
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 |