CN111090069B - 定量磁共振成像参数确定方法、装置、设备及存储介质 - Google Patents
定量磁共振成像参数确定方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN111090069B CN111090069B CN201911148335.7A CN201911148335A CN111090069B CN 111090069 B CN111090069 B CN 111090069B CN 201911148335 A CN201911148335 A CN 201911148335A CN 111090069 B CN111090069 B CN 111090069B
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- signal
- radio frequency
- voxel
- mrf
- 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
- 238000002595 magnetic resonance imaging Methods 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims abstract description 82
- 239000013598 vector Substances 0.000 claims abstract description 72
- 230000005415 magnetization Effects 0.000 claims abstract description 71
- 230000005284 excitation Effects 0.000 claims description 61
- 238000004590 computer program Methods 0.000 claims description 24
- 238000005070 sampling Methods 0.000 claims description 23
- 238000002075 inversion recovery Methods 0.000 claims description 20
- 238000004088 simulation Methods 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 claims description 7
- 230000035945 sensitivity Effects 0.000 abstract description 8
- 238000001514 detection method Methods 0.000 abstract description 7
- 210000001519 tissue Anatomy 0.000 description 60
- 238000010586 diagram Methods 0.000 description 13
- 238000002360 preparation method Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 230000006870 function Effects 0.000 description 8
- 230000008520 organization Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 6
- 229910052739 hydrogen Inorganic materials 0.000 description 6
- 239000001257 hydrogen Substances 0.000 description 6
- 230000009471 action Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 3
- 238000001727 in vivo Methods 0.000 description 3
- 230000001575 pathological effect Effects 0.000 description 3
- 238000011002 quantification Methods 0.000 description 3
- 230000002238 attenuated effect Effects 0.000 description 2
- 238000010304 firing Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005281 excited state Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本申请属于磁共振成像技术领域，提供了一种定量磁共振成像参数确定方法、装置、设备及存储介质。方法包括基于预设的MRF脉冲序列获取被试对象的多个磁共振图像，从多个磁共振图像中提取每个体素的演变信号；根据预设的字典中与每个体素的演变信号匹配的模拟信号，确定每个体素的横向弛豫时间；其中，MRF脉冲序列包括顺序激发的多个射频脉冲串，射频脉冲串用于激发被试对象的横向磁化矢量；演变信号用于表征每个体素在多个磁共振图像中信号强度的变化。本申请实施例提供的定量磁共振成像参数确定方法增加了对横向弛豫时间检测的灵敏度，提高了基于MRF检测组织横向弛豫时间的准确度。
Description
技术领域
本申请涉及磁共振成像技术领域，尤其涉及一种定量磁共振成像参数确定方法、装置、设备及存储介质。
背景技术
定量磁共振成像对于人体病理组织的检测具有组织特异性，可以帮助医生更好地区分健康组织和病理组织，具有非常重要的临床诊断和研究价值。
磁共振指纹成像(Magnetic Resonance Fingerprinting，MRF)是一种新型的快速定量磁共振成像参数确定技术。MRF采用一组伪随机变化的脉冲序列扫描被试对象，使得来自被试对象不同的组织类型(以体素为载体)呈现特有的信号演变(即指纹)。然后将得到的指纹与预先获取的字典中的多个字典条目进行匹配，以得到各组织类型多个组织属性参数的定量值。组织属性参数包括纵向弛豫时间(以下简称T1)和横向弛豫时间(以下简称T2)。
目前，MRF的脉冲序列主要是基于反转恢复的稳态自由进动序列，反转恢复序列对T1特别敏感，而稳态自由进动序列具有相对高的信噪比，因此基于该脉冲序列对获取不同被试对象的组织属性参数具有一定的优势。但是由于磁场的不均匀性以及T2持续时间较短，基于当前的MRF脉冲序列获得的磁共振图像中，对各组织类型T2差异不敏感，导致获得的T2准确度较低。
发明内容
有鉴于此，本申请实施例提供了一种定量磁共振成像参数确定方法、装置、设备及存储介质，以解决现有技术中定量磁共振成像中横向弛豫时间不准确的技术问题。
第一方面，本申请实施例提供了一种定量磁共振成像参数确定方法，包括：
基于预设的MRF脉冲序列获取被试对象的多个磁共振图像；其中，MRF 脉冲序列包括顺序激发的多个射频脉冲串，射频脉冲串用于激发被试对象的磁化矢量；
从多个磁共振图像中提取每个体素的演变信号；其中，演变信号用于表征每个体素在多个磁共振图像中信号强度的变化；
根据预设的字典中与每个体素的演变信号最匹配的模拟信号，确定每个体素所代表组织的横向弛豫时间；其中，字典包括对MRF脉冲序列进行信号模拟演变生成的多个模拟信号，以及与每个模拟信号一一对应的横向弛豫时间。
在第一方面的一种可能的实现方式中，射频脉冲串包括顺序激发的90°激发脉冲、至少一个180°激发脉冲和-90°激发脉冲。
在第一方面的一种可能的实现方式中，射频脉冲串中的相邻两个激发脉冲之间的间隔时间相同。
在第一方面的一种可能的实现方式中，每个射频脉冲串的持续时间为磁化矢量翻转至横向平面后的停留时间，多个所述射频脉冲串的持续时间依次递增。
在第一方面的一种可能的实现方式中，每个射频脉冲串之后跟随有一个用于采样的射频脉冲，每个射频脉冲的翻转角均小于等于10°。
在第一方面的一种可能的实现方式中，根据预设的字典中与每个体素的演变信号最匹配的模拟信号，确定每个体素所代表组织的横向弛豫时间，包括：
将每个体素的时间信号与预设的字典中所有的模拟信号进行比较，获得字典中与各体素的磁共振信号匹配度最高的模拟信号；
针对每个体素，将与体素的演变信号匹配度最高的模拟信号对应的横向弛豫时间作为体素的横向弛豫时间。
在第一方面的一种可能的实现方式中，MRF脉冲序列还包括基于反转恢复的梯度回波序列。
第二方面，本申请实施例提供了一种定量磁共振成像参数确定装置，包括：
预成像模块，用于基于预设的MRF脉冲序列获取被试对象的多个磁共振图像；其中，MRF脉冲序列包括顺序激发的多个射频脉冲串，射频脉冲串用于激发被试对象的磁化矢量；
信号获取模块，用于从多个磁共振图像中提取每个体素的演变信号；其中，演变信号用于表征每个体素在多个磁共振图像中信号强度的变化；
参数确定模块，用于根据预设的字典中与每个体素的演变信号最匹配的模拟信号，确定每个体素所代表组织的横向弛豫时间；其中，字典包括对MRF 脉冲序列进行信号模拟演变生成的多个模拟信号，以及与每个模拟信号一一对应的横向弛豫时间。
第三方面，本申请实施例提供了一种定量磁共振成像参数确定设备，包括存储器、处理器以及存储在存储器中并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述第一方面任一项方法的步骤。
第四方面，本申请实施例提供了一种计算机可读存储介质，计算机可读存储介质存储有计算机程序，计算机程序被处理器执行时实现上述第一方面任一项方法的步骤。
第五方面，本申请实施例提供了一种计算机程序产品，当计算机程序产品在终端设备上运行时，使得终端设备执行上述第一方面中任一项的方法。
本申请实施例提供的定量磁共振成像参数确定方法，基于预设的MRF脉冲序列获取被试对象的多个磁共振图像；从多个磁共振图像中提取每个体素的演变信号，根据预设的字典中与每个体素的演变信号匹配的模拟信号，确定每个体素所代表组织的横向弛豫时间，实现了磁共振成像中横向弛豫时间的定量。与现有技术中的MRF脉冲序列相比，本申请实施例中的MRF序列包括顺序激发的多个射频脉冲串，在MRF序列被执行以进行磁共振成像过程中，每个射频脉冲串中的射频脉冲使得被试对象的体内组织的氢核翻转产生横向磁化矢量，并在射频脉冲串的持续时间内，对横向磁化矢量进行再次激发，实现了T2驰豫的加权，因此本申请实施例提供的定量磁共振成像参数确定方法，提高了横向弛豫时间检测的灵敏度，可获得准确的横向弛豫时间。
可以理解的是，上述第二方面至第五方面的有益效果可以参见上述第一方面中的相关描述，在此不再赘述。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本申请一实施例提供的定量磁共振成像系统的架构示意图；
图2是本申请一实施例提供的定量磁共振成像参数确定方法的流程示意图；
图3是本申请一实施例提供的射频脉冲串的示意图；
图4是本申请一实施例提供的单个体素的演变信号；
图5是本申请一实施例提供的确定体素横向弛豫时间的流程示意图；
图6是本申请另一实施例提供的脉冲序列的示意图；
图7是本申请一实施例提供的伪随机变化的翻转角和重复时间示意图；
图8是本申请一实施例提供的字典示意图；
图9是本申请一实施例提供的定量磁共振成像参数确定装置的组成示意图；
图10是本申请一实施例提供的定量磁共振成像参数确定设备的结构示意图。
具体实施方式
以下描述中，为了说明而不是为了限定，提出了诸如特定系统结构、技术之类的具体细节，以便透彻理解本申请实施例。然而，本领域的技术人员应当清楚，在没有这些具体细节的其它实施例中也可以实现本申请。在其它情况中，省略对众所周知的系统、装置、电路以及方法的详细说明，以免不必要的细节妨碍本申请的描述。
应当理解，当在本申请说明书和所附权利要求书中使用时，术语“包括”指示所描述特征、整体、步骤、操作、元素和/或组件的存在，但并不排除一个或多个其它特征、整体、步骤、操作、元素、组件和/或其集合的存在或添加。
另外，在本申请说明书和所附权利要求书的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
首先，对本申请中的涉及到的概念进行解释如下：
磁共振成像(Magnetic Resonance Imaging，MRI)基本原理为：组织的氢原子自旋运动产生磁矩。在强均匀的主磁场作用下，原本无规律排列的氢原子自旋磁矩会沿主磁场方向排列，形成宏观磁化矢量(包括纵向磁化矢量和横向磁化矢量)。在射频脉冲激励下，宏观磁化矢量将翻转到与主磁场垂直的方向，氢原子进动旋转中磁化矢量切割射频线圈，产生电磁感应信号，经过数据重建形成磁共振图像。
磁共振指纹成像(Magnetic Resonance Fingerprinting，MRF)：是一种快速、可同时获得多个组织属性参数的定量磁共振成像方法。
脉冲序列(Pulse Sequence)是指射频脉冲、空间编码、回波采样在一定时序上的组合。
重复时间(Repetition Time，TR)指脉冲序列中两次相邻激发脉冲之间的时间间隔。
回波时间(Echedelay Time，TE)指激发脉冲激发后到采集回波信号之间的时间。
弛豫：在激发脉冲的激发下，被试对象内的质子吸收能量处于激发状态；激发脉冲终止后，处于激发状态的质子恢复到其原始状态的过程称为弛豫。
T1：纵向弛豫时间，纵向磁化矢量从零恢复至总信号强度的三分之二所需的时间。
T2：横向弛豫时间，横向磁化矢量从百分百衰减至三分之一所需的时间。
翻转角(Flip Angle,FA)：激发脉冲使磁化矢量沿主磁场方向偏离的角度。 FA为90°时，磁化矢量垂直于主磁场方向，翻转至横向平面。
定量磁共振成像参数确定是指确定被试对象的组织属性参数的值，其中被试对象的组织属性参数包括T1、T2、质子密度等。定量组织属性参数在绝对意义上可以帮助医生更好地区分健康组织和病理组织，也便于医生在后续研究中更客观地比较不同的检查。目前应用最广泛的组织属性参数包括T1和T2。
磁共振指纹成像(以下简称MRF)是一种快速定量磁共振成像参数确定方法。相比于传统的T1、T2定量方法，MRF可以同时获得多个组织属性参数的参数值。MRF主要包括以下步骤：在脉冲序列的每一次激发中采用不同的TR 和FA进行被试对象感兴趣区域的数据采集；然后根据脉冲序列的参数，基于经典的一阶梯布洛赫模型生成字典，其中，字典中包括多个模拟时间信号和各个模拟时间信号对应的组织属性参数值；最后将采集到的数据重建后，获取单个体素的随时间变化的演变信号(也即“指纹”)；将每个体素的演变信号与字典中包含的模拟时间信号进行一一比较，从字典中查找获得各体素的最佳匹配模拟时间信号，将最佳匹配模拟时间信号对应的组织属性参数值作为体素的组织属性参数值。由于被试对象不同组织类型的组织属性参数不同，基于相同的脉冲序列可以获得的不同组织类型的不同演变信号，重复上述过程可获得被试对象不同组织类型(所有体素)的组织属性参数值。
图1为本申请一实施例提供的定量磁共振成像系统的架构示意图，如图1 所示，定量磁共振成像系统包括磁共振设备10和图像重建设备20。
其中磁共振设备10用于执行用户输入的扫描指令并获取磁共振数据，磁共振设备10包括显示器、一个或多个输入设备以及处理器。显示器提供用户输入界面，用于接收用户输入的扫描指令，该扫描指令用于生成脉冲序列。
磁共振设备10可以通过通信系统与图像重建设备20连接，用于将用户输入的处理指令以及磁共振数据发送至图像重建设备20，通信系统包括有线或无线的网络连接。
图像重建设备20为服务器，其接收磁共振设备10发送的磁共振数据和处理指令，并根据处理指令来处理磁共振数据，以生成重建的磁共振图像，并根据该磁共振图像确定组织属性参数的定量值。处理指令包括磁共振数据的傅里叶变换、图像重建、模式匹配以及字典生成等。
图像重建设备20也可以包括显示器、一个或多个输入设备，用于接收用户输入的扫描指令或处理指令。
定量磁共振成像参数确定系统还可以包括一个或多个联网用户终端，例如，移动手机、智能计算机等。联网用户终端可以通过通信系统与磁共振设备10 连接，以实现磁共振设备10的远程访问和控制。
实际应用中，磁共振设备10接收用户输入的扫描指令，并生成相应的脉冲信号，具体包括产生用于执行规定扫描的梯度波形以及用于控制射频发生器产生期望的频率、相位和脉冲振幅波形的射频脉冲。
当时激发脉冲终止后，磁共振设备10采集获得磁共振数据并将该磁共振数据发送图像重建设备20,图像重建设备20执行磁共振指纹成像(MRF)中的图像重建以及模式匹配工作，确定组织属性参数的定量值。
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行示例性地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于下文中列举的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图2为本申请一实施例提供的定量磁共振成像参数确定方法的流程示意图，主要涉及在MRF中如何准确定量T2的过程。如图1所示，该定量磁共振成像参数确定方法包括：
S201、基于预设的MRF脉冲序列获取被试对象的多个磁共振图像；其中， MRF脉冲序列包括顺序激发的多个射频脉冲串，射频脉冲串用于激发被试对象的磁化矢量。
被试对象包括多个不同组织类型，基于相同的脉冲序列扫描被试对象，不同组织类型可以获得不同演变信号。示例性的，被试对象可以包括仿体、活体 (动物或人体)，离体器官或组织等。
其中，不同组织类型的共振特性不同。示例性的，组织类型可以包括水、脂肪、骨头、肌肉、血液等。
MRF脉冲序列为基于MRF框架生成的脉冲序列。不同于常规的磁共振成像，MRF框架采用一系列变化的序列块来生成目标的脉冲序列，这些序列块产生射频脉冲顺序激发并施加到不同组织类型中以产生不同的信号演变。
在本实施例中，MRF脉冲序列包括顺序激发的多个射频脉冲串，多个射频脉冲串用于激发被试对象的磁化矢量以进行T2磁化准备。
其中，T2磁化准备是指对被试对象的氢核进行横向磁化准备，以提高T2 检测的灵敏度。由于T2用于表征横向磁化矢量的衰减特性，该磁化准备可以为对横向磁化矢量进行多次聚焦，以对被试对象的横向弛豫特异性进行多次编码。
其中，横向磁化矢量为被试对象的多个氢核的磁矩在横向平面的宏观表现。当被试对象进入主磁场后，被试对象的多个氢核的磁矩绕着主磁场的方向章动。每个磁矩都是矢量，所有磁矩的矢量和即为磁化矢量。磁化矢量包括纵向磁化矢量和横向磁化矢量。在没有激励脉冲的情况下，磁化矢量的方向是平行于主磁场的，平行于主磁场方向的磁化矢量为纵向磁化矢量，而垂直于磁场方向的磁化矢量为横向磁化矢量。
在本实施例中，每个射频脉冲串包括顺序激发的激励脉冲、180°翻转脉冲以及90°翻转脉冲，其中激励脉冲用于将磁化矢量翻转至横向平面，获得较大的横向磁化矢量；180°翻转脉冲用于将构成横向磁化矢量的所有磁矩的方向进行180°翻转；翻转脉冲用于将弛豫后的横向磁化矢量翻转回纵向平面。在射频脉冲的激励脉冲和90°翻转脉冲之间，180°翻转脉冲每激发一次，即对当前横向磁化矢量的矢量方向进行180°翻转，可以实现对横向弛豫特异性的多次编码，提高T2检测的灵敏度。其中，180°翻转脉冲可以为一个或多个。
在一种实施方式中，射频脉冲串包括顺序激发的90°激发脉冲、至少一个 180°激发脉冲和-90°激发脉冲；其中，90°激发脉冲为上述的激励脉冲，180°激发脉冲为上述的180°翻转脉冲，-90°激发脉冲为上述的90°翻转脉冲。射频脉冲串中的相邻两个激发脉冲之间的间隔时间相同。其中，180°激发脉冲可以有多个，每一个180°激发脉冲进行一次横向磁化矢量在横向平面的翻转，通过设置多个180°激发脉冲，可以对横向弛豫特异性进行多次编码。。
请一并参阅图3，图3为本实施例提供的一种射频脉冲串的示意图，如图3 所示，脉冲序列包括两个射频脉冲串，每个射频脉冲串均包括顺序激发的90°激发脉冲、一个180°激发脉冲和一个-90°激发脉冲。90°激发脉冲和180°激发脉冲之间的间隔时间，与180°激发脉冲和-90°激发脉冲之间的间隔时间相同。基于该脉冲序列扫描被试对象时，90°激发脉冲直接将磁化矢量完全翻转至横向平面，获得了最大的横向磁化矢量，90激发脉冲激发完毕后，横向磁化矢量进入弛豫状态，开始衰减，直至180°激发脉冲将当前横向磁化矢量中所有磁矩的矢量方向进行180°翻转，获得翻转后横向磁化矢量，翻转后的横向磁化矢量继续衰减，直至-90°激发脉冲将衰减后的横向磁化矢量翻转回纵向平面，至此完成被试对象的T2磁化准备。
每个射频脉冲串的持续时间为磁化矢量翻转至横向平面后的停留时间，具体体现为该射频脉冲串的激励脉冲和翻转脉冲之间的间隔时间；示例性的，如图3所示，每个射频脉冲串的持续时间为90°激发脉冲和-90°激发脉冲之间的间隔时间。
每个射频脉冲串的持续时间可以不同，由于不同体素所代表的组织的T2 不同，对不同体素的T2的敏感的准备时间也不同。将每个射频脉冲串的持续时间设置为对不同体素所代表的组织的T2敏感的准备时间，多个射频脉冲串的持续时间依次递增。示例性的，射频激发脉冲串包括3个，3个射频激发脉冲串的持续时间分别为5毫秒，25毫秒以及35毫秒。射频激发脉冲串的持续时间一般不超过100毫秒。
在本实施例中，视频脉冲串用于进行被试对象的T2准备，因此每个射频脉冲串之后，可以设置一个用于采样的射频脉冲，其中每个射频脉冲的翻转角均小于等于10°。实际应用中，每一个射频激发脉冲串结束后，施加一个射频脉冲，此时每个体素产生一个回波，通过一个采样窗(模数转换装置A/D)进行回波信号的采集。
本实施中，基于上述MRF脉冲序列，进行被试对象的扫描，获取被试对象的磁共振图像，具体包括S2011至S2013：
S2011、获取MRF脉冲序列的脉冲序列参数；
脉冲序列参数用于表征MRF脉冲序列的射频脉冲串、用于采样的射频脉冲在时序上的组合。同时，脉冲序列参数还包括了各脉冲的翻转角FA、重复时间TR、回波时间TE以及射频脉冲在每一次激发中的相位变化。
S2012、将脉冲序列参数发送至磁共振设备，以使磁共振设备基于脉冲序列参数生成脉冲信号，并获取被试对象的多个采样数据。
磁共振设备基于脉冲序列参数获取采样数据的方式可以基于现有技术实现，在此不再赘述。
S2013、接收多个采样数据，并对每个采样数据进行傅里叶变换，生成对应每个采样数据的被试对象的磁共振图像。
磁共振设备将采样得到的采样数据按一定的编排方式进行存放，构成采样数据空间。根据采样数据空间每一个采样点的信号强度及其所在的位置，生成磁共振图像。具体地，对于直线式采样得到的数据，其数据均匀分布网格点上，经过傅里叶逆变换就可以得到二维平面内上，根据每个体素的信号强度，将其转换为相应的灰度值，就得到磁共振图像。
S202、从多个磁共振图像中提取每个体素的演变信号；其中，演变信号用于表征每个体素在多个磁共振图像中信号强度的变化。
每个体素在不同的采样时间信号强度不同，体现到磁共振图像上体现为灰度值大小不同。从多个磁共振图像中获取每个体素点的信号强度，并按照磁共振图像的采样时间进行排序组合，即获得每个体素的演变信号。
请一并参阅图4，图4为本申请一实施例提供的单个体素的演变信号。实际应用中，预先设置采样时间点为1000，在每一个采样时间点执行该MRF脉冲序列，即得到1000张欠采样的磁共振图像。将1000张磁共振图像中每个像素点的信号强度按照时间进行组合，获得每个像素点的演变信号曲线，即各体素的演变信号。应理解的是，磁共振图像中的像素点与被试对象的体素一一对应。
S203、根据预设的字典中与每个体素的演变信号最匹配的模拟信号，确定每个体素所代表组织的横向弛豫时间；其中，字典包括对所述MRF脉冲序列进行信号模拟演变生成的多个模拟信号，以及与每个模拟信号一一对应的横向弛豫时间。
字典可以预先获得，字典包括基于上述MRF脉冲序列进行信号模拟演变生成的多个模拟信号，以及与每个模拟信号一一对应的横向弛豫时间。
信号模拟演变可以通过布洛赫模型模拟计算获得。每一个模拟信号为字典中的一个字典条目，每条字典条目对应一个T2。
示例性的，首先根据被试对象生理特性，将T2不均匀划分。例如，为了覆盖被试对象，T2的范围设置为0-500毫秒。将小于100毫秒的T2以2毫秒的梯度增加，大于100毫秒小于200毫秒的部分以5毫秒的梯度增加，大于200 毫秒部分以50毫秒的梯度增加，生成多组T2。然后根据步骤S201中的MRF 脉冲序列的参数集合，用布洛赫模型模拟计算获得模拟时间信号，每变换一组 T2，则获得一个字典条目，直至所有组T2都进行模拟计算，生成包含多个字典条目的字典。其中，布洛赫模型表示的是质子在磁场作用下磁化矢量随弛豫时间的变化。
请一并参阅图5，图5主要涉及如何确定每个体素的T2，如图5所示，根据预设的字典中与每个体素的演变信号匹配的模拟信号，确定每个体素的横向弛豫时间，具体包括S501和S502：
S501、将每个体素的时间信号与预设的字典中所有的模拟信号进行比较，获得所述字典中与各体素的磁共振信号匹配度最高的模拟信号。
S502、针对每个体素，将与所述体素的演变信号匹配度最高的模拟信号对应的横向弛豫时间作为所述体素的横向弛豫时间。
首先，对字典中的每一个字典条目进行归一化处理。然后基于点乘操作进行比较。具体地，在获得字典以及每个体素的演变信号后，针对每个体素，将该体素的演变信号与预设字典中每一个字典条目一一进行比较，获得字典中与该体素的演变信号匹配度最高的字典条件；将该字典条目对应的T2作为体素的T2。字典的引入，使得MRF可以融入一些磁共振设备的系统参数，如主磁场的不均匀性以及射频脉冲的不均匀性的，提高了MRF的鲁棒性，使其适用于不同组织位置和不同时间段的定量磁共振成像参数确定。
本申请实施例提供的定量磁共振成像参数确定方法，基于预设的MRF脉冲序列获取被试对象的多个磁共振图像；从多个磁共振图像中提取每个体素的演变信号，根据预设的字典中与每个体素的演变信号匹配的模拟信号，确定每个体素所代表组织的横向弛豫时间，实现了磁共振成像中横向弛豫时间的定量。与现有技术中的MRF脉冲序列相比，本申请实施例中的MRF序列包括顺序激发的多个射频脉冲串，在MRF序列被执行以进行磁共振成像过程中，每个射频脉冲串中的射频脉冲使得被试对象的体内组织的氢核翻转产生横向磁化矢量，并在射频脉冲串的持续时间内，对横向磁化矢量进行再次激发，实现了T2驰豫的加权，因此本申请实施例提供的定量磁共振成像参数确定方法，提高了横向弛豫时间检测的灵敏度，可获得准确的横向弛豫时间。
MRF脉冲序列为基于MRF框架生成的脉冲序列。不同于常规的磁共振成像，MRF框架采用一系列变化的序列块来生成目标的脉冲序列，这些序列块产生射频脉冲顺序激发并施加到不同组织类型中以产生不同的信号演变。并根据该信号演变，可以同时获得多个组织属性参数的参数值。其中，组织属性参数可以包括纵向弛豫时间T1、主磁场强度B0等。
具体地，字典可以包括基于MRF脉冲序列进行信号模拟演变生成的多个模拟信号，以及与每个模拟信号一一对应的特征物理参数。其中特征物理参数除了包括横向弛豫时间之外，还可以包括纵向弛豫时间T1，主磁场强度B0,等。然后根据预设的字典中与每个体素的演变信号最匹配的模拟信号，确定每个体素所代表组织的特征物理参数。
请一并参阅图6，图6为本申请另一实施例提供的脉冲序列的示意图，如图6所示，MRF脉冲序列包括两个顺序激发的射频脉冲串、反转恢复脉冲和平衡稳态自动进动序列。
本实施例中，MRF可以包括T2准备序列块，T2准备序列块用于产生顺序激发的多个射频脉冲串，MRF还可以包括反转恢复脉冲序列块，用于产生反转恢复脉冲。MRF通过序列块的框架设计，实现了多个脉冲序列块的快速添加。
在本实施例中，MRF脉冲序列还可以包括基于反转恢复的梯度回波序列。在MRF脉冲序列每一次激发中，基于反转恢复的梯度回波序列的重复时间和翻转角均不相同。其中，基于反转恢复的梯度回波序列包括顺序激发的反转恢复脉冲和梯度回波序列脉冲，可以通过在MRF框架中添加一个反转恢复序列块和一个梯度回波采集序列块实现。反转恢复序列对T1敏感，通过设置反转恢复序列可以提高磁共振信号对T1的灵敏度，获得准确度T1值。
反转恢复脉冲(inversion recovery，IR)包括顺序激发180°射频脉冲和 90°射频脉冲。两个射频脉冲之间的间隔为反转等待时间。使用时，首先基于 180°射频脉冲激发磁化矢量，然后等待一段时间(反转时间TI)后使用90°射频脉冲在此进行激发，两个射频脉冲形成回波。
示例性的，首先通过反转恢复180°脉冲激发磁化矢量，使其反转180°，之后等待第一个T1的时间，此时纵向磁化矢量将部分恢复，然后施加多个小角度射频脉冲(具有90°脉冲分量)将目前大小的纵向磁化矢量反转到横向坐标，以进行数据采集。待纵向磁化矢量完全恢复稳态之后，改变反转时间TI 的值的大小重复以上操作，选择不同反转时间，可以降低化学移位伪影，获得较准确的T1值。
梯度回波序列用于进行回波采样，梯度回波序列可以为平衡稳态自由进动序列。平衡稳态自由进动序列的特征是横向磁化矢量和纵向磁化矢量均在多个射频激发之后达到稳态。由于激发的射频角度可以很小，TR可以设置的很短，因此平衡稳态自由进动序列的扫描时间相对于普通的梯度回波序列的扫描时间大大减少。
实际应用中，MRF可以通过在每个采集时间点改变脉冲序列参数获得来在空间上不相干和或时间上不相干的演变信号。能够改变的脉冲序列参数包括各脉冲的反转角FA、脉冲相位、重复时间TR、回波时间TE以及采样模式。示例性的，在本实施例中，在MRF脉冲序列的每一次激发中，平衡稳态自动进动序列的重复时间TR和翻转角FA位随机变化，以获得在空间上不相干的演变信号，即各体素的演变信号。平衡稳态自动进动序列的重复时间TR和翻转角 FA位随机变化曲线可以参阅图7。
基于本实施例的MRF脉冲序列进行被试对象的扫描，获取被试对象的磁共振图像的步骤与S2011至S2013相同，仅脉冲序列参数不同。具体地，脉冲序列参数用于表征MRF脉冲序列的射频脉冲串、用于采样的射频脉冲、反转恢复脉冲以及平衡稳态自动进动脉冲在时序上的组合。同时，脉冲序列参数还包括了各脉冲的翻转角FA、重复时间TR、回波时间TE以及射频脉冲在每一次激发中的相位变化。示例性的，平衡稳态自动进动脉冲的翻转角FA和重复时间TR随时间伪随机变化。
当获取被试对象的多个磁共振图像后，从多个磁共振图像中提取被试对象的每个体素的演变信号；其中，演变信号用于表征每个体素在多个磁共振图像中信号强度的变化。然后根据预设的字典中与每个体素的演变信号匹配的模拟信号，确定每个体素的组织属性参数值；其中，字典包括对图6所示MRF脉冲序列进行信号模拟演变生成的多个模拟信号，以及与每个模拟信号一一对应的组织属性参数值。组织属性参数值包括T1和T2。
其中，字典可以通过布洛赫模型模拟预先计算获得。
示例性的，将T1的范围设置为0-5000毫秒，T2的范围设置为0-500毫秒。首先对组织参数值进行排列组合，生成多组组织参数值，其中每组组织参数包括T1和T2，且每组组织参数的T2和T1不完全相同；然后根据步骤S201中的MRF脉冲序列的参数集合，用布洛赫模型模拟计算获得模拟时间信号，每变换一组组织参数值，则获得一个字典条目，直至所有组的组织参数值都进行模拟计算，生成包含多个字典条目的字典。其中，布洛赫模型表示的是质子在磁场作用下磁化矢量随弛豫时间的变化。
请一并参阅图8，图8为本申请实施例提供的字典示意图，如图8所示，每一组组织属性参数值对应一个模拟信号。组织属性参数包括T1和T2，字典的字典条目数目由T2和T1的排列组合数确定。应理解的是，组织属性参数还可以包括质子密度等其他组织属性参数。
本实施例提供的定量磁共振成像方法中，MRF脉冲序列不仅包括用于进行激发磁化矢量的射频脉冲串，还包括基于反转恢复的梯度回波序列，在MRF 脉冲序列的每一次激发中，该梯度回波序列的FA和TR伪随机变化，故可以实现被试对象的各体素的信号在空间上编码。一方面，该MRF脉冲序列通过射频脉冲串和反转恢复脉冲的设置，能够同时获得准确的T1和T2。另一方面， MRF方案中通过字典的设计，可以融入一些磁共振设备的系统参数，如主磁场的不均匀性以及射频脉冲的不均匀性等，提高了MRF的鲁棒性，使其适用于不同组织类型和不同时间段的定量磁共振成像。
应理解，上述实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本申请实施例的实施过程构成任何限定。
基于上述实施例所提供的定量磁共振成像参数确定，本发明实施例进一步给出实现上述方法实施例的装置实施例。
图9为本申请一实施例提供的定量磁共振成像参数确定装置的组成示意图。如图9所示，定量磁共振成像参数确定装置90包括：预成像模块901、信号获取模块902以及参数确定模块903。
预成像模块901，用于基于预设的MRF脉冲序列获取被试对象的多个磁共振图像；其中，MRF脉冲序列包括顺序激发的多个射频脉冲串，射频脉冲串用于激发被试对象的磁化矢量；
信号获取模块902，用于从多个磁共振图像中提取每个体素的演变信号；其中，演变信号用于表征每个体素在多个磁共振图像中信号强度的变化；
参数确定模块903，用于根据预设的字典中与每个体素的演变信号最匹配的模拟信号，确定每个体素所代表组织的横向弛豫时间；其中，字典包括对 MRF脉冲序列进行信号模拟演变生成的多个模拟信号，以及与每个模拟信号一一对应的横向弛豫时间。
可选地，预成像模块901中，射频脉冲串包括顺序激发的90°激发脉冲、至少一个180°激发脉冲和-90°激发脉冲。
射频脉冲串中的相邻两个激发脉冲之间的间隔时间相同。
每个射频脉冲串的持续时间为磁化矢量翻转至横向平面后的停留时间，多个射频脉冲串的持续时间依次递增。
每个射频脉冲串之后跟随有一个用于采样的射频脉冲，每个射频脉冲的翻转角均小于等于10°。
MRF脉冲序列还包括基于反转恢复的梯度回波序列。
可选地，参数确定模块903，具体用于：
将每个体素的时间信号与预设的字典中所有的模拟信号进行比较，获得字典中与各体素的磁共振信号匹配度最高的模拟信号；其中，字典包括对MRF 脉冲序列进行信号模拟演变生成的多个模拟信号，以及每个模拟信号与横向弛豫时间之间的对应关系。
针对每个体素，将与体素的演变信号匹配度最高的模拟信号对应的横向弛豫时间作为体素的横向弛豫时间。
本实施例提供的定量磁共振成像参数确定装置，基于预设的MRF脉冲序列获取被试对象的多个磁共振图像；从多个磁共振图像中提取每个体素的演变信号，根据预设的字典中与每个体素的演变信号匹配的模拟信号，确定每个体素所代表组织的横向弛豫时间，实现了磁共振成像中横向弛豫时间的定量。与现有技术中的MRF脉冲序列相比，本申请实施例中的MRF序列包括顺序激发的多个射频脉冲串，在MRF序列被执行以进行磁共振成像过程中，每个射频脉冲串中的射频脉冲使得被试对象的体内组织的氢核翻转产生横向磁化矢量，并在射频脉冲串的持续时间内，对横向磁化矢量进行再次激发，实现了T2驰豫的加权，因此本申请实施例提供的定量磁共振成像参数确定方法，提高了横向弛豫时间检测的灵敏度，可获得准确的横向弛豫时间。
图9所示实施例提供的定量磁共振成像参数确定装置，可用于执行上述方法实施例中的技术方案，其实现原理和技术效果类似，本实施例此处不再赘述。
所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，仅以上述各功能单元、模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能单元、模块完成，即将所述装置的内部结构划分成不同的功能单元或模块，以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中，上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。另外，各功能单元、模块的具体名称也只是为了便于相互区分，并不用于限制本申请的保护范围。上述系统中单元、模块的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
图10是本申请一实施例提供的定量磁共振成像参数确定设备的示意图。如图10所示，该实施例的定量磁共振成像参数确定终端设备10包括：至少一个处理器1001、存储器1002以及存储在存储器1002中并可在处理器1001上运行的计算机程序。定量磁共振成像参数确定设备还包括通信部件1003，其中，处理器1001、存储器1002以及通信部件1003通过总线1004连接。
处理器1001执行计算机程序时实现上述各个定量磁共振成像参数确定方法实施例中的步骤，例如图2所示实施例中的步骤S201至步骤S203。或者，处理器1001执行所述计算机程序时实现上述各装置实施例中各模块/单元的功能，例如图9所示模块901至903的功能。
示例性的，计算机程序可以被分割成一个或多个模块/单元，一个或者多个模块/单元被存储在存储器1002中，并由处理器1001执行，以完成本申请。所述一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段用于描述计算机程序在定量磁共振成像图像重建设备100中的执行过程。
在本实施例中，定量磁共振成像参数确定设备可以为云端服务器或用户终端。用户终端可以是但不限于各种能运行应用的个人计算机、笔记本电脑、智能手机等。云端服务器可以是实现单一功能的服务器，也可以是实现多种功能的服务器，具体可以是独立的物理服务器，也可以是物理服务器集群。
本领域技术人员可以理解，图10仅仅是定量磁共振成像参数确定设备的示例，并不构成对定量磁共振成像参数确定设备的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如输入输出设备、网络接入设备、总线等。
处理器1001可以是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列 (Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
存储器1002可以是定量磁共振成像参数确定设备的内部存储单元，也可以是定量磁共振成像图像参数确定的外部存储设备，例如插接式硬盘，智能存储卡(Smart MediaCard，SMC)，安全数字(Secure Digital，SD)卡，闪存卡(Flash Card)等。所述存储器1002用于存储所述计算机程序以及定量磁共振成像图像重建设备所需的其他程序和数据。所述存储器1002还可以用于暂时地存储已经输出或者将要输出的数据。
总线可以是工业标准体系结构(Industry Standard Architecture，ISA)总线、外部设备互连(Peripheral Component，PCI)总线或扩展工业标准体系结构 (ExtendedIndustry Standard Architecture，EISA)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示，本申请附图中的总线并不限定仅有一根总线或一种类型的总线。
本申请实施例还提供了一种计算机可读存储介质，计算机可读存储介质存储有计算机程序，计算机程序被处理器执行时实现可实现上述各个方法实施例中的步骤。
本申请实施例提供了一种计算机程序产品，当计算机程序产品在定量磁共振成像参数确定设备上运行时，使得定量磁共振成像参数确定设备执行时实现可实现上述各个方法实施例中的步骤。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请实现上述实施例方法中的全部或部分流程，可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质至少可以包括：能够将计算机程序代码携带到定量磁共振成像参数确定设备的任何实体或装置、记录介质、计算机存储器、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、电载波信号、电信信号以及软件分发介质。例如U盘、移动硬盘、磁碟或者光盘等。在某些司法管辖区，根据立法和专利实践，计算机可读介质不可以是电载波信号和电信信号。
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述或记载的部分，可以参见其它实施例的相关描述。
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本申请的范围。
在本申请所提供的实施例中，应该理解到，所揭露的装置/网络设备和方法，可以通过其它的方式实现。例如，以上所描述的装置/网络设备实施例仅仅是示意性的，例如，所述模块或单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口，装置或单元的间接耦合或通讯连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
以上所述实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围，均应包含在本申请的保护范围之内。
Claims (7)
1.一种定量磁共振成像参数确定方法，其特征在于，包括：
基于预设的MRF脉冲序列获取被试对象的多个磁共振图像；其中，所述MRF脉冲序列包括顺序激发的多个射频脉冲串，所述射频脉冲串用于激发所述被试对象的磁化矢量；
从所述多个磁共振图像中提取每个体素的演变信号；其中，所述演变信号用于表征每个体素在所述多个磁共振图像中信号强度的变化；
根据预设的字典中与每个体素的演变信号最匹配的模拟信号，确定每个体素所代表组织的横向弛豫时间；其中，所述字典包括对所述MRF脉冲序列进行信号模拟演变生成的多个模拟信号，以及与每个模拟信号一一对应的横向弛豫时间；
其中，所述射频脉冲串包括顺序激发的90°激发脉冲、至少一个180°翻转脉冲和-90°翻转脉冲；所述射频脉冲串中的相邻两个激发脉冲之间的间隔时间相同；每一个180°翻转脉冲进行一次横向磁化矢量在横向平面的翻转，对当前横向磁化矢量的矢量方向进行180°翻转，以实现对横向弛豫特异性的多次编码；
每个所述射频脉冲串的持续时间为所述磁化矢量翻转至横向平面后的停留时间，多个所述射频脉冲串的持续时间依次递增，且所述射频脉冲串的持续时间不超过100毫秒。
2.如权利要求1所述的定量磁共振成像参数确定方法，其特征在于，每个射频脉冲串之后跟随有一个用于采样的射频脉冲，每个所述射频脉冲的翻转角均小于等于10°。
3.如权利要求1或2所述的定量磁共振成像参数确定方法，其特征在于，所述根据预设的字典中与每个体素的演变信号最匹配的模拟信号，确定每个体素所代表组织的横向弛豫时间，包括：
将每个体素的时间信号与预设的字典中所有的模拟信号进行比较，获得所述字典中与各体素的磁共振信号匹配度最高的模拟信号；
针对每个体素，将与所述体素的演变信号匹配度最高的模拟信号对应的横向弛豫时间作为所述体素的横向弛豫时间。
4.如权利要求1所述的定量磁共振成像参数确定方法，其特征在于，所述MRF脉冲序列还包括基于反转恢复的梯度回波序列。
5.一种定量磁共振成像参数确定装置，其特征在于，包括：
预成像模块，用于基于预设的MRF脉冲序列获取被试对象的多个磁共振图像；其中，MRF脉冲序列包括顺序激发的多个射频脉冲串，射频脉冲串用于激发被试对象的磁化矢量；
信号获取模块，用于从多个磁共振图像中提取每个体素的演变信号；其中，演变信号用于表征每个体素在多个磁共振图像中信号强度的变化；
参数确定模块，用于根据预设的字典中与每个体素的演变信号最匹配的模拟信号，确定每个体素所代表组织的横向弛豫时间；其中，字典包括对MRF脉冲序列进行信号模拟演变生成的多个模拟信号，以及与每个模拟信号一一对应的横向弛豫时间；
其中，所述射频脉冲串包括顺序激发的90°激发脉冲、至少一个180°翻转脉冲和-90°翻转脉冲；所述射频脉冲串中的相邻两个激发脉冲之间的间隔时间相同；每一个180°翻转脉冲进行一次横向磁化矢量在横向平面的翻转；对当前横向磁化矢量的矢量方向进行180°翻转，以实现对横向弛豫特异性的多次编码；
每个所述射频脉冲串的持续时间为所述磁化矢量翻转至横向平面后的停留时间，多个所述射频脉冲串的持续时间依次递增，且所述射频脉冲串的持续时间不超过100毫秒。
6.一种定量磁共振成像参数确定设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至4任一项所述方法的步骤。
7.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至4任一项所述方法的步骤。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911148335.7A CN111090069B (zh) | 2019-11-21 | 2019-11-21 | 定量磁共振成像参数确定方法、装置、设备及存储介质 |
| PCT/CN2020/129571 WO2021098690A1 (zh) | 2019-11-21 | 2020-11-17 | 定量磁共振成像参数确定方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911148335.7A CN111090069B (zh) | 2019-11-21 | 2019-11-21 | 定量磁共振成像参数确定方法、装置、设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111090069A CN111090069A (zh) | 2020-05-01 |
| CN111090069B true CN111090069B (zh) | 2022-03-29 |
Family
ID=70394042
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911148335.7A Active CN111090069B (zh) | 2019-11-21 | 2019-11-21 | 定量磁共振成像参数确定方法、装置、设备及存储介质 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN111090069B (zh) |
| WO (1) | WO2021098690A1 (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111090069B (zh) * | 2019-11-21 | 2022-03-29 | 深圳先进技术研究院 | 定量磁共振成像参数确定方法、装置、设备及存储介质 |
| CN111685764B (zh) * | 2020-05-13 | 2021-08-20 | 复旦大学 | 一种基于磁共振指纹成像的快速定量T1ρ的方法 |
| CN112617796B (zh) * | 2020-12-30 | 2024-10-15 | 深圳市联影高端医疗装备创新研究院 | 磁共振功能成像方法、装置、设备、存储介质及系统 |
| CN113030817B (zh) * | 2021-03-02 | 2021-12-28 | 深圳市儿童医院 | 一种磁共振成像方法、设备及存储介质 |
| CN113805133B (zh) * | 2021-09-13 | 2022-05-20 | 上海交通大学 | 一种减少图像失真的磁共振平面回波成像方法 |
| CN114089242B (zh) * | 2021-11-03 | 2024-08-23 | 深圳市联影高端医疗装备创新研究院 | 激发脉冲的设计方法、装置、计算机设备和存储介质 |
| CN114217255B (zh) * | 2021-11-29 | 2022-09-20 | 浙江大学 | 一种快速的肝脏多参数定量成像方法 |
Family Cites Families (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5329818B2 (ja) * | 2007-03-26 | 2013-10-30 | 株式会社東芝 | 磁気共鳴イメージング装置および磁気共鳴イメージング装置における撮影条件設定方法 |
| JP5624273B2 (ja) * | 2007-11-22 | 2014-11-12 | 株式会社東芝 | 磁気共鳴イメージング装置 |
| US8797031B2 (en) * | 2011-01-27 | 2014-08-05 | Siemens Medical Solutions Usa, Inc. | MR imaging system for discriminating between imaged tissue types |
| US9618596B2 (en) * | 2012-07-02 | 2017-04-11 | Syntheticmr Ab | Methods and systems for improved magnetic resonance acquisition using a single acquisition type |
| RU2693837C2 (ru) * | 2014-11-14 | 2019-07-05 | Конинклейке Филипс Н.В. | Магнитно-резонансный метод пальцевых отпечатков |
| CN105988098B (zh) * | 2015-01-30 | 2021-07-27 | Ge医疗系统环球技术有限公司 | 磁共振信号采集系统及方法 |
| US10241178B2 (en) * | 2015-05-29 | 2019-03-26 | Case Western Reserve University | System and method for magnetic resonance fingerprinting at high field strengths |
| JP5943159B2 (ja) * | 2016-01-04 | 2016-06-29 | 株式会社日立製作所 | 磁気共鳴イメージング装置 |
| US10359486B2 (en) * | 2016-04-03 | 2019-07-23 | Q Bio, Inc. | Rapid determination of a relaxation time |
| US10495709B2 (en) * | 2016-07-25 | 2019-12-03 | Siemens Healthcare Gmbh | Multi-contrast simultaneous multislice magnetic resonance imaging with binomial radio-frequency pulses |
| WO2018108643A1 (en) * | 2016-12-15 | 2018-06-21 | Koninklijke Philips N.V. | Multi-state magnetic resonance fingerprinting |
| CN108387857B (zh) * | 2017-12-25 | 2020-11-10 | 深圳先进技术研究院 | 一种用于磁共振成像的局部匀场系统及匀场方法 |
| CN108175409B (zh) * | 2018-01-05 | 2021-03-23 | 郜发宝 | 一种定量快速锁频磁共振成像方法 |
| US11373392B2 (en) * | 2018-04-06 | 2022-06-28 | Case Western Reserve University | System and method for magnetic resonance fingerprinting using a plurality of pulse sequence types |
| CN108872897B (zh) * | 2018-04-19 | 2021-05-21 | 上海市东方医院 | 核磁共振t2图像成像方法 |
| CN108937935A (zh) * | 2018-05-25 | 2018-12-07 | 上海东软医疗科技有限公司 | 磁共振成像方法和装置 |
| CN111090069B (zh) * | 2019-11-21 | 2022-03-29 | 深圳先进技术研究院 | 定量磁共振成像参数确定方法、装置、设备及存储介质 |
-
2019
- 2019-11-21 CN CN201911148335.7A patent/CN111090069B/zh active Active
-
2020
- 2020-11-17 WO PCT/CN2020/129571 patent/WO2021098690A1/zh active Application Filing
Also Published As
| Publication number | Publication date |
|---|---|
| CN111090069A (zh) | 2020-05-01 |
| WO2021098690A1 (zh) | 2021-05-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111090069B (zh) | 定量磁共振成像参数确定方法、装置、设备及存储介质 | |
| US10663545B2 (en) | Method and apparatus for low-artifact magnetic resonance fingerprinting scan | |
| CN110807492B (zh) | 一种磁共振多参数同时定量成像方法及系统 | |
| US7439737B2 (en) | Contrast prepared MRI involving non-cartesian trajectories with oversampling of the center of k-space | |
| US20170003365A1 (en) | System and method for magnetic resonance fingerprinting | |
| CN110889897B (zh) | 一种体素内不相干运动磁共振成像参数重建方法及系统 | |
| CN107126211A (zh) | 磁共振成像方法和系统 | |
| CN109613460B (zh) | 磁共振成像方法、装置、计算机设备和存储介质 | |
| CN114167334B (zh) | 磁共振图像的重建方法、装置及电子设备 | |
| US10241173B2 (en) | Systems and methods for designing magnetic resonance fingerprinting imaging parameters | |
| CN111062994B (zh) | 一种化学交换饱和转移对比图像的重建方法及系统 | |
| JPH10502858A (ja) | 対象の磁気共鳴画像化方法及び装置 | |
| CN114217255B (zh) | 一种快速的肝脏多参数定量成像方法 | |
| CN110095742A (zh) | 一种基于神经网络的平面回波成像方法和装置 | |
| WO2016050938A1 (en) | Magnetic resonance imaging with enhanced bone visualization | |
| US8781554B2 (en) | Method and system of simulating magnetic resonance imaging signals | |
| US9772390B2 (en) | Magnetic resonance imaging device and method for generating image using same | |
| CN110992435B (zh) | 图像重建方法及设备、成像数据的处理方法及装置 | |
| WO2009081786A1 (ja) | 磁気共鳴イメージング装置及び磁化率強調画像撮影方法 | |
| JP2023541775A (ja) | ボリュームを網羅する拡散強調磁気共鳴画像のシーケンスを取得し再構成する方法及び装置 | |
| US11921179B2 (en) | Methods and systems for spin-echo train imaging using spiral rings with retraced trajectories | |
| WO2021012972A1 (zh) | 核磁共振成像控制方法、装置及计算机可读存储介质 | |
| Netreba et al. | Correlation reconstructed spine and time relaxation spatial distribution of atomic systems in MRI | |
| Walker-Samuel | Control of a simulated MRI scanner with deep reinforcement learning | |
| CN118135096B (zh) | 基于磁共振指纹超高场T1，T2和T1ρ定量测量方法 |
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 |