CN110174632B - 磁共振成像方法、装置、成像设备及磁共振成像*** - Google Patents
磁共振成像方法、装置、成像设备及磁共振成像*** Download PDFInfo
- Publication number
- CN110174632B CN110174632B CN201910496273.2A CN201910496273A CN110174632B CN 110174632 B CN110174632 B CN 110174632B CN 201910496273 A CN201910496273 A CN 201910496273A CN 110174632 B CN110174632 B CN 110174632B
- Authority
- CN
- China
- Prior art keywords
- echo
- echoes
- magnetic
- magnetic resonance
- phases
- 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
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000002595 magnetic resonance imaging Methods 0.000 title claims abstract description 41
- 238000003384 imaging method Methods 0.000 title claims abstract description 35
- 238000002592 echocardiography Methods 0.000 claims abstract description 77
- 230000035945 sensitivity Effects 0.000 claims abstract description 34
- 238000001514 detection method Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 238000007781 pre-processing Methods 0.000 claims description 4
- 230000014509 gene expression Effects 0.000 claims description 3
- 238000011282 treatment Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 9
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 5
- 238000003745 diagnosis Methods 0.000 description 4
- 208000032843 Hemorrhage Diseases 0.000 description 3
- 230000000740 bleeding effect Effects 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000000306 component Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 208000021245 head disease Diseases 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 206010008111 Cerebral haemorrhage Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000012503 blood component Substances 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000002308 calcification Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000005786 degenerative changes Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 210000000653 nervous system Anatomy 0.000 description 1
- 230000004770 neurodegeneration Effects 0.000 description 1
- 208000015122 neurodegenerative disease Diseases 0.000 description 1
- 230000000626 neurodegenerative effect Effects 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 238000010587 phase diagram Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 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
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/5602—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by filtering or weighting based on different relaxation times within the sample, e.g. T1 weighting using an inversion pulse
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
本公开提出一种磁共振成像方法、装置、成像设备及磁共振成像***,该方法包括:采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位的信息;利用所述N个回波的幅度和相位的信息,获得不同类型的磁共振图像,其中,获得所述磁敏感加权图像包括:对所述N个回波的相位进行拟合,获得每个回波的有效相位;基于所述有效相位对所述N个回波的幅度进行加权,生成磁敏感加权图像。通过本实施例,在获得磁敏感加权图像时,通过对所采集的N个回波的相位进行拟合,获得每个回波的具有高信噪比的有效相位,并基于该有效相位对N个回波的幅度进行加权,提高了图像的信噪比,有效地降低强磁敏感的组织带来的磁化率伪影。
Description
技术领域
本申请涉及磁共振成像领域,尤其涉及一种磁共振成像方法、装置、成像设备及磁共振成像***。
背景技术
磁共振成像(Magnetic Resonance Imaging,MRI)作为一种多参数、多对比度的成像技术,是现代医疗影像学中主要的成像方式之一。MRI***的基本工作原理是根据磁共振现象,采用射频发射线圈激励受检体中的氢质子,并运用梯度场进行空间编码,随后采用射频接收线圈接收带位置信息的电磁信号,最终利用傅里叶变换重建出被检体的图像信息。
磁敏感加权成像(Susceptibility weighting imaging,SWI)是MRI中通过组合相位信息和幅度信息来增强图像中呈现的不同组织的对比度的成像方法。在传统的SWI中,使用单个梯度回波序列,并且以长回波时间(TE)采集单回波信号。然而,对于一些具有强磁敏感的组织,例如空腔或出血部位,可能产生高空间频率的卷绕/缠绕相位,可能导致图像伪影和成像对比失真。
发明内容
为克服相关技术中存在的问题,本公开提供了一种磁共振成像方法、装置、成像设备及磁共振成像***。
具体地,本公开是通过如下技术方案实现的:
第一方面,提供一种磁共振成像方法,所述方法包括:
采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位的信息,其中,所述N为大于1的自然数;
利用所述N个回波的幅度和相位的信息,获得不同类型的磁共振图像,所述不同类型的磁共振图像包括至少一个下述图像:磁敏感加权图像、磁化率空间分布定量图、增强磁敏感加权图像、T2*定量拓扑图,
其中,获得所述磁敏感加权图像包括:
对所述N个回波的相位进行拟合,获得每个回波的有效相位;
基于所述有效相位对所述N个回波的幅度进行加权,生成磁敏感加权图像。
第二方面,提供一种磁共振成像装置,所述装置包括:
采集单元,用于采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位,其中,所述N为大于1的自然数;
获得单元,用于利用所述N个回波的幅度信息和相位信息,获得不同类型的磁共振图像,所述不同类型的磁共振图像包括至少一个下述图像:磁敏感加权图像、磁化率空间分布定量图、增强磁敏感加权图像、T2*定量拓扑图;
所述获得单元包括:
拟合子单元,用于对所述N个回波的相位进行拟合,获得每个回波的有效相位;
生成子单元,用于基于所述有效相位对所述N个回波的幅度进行加权,生成磁敏感加权图像。
第三方面,提供一种成像设备,所述设备包括:内部总线,以及通过内部总线连接的存储器、处理器和外部接口;其中,
所述外部接口,用于获得回波信号;
所述存储器,用于存储磁共振成像对应的机器可读指令;
所述处理器,用于读取所述存储器上的所述机器可读指令,并执行所述指令以实现如下操作:
采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位,其中,所述N为大于1的自然数;
利用所述N个回波的幅度信息和相位信息,获得不同类型的磁共振图像,所述不同类型的磁共振图像包括至少一个下述图像:磁敏感加权图像、磁化率空间分布定量图、增强磁敏感加权图像、T2*定量拓扑图,
其中,获得所述磁敏感加权图像包括:
对所述N个回波的相位进行拟合,获得每个回波的有效相位;
基于所述有效相位对所述N个回波的幅度进行加权,生成磁敏感加权图像。
第四方面,提供一种磁共振成像***,所述***包括:磁场发生器、信号探测器、扫描床和成像设备;其中,
所述磁场发生器,用于产生均匀磁场;
所述信号探测器,用于向探测区域发送信号和从所述探测区域接收信号并成生成像数据;
所述扫描床,用于承载待检测对象;
所述成像设备,用于:
采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位,其中,所述N为大于1的自然数;
利用所述N个回波的幅度信息和相位信息,获得不同类型的磁共振图像,所述不同类型的磁共振图像包括至少一个下述图像:磁敏感加权图像、磁化率空间分布定量图、增强磁敏感加权图像、T2*定量拓扑图,
其中,获得所述磁敏感加权图像包括:
对所述N个回波的相位进行拟合,获得每个回波的有效相位;
基于所述有效相位对所述N个回波的幅度进行加权,生成磁敏感加权图像。
本公开的实施例提供的技术方案可以包括以下有益效果:
通过利用所采集的N个回波的幅度信息和相位信息,获得不同类型的磁共振图像。在获得磁敏感加权图像时,通过对所采集的N个回波的相位进行拟合,获得每个回波的具有高信噪比的有效相位,并基于该有效相位对N个回波的幅度进行加权,提高了图像的信噪比,有效地降低强磁敏感的组织,例如空腔或出血部位,带来的磁化率伪影。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本说明书的实施例,并与说明书一起用于解释本说明书的原理。
图1示出MRI***的应用场景示意图;
图2示出本公开根据一示例性实施例示出的一种磁共振成像方法的流程图;
图3示出本公开根据一示例性实施例示出的回波示意图;
图4示出本公开根据一示例性实施例示出的另一种磁共振成像方法的流程图;
图5示出本公开根据一示例性实施例示出的另一种磁共振成像方法的流程图;
图6示出本公开根据一示例性实施例示出的另一种磁共振成像方法的流程图;
图7A和图7B示出通过本公开的磁共振成像方法进行成像的测试例;
图8示出本公开根据一示例性实施例示出的一种磁共振成像装置的示意图;
图9示出本公开根据一示例性实施例示出的一种成像设备的结构示意图;
图10示出本公开根据一示例性实施例示出的一种磁共振成像***的结构示意图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。
为了更好的理解本申请,首先介绍MRI***。如图1所示,简单示意了MRI***的组成,主要包括检查床110、磁体120、梯度线圈131-133、射频线圈140、主计算机150、梯度放大器160、射频控制器170和控制台180。磁体120是产生主磁场的装置。梯度线圈131-133、梯度放大器160等组成梯度***,主要用于产生梯度磁场以能够进行磁共振信号的空间定位编码。其中,梯度线圈由三组独立的线圈构成,分别为X轴梯度线圈131,Y轴梯度线圈132和Z轴梯度线圈133。射频线圈140、射频控制器170等组成射频***,主要用于向被检体发射射频信号,再从被检体接收射频信号,因此射频线圈140有射频发射线圈和射频接收线圈之分,射频发射线圈和射频接收线圈可以为同一个线圈,也可以分为不同的线圈。为了简单起见,图1中仅示出一个线圈140。主计算机150负责MR成像序列的发送、采集的射频数据的运算、磁共振图像重建和显示等。
磁共振成像基于不同对比度性质生成的定性图像包括T1加权,T2加权,质子密度加权,弥散加权,磁敏感加权等等,这些图像都以幅度图像的形式展现。磁共振成像是一种多参数的成像方式,能提供的远远不止这些定性的幅度信息,磁共振数据中还包含丰富的相位信息和定量信息。尤其在磁敏感加权成像方面,相位信息和定量信息对于疾病诊断尤其是在脑神经科学研究以及神经退行性疾病临床应用方面也很重要。磁敏感加权成像(Susceptibility weighting imaging,SWI)和增强磁敏感加权成像(enhancedSusceptibility weighting imaging eSWI)对于显示静脉血管、血液成份等非常敏感;T2*和R2*定量信息图以及相位信息输出图对组织含铁元素变化等导致磁化率变化的因素非常敏感,在铁质沉积以及多种神经退行性疾病中有很好的应用价值;磁化率空间分布定量图(Quantitative Susceptibility Mapping,QSM)是近年来比较新兴的技术,可用来定量组织的铁含量,并且可以帮助医生快速地区分脑出血(高信号)和钙化(低信号)。
为了能够综合运用这几个定量信息,以对组织病理尤其是神经***退行性病变进行精准诊断,在本申请实施例中,根据多个回波的幅度信息和相位信息获得不同类型的磁共振图像。
参见图2,为本申请磁共振成像方法的一个实施例流程图,该实施例可以包括以下步骤:
在步骤201中,采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位,其中,所述N为大于1的自然数。
多回波序列,是指一次射频脉冲激发后会通过射频接收线圈采集多个回波。也即在一个TR(Repetition Time,重复时间,是指相同切片的两次连续激励之间的时间)内,有N个回波,这N个回波分别对应回波时间TE1、TE2、…、TEN。
在一个示例中,可以利用三维梯度多回波序列采集不同回波时间的N个磁共振回波信号,N优选大于等于5。
采集序列例如图3所示。该梯度回波序列使用负极性的频率编码(读出)梯度场,自旋在小于90°的激励脉冲激发后立即开始散相,然后通过梯度场翻转,在TE时刻读出场中央自旋重聚得到回波。
图3中的横轴表示时间,纵轴简略地表明相对应的线圈的设置情况。此外,图3中的Gz表示Z方向的梯度磁场强度,是由Z轴梯度线圈133发出的,一般用于对人体进行选层;Gy表示Y方向的梯度磁场强度,是由Y轴梯度线圈132发出的,一般用于对人体进行相位编码;Gx表示X方向的梯度磁场强度,是由X轴梯度线圈131发出的,一般用于对人体进行频率编码。
虚线框中所示为所采集回波中的一个时序,TE和TR的设置是根据被检体的诊断需求由操作员确定。
采用梯度多回波序列进行磁共振回波信号的采集能够缩短梯度回波脉冲序列的TR,从而减少扫描时间。
在一个示例中,N个回波信号对应的翻转角θ为恩斯特角度(Ernst Angle)。
翻转角是指射频脉冲作用后,能够把宏观纵向磁化矢量翻转到偏离原来纵向方向的角度。在梯度回波序列中,由于TR比较短,因此存在一个最佳的角度,即恩斯特角度,使得一组给定组织在特定的TR中能够达到最佳信噪比。
在步骤202中,对所述N个回波的相位进行拟合,获得每个回波的有效相位。
在重建磁敏感加权图像SWI图像时,对各个回波数据对应的相位先进行拟合,得到一个高信噪比的有效相位图,并确定每个回波的有效相位。
在一个示例中,可以采用最小二乘法进行相位拟合。本领域技术人员应当了解,也可以采用其他适当的方法进行相位拟合。
在对N个回波的相位进行拟合前,所述方法还包括:对所述N个回波进行预处理;
所述预处理可以包括至少一种下述处理:
对所述N个回波的相位进行去卷绕和去背景场。
通过执行相位去卷绕,能够减少或移除原始相位图像中的相位卷绕伪影;通过执行相位去背景场,能够减小或移除背景场对相位信息的影响。因此,通过对采集的所有回波的对应相位进行预处理,有利于得到真实的组织相位。
在步骤203中,基于所述有效相位对所述N个回波的幅度进行加权,生成磁敏感加权图像。
在一个示例中,可以利用有效相位生成相位蒙版,利用该相位蒙版与N个回波的幅度进行加权,生成磁敏感加权图像。
相位蒙版用于抑制或增强特定范围内的相位,并且可以增强图像中呈现的感兴趣组织的对比度。相位蒙版的振幅例如可以在0与1之间变化,当相位蒙版与N个回波的幅度进行加权,如果相位蒙版的值为1,则相位蒙版对回波幅度不产生影响;如果相位蒙版的值小于1,则相位蒙版可以抑制相应像素的幅度值并增强原始幅度图像的对比度。
在一个示例中,对步骤202中得到的有效相位进行高通滤波,利用经滤波的有效相位生成相位蒙版,对N个回波的幅度进行加权。
通过对有效相位进行高通滤波,能够有效地降低空腔带来的磁化率伪影。
在本实施例中,通过对所采集的N个回波的相位进行拟合,获得每个回波的具有高信噪比的有效相位,并基于该有效相位对N个回波的幅度进行加权。由于TE短的回波受磁化率的影响小,因此通过于综合利用受磁化率影响小的回波相位值,提高了图像的信噪比,有效地降低强磁敏感的组织,例如空腔或出血部位,带来的磁化率伪影。
图4示出本申请磁共振成像方法的另一个实施例流程图,该实施例用于生成磁化率空间分布定量图,可以包括以下步骤:
在步骤401中,采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位的信息,其中,所述N为大于1的自然数;
在步骤402中,基于所述N个回波信号的相位获得每个回波对应的磁化率空间分布图;
磁化率空间分布QSM表征的是不同物质在主磁场中引起的局部磁场的偏移,其单位是偏离主磁场的幅度,用ppm(百万分之)或ppb(十亿分之)表示。
在一个示例中,可以利用磁化率反演算法,根据每个回波的相位计算N个回波下对应的磁化率空间分布QSM1、QSM2、…、QSMN。
磁化率反演公式如下所示:
其中,χn(r)为第n个回波中空间位置r处组织对应的磁化率,g(k)为格林函数,kx、ky、kz分别表示k空间中三个方向上的空间分量,FT为傅立叶变换,FT-1为傅立叶逆变换,φn(r)为第n个回波的经过相位解卷绕和背景去除处理后的相位,B0为主磁场强度,TE为回波时间,γ为旋磁比。
根据公式(1),对各个回波分别做磁化率反演计算求得各自对应的空间磁敏感分布定量图,按照TE时间对各个回波的空间磁敏感分布定量图进行加权合成得到最终的QSM输出。
在本步骤中,优选对N个回波的相位进行预处理后再进行磁化率空间分布的计算。预处理可以包括至少一种下述处理:
对所述N个回波的相位进行去卷绕和去背景场。
在步骤403中,利用每个回波对应的回波时间TE对每个回波对应的磁化率空间分布图进行加权,获得最终的磁化率空间分布定量图。
在本实施例中,通过对N个回波的相位进行磁化率反演获得各自对应的空间磁敏感分布定量图,并按照TE时间进行加权,得到最终的QSM输出。本实施例通过从一次数据采集中获取完全配准的多个磁共振定量拓扑信息和磁化率敏感加权图像以及对应的相位图像,在磁化率敏感相关的头部疾病精确诊断中具有重要价值。
图5示出本申请磁共振成像方法的另一个实施例流程图,该实施例用于生成增强磁敏感加权图像,可以包括以下步骤:
在步骤501中,采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位的信息,其中,所述N为大于1的自然数;
在步骤502中,基于所述N个回波信号的相位获得每个回波对应的磁化率空间分布图;
在步骤503中,利用第N个和第N-1个回波对应的磁化率空间分布图对第N个和第N-1个回波的幅值图像进行加权,获得各自对应的增强磁敏感加权图像;
在步骤504中,对第N个和第N-1个回波的增强磁敏感加权图像进行均值处理,得到最终的增强磁敏感加权图像。
在本实施例中,基于N个回波中,倒数两个回波对应的磁化率空间分布图,通过将其作为权重因子对对应的回波幅值图像进行加权,获得增强磁敏感加权图像,能够从一次数据采集中获取完全配准的磁共振定量拓扑信息和磁化率敏感加权图像,在磁化率敏感相关的头部疾病精确诊断中具有重要价值。
图6示出本申请磁共振成像方法的另一个实施例流程图,该实施例用于生成T2*定量拓扑图,可以包括以下步骤:
在步骤601中,采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位的信息,其中,所述N为大于1的自然数。
在本步骤中,回波信号例如可以描述为:
其中,ρ0为组织质子密度,θ为翻转角,TR为重复时间,TE为回波时间,T1为组织的T1弛豫时间,T2*为组织的横向弛豫衰减时间。
在步骤602中,基于N个回波信号获取对应的T2*弛豫时间。
T2*弛豫时间是本实施例中需要求解的定量值,根据公式(2)进行T2*弛豫时间的求解。
在步骤603中,将每个回波信号的幅度与T2*弛豫时间的非线性关系式转换为线性关系式。
在一个示例中,将公式(2)中指数项exp(-TE/T2*)之前的部分表示为:
可得:
S=ρθ*exp(-TE/T2*) (4)
对公式(4)两边求对数,可得:
logS=log(ρθ)+(-TE)*R2* (5)
通过以上过程,将为非线性表达式(2)转换为了线性表达式(5)。
在步骤604中,对N个回波对应的N组线性关系式(5)联立线性方程组,并利用最小二乘法获得T2*定量拓扑图和R2*定量拓扑图。
在本实施例中,通过将每个回波信号的幅度与T2*弛豫时间的非线性关系式转换为线性关系式,通过联立线性方程组获得T2*定量拓扑图和R2*定量拓扑图,充分利用了所采集的回波信息的幅度信息。
图7A和图7B示出通过本公开的磁共振成像方法进行成像的测试例,其中,图7A示出根据所采集的多个回波的相位信息,一次性获得的配准的SWI成像、QSM成像、eSWI成像以及高通相位图;图7B示出根据所采集的回波的幅度信息(模图),同时获得的R2*定量拓扑图和T2*定量拓扑图。
与前述方法的实施例相对应,本公开还提供了装置、终端设备和***的实施例。
参见图8,为本公开磁共振成像装置的一个实施例框图。该装置包括:
采集单元810,用于采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位,其中,所述N为大于1的自然数;
获得单元820用于利用所述N个回波的幅度信息和相位信息,获得不同类型的磁共振图像,所述不同类型的磁共振图像包括至少一个下述图像:磁敏感加权图像、磁化率空间分布定量图、增强磁敏感加权图像、T2*定量拓扑图;
其中,获得单元820包括:
拟合子单元821,用于对所述N个回波的相位进行拟合,获得每个回波的有效相位;
生成子单元822,用于基于所述有效相位对所述N个回波的幅度进行加权,生成磁敏感加权图像。
上述装置中各个单元的功能和作用的实现过程具体详见上述方法中对应步骤的实现过程,在此不再赘述。
对于装置实施例而言,由于其基本对应于方法实施例,所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本申请方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
参见图9,为本公开成像设备的一个实施例框图。该设备可以包括:内部总线910,以及通过内部总线连接的存储器920、处理器930和外部接口940。
其中,外部接口940,用于获得回波信号;
存储器920,用于存储磁共振成像对应的机器可读指令;
处理器930,用于读取存储器上的所述机器可读指令,并执行所述指令实现如下操作:
采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位,其中,所述N为大于1的自然数;
利用所述N个回波的幅度信息和相位信息,获得不同类型的磁共振图像,所述不同类型的磁共振图像包括至少一个下述图像:磁敏感加权图像、磁化率空间分布定量图、增强磁敏感加权图像、T2*定量拓扑图,
其中,获得所述磁敏感加权图像包括:
对所述N个回波的相位进行拟合,获得每个回波的有效相位;
基于所述有效相位对所述N个回波的幅度进行加权,生成磁敏感加权图像。
参见图10,为本公开磁共振成像***的一个实施例框图。该***可以包括:成像设备1010、扫描床1020、磁场发生器1030、信号探测器1040以及成像设备1050。
其中,信号探测器1040,用于向探测区域发送信号和从所述探测区域接收信号并成生成像数据;
磁场发生器1030,用于产生均匀磁场;
扫描床1020,用于承载待检测对象;
成像设备1010,用于采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位,其中,所述N为大于1的自然数;
利用所述N个回波的幅度信息和相位信息,获得不同类型的磁共振图像,所述不同类型的磁共振图像包括至少一个下述图像:磁敏感加权图像、磁化率空间分布定量图、增强磁敏感加权图像、T2*定量拓扑图,
其中,获得所述磁敏感加权图像包括:
对所述N个回波的相位进行拟合,获得每个回波的有效相位;
基于所述有效相位对所述N个回波的幅度进行加权,生成磁敏感加权图像。
在本公开实施例中,计算机可读存储介质可以是多种形式,比如,在不同的例子中,所述机器可读存储介质可以是:RAM(Radom Access Memory,随机存取存储器)、易失存储器、非易失性存储器、闪存、存储驱动器(如硬盘驱动器)、固态硬盘、任何类型的存储盘(如光盘、dvd等),或者类似的存储介质,或者它们的组合。特殊的,所述的计算机可读介质还可以是纸张或者其他合适的能够打印程序的介质。使用这些介质,这些程序可以被通过电学的方式获取到(例如,光学扫描)、可以被以合适的方式编译、解释和处理,然后可以被存储到计算机介质中。
以上所述仅为本公开的较佳实施例而已,并不用以限制本公开,凡在本公开的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本公开保护的范围之内。
Claims (11)
1.一种磁共振成像方法,其特征在于,包括:
采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位的信息,其中,所述N为大于1的自然数;
利用所述N个回波的幅度和相位的信息,获得不同类型的磁共振图像,所述不同类型的磁共振图像包括磁敏感加权图像和至少一个下述图像:磁化率空间分布定量图、增强磁敏感加权图像、T2*定量拓扑图,
其中,获得所述磁敏感加权图像包括:
对所述N个回波的相位进行拟合,获得每个回波的有效相位;
基于所述有效相位对所述N个回波的幅度进行加权,生成磁敏感加权图像。
2.根据权利要求1所述的方法,其特征在于,所述采集不同回波时间的N 个回波信号,获取所述N个回波信号的幅度和相位包括:
利用三维梯度多回波序列,采集不同回波时间的N 个回波信号。
3.根据权利要求1所述的方法,其特征在于,所述N个回波信号对应的翻转角为恩斯特角度。
4.根据权利要求1所述的方法,其特征在于,在对所述N个回波的相位进行拟合前,所述方法还包括:对所述N个回波进行预处理;
所述预处理包括至少一种下述处理:
对所述N个回波的相位进行去卷绕和去背景场。
5.根据权利要求1所述的方法,其特征在于,所述基于所述有效相位对所述N个回波的幅度进行加权包括:
对所述有效相位进行高通滤波,利用经滤波的有效相位对所述N个回波的幅度进行加权。
6.根据权利要求1所述的方法,其特征在于,获得所述磁化率空间分布定量图包括:
基于所述N个回波信号的相位获得每个回波对应的磁化率空间分布图;
利用每个回波对应的回波时间对每个回波对应的磁化率空间分布图进行加权,获得磁化率空间分布定量图。
7.根据权利要求6所述的方法,其特征在于,获得增强磁敏感加权图像包括:
利用第N个和第N-1个回波对应的磁化率空间分布图对第N个和第N-1个回波的幅值图像进行加权,获得各自对应的增强磁敏感加权图像;
对第N个和第N-1个回波的增强磁敏感加权图像进行均值处理,得到最终的增强磁敏感加权图像。
8.根据权利要求1所述的方法,其特征在于,获得T2*定量拓扑图包括:
基于N个回波信号获取对应的T2*弛豫时间;
将每个回波信号的幅度与T2*弛豫时间的非线性关系式转换为线性关系式;
对所述N个回波对应的N组线性关系式联立线性方程组,并利用最小二乘法获得T2*定量拓扑图和R2*定量拓扑图,其中,R2*=1/ T2*。
9.一种磁共振成像装置,其特征在于,包括:
采集单元,用于采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位,其中,所述N为大于1的自然数;
获得单元,用于利用所述N个回波的幅度信息和相位信息,获得不同类型的磁共振图像,所述不同类型的磁共振图像包括磁敏感加权图像和至少一个下述图像:磁化率空间分布定量图、增强磁敏感加权图像、T2*定量拓扑图;
所述获得单元包括:
拟合子单元,用于对所述N个回波的相位进行拟合,获得每个回波的有效相位;
生成子单元,用于基于所述有效相位对所述N个回波的幅度进行加权,生成磁敏感加权图像。
10.一种成像设备,其特征在于,所述设备包括:内部总线,以及通过内部总线连接的存储器、处理器和外部接口;其中,
所述外部接口,用于获得回波信号;
所述存储器,用于存储磁共振成像对应的机器可读指令;
所述处理器,用于读取所述存储器上的所述机器可读指令,并执行所述指令以实现如下操作:
采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位,其中,所述N为大于1的自然数;
利用所述N个回波的幅度信息和相位信息,获得不同类型的磁共振图像,所述不同类型的磁共振图像包括磁敏感加权图像和至少一个下述图像:磁化率空间分布定量图、增强磁敏感加权图像、T2*定量拓扑图,
其中,获得所述磁敏感加权图像包括:
对所述N个回波的相位进行拟合,获得每个回波的有效相位;
基于所述有效相位对所述N个回波的幅度进行加权,生成磁敏感加权图像。
11.一种磁共振成像***,其特征在于,所述***包括磁场发生器、信号探测器、扫描床和成像设备;其中,
所述磁场发生器,用于产生均匀磁场;
所述信号探测器,用于向探测区域发送信号和从所述探测区域接收信号并成生成像数据;
所述扫描床,用于承载待检测对象;
所述成像设备,用于:
采集不同回波时间的N个回波信号,获取所述N个回波信号的幅度和相位,其中,所述N为大于1的自然数;
利用所述N个回波的幅度信息和相位信息,获得不同类型的磁共振图像,所述不同类型的磁共振图像包括磁敏感加权图像和至少一个下述图像:磁化率空间分布定量图、增强磁敏感加权图像、T2*定量拓扑图,
其中,获得所述磁敏感加权图像包括:
对所述N个回波的相位进行拟合,获得每个回波的有效相位;
基于所述有效相位对所述N个回波的幅度进行加权,生成磁敏感加权图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910496273.2A CN110174632B (zh) | 2019-06-10 | 2019-06-10 | 磁共振成像方法、装置、成像设备及磁共振成像*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910496273.2A CN110174632B (zh) | 2019-06-10 | 2019-06-10 | 磁共振成像方法、装置、成像设备及磁共振成像*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110174632A CN110174632A (zh) | 2019-08-27 |
CN110174632B true CN110174632B (zh) | 2021-06-01 |
Family
ID=67698166
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910496273.2A Active CN110174632B (zh) | 2019-06-10 | 2019-06-10 | 磁共振成像方法、装置、成像设备及磁共振成像*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110174632B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110992435B (zh) * | 2019-11-06 | 2023-10-20 | 上海东软医疗科技有限公司 | 图像重建方法及设备、成像数据的处理方法及装置 |
CN111025210B (zh) * | 2019-12-20 | 2022-02-01 | 东软医疗***股份有限公司 | 磁共振成像方法、装置、电子设备、存储介质 |
CN114910850B (zh) * | 2021-02-10 | 2024-04-23 | 清华大学 | 一种双核mri的图像增强超构表面器件 |
CN114137462B (zh) * | 2021-11-19 | 2023-10-20 | 中国科学院深圳先进技术研究院 | 一种用于低场磁共振的多对比度成像方法和设备 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102116856A (zh) * | 2010-12-30 | 2011-07-06 | 中国科学院深圳先进技术研究院 | 横向弛豫时间测量方法及*** |
CN103871025A (zh) * | 2012-12-07 | 2014-06-18 | 深圳先进技术研究院 | 医学影像增强方法及其*** |
CN104267361A (zh) * | 2014-10-13 | 2015-01-07 | 厦门大学 | 基于结构特征的自适应定量磁化率分布图复合重建的方法 |
CN105997074A (zh) * | 2016-04-28 | 2016-10-12 | 华东师范大学 | 一种磁共振定量磁化率成像多回波相位拟合方法 |
CN107126211A (zh) * | 2017-06-07 | 2017-09-05 | 上海联影医疗科技有限公司 | 磁共振成像方法和*** |
CN107209243A (zh) * | 2015-02-03 | 2017-09-26 | 皇家飞利浦有限公司 | 用于磁敏感加权磁共振成像的方法和*** |
CN108375746A (zh) * | 2017-01-18 | 2018-08-07 | 上海联影医疗科技有限公司 | 一种相位反卷折方法与设备 |
CN108714028A (zh) * | 2018-04-11 | 2018-10-30 | 上海联影医疗科技有限公司 | 一种磁共振成像方法、装置及医学成像*** |
CN108956678A (zh) * | 2018-06-11 | 2018-12-07 | 西南石油大学 | 一种基于核磁共振测井的t2谱敏感参数提取方法 |
CN109100669A (zh) * | 2018-07-12 | 2018-12-28 | 厦门大学 | 基于重叠回波的单扫描同步磁共振扩散及t2成像方法 |
CN109477877A (zh) * | 2016-07-14 | 2019-03-15 | 上海联影医疗科技有限公司 | 磁共振成像***和方法 |
CN109564268A (zh) * | 2016-08-09 | 2019-04-02 | 皇家飞利浦有限公司 | 在多梯度回波mri中对场波动的回顾性校正 |
WO2019098714A1 (ko) * | 2017-11-17 | 2019-05-23 | 울산과학기술원 | 뇌에 존재하는 철의 시각화를 위한 방법 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10534057B2 (en) * | 2017-03-10 | 2020-01-14 | Maria Drangova | Method for dixon MRI, multi-contrast imaging and multi-parametric mapping with a single multi-echo gradient-recalled echo acquisition |
-
2019
- 2019-06-10 CN CN201910496273.2A patent/CN110174632B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102116856A (zh) * | 2010-12-30 | 2011-07-06 | 中国科学院深圳先进技术研究院 | 横向弛豫时间测量方法及*** |
CN103871025A (zh) * | 2012-12-07 | 2014-06-18 | 深圳先进技术研究院 | 医学影像增强方法及其*** |
CN104267361A (zh) * | 2014-10-13 | 2015-01-07 | 厦门大学 | 基于结构特征的自适应定量磁化率分布图复合重建的方法 |
CN107209243A (zh) * | 2015-02-03 | 2017-09-26 | 皇家飞利浦有限公司 | 用于磁敏感加权磁共振成像的方法和*** |
CN105997074A (zh) * | 2016-04-28 | 2016-10-12 | 华东师范大学 | 一种磁共振定量磁化率成像多回波相位拟合方法 |
CN109477877A (zh) * | 2016-07-14 | 2019-03-15 | 上海联影医疗科技有限公司 | 磁共振成像***和方法 |
CN109564268A (zh) * | 2016-08-09 | 2019-04-02 | 皇家飞利浦有限公司 | 在多梯度回波mri中对场波动的回顾性校正 |
CN108375746A (zh) * | 2017-01-18 | 2018-08-07 | 上海联影医疗科技有限公司 | 一种相位反卷折方法与设备 |
CN107126211A (zh) * | 2017-06-07 | 2017-09-05 | 上海联影医疗科技有限公司 | 磁共振成像方法和*** |
WO2019098714A1 (ko) * | 2017-11-17 | 2019-05-23 | 울산과학기술원 | 뇌에 존재하는 철의 시각화를 위한 방법 |
CN108714028A (zh) * | 2018-04-11 | 2018-10-30 | 上海联影医疗科技有限公司 | 一种磁共振成像方法、装置及医学成像*** |
CN108956678A (zh) * | 2018-06-11 | 2018-12-07 | 西南石油大学 | 一种基于核磁共振测井的t2谱敏感参数提取方法 |
CN109100669A (zh) * | 2018-07-12 | 2018-12-28 | 厦门大学 | 基于重叠回波的单扫描同步磁共振扩散及t2成像方法 |
Non-Patent Citations (2)
Title |
---|
Multi-atlas tool for automated segmentation of brain gray matter nuclei and quantification of their magnetic susceptibility;Xu Li;《NeuroImage》;20190207;全文 * |
三维多回波流动补偿定量磁化率成像方法研究;翁鹏飞;《中国优秀硕士学位论文全文数据库 信息科技辑》;20151015;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110174632A (zh) | 2019-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110174632B (zh) | 磁共振成像方法、装置、成像设备及磁共振成像*** | |
US8842896B2 (en) | Imaging of moving objects | |
Jiang et al. | Diffusion tensor imaging (DTI) of the brain in moving subjects: application to in‐utero fetal and ex‐utero studies | |
US8116541B2 (en) | Method and apparatus for multi-coil magnetic resonance imaging | |
US9766316B2 (en) | Magnetic resonance imaging device and quantitative susceptibility mapping method | |
US7800367B2 (en) | Method and apparatus for generating T2* weighted magnetic resonance images | |
US20080150532A1 (en) | Method and apparatus for measuring t1 relaxation | |
US8085041B2 (en) | Three-point method and system for fast and robust field mapping for EPI geometric distortion correction | |
US7999543B2 (en) | MR method for spatially-resolved determination of relaxation parameters | |
US10203387B2 (en) | MR imaging with enhanced susceptibility contrast | |
CN105308469A (zh) | 具有B0失真校正的并行MR成像以及使用规则化的SENSE重建的多点Dixon水-脂肪分离 | |
WO2012077543A1 (ja) | 磁気共鳴イメージング装置及びコントラスト強調画像取得方法 | |
RU2015116879A (ru) | Опорное сканирование при устойчивой к металлам мр визуализации | |
Kleban et al. | Strong diffusion gradients allow the separation of intra-and extra-axonal gradient-echo signals in the human brain | |
Cao et al. | Optimized multi‐axis spiral projection MR fingerprinting with subspace reconstruction for rapid whole‐brain high‐isotropic‐resolution quantitative imaging | |
Lee et al. | Rapid dual‐RF, dual‐echo, 3D ultrashort echo time craniofacial imaging: A feasibility study | |
JP2793125B2 (ja) | 核スピン断層撮影装置 | |
KR102163337B1 (ko) | 언더샘플링 방향 변화를 통한 다중영상획득 mri의 고속화 방법 및 이를 위한 장치 | |
US11815582B2 (en) | Dual echo steady state MR imaging using bipolar diffusion gradients | |
Naeyaert et al. | Accelerating in vivo fast spin echo high angular resolution diffusion imaging with an isotropic resolution in mice through compressed sensing | |
US8854039B2 (en) | System and method for quantifying MR data acquired during a UTE sequence | |
Li et al. | Comparison of uniform‐density, variable‐density, and dual‐density spiral samplings for multi‐shot DWI | |
Geraghty et al. | Partial Fourier reconstruction for improved resolution in 3D hyperpolarized 13C EPI | |
Yang et al. | Acquisition of Diffusion MRI Data | |
Kadalie et al. | Rapid whole brain 3D T2 mapping respiratory‐resolved Double‐Echo Steady State (DESS) sequence with improved repeatability |
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 |