CN111077486A - 三维正对比磁共振成像方法、装置、设备和存储介质 - Google Patents

三维正对比磁共振成像方法、装置、设备和存储介质 Download PDF

Info

Publication number
CN111077486A
CN111077486A CN201811224423.6A CN201811224423A CN111077486A CN 111077486 A CN111077486 A CN 111077486A CN 201811224423 A CN201811224423 A CN 201811224423A CN 111077486 A CN111077486 A CN 111077486A
Authority
CN
China
Prior art keywords
magnetic resonance
magnetic field
data
body layer
shift
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.)
Pending
Application number
CN201811224423.6A
Other languages
English (en)
Inventor
史彩云
梁栋
王海峰
谢国喜
陈汉威
程静
刘新
郑海荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Institute of Advanced Technology of CAS
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN201811224423.6A priority Critical patent/CN111077486A/zh
Publication of CN111077486A publication Critical patent/CN111077486A/zh
Pending legal-status Critical Current

Links

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
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (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

本发明实施例公开了一种三维正对比磁共振成像方法、装置、设备和存储介质，该方法包括：基于可变翻转角的SPACE序列，获取包含金属异物的目标对象的每一体层的回波读出梯度所对应的偏移后数据和偏移前数据；求取每一体层的所述偏移后数据与所述偏移前数据所对应的真实相位的相位差；根据每一体层的相位差，确定由所述金属异物的磁场强度在每一体层所形成的局部磁场；根据所述局部磁场求取目标对象的三维正对比磁共振图像。解决了现有基于快速自旋回波序列的三维正对比磁共振成像耗时较长的技术问题，达到了缩短三维正对比磁共振成像方法的技术效果。

Description

三维正对比磁共振成像方法、装置、设备和存储介质
技术领域
本发明实施例涉及图像处理领域，尤其涉及一种三维正对比磁共振成像方法、装置、设备和存储介质。
背景技术
磁共振成像技术(Magnetic Resonance Imaging，简称MRI)是以人体组织内存在的氢质子为成像基础，通过外加一个与氢质子自旋相同频率的射频脉冲信号来产生共振信号。由于磁兼容的介入/植入装置本身是不存在氢质子的，从而在传统MRI图像的介入/植入装置的实际位置区域表现出信号缺失，而且介入/植入装置通常为金属材质，在MRI磁场中会被磁化，产生局部磁场，从而对介入/植入装置周围组织区域产生干扰，即磁化率伪影，因此在介入/植入装置的周围区域表现出一个较大的黑洞，该区域远远大于介入/植入装置本身的大小，即负对比图像。图像的负对比使诊断者很难将前述黑洞与组织空洞区分开来，从而无法进行装置的精准定位和评估。
为了清晰地显示出介入/植入装置在磁共振图像中的位置，通常通过正(亮)对比方式，即通过正对比磁共振图像，显示介入/植入装置的位置。但现有的三维正对比磁共振成像方法是以快速自旋回波(GRE)作为数据采集基础，对相位信息采取去除背景场处理后得到一张局部场图，再结合定量磁化率成像方法重建出磁化率图像，最后正对比的显示出金属粒子位置。但以快速自旋回波去实现三维数据采集，通常需要很长的时间才能完成数据采集，从而限制了定量磁化率成像在临床上的应用，进而限制了三维正对比磁共振成像在临床上的应用。
发明内容
本发明实施例提供了一种三维正对比磁共振成像方法、装置、设备和存储介质，以解决现有技术的三维正对比磁共振成像方法耗时较长的技术问题。
第一方面，本发明实施例提供了一种三维正对比磁共振成像方法，包括：
基于可变翻转角的SPACE序列，获取包含金属异物的目标对象的每一体层的回波读出梯度所对应的偏移后数据和偏移前数据；
求取每一体层的所述偏移后数据与所述偏移前数据所对应的真实相位的相位差；
根据每一体层的相位差，确定由所述金属异物的磁场强度在每一体层所形成的局部磁场；
根据所述局部磁场求取目标对象的三维正对比磁共振图像。
第二方面，本发明实施例还提供了一种三维正对比磁共振成像装置，包括：
数据获取模块，用于基于可变翻转角的SPACE序列，获取包含金属异物的目标对象的每一体层的回波读出梯度所对应的偏移后数据和偏移前数据；
相位差确定模块，用于求取每一体层的所述偏移后数据与所述偏移前数据所对应的真实相位的相位差；
局部磁场确定模块，用于根据每一体层的相位差，确定由所述金属异物的磁场强度在每一体层所形成的局部磁场；
正对比磁共振图像模块，用于根据所述局部磁场求取目标对象的三维正对比磁共振图像。
第三方面，本发明实施例还提供了一种计算机设备，所述计算机设备包括：
一个或多个处理器；
存储装置，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如第一方面所述的三维正对比磁共振成像方法。
第四方面，本发明实施例还提供了一种包含计算机可执行指令的存储介质，所述计算机可执行指令在由计算机处理器执行时用于执行如第一方面所述的三维正对比磁共振成像方法。
本实施例提供的三维正对比磁共振成像方法的技术方案，基于可变翻转角的SPACE序列，获取包含金属异物的目标对象的每一体层的回波读出梯度所对应的偏移后数据和偏移前数据；求取每一体层的偏移后数据与偏移前数据所对应的真实相位的相位差；根据每一体层的相位差，确定由金属异物的磁场强度在每一体层所形成的局部磁场；根据局部磁场求取目标对象的三维正对比磁共振图像。通过可变翻转角的SPACE序列可使组织信号在回波链的大部分时间内保持稳态，且回波链长度达到100-200，极大地提高了T2加权图像的数据采集效率，通过使回波读出梯度的相位发生偏移，并通过回波读出梯度相位偏移前数据和相位偏移后数据来确定二者的相位差，然后根据相位差确定金属对周围组织产生的局部磁场图，进而确定三维正对比磁共振图像，由于基于SPACE的三维数据采集效率的大幅度提高，大大缩短了三维正对比磁共振成像的时间，有利于三维正对比磁共振成像的临床推广。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图做一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例一提供的三维正对比磁共振成像方法的流程图；
图2是本发明实施例一提供的SPACE序列示意图；
图3A是本发明实施例一提供的基于SPACE序列采集的偏移前数据重建的幅值图；
图3B是本发明实施例一提供的三维正对比磁共振图像；
图4是本发明实施例二提供的三维正对比磁共振成像装置的结构框图；
图5是本发明实施例三提供的计算机设备的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，以下将参照本发明实施例中的附图，通过实施方式清楚、完整地描述本发明的技术方案，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例一
图1是本发明实施例一提供的三维正对比磁共振成像方法的流程图。本实施例的技术方案适用于定位植入目标对象体内的金属异物位置的情况。该方法可以由本发明实施例提供的三维正对比磁共振成像装置来执行，该装置可以采用软件和/或硬件的方式实现，并配置在处理器中应用。该方法具体包括如下步骤：
S101、基于可变翻转角的SPACE序列，获取包含金属异物的目标对象的每一体层的回波读出梯度所对应的偏移后数据和偏移前数据。
SPACE序列是在快速自旋回波序列的基础上发展而来的，而快速自旋回波序列是在自旋回波序列的基础上发展起来的，传统的快速自旋回波序列的翻转角通常为一常数，比如180度，组织信号以指数衰减，且衰减速度较快，而当回波链到达20左右时，其回波信号基本不能用于成像。但采用可变翻转角后，选择合适的翻转角可使组织信号在回波链的大部分时间内都保持稳态，回波链长度可以达到100～200，极大地提高了T2加权图像的采集效率。将可变翻转角的快速自旋回波序列与选层激发脉冲结合得到可变翻转角的SPACE序列，该序列使得FSE可以精确地选层，并能应用于不同的身体部位，基于此，本实施例采用SPACE序列作为磁共振测量数据的采集序列。
如图2所示，由SPACE的时序图可以看出：射频脉冲的第一个脉冲为90度的选层激发脉冲，第二个脉冲为180度的回聚脉冲，随后是可变翻转角回聚脉冲。其中，180度和可变翻转角均采用非选择性回聚脉冲，以缩短回波信号的间距，从而提高磁共振数据的采集效率。
为了适用90度选层激发脉冲较长的持续时间，180度回聚脉冲的回波间隔(EchoSpacing，简称ESP)被延长，称为ESP1，其后可变翻转角的回波间隔称为EP2，EP2越短越好，并由梯度***、分辨率和带宽所决定。在180度回聚脉冲之后，经过ESP1/2时间在横向平面上产生了完全回波，可变翻转角延续了完全回波在横向平面上的时间，使回波链延长，提高了磁共振数据的采集效率。
需要说明的是，由于可变翻转角之间的脉冲回波间隔是ESP2，所以第一个可变翻转角与180度脉冲间隔的时间是ESP1/2+ESP2/2。当然，SPACE序列的可变翻转角还需要满足CPMG(Carr-Purcell-Meiboom-Gill)的条件。
从图2还可看出，本实施例的SPACE序列还对回波读出梯度的位置进行了偏移，该偏移时间Tshift很短，优选但不限于0.3～0.8ms，从而使得有效回波时间很短。为了消除由磁共振***本身磁化率的差异所造成的相位变换，本实施例通过该SPACE序列采集包含金属异物的目标对象的每一体层的回波读出梯度所对应的偏移后数据和偏移前数据。
S102、求取每一体层的偏移后数据与偏移前数据所对应的真实相位的相位差。
由于目标对象的金属异物会对目标对象原有组织的磁场信息产生影响，使其发生变化，而且磁共振测量数据的相位信息包含了丰富的组织磁化率变化信息，因此可以采用包含相位信息的相位图像来确定金属异物的位置，基于此，本实施例需要先确定每一体层的偏移后数据与偏移前数据所对应的真实相位的相位差，以得到由组织自身磁化率差异(局部磁场)引起的相位变化。
在确定每一体层的偏移后数据与偏移前数据所对应的真实相位时，考虑到磁共振***所产生的相位值被限制在(-π，π)范围内，该相位值与真实的相位相差2kπ,因此实际得到的相位值与真实的相位值存在着不同程度的相位周期性模糊。根据缠绕相位，也就是实际得到的相位值，得到真实相位的处理过程叫做相位解缠绕。
本实施例先对每一体层的偏移后数据和偏移前数据在空间上进行解缠绕，以得到每一体层的偏移后数据的真实相位和偏移前数据的真实相位；将每一体层的偏移后数据的真实相位与偏移前数据的真实相位的差值作为每一体层的相位差，以得到由组织本身磁化率差异(局部磁场)而产生的相位变化。
S103、根据每一体层的相位差，确定由金属异物的磁场强度在每一体层所形成的局部磁场。
其中，组织本身磁化率差异，即局部磁场表示为：ΔB＝Δθ/γB0Tshift。其中，ΔB表示金属异物在每一体层所形成的局部磁场，也就是局部场图，Δθ表示相位差，γ表示旋磁比，为常数；B0表示磁共振成像的主磁场，Tshift表示回波读出梯度的偏移时间。
S104、根据局部磁场求取目标对象的三维正对比磁共振图像。
由于正对比磁共振图像是基于磁化率图像生成的，因此本实施例先通过定量磁敏感成像方法重建出磁化率图像，然后基于磁化率图像确定正对比度磁共振图像。
为了重建出磁化率图像，本实施例先确定磁化率与磁场在傅里叶域的关系表达式，具体为：
Figure BDA0001835582810000071
其中，k为傅氏域坐标，
Figure BDA0001835582810000072
kx，ky和kz分别是k在x、y和z轴上的映射，B(k)是傅里叶域的磁场，由回波时间内的相位偏移决定，反映了由磁化率分布决定的相位分布，属于相位的全局特征，X(k)是傅里叶域的磁化率分布。
由前述公式可知，要求磁化率分布X(k)，需要避免
Figure BDA0001835582810000081
因此用通常得到的磁场来求解任意磁化率分布式是不适定的。为此，本实施例引入L1范数约束来求解该不适定问题的最优解，同时，考虑到金属异物的磁化率远大于人体组织的磁化率的特点，增加两个约束矩阵，即加权矩阵W和掩膜矩阵M以提高图像重建的质量。重建公式为：
其中，λ为正则化参数，以保证数据的一致性和稀疏性，χ为磁化率，ΔB表示目标区域内金属异物对周围组织产生的局部磁场，即局部图，Δθ表示相位差，γ表示旋磁比，且为常数，B0表示磁共振成像的主磁场，Tshift表示回波读出梯度偏移时间，W为加权矩阵，M为掩膜矩阵，C为由离散化卷积核得到的极化核矩阵，G表示在三维空间中三个方向的梯度算子。
为了提高f(χ,λ)的求解速度和准确率，本实施例利用原始对偶算法迭代求解f(χ,λ)，该算法用来解决如下原始最小化问题：
Figure BDA0001835582810000083
该原始最小化问题所对应的对偶问题为：
Figure BDA0001835582810000084
利用勒让德变换，得到
Figure BDA0001835582810000085
带入公式前述原始最小化问题表达式，于是最小化问题转变为鞍点问题
Figure BDA0001835582810000086
然后采用近似映射(Proximal Mapping)分别求解原始问题和对偶问题，以确定矩阵C和矩阵G。
矩阵C和矩阵G确定后，采用一阶原始对偶算法，比如Chambolle-Pock一阶原始对偶算法进行磁化率图像的重建，具体过程如下：
1、
Figure BDA0001835582810000091
2、χ0，p0，q0均初始化为0。
3、
Figure BDA0001835582810000094
4、按照如下公式顺序迭代，直至n≥N。
Figure BDA0001835582810000092
Figure BDA0001835582810000093
χn+1＝χn-τCTWTpn+1+τGTMTqn+1
Figure BDA0001835582810000095
n＝n+1
上述迭代完成后，得到磁化率图像，该磁化率图像便是最终的三维正对比磁共振图像。由于一阶原始对偶算法具有较高的运算效率，因此可以提高磁化率图像的重建速度，进而提高三维正对比磁共振图像的重建速度，便于临床推广。需要说明的是，本实施例采用现有技术将磁化率图像转换成三维正对比磁共振图像即可。
为了更好地说明本实施例所述的三维正对比磁共振成像方法所成图像具有较高的图像质量，本实施例还给出了一模拟实验以及实验结果，具体为：将由钛金属制成的直径为2mm的穿刺针置于1mg/ml的硫酸铜溶液的实验瓶中，然后基于SPACE序列采集该实验瓶的梯度回波偏移前数据和偏移后数据。通过现有技术的磁共振成像方法对偏移前数据进行磁共振图像重建，以得到基于幅值信息的磁共振图像，如图3A所示。通过本实施例所述的三维正对比磁共振成像方法对偏移前数据和偏移后数据进行图像重建，以得到基于相位信息的三维正对比磁共振图像，如图3B所示。
对比图3A和图3B可以看出，图3A的穿刺针产生了很大的伪影(黑洞)，具***置参见箭头位置；而图3B的正对比磁共振图像显示出了穿刺针的具***置，因此，基于本实施例所述的三维正对比磁共振成像方法重建的正对比磁共振图像，可以清晰地显示出目标对象的金属异物的准确位置，有利于临床推广，以及提高临床诊断的准确性。
本实施例提供的三维正对比磁共振成像方法的技术方案，基于可变翻转角的SPACE序列，获取包含金属异物的目标对象的每一体层的回波读出梯度所对应的偏移后数据和偏移前数据；求取每一体层的偏移后数据与偏移前数据所对应的真实相位的相位差；根据每一体层的相位差，确定由金属异物的磁场强度在每一体层所形成的局部磁场；根据局部磁场求取目标对象的三维正对比磁共振图像。通过可变翻转角的SPACE序列可使组织信号在回波链的大部分时间内保持稳态，且回波链长度达到100-200，极大地提高了T2加权图像的数据采集效率，通过使回波读出梯度的相位发生偏移，并通过回波读出梯度相位偏移前数据和相位偏移后数据来确定二者的相位差，然后根据相位差确定局部磁场，进而确定三维正对比磁共振图像，由于数据采集效率的大幅度提高，大大缩短了三维正对比磁共振成像的时间，有利于三维正对比磁共振成像的临床推广。
实施例二
图4是本发明实施例二提供的三维正对比磁共振成像装置的结构框图。该装置用于执行上述任意实施例所提供的三维正对比磁共振成像方法，该装置可选为软件或硬件实现。该装置包括：
数据获取模块11，用于基于可变翻转角的SPACE序列，获取包含金属异物的目标对象的每一体层的回波读出梯度所对应的偏移后数据和偏移前数据；
相位差确定模块12，用于求取每一体层的所述偏移后数据与所述偏移前数据所对应的真实相位的相位差；
局部磁场确定模块13，用于根据每一体层的相位差，确定由所述金属异物的磁场强度在每一体层所形成的局部磁场；
正对比磁共振图像模块14，用于根据所述局部磁场求取目标对象的三维正对比磁共振图像。
本实施例提供的三维正对比磁共振成像装置的技术方案，通过数据获取模块，基于可变翻转角的SPACE序列，获取包含金属异物的目标对象的每一体层的回波读出梯度所对应的偏移后数据和偏移前数据；通过相位差确定模块求取每一体层的偏移后数据与偏移前数据所对应的真实相位的相位差；通过局部磁场确定模块根据每一体层的相位差，确定由金属异物的磁场强度在每一体层所形成的局部磁场；通过正对比磁共振图像模块根据局部磁场求取目标对象的三维正对比磁共振图像。通过可变翻转角的SPACE序列可使组织信号在回波链的大部分时间内保持稳态，且回波链长度达到100-200，极大地提高了T2加权图像的数据采集效率，通过使回波读出梯度的相位发生偏移，并通过回波读出梯度相位偏移前数据和相位偏移后数据来确定二者的相位差，然后根据相位差确定局部磁场，进而确定正对比磁共振图像，由于数据采集效率的大幅度提高，大大缩短了三维正对比磁共振成像的时间，有利于三维正对比磁共振成像的临床推广。
本发明实施例所提供的三维正对比磁共振成像装置可执行本发明任意实施例所提供的三维正对比磁共振成像方法，具备执行方法相应的功能模块和有益效果。
实施例三
图5为本发明实施例三提供的计算机设备的结构框图，如图5所示，该计算机设备包括处理器201、存储器202、输入装置203以及输出装置204；计算机设备中处理器201的数量可以是一个或多个，图5中以一个处理器201为例；设备中的处理器201、存储器202、输入装置203以及输出装置204可以通过总线或其他方式连接，图5中以通过总线连接为例。
存储器202作为一种计算机可读存储介质，可用于存储软件程序、计算机可执行程序以及模块，如本发明实施例中的三维正对比磁共振成像方法对应的程序指令/模块(例如，数据获取模块11、相位差确定模块12、局部磁场确定模块13以及正对比磁共振图像模块14)。处理器201通过运行存储在存储器202中的软件程序、指令以及模块，从而执行设备的各种功能应用以及数据处理，即实现上述的三维正对比磁共振成像方法。
存储器202可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作***、至少一个功能所需的应用程序；存储数据区可存储根据终端的使用所创建的数据等。此外，存储器202可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中，存储器202可进一步包括相对于处理器201远程设置的存储器，这些远程存储器可以通过网络连接至设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
输入装置203可用于接收输入的数字或字符信息，以及产生与设备的用户设置以及功能控制有关的键信号输入。
输出装置204可包括显示屏等显示设备，例如，用户终端的显示屏。
实施例四
本发明实施例四还提供一种包含计算机可执行指令的存储介质，所述计算机可执行指令在由计算机处理器执行时用于执行一种三维正对比磁共振成像方法，该方法包括：
基于可变翻转角的SPACE序列，获取包含金属异物的目标对象的每一体层的回波读出梯度所对应的偏移后数据和偏移前数据；
求取每一体层的所述偏移后数据与所述偏移前数据所对应的真实相位的相位差；
根据每一体层的相位差，确定由所述金属异物的磁场强度在每一体层所形成的局部磁场；
根据所述局部磁场求取目标对象的三维正对比磁共振图像。
当然，本发明实施例所提供的一种包含计算机可执行指令的存储介质，其计算机可执行指令不限于如上所述的方法操作，还可以执行本发明任意实施例所提供的三维正对比磁共振成像方法中的相关操作。
通过以上关于实施方式的描述，所属领域的技术人员可以清楚地了解到，本发明可借助软件及必需的通用硬件来实现，当然也可以通过硬件实现，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如计算机的软盘、只读存储器(Read-Only Memory，简称ROM)、随机存取存储器(RandomAccess Memory，简称RAM)、闪存(FLASH)、硬盘或光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的三维正对比磁共振成像方法。
值得注意的是，上述三维正对比磁共振成像装置的实施例中，所包括的各个单元和模块只是按照功能逻辑进行划分的，但并不局限于上述的划分，只要能够实现相应的功能即可；另外，各功能单元的具体名称也只是为了便于相互区分，并不用于限制本发明的保护范围。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。

Claims (10)

1.一种三维正对比磁共振成像方法，其特征在于，包括：
基于可变翻转角的SPACE序列，获取包含金属异物的目标对象的每一体层的回波读出梯度所对应的偏移后数据和偏移前数据；
求取每一体层的所述偏移后数据与所述偏移前数据所对应的真实相位的相位差；
根据每一体层的相位差，确定由所述金属异物的磁场强度在每一体层所形成的局部磁场；
根据所述局部磁场求取目标对象的三维正对比磁共振图像。
2.根据权利要求1所述的方法，其特征在于，所述求取每一体层的所述偏移后数据与所述偏移前数据所对应的真实相位的相位差，包括：
对每一体层的所述偏移后数据和所述偏移前数据在空间上进行解缠绕，以得到每一体层的所述偏移后数据的真实相位和所述偏移前数据的真实相位；
将每一体层的所述偏移后数据的真实相位与所述偏移前数据的真实相位的差值作为每一体层的相位差。
3.根据权利要求1所述的方法，其特征在于，所述局部磁场的表达式为B＝Δθ/γB0Tshift，其中，ΔB表示金属异物在每一体层所形成的局部磁场，Δθ表示相位差，γ表示旋磁比，为常数；B0表示磁共振成像的主磁场，Tshift表示回波读出梯度的偏移时间。
4.根据权利要求3所述的方法，其特征在于，所述根据所述局部磁场求取目标对象的三维正对比磁共振图像，包括：
根据所述局部磁场，通过定量磁敏感成像方法进行磁共振图像重建，以得到三维正对比度磁共振图像。
5.根据权利要求4所述的方法，其特征在于，所述根据所述局部磁场，通过定量磁敏感成像方法进行磁共振图像重建，以得到三维正对比度磁共振图像，包括：
确定磁化率与磁场在傅里叶域的关系表达式；
利用局部磁场对所述关系表达式进行正则化约束重建以获得目标对象的磁化率图像；
根据所述磁化率图像确定三维正对比磁共振图像。
6.根据权利要求5所述的方法，其特征在于，所述利用局部磁场对所述关系表达式进行正则化约束重建以获得目标对象的磁化率图像，包括：
利用局部磁场对所述关系表达式进行正则化约束重建的表达式为：
Figure FDA0001835582800000021
其中，λ为正则化参数，χ为磁化率，ΔB表示目标区域内金属异物对周围组织产生的局部磁场，即局部场图，θ表示相位差，γ表示旋磁比，且为常数，B0表示磁共振成像的主磁场，Tshift表示回波读出梯度偏移时间，W为加权矩阵，M为掩膜矩阵，C为由离散化卷积核得到的极化核矩阵，G表示在三维空间中三个方向的梯度算子；
通过一阶原始对偶算法求解f(χ,λ)以得到磁化率图像。
7.根据权利要求1所述的方法，其特征在于，所述回波读出梯度的偏移时间Tshift的范围为0.3-0.8ms。
8.一种三维正对比磁共振成像装置，其特征在于，包括：
数据获取模块，用基于可变翻转角的SPACE序列，获取包含金属异物的目标对象的每一体层的回波读出梯度所对应的偏移后数据和偏移前数据；
相位差确定模块，用于求取每一体层的所述偏移后数据与所述偏移前数据所对应的真实相位的相位差；
局部磁场确定模块，用于根据每一体层的相位差，确定由所述金属异物的磁场强度在每一体层所形成的局部磁场；
正对比磁共振图像模块，用于根据所述局部磁场求取目标对象的三维正对比磁共振图像。
9.一种计算机设备，其特征在于，所述计算机设备包括：
一个或多个处理器；
存储装置，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如权利要求1-7任一所述的三维正对比磁共振成像方法。
10.一种包含计算机可执行指令的存储介质，其特征在于，所述计算机可执行指令在由计算机处理器执行时用于执行如权利要求1-7任一所述的三维正对比磁共振成像方法。
CN201811224423.6A 2018-10-19 2018-10-19 三维正对比磁共振成像方法、装置、设备和存储介质 Pending CN111077486A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811224423.6A CN111077486A (zh) 2018-10-19 2018-10-19 三维正对比磁共振成像方法、装置、设备和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811224423.6A CN111077486A (zh) 2018-10-19 2018-10-19 三维正对比磁共振成像方法、装置、设备和存储介质

Publications (1)

Publication Number Publication Date
CN111077486A true CN111077486A (zh) 2020-04-28

Family

ID=70309451

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811224423.6A Pending CN111077486A (zh) 2018-10-19 2018-10-19 三维正对比磁共振成像方法、装置、设备和存储介质

Country Status (1)

Country Link
CN (1) CN111077486A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113866694A (zh) * 2021-09-26 2021-12-31 上海交通大学 一种快速三维磁共振t1定量成像方法、***及介质
WO2024139909A1 (zh) * 2022-12-26 2024-07-04 中国科学院深圳先进技术研究院 基于深度学习的正对比磁共振三维成像方法及装置

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102435966A (zh) * 2011-09-02 2012-05-02 中国科学院深圳先进技术研究院 三维磁共振成像方法及***
CN102631196A (zh) * 2011-12-08 2012-08-15 中国科学院深圳先进技术研究院 磁共振弹性成像三维可视化方法及***
CN102707251A (zh) * 2011-12-12 2012-10-03 中国科学院深圳先进技术研究院 计算space序列信号的方法和***及主动脉信号的采集方法
CN102727206A (zh) * 2011-12-08 2012-10-17 中国科学院深圳先进技术研究院 颅脑、颈动脉和主动脉三位一体的扫描方法及扫描***
CN102914752A (zh) * 2011-08-01 2013-02-06 西门子公司 在预定区域内采集磁共振数据的方法和相应的磁共振设备
CN103006216A (zh) * 2012-12-17 2013-04-03 中国科学院深圳先进技术研究院 磁共振弹性成像重建方法及成像***
CN103076582A (zh) * 2011-10-25 2013-05-01 西门子公司 采集三维体积片段中的磁共振数据的方法以及磁共振装置
CN103576114A (zh) * 2012-07-24 2014-02-12 西门子公司 获取磁共振数据的方法、磁共振设备、程序和数据载体
US8886283B1 (en) * 2011-06-21 2014-11-11 Stc.Unm 3D and 4D magnetic susceptibility tomography based on complex MR images
CN105249967A (zh) * 2015-11-18 2016-01-20 深圳先进技术研究院 一种正对比磁共振成像方法及装置
CN106264530A (zh) * 2016-11-04 2017-01-04 深圳先进技术研究院 用于三维动态磁共振成像的笛卡尔k空间采集方法和***
CN107205718A (zh) * 2015-12-30 2017-09-26 深圳先进技术研究院 一种聚焦超声位移成像方法及装置

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8886283B1 (en) * 2011-06-21 2014-11-11 Stc.Unm 3D and 4D magnetic susceptibility tomography based on complex MR images
CN102914752A (zh) * 2011-08-01 2013-02-06 西门子公司 在预定区域内采集磁共振数据的方法和相应的磁共振设备
CN102435966A (zh) * 2011-09-02 2012-05-02 中国科学院深圳先进技术研究院 三维磁共振成像方法及***
CN103076582A (zh) * 2011-10-25 2013-05-01 西门子公司 采集三维体积片段中的磁共振数据的方法以及磁共振装置
CN102631196A (zh) * 2011-12-08 2012-08-15 中国科学院深圳先进技术研究院 磁共振弹性成像三维可视化方法及***
CN102727206A (zh) * 2011-12-08 2012-10-17 中国科学院深圳先进技术研究院 颅脑、颈动脉和主动脉三位一体的扫描方法及扫描***
CN102707251A (zh) * 2011-12-12 2012-10-03 中国科学院深圳先进技术研究院 计算space序列信号的方法和***及主动脉信号的采集方法
CN103576114A (zh) * 2012-07-24 2014-02-12 西门子公司 获取磁共振数据的方法、磁共振设备、程序和数据载体
CN103006216A (zh) * 2012-12-17 2013-04-03 中国科学院深圳先进技术研究院 磁共振弹性成像重建方法及成像***
CN105249967A (zh) * 2015-11-18 2016-01-20 深圳先进技术研究院 一种正对比磁共振成像方法及装置
CN107205718A (zh) * 2015-12-30 2017-09-26 深圳先进技术研究院 一种聚焦超声位移成像方法及装置
CN106264530A (zh) * 2016-11-04 2017-01-04 深圳先进技术研究院 用于三维动态磁共振成像的笛卡尔k空间采集方法和***

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈敏等: "基于磁敏感的磁兼容金属植入物的正对比磁共振成像", 《集成技术》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113866694A (zh) * 2021-09-26 2021-12-31 上海交通大学 一种快速三维磁共振t1定量成像方法、***及介质
CN113866694B (zh) * 2021-09-26 2022-12-09 上海交通大学 一种快速三维磁共振t1定量成像方法、***及介质
WO2024139909A1 (zh) * 2022-12-26 2024-07-04 中国科学院深圳先进技术研究院 基于深度学习的正对比磁共振三维成像方法及装置

Similar Documents

Publication Publication Date Title
Cole et al. Analysis of deep complex‐valued convolutional neural networks for MRI reconstruction and phase‐focused applications
CN108957375B (zh) 磁共振成像方法、装置、设备及存储介质
Buonincontri et al. MR fingerprinting with simultaneous B1 estimation
Doneva et al. Matrix completion-based reconstruction for undersampled magnetic resonance fingerprinting data
US10768260B2 (en) System and method for controlling noise in magnetic resonance imaging using a local low rank technique
Ben‐Eliezer et al. Rapid and accurate T2 mapping from multi–spin‐echo data using Bloch‐simulation‐based reconstruction
US10175331B2 (en) Propeller with Dixon water fat separation
Hamilton et al. Deep learning reconstruction for cardiac magnetic resonance fingerprinting T1 and T2 mapping
US8379951B2 (en) Auto calibration parallel imaging reconstruction method from arbitrary k-space sampling
US10895622B2 (en) Noise suppression for wave-CAIPI
US11181598B2 (en) Multi-contrast MRI image reconstruction using machine learning
Cole et al. Analysis of deep complex-valued convolutional neural networks for MRI reconstruction
CN113994225A (zh) 用于磁共振成像中的图像重建的***和方法
van der Heide et al. High‐resolution in vivo MR‐STAT using a matrix‐free and parallelized reconstruction algorithm
Maier et al. PyQMRI: an accelerated Python based quantitative MRI toolbox
Barbieri et al. Circumventing the curse of dimensionality in magnetic resonance fingerprinting through a deep learning approach
CN111077486A (zh) 三维正对比磁共振成像方法、装置、设备和存储介质
Benjamin et al. Multi-shot Echo Planar Imaging for accelerated Cartesian MR Fingerprinting: an alternative to conventional spiral MR Fingerprinting
Tu et al. K-space and image domain collaborative energy-based model for parallel MRI reconstruction
Gan et al. SS-JIRCS: Self-supervised joint image reconstruction and coil sensitivity calibration in parallel MRI without ground truth
Licht et al. Multidimensional compressed sensing to advance 23Na multi‐quantum coherences MRI
WO2021129235A1 (zh) 三维磁共振快速参数成像方法和装置
Hu et al. SPICE: Self-supervised learning for MRI with automatic coil sensitivity estimation
CN114469047A (zh) 基于可变翻转角的快速正对比磁共振成像方法和装置
US20200341089A1 (en) System and method for improved magnetic resonance fingerprinting using inner product space

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20200428

RJ01 Rejection of invention patent application after publication