CN102525460B - 一种磁共振成像水脂图像辨析方法及装置 - Google Patents
一种磁共振成像水脂图像辨析方法及装置 Download PDFInfo
- Publication number
- CN102525460B CN102525460B CN 201010610002 CN201010610002A CN102525460B CN 102525460 B CN102525460 B CN 102525460B CN 201010610002 CN201010610002 CN 201010610002 CN 201010610002 A CN201010610002 A CN 201010610002A CN 102525460 B CN102525460 B CN 102525460B
- Authority
- CN
- China
- Prior art keywords
- projection
- water
- fat image
- phase
- fat
- 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
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4869—Determining body composition
- A61B5/4872—Body fat
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4869—Determining body composition
- A61B5/4875—Hydration status, fluid retention of the body
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Radiology & Medical Imaging (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种磁共振成像水脂图像辨析方法和装置。方法包括:获取没有相位编码的三个回波作为参考扫描;利用参考扫描,计算相位编码方向的参考水脂图像投影;根据三点Dixon法得到的水图像和脂肪图像,计算相位编码方向的完全水脂图像投影;计算参考水脂图像投影与完全水脂图像投影之间的相关性,得到至少两个相关值,并从中获取最大相关值,确定所述用三点Dixon法计算出的图像的类型。本发明中的技术方案能够解决现有技术中的水脂调换问题,确定计算出的水脂图像的准确类型。
Description
技术领域
本发明涉及磁共振成像技术领域,尤其涉及一种磁共振成像水脂图像辨析方法及装置。
背景技术
磁共振成像是利用磁共振现象进行成像的一种技术。磁共振现象的原理主要包括:包含单数质子的原子核,例如人体内广泛存在的氢原子核,其质子具有自旋运动,犹如一个小磁体,并且这些小磁体的自旋轴无一定的规律,如果施加外在磁场,这些小磁体将按外在磁场的磁力线重新排列,具体为在平行于或反平行于外在磁场磁力线的两个方向排列,将上述平行于外在磁场磁力线的方向称为正纵向轴,将上述反平行于外在磁场磁力线的方向称为负纵向轴,原子核只具有纵向磁化分量,该纵向磁化分量既具有方向又具有幅度。用特定频率的射频(RF,Radio Frequency)脉冲激发处于外在磁场中的原子核,使这些原子核的自旋轴偏离正纵向轴或负纵向轴,产生共振,这就是磁共振现象。上述被激发原子核的自旋轴偏离正纵向轴或负纵向轴之后,原子核具有了横向磁化分量。
停止发射射频脉冲后,被激发的原子核发射回波信号,将吸收的能量逐步以电磁波的形式释放出来,其相位和能级都恢复到激发前的状态,将原子核发射的回波信号经过空间编码等进一步处理即可重建图像。上述被激发原子核向激发前状态的恢复过程称为驰豫过程,恢复到平衡状态所需的时间称为驰豫时间。
由于人体内部脂肪中的氢原子核和水中的氢原子核所处的分子环境不相同,在使用相同射频脉冲进行激发时,它们的共振频率不相同,驰豫时间也不相同。在不同的回波时间采集信号,脂肪组织和水表现出不同的相位以及信号强度。
狄克逊(Dixon)法是在磁共振成像中用以产生纯水质子图像的方法,其基本原理是分别采集水和脂肪质子的同相位(In Phase)和反相位(Out phase)两种回波信号,两种不同相位的信号通过运算,去除脂肪信号,产生一幅纯水质子的影像,从而达到脂肪抑制的目的。Dixon法的缺点是受磁场不均匀性影响较大,易受呼吸运动影响,计算方法复杂并容易出现错误。
为了同时得到分离的纯水和纯脂肪图像,一种改进的三点Dixon方法被广泛应用,该方法的原理是同时取得一幅同相位(或反相位)图像和两幅反相位(或同相位)图像,根据两幅反相位(或同相位)图像,求得磁场不均匀导致的附加相位,对两幅反相位(或同相位)图像进行相位校正,然后与同相位(或反相位)图像一起求得水的图像及脂肪的图像。
三点Dixon法的具体操作过程有多种,例如,以利用一幅同相位图像和两幅反相位图像的情况为例,发明人为汪坚敏和翁得河、申请号为200510008973.0的中国专利申请中公开的一种水脂分离图像重建方法,该方法包括以下步骤:取得一幅同相位图像和两幅反相位图像;求各通道的数据线圈灵敏度分布;合成各通道图像;求两幅反相位图像的相位差;检测同相位图像中的一些特征区域以作为修正相位的判据;以及修正反相位图像的相位,计算出水和脂肪的图像。
此外发明人为吴东梅和王振奎、申请号为200910119608.5的中国专利申请中也公开的一种水脂分离实现方法,该方法包括以下步骤:获取一幅同相位图像以及两幅反相位图像;计算所述两幅反相位图像的相位差,并利用该相位差对所述两幅反相位图像中由磁场不均匀性所引起的相位进行修正;利用修正后的两幅反相位图像,对所述同相位图像中由涡电流所引起的线性相位差进行修正;根据修正后的三幅图像,计算得到水和脂肪图像。
图1给出了基于快速自旋回波(TSE)序列的三点Dixon法获取两个反相位图像和一个同相位图像的数据采集及编码序列的示意图。其中,RF、RO和PE分别对应射频脉冲、读出梯度和编码梯度。
如图1所示,磁共振成像设备首先发射一个90度射频脉冲RF_0,然后再发射一个180度重聚相射频脉冲RF_1。在距离90度射频脉冲RF_0的回波时间(TE)之后,磁共振成像设备在读出梯度方向上施加读出梯度,分别读取三个回波1、2、3。接着再发射一个180度重聚相射频脉冲RF_2,并在读出梯度方向上施加读出梯度,分别读出三个回波4、5、6;上述操作可重复多次。图1序列在改变相位编码梯度后,不断地被重复,直至读取k空间上的所有回波。其中,回波1、4、7……等对应的数据线构成一幅反相位图像的原始数据,回波2、5、8……等对应的数据线构成一幅正相位图像的原始数据,回波3、6、9……等对应的数据线构成另一幅反相位图像的原始数据。这些回波通过用在相位编码方向施加的相位编码梯度的适当调制而得到k空间上不同行的数据,当k空间被填满后,经过傅立叶变换,构成对应的一幅同相位图像和两幅反相位图像。
但现有三点Dixon方法的不足在于:由于两幅反相位(同相位)图像求得的附加相位不能直接用来对其相位进行纠正,而必须对相位进行反缠绕,而由于相位反缠绕固有的不稳定性,使得有时计算出来的水和脂肪的图像可能调换,即实际上是水的图像,计算出来却可能是脂肪的图像;而实际上是脂肪的图像,计算出来却可能是水的图像。
基于这种水脂图像调换的问题,目前尚没有相应的解决方案,为此如何在磁共振成像设备中确定当前计算出的水(脂肪)图像实际上究竟是水的图像还是脂肪的图像,已成为磁共振成像技术中亟待解决的问题。
发明内容
有鉴于此,本发明提出了一种磁共振成像水脂图像辨析方法及装置,用以确定计算出的水脂图像的准确类型。
本发明提供的磁共振成像水脂图像辨析方法,包括:
获取没有相位编码的三个回波作为参考扫描;
利用所述参考扫描,计算相位编码方向的参考水图像投影和/或参考脂肪图像投影,得到参考水脂图像投影;
根据三点Dixon法计算得到的水图像和脂肪图像,计算相位编码方向的完全水图像投影和/或完全脂肪图像投影,得到完全水脂图像投影;
计算参考水脂图像投影与完全水脂图像投影之间的相关性,得到至少两个相关值;
从计算得到的相关值中获取最大相关值,将所述最大相关值对应的完全图像投影的类型确定为所述最大相关值对应的参考图像投影的类型,根据所述完全图像投影的类型确定所述用三点Dixon法计算出的图像的类型。
本发明提供的磁共振成像水脂图像辨析装置,包括:
参考扫描获取单元,用于获取没有相位编码的三个回波作为参考扫描;
参考图像投影计算单元,用于利用所述参考扫描获取单元获取的参考扫描,计算相位编码方向的参考水图像投影和/或参考脂肪图像投影,得到参考水脂图像投影;
完全图像投影计算单元,用于根据三点Dixon法计算得到的水图像和脂肪图像,计算相位编码方向的完全水图像投影和/或完全脂肪图像投影,得到完全水脂图像投影;
相关性计算单元,计算参考水脂图像投影与完全水脂图像投影之间的相关性,得到至少两个相关值;
结果判定单元,用于从所述相关性计算单元计算得到的相关值中获取最大相关值,将所述最大相关值对应的完全图像投影的类型确定为所述最大相关值对应的参考图像投影的类型,根据所述完全图像投影的类型确定对应的用三点Dixon法计算出的图像的类型。
从上述方案中可以看出,本发明实施例中利用没有相位编码的三个回波作为参考扫描,并利用该参考扫描计算相位编码方向的参考水图像投影和参考脂肪图像投影,即得到水脂类型准确的图像投影;之后以该参考水图像投影和参考脂肪图像投影作为标准,对根据三点Dixon法得到的水图像和脂肪图像计算得到的相位编码方向的完全水图像投影和完全脂肪图像投影进行相关性计算,根据相关值的大小,对三点Dixon法得到的水脂图像进行辨析,从而解决了现有技术中水脂图像调换的问题。
进一步地,通过在进行相关性计算之前,对参考水图像投影和参考脂肪图像投影进行归一化处理,使得参考水图像投影和参考脂肪图像投影与完全水图像投影和完全脂肪图像投影具有相同的功率谱,提高了相关性的准确度,进而增强了对三点Dixon法得到的水脂图像辨析的准确率。
附图说明
下面将通过参照附图详细描述本发明的优选实施例,使本领域的普通技术人员更清楚本发明的上述及其它特征和优点,附图中:
图1为现有技术中基于快速自旋回波序列的三点Dixon法获取两个反相位图像和一个同相位图像的数据采集及编码序列的示意图。
图2为本发明实施例中磁共振成像水脂图像辨析方法的示例性流程示意图。
图3为本发明实施例中磁共振成像水脂图像辨析装置的示例性结构示意图。
图4为本发明实施例中参考图像投影计算单元的结构示意图。
图5为本发明实施例中相关性计算单元的结构示意图。
图6为在0.35T磁共振成像设备上利用本发明技术方案得到的参考水脂图像投影和完全水脂图像投影。
图7为在1.5T磁共振成像设备上利用本发明技术方案得到的参考水脂图像投影和完全水脂图像投影。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,以下举实施例对本发明进一步详细说明。
本发明实施例中,首先利用没有相位编码的三个回波作为参考扫描,并利用参考扫描,计算得到参考水脂图像投影,该参考水脂图像投影在计算过程中由于不涉及相位编码及相位反缠绕,因此不会出现图像的调换,即参考水图像投影对应的一定是水的图像投影,参考脂肪图像投影对应的一定是脂肪的图像投影。之后,根据三点Dixon法计算的水脂图像得到完全水脂图像投影,再计算参考水脂图像投影与完全水脂图像投影之间的相关性,根据相关值的大小,确定完全水脂图像投影的类别,进而确定三点Dixon法计算的水脂图像的类别,解决图像调换的问题。
图2给出了本发明实施例中磁共振成像水脂图像辨析方法的示例性流程示意图,如图2所示,该方法包括如下步骤:
步骤201,获取没有相位编码的三个回波作为参考扫描,所述三个回波包括两个反相位回波和一个同相位回波,或者所述三个回波包括两个同相位回波和一个反相位回波。
本步骤中,没有相位编码的三个回波指的是没有通过在相位编码方向施加的相位编码梯度进行调制的三个回波。并且三个回波可以基于各种序列进行采集,且可以有不同的采集顺序,以三个回波包括两个反相位回波和一个同相位回波的情况为例,则三个回波可以是按照图1所示方法读出的某三个回波,也可以是基于其它序列,如回波序列(SE)或梯度回波序列(GRE)等的三点Dixon方法读出的某三个回波,并且三个回波可以是如图1所示的按照反相位-同相位-反相位的顺序依次采集的三个回波,也可以是先采集一个同相位的回波,然后采集两个反相位的回波,还可以是先采集两个反相位的回波,然后采集一个同相位的回波。本实施例中,假设获取的没有相位编码的三个回波为图1中的回波1、2、3。当然,实际应用中,没有相位编码的三个回波也可以是其它位置处的三个回波。
对于三个回波包括两个同相位回波和一个反相位回波的情况,采集方法与上述方法类似。同样可以基于各种序列进行采集,并且同样可以有不同的采集顺序,例如,三个回波可以是如按照同相位-反相位-同相位的顺序依次采集的三个回波,也可以是先采集一个反相位的回波,然后采集两个同相位的回波,还可以是先采集两个同相位的回波,然后采集一个反相位的回波。
步骤202,利用所述参考扫描,计算相位编码方向的参考水图像投影和参考脂肪图像投影,得到参考水脂图像投影。
本步骤中,对于三个回波包括两个反相位回波和一个同相位回波的情况,可首先利用两个反相位回波计算两个反相位之间的相位差,之后利用得到的相位差在相位编码方向计算参考水图像投影和参考脂肪图像投影。
具体实现时,可根据三个回波水和脂肪的化学位移相位情况进行计算,例如,若两个反相位回波和一个同相位回波水和脂肪的化学位移相位为(-π,0,π),则可根据如下公式(1)计算参考水脂图像投影。
Δφ=ang(l3*conj(l1))/2
wpr=l2+0.5·(l3·e-j·Δφ+l1·ej·Δφ) (1)
fpr=l2-0.5·(l3·e-j·Δφ+l1·ej·Δφ)
又如,若两个反相位和一个同相位回波水和脂肪的化学位移相位为(-3π,0,π),则可按照如下公式(2)计算参考水脂图像投影。
Δφ=ang(l3*conj(l1))/4
wpr=l2+0.5·(l3·e-j·Δφ+l1·ej·3·Δφ) (2)
fpr=l2-0.5·(l3·e-j·Δφ+l1·ej·3·Δφ)
其中,l1、l3为两个反相位回波,l2为一个同相位回波,若步骤201中获取的没有相位编码的三个回波为图1中的回波1、2、3,则l1、l2、l3分别为回波1、2、3;ang()为相位获取函数,Δφ为两个反相位之间的相位差;conj()为共轭函数;wpr和fpr分别为直接利用Δφ进行相位纠正后计算得到的参考水图像投影和参考脂肪图像投影。
上述过程中,在进行相位修正时直接利用Δφ进行相位纠正,而无须对其进行反缠绕处理后再进行相位修正,因此不会发生图像的调换问题,即参考水图像投影一定会是对应水的参考图像投影,而不会是对应脂肪的参考图像投影;反之亦然。
对于三个回波包括两个同相位回波和一个反相位回波的情况,计算方法与上述方法类似。同样可首先利用两个同相位回波计算两个同相位之间的相位差,之后利用得到的相位差在相位编码方向计算参考水图像投影和参考脂肪图像投影。
具体实现时,同样可以根据三个回波水和脂肪的化学位移相位情况进行计算,例如,若两个同相位回波和一个反相位回波水和脂肪的化学位移相位为(0,π,2π),则可按照如下公式(3)计算参考水脂图像投影。
Δφ=ang(l3*conj(l1))/2
wpr=0.5·(l3·e-j·2·Δφ+l1)+l2·e-j·Δφ (3)
fpr=0.5·(l3·e-j·2·Δφ+l1)-l2·e-j·Δφ
又如,若两个同相位和一个反相位回波水和脂肪的化学位移相位为(0,π,4π),则可按照如下公式(4)计算参考水脂图像投影。
Δφ=ang(l3*conj(l1))/4
wpr=0.5·(l3·e-j·4·Δφ+l1)+l2·e-j·Δφ (4)
fpr=0.5·(l3·3-j·4·Δφ+l1)-l2·e-j·Δφ
其中,l1、l3为两个同相位回波,l2为一个反相位回波;ang()为相位获取函数,Δφ为两个同相位之间的相位差;conj()为共轭函数;wpr和fpr分别为直接利用Δφ进行相位纠正后计算得到的参考水图像投影和参考脂肪图像投影。
同样,上述过程中,在进行相位修正时也是直接利用Δφ进行相位纠正,而无须对其进行反缠绕处理后再进行相位修正,因此不会发生图像的调换问题,即参考水图像投影一定会是对应水的参考图像投影,而不会是对应脂肪的参考图像投影;反之亦然。
实际应用中,还可以采用其他计算方法,只要能够利用所述参考扫描,计算得到相位编码方向的参考水图像投影和参考脂肪图像投影即可,此处不再一一列举。
步骤203,利用三点Dixon法计算得到水图像和脂肪图像。
本步骤中,可利用现有的及各种新增的三点Dixon法。
例如,可以利用发明人为汪坚敏和翁得河的中国专利申请200510008973.0中公开的水脂分离图像重建方法,也可以利用发明人为吴东梅和王振奎的中国专利申请200910119608.5中公开的水脂分离实现方法。此外,也可以采用利用两幅同相位图像和一幅反相位图像进行水脂图像计算的三点Dixon法,此处不再一一列举。
步骤204,利用所述水图像和脂肪图像,计算相位编码方向的水图像投影和脂肪图像投影。为与步骤202中的参考水图像投影和参考脂肪图像投影相区别,此处分别称为完全水图像投影和完全脂肪图像投影,从而得到完全水脂图像投影。
本步骤中,计算相位编码方向的完全水图像投影和完全脂肪图像投影的方法可有多种,如下的公式(5)列出了其中的一种计算方法。
其中,wpc和fpc分别为完全水图像投影和完全脂肪图像投影,wi和fi分别为沿相位编码方向的体素点数据,n为相位编码方向的分辨率。
具体实现时,也可以采用除公式(5)之外的其他计算方法,只要能够利用所述水图像和脂肪图像,计算得到相位编码方向的完全水图像投影和完全脂肪图像投影即可,此处不再一一列举。
步骤205,计算参考水脂图像投影与完全水脂图像投影之间的相关性,得到至少两个相关值。
本步骤中,计算参考水脂图像投影与完全水脂图像投影之间的相关性时,可有多种组合方式。例如,可计算参考水图像投影与完全水图像投影的相关性,参考水图像投影与完全脂肪图像投影的相关性,参考脂肪图像投影与完全水图像投影的相关性,以及参考脂肪图像投影与完全脂肪图像投影的相关性,得到四个相关值。
或者,也可仅计算其中的两个具有可比性的相关值。如,可只计算参考水图像投影与完全水图像投影的相关性,以及参考水图像投影与完全脂肪图像投影的相关性,得到两个相关值。这种情况下,步骤202中可仅计算参考水图像投影,即参考脂肪图像投影可不计算,此时,参考水脂图像投影中可仅包括参考水图像投影。当然,也可计算参考脂肪图像投影。
又如,也可只计算参考脂肪图像投影与完全水图像投影的相关性,以及参考脂肪图像投影与完全脂肪图像投影的相关性,得到两个相关值。这种情况下,步骤202中可仅计算参考脂肪图像投影,即参考水图像投影可不计算,此时,参考水脂图像投影中可仅包括参考脂肪图像投影。当然,也可计算参考水图像投影。
再如,还可只计算参考水图像投影与完全水图像投影的相关性,以及参考脂肪图像投影与完全水图像投影的相关性,得到两个相关值。这种情况下,步骤204中可仅计算完全水图像投影,即完全脂肪图像投影可不计算,此时,完全水脂图像投影中可仅包括完全水图像投影。当然,也可计算完全脂肪图像投影。
再如,还可只计算参考水图像投影与完全脂肪图像投影的相关性,以及参考脂肪图像投影与完全脂肪图像投影的相关性,得到两个相关值。这种情况下,步骤204中可仅计算完全脂肪图像投影,即完全水图像投影可不计算,此时,完全水脂图像投影中可仅包括完全脂肪图像投影。当然,也可计算完全水图像投影。
本步骤中,对于分别表示参考水图像投影和参考脂肪图像投影的wpr和fpr,以及分别表示完全水图像投影和完全脂肪图像投影的wpc和fpc,可分别计算wpr与wpc的相关性,wpr与fpc的相关性,fpr与wpc的相关性,fpr与fpc的相关性,共得到四个相关值。或者,也可仅计算其中具有可比性的两个相关值,例如,仅计算wpr与wpc的相关性,wpr与fpc的相关性;或fpr与wpc的相关性,fpr与fpc的相关性;或wpr与wpc的相关性,fpr与wpc的相关性;或wpr与fpc的相关性,fpr与fpc的相关性。
具体实现时,可利用公式(6)计算上述相关性。
Crc=max(FT(FT(Pr)·conj(FT(Pc)))) (6)
其中,Crc为相关值,Pr为wpr或fpr,Pc为wpc或fpc,FT()为傅立叶变换函数,conj()为共轭函数;max()为求最大值的函数。
实际应用时,由于步骤202中得到的参考水脂图像投影与步骤204中得到的完全水脂图像投影,通常具有不同的功率谱幅度,因此为了获取与步骤204中得到的完全水脂图像投影相同的功率谱,在相关性计算之前,可先对参考水脂图像投影进行归一化,之后,利用归一化后的参考水脂图像投影进行上述的相关性计算。当然,也可以不进行归一化。
具体实现时,也可以采用除公式(6)之外的其他计算方法,例如,采用如下所示公式(7)计算相关性。
Crc=sum(abs(Pr-Pc)) (7)
其中,Crc为相关值,Pr为wpr或fpr,Pc为wpc或fpc,sum为求和算子,abs为求绝对值算子。
此外,还可以采用其他的计算方法,只要能够计算得到参考水图像投影与完全水图像投影及与完全脂肪图像投影的相关性,以及参考脂肪图像投影与完全水图像投影及与完全脂肪图像投影的相关性即可,此处不再一一列举。
步骤206,从计算得到的相关值中获取最大相关值,将所述最大相关值对应的完全图像投影的类型确定为所述最大相关值对应的参考图像投影的类型,根据所述完全图像投影的类型确定对应的用三点Dixon法计算出的图像的类型,即步骤203中计算出的图像的类型。
通过上述方法,可以验证步骤203中利用三点Dixon法计算出的水图像和脂肪图像是否发生了调换。
具体实现时,步骤201-步骤202和步骤203-步骤204之间没有绝对的先后顺序,即也可以先执行步骤203-步骤204,再执行步骤201-步骤202,或者步骤201-步骤202和步骤203-步骤204也可以并行或交叉执行。
以上对本发明实施例中的磁共振成像水脂图像辨析方法进行了详细描述,下面再对本发明实施例中磁共振成像水脂图像辨析装置进行详细描述。
图3给出了本发明实施例中磁共振成像水脂图像辨析装置的示例性结构示意图,如图3所示,该装置包括:一个参考扫描获取单元310、一个参考图像投影计算单元320、一个完全图像投影计算单元330、一个相关性计算单元340和一个结果判定单元350。
其中,参考扫描获取单元310用于获取没有相位编码的三个回波作为参考扫描。其中,所述三个回波包括两个反相位回波和一个同相位回波,或者所述三个回波包括两个同相位回波和一个反相位回波。具体实现时,参考扫描获取单元310的具体操作过程可与步骤201中描述的具体操作过程一致。
参考图像投影计算单元320用于利用所述参考扫描获取单元310获取的参考扫描,计算相位编码方向的参考水图像投影和参考脂肪图像投影,得到参考水脂图像投影。具体实现时,参考图像投影计算单元320的具体操作过程可与步骤202中描述的具体操作过程一致。
完全图像投影计算单元330用于根据三点Dixon法计算得到的水图像和脂肪图像,计算相位编码方向的完全水图像投影和完全脂肪图像投影,得到完全水脂图像投影。具体实现时,完全图像投影计算单元330的具体操作过程可与步骤204中描述的具体操作过程一致。
这里,三点Dixon法计算得到的水图像和脂肪图像可以是利用本实施例中的装置之外的其他装置计算得到的,此时,本实施例中的装置不包括水脂图像计算单元;或者,也可以是利用本实施例中的装置计算得到的,此时,本实施例中的装置中还包括水脂图像计算单元360(如图3中的虚线部分所示),用于利用三点Dixon法计算得到水图像和脂肪图像,将所述水图像和脂肪图像提供给所述完全图像投影计算单元330进行完全水脂图像投影的计算。
相关性计算单元340用于利用所述参考图像投影计算单元320得到的参考水脂图像投影,以及所述完全图像投影计算单元330得到的完全水脂图像投影,计算参考水脂图像投影与完全水脂图像投影之间的相关性,得到至少两个相关性值。具体实现时,相关性计算单元340的具体操作过程可与步骤205中描述的具体操作过程一致。即同样存在与步骤205中所描述的各种组合方式的相关性计算。相应地,当只需计算具有可比性的两个相关值时,参考图像投影计算单元320或完全图像投影计算单元330可只计算进行相关性计算时所需的一个图像投影,此时参考水脂图像投影或完全水脂图像投影中可仅包括一个对应的图像投影。当然,如果都计算也可以。
结果判定单元350用于从所述相关性计算单元340计算得到的相关值中获取最大相关值,将所述最大相关值对应的完全图像投影的类型确定为所述最大相关值对应的参考图像投影的类型,根据所述完全图像投影的类型确定对应的用三点Dixon法计算出的图像的类型。
具体实现时,若为了使参考水脂图像投影获取与完全水脂图像投影相同的功率谱,在相关性计算之前,可先对参考水脂图像投影进行归一化,之后,利用归一化后的参考水脂图像投影进行上述的相关性计算。此时,本实施例的装置中可进一步包括归一化单元370(如图3中的虚线部分所示),用于对所述参考图像投影计算单元320得到的参考水脂图像投影进行归一化。相应地,相关性计算单元340用于利用所述归一化单元370归一化后的参考水脂图像投影,以及所述完全图像投影计算单元330得到的完全水脂图像投影,计算参考水脂图像投影与完全水脂图像投影之间的相关性,得到至少两个相关值。当然,也可以不进行归一化,此时,本实施例的装置中可无需包括归一化单元370。
具体实现时,基于不同的计算方法,参考图像投影计算单元320可具有不同的内部结构,图4示出了参考图像投影计算单元320的一种内部结构示意图。如图4中的实线部分所示,该参考图像投影计算单元320可具体包括:
一个第一计算子单元321,用于利用两个反相位回波计算两个反相位之间的相位差。
一个第二计算子单元322,用于利用所述第一计算子单元321得到的相位差在相位编码方向计算参考水图像投影和/或参考脂肪图像投影。
或者,该参考图像投影计算单元320也可如图4中的虚线部分所示,包括:
一个第三计算子单元323,用于利用两个同相位回波计算两个同相位之间的相位差。
一个第四计算子单元324,用于利用所述第三计算子单元323得到的相位差在相位编码方向计算参考水图像投影和/或参考脂肪图像投影。
又或者,该参考图像投影计算单元320也可同时包括上述四个计算子单元,具体进行计算时,该参考图像投影计算单元320可根据参考扫描的类型进行选用,即由选择单元判断参考扫描的类型,如果参考扫描包括两个反相位回波和一个同相位回波,则该参考图像投影计算单元320选用第一计算子单元321和第二计算子单元322;如果参考扫描包括两个同相位回波和一个反相位回波,则该参考图像投影计算单元320选用第三计算子单元323和第四计算子单元324。
具体实现时,基于不同的相关性组合方式,相关性计算单元340也可具有不同的内部结构,图5示出了相关性计算单元340的一种内部结构示意图。如图5所示,该相关性计算单元340包括:
一个第一相关性计算子单元341,用于计算参考水图像投影与完全水图像投影的相关性;
一个第二相关性计算子单元342,用于计算参考水图像投影与完全脂肪图像投影的相关性;
一个第三相关性计算子单元343,用于计算参考脂肪图像投影与完全水图像投影的相关性;
一个第四相关性计算子单元344,用于计算参考脂肪图像投影与完全脂肪图像投影的相关性。
此外,基于其它相关性组合,相关性计算单元340也可只包括上述第一相关性计算子单元341和第二相关性计算子单元342;或者只包括上述第三相关性计算子单元341和第四相关性计算子单元342;或者只包括上述第一相关性计算子单元341和第三相关性计算子单元342;或者只包括上述第二相关性计算子单元341和第四相关性计算子单元342。
上面对本发明一实施例中的磁共振成像水脂图像辨析装置进行了详细描述。具体实现时,磁共振成像水脂图像辨析装置也可以有其他形式的内部结构,例如,可对上述各功能单元进行重组或进一步细化。并且各功能单元可以位于不同的实体中,也可以全部位于同一实体中,具体可根据实际需要进行配置。
本申请的发明人利用一台0.35T的磁共振成像设备对头部进行了水脂分离成像。所用参数包括:视野(FOV)=24×24cm,回波时间/重复时间(TE/TR)=34/695ms,矩阵(matrix)=256×256。图6为本申请的发明人利用本发明的方法计算得到的参考水图像投影、参考脂肪图像投影、完全水图像投影和完全脂肪图像投影。图中横坐标为像素点,纵坐标为投影的值。通过图6可以看出,wpr与fpc有较好的相关性,fpr与wpc有较好的相关性,因此可以确定脂肪图像和水图像。
本申请的发明人还利用一台1.5T的磁共振成像设备对膝盖进行了水脂分离成像。所用参数包括:视野(FOV)=26.3×35cm,回波时间/重复时间(TE/TR)=13/3500ms,矩阵(matrix)=231×256。图7为本申请的发明人利用本发明的方法计算得到的参考水图像投影、参考脂肪图像投影、完全水图像投影和完全脂肪图像投影。图中横坐标为像素点,纵坐标为投影的值。通过图7可以看出,wpr与fpc有较好的相关性,fpr与wpc有较好的相关性,因此可以确定脂肪图像和水图像。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (12)
1.一种磁共振成像水脂图像辨析方法,该方法包括:
获取没有相位编码的三个回波作为参考扫描;
利用所述参考扫描,计算相位编码方向的参考水图像投影和/或参考脂肪图像投影,得到参考水和/或脂图像投影;
根据三点Dixon法计算得到的水图像和脂肪图像,计算相位编码方向的完全水图像投影和/或完全脂肪图像投影,得到完全水和/或脂图像投影;
计算参考水脂图像投影与完全水脂图像投影之间的相关性,得到至少两个相关值;
从计算得到的相关值中获取最大相关值,并确定用三点Dixon法计算出的图像的类型。
2.根据权利要求1所述的方法,其特征在于,该方法进一步包括:在计算相关性之前,对参考水和/或脂图像投影进行归一化。
3.根据权利要求1或2所述的方法,其特征在于,计算相关性的步骤包括:
计算参考水图像投影与完全水图像投影之间的相关性、参考水图像投影与完全脂肪图像投影之间的相关性、参考脂肪图像投影与完全水图像投影之间的相关性、以及参考脂肪图像投影与完全脂肪图像投影之间的相关性中的至少两者。
4.根据权利要求3所述的方法,其特征在于,
按照如下公式计算相关性
Crc=max(FT(FT(Pr)·conj(FT(Pc))))
或者,按照如下公式计算相关性
Crc=sum(abs(Pr-Pc))
其中,Crc为相关值,Pr为wpr或fpr,Pc为wpc或fpc,FT()为傅立叶变换函数,conj()为共轭函数,max()为求最大值的函数,sum为求和算子,abs为求绝对值算子;其中,wpr和fpr分别为参考水图像投影和参考脂肪图像投影,wpc和fpc分别为完全水图像投影和完全脂肪图像投影。
5.根据权利要求1所述的方法,其特征在于,所述三个回波包括两个反相位回波和一个同相位回波;利用两个反相位回波计算两个反相位之间的相位差,利用得到的相位差在相位编码方向计算参考水图像投影和/或参考脂肪图像投影;
或者,
所述三个回波包括两个同相位回波和一个反相位回波;利用两个同相位回波计算两个同相位之间的相位差,利用得到的相位差在相位编码方向计算参考水图像投影和/或参考脂肪图像投影。
6.根据权利要求5所述的方法,其特征在于,
在两个反相位回波和一个同相位回波水和脂肪的化学位移相位为(-π,0,π)时,按照如下公式计算参考水图像投影和参考脂肪图像投影
Δφ=ang(l3*conj(l1))/2
wpr=l2+0.5·(l3·e-j·Δφ+l1·ej·Δφ)
fpr=l2-0.5·(l3·e-j·Δφ+l1·ej·Δφ)
在两个反相位和一个同相位回波水和脂肪的化学位移相位为(-3π,0,π)时,按照如下公式计算参考水图像投影和参考脂肪图像投影
Δφ=ang(l3*conj(l1))/4
wpr=l2+0.5·(l3·e-j·Δφ+l1·ej·3·Δφ)
fpr=l2-0.5(l3·e-j·Δφ+l1·ej·3·Δφ)
其中,l1、l3为两个反相位回波,l2为一个同相位回波;ang()为相位获取函数,Δφ为两个反相位之间的相位差;conj()为共轭函数;wpr和fpr分别为直接利用Δφ进行相位纠正后计算得到的参考水图像投影和参考脂肪图像投影。
7.根据权利要求5所述的方法,其特征在于,
在两个同相位回波和一个反相位回波水和脂肪的化学位移相位为(0,π,2π)时,按照如下公式计算参考水图像投影和参考脂肪图像投影
Δφ=ang(l3*conj(l1))/2
wpr=0.5·(l3·e-j·2·Δφ+l1)+l2·e-j·Δφ
fpr=0.5·(l3·e-j2·Δφ+l1)-l2·e-j·Δφ
在两个同相位和一个反相位回波水和脂肪的化学位移相位为(0,π,4π)时,按照如下公式计算参考水图像投影和参考脂肪图像投影
Δφ=ang(l3*conj(l1))/4
wpr=0.5·(l3·e-j·4·Δφ+l1)+l2·e-j·Δφ
fpr=0.5·(l3·e-j·4·Δφ+l1)-l2·e-j·Δφ
其中,l1、l3为两个同相位回波,l2为一个反相位回波;ang()为相位获取函数,Δφ为两个同相位之间的相位差;conj()为共轭函数;wpr和fpr分别为直接利用Δφ进行相位纠正后计算得到的参考水图像投影和参考脂肪图像投影。
8.一种磁共振成像水脂图像辨析装置,该装置包括:
一个参考扫描获取单元(310),用于获取没有相位编码的三个回波作为参考扫描;
一个参考图像投影计算单元(320),用于利用所述参考扫描,计算相位编码方向的参考水图像投影和/或参考脂肪图像投影,得到参考水和/或脂图像投影;
一个完全图像投影计算单元(330),用于根据三点Dixon法计算得到的水图像和脂肪图像,计算相位编码方向的完全水图像投影和/或完全脂肪图像投影,得到完全水和/或脂图像投影;
相关性计算单元(340),计算参考水脂图像投影与完全水脂图像投影之间的相关性,得到至少两个相关值;
结果判定单元(350),用于从所述相关值中获取最大相关值,并确定对应的用三点Dixon法计算出的图像的类型。
9.根据权利要求8所述的装置,其特征在于,该装置进一步包括:
一个归一化单元(370),用于对所述参考图像投影计算单元(320)得到的参考水脂图像投影进行归一化,将归一化后的参考水和/或脂图像投影提供给所述相关性计算单元(340)。
10.根据权利要求8所述的装置,其特征在于,该装置进一步包括:
一个水脂图像计算单元(360),用于利用三点Dixon法计算得到水图像和脂肪图像,将所述水图像和脂肪图像提供给所述完全图像投影计算单元(330)。
11.根据权利要求8所述的装置,其特征在于,所述参考图像投影计算单元(320)包括:一个第一计算子单元(321),用于利用两个反相位回波计算两个反相位之间的相位差;一个第二计算子单元(322),用于利用所述第一计算子单元(321)得到的相位差在相位编码方向计算参考水图像投影和/或参考脂肪图像投影;
和/或,
一个第三计算子单元(323),用于利用两个同相位回波计算两个同相位之间的相位差;一个第四计算子单元(324),用于利用所述第三计算子单元(323)得到的相位差在相位编码方向计算参考水图像投影和/或参考脂肪图像投影。
12.根据权利要求8所述的装置,其特征在于,所述相关性计算单元(340)包括:第一相关性计算子单元(341)、第二相关性计算子单元(342)、第三相关性计算子单元(343)和第四相关性计算子单元(344)中的至少两者;
其中,所述第一相关性计算子单元(341)用于计算参考水图像投影与完全水图像投影的相关性;
所述第二相关性计算子单元(342)用于计算参考水图像投影与完全脂肪图像投影的相关性;
所述第三相关性计算子单元(343)用于计算参考脂肪图像投影与完全水图像投影的相关性;
所述第四相关性计算子单元(344)用于计算参考脂肪图像投影与完全脂肪图像投影的相关性。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010610002 CN102525460B (zh) | 2010-12-29 | 2010-12-29 | 一种磁共振成像水脂图像辨析方法及装置 |
US13/339,660 US8890526B2 (en) | 2010-12-29 | 2011-12-29 | Method and apparatus for making distinction in a magnetic resonance imaging water-fat image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010610002 CN102525460B (zh) | 2010-12-29 | 2010-12-29 | 一种磁共振成像水脂图像辨析方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102525460A CN102525460A (zh) | 2012-07-04 |
CN102525460B true CN102525460B (zh) | 2013-11-06 |
Family
ID=46334426
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010610002 Active CN102525460B (zh) | 2010-12-29 | 2010-12-29 | 一种磁共振成像水脂图像辨析方法及装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US8890526B2 (zh) |
CN (1) | CN102525460B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6446027B2 (ja) * | 2013-04-03 | 2018-12-26 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 高snrの同相画像及び低snrの少なくとも部分的に異相の画像を使用するディクソン型の水/脂肪分離mri |
WO2014195384A1 (en) * | 2013-06-06 | 2014-12-11 | Koninklijke Philips N.V. | Parallel mri with b0 distortion correction and multi-echo dixon water-fat separation using regularised sense reconstruction |
JP6014266B2 (ja) * | 2013-08-07 | 2016-10-25 | 株式会社日立製作所 | 磁気共鳴撮影装置および水脂肪分離方法 |
JP6626507B2 (ja) * | 2014-12-04 | 2019-12-25 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 事前の知識を使うディクソン磁気共鳴撮像 |
DE102015202062A1 (de) * | 2015-02-05 | 2016-08-11 | Siemens Healthcare Gmbh | Rekonstruktion von Magnetresonanzbilddaten für mehrere chemische Substanzarten bei Multi-Echo-Bildgebungsverfahren |
CN107167752B (zh) * | 2017-07-04 | 2020-11-24 | 南京拓谱医疗科技有限公司 | 一种超快速磁共振水脂分离成像方法 |
EP3680680A1 (en) * | 2019-01-11 | 2020-07-15 | Siemens Healthcare GmbH | Method for obtaining an optimised operating parameter, storage medium and magnetic resonance apparatus |
CN113425254B (zh) * | 2021-05-10 | 2022-10-18 | 复旦大学 | 基于混合数据输入的体脂率预测模型的青年男性体脂率预测方法 |
CN114159027B (zh) * | 2021-11-26 | 2024-04-12 | 浙江大学 | 一种变回波个数的磁共振指纹成像方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5909119A (en) * | 1995-08-18 | 1999-06-01 | Toshiba America Mri, Inc. | Method and apparatus for providing separate fat and water MRI images in a single acquisition scan |
US6459922B1 (en) * | 1999-03-30 | 2002-10-01 | Toshiba America Mri, Inc. | Post data-acquisition method for generating water/fat separated MR images having adjustable relaxation contrast |
US6583623B1 (en) * | 2000-03-31 | 2003-06-24 | University Of Rochester | Interleaved water and fat dual-echo spin echo magnetic resonance imaging with intrinsic chemical shift elimination |
US7084626B2 (en) * | 2003-06-30 | 2006-08-01 | Board Of Regents, The University Of Texas System | Methods and apparatuses for fast chemical shift magnetic resonance imaging |
CN100396239C (zh) * | 2005-02-28 | 2008-06-25 | 西门子(中国)有限公司 | 磁共振多通道成像水脂分离重建算法 |
US7777486B2 (en) * | 2007-09-13 | 2010-08-17 | The Board Of Trustees Of The Leland Stanford Junior University | Magnetic resonance imaging with bipolar multi-echo sequences |
CN101843483B (zh) | 2009-03-23 | 2012-07-18 | 西门子(深圳)磁共振有限公司 | 一种水脂分离实现方法和装置 |
-
2010
- 2010-12-29 CN CN 201010610002 patent/CN102525460B/zh active Active
-
2011
- 2011-12-29 US US13/339,660 patent/US8890526B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
US8890526B2 (en) | 2014-11-18 |
US20120169337A1 (en) | 2012-07-05 |
CN102525460A (zh) | 2012-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102525460B (zh) | 一种磁共振成像水脂图像辨析方法及装置 | |
Mori et al. | Introduction to diffusion tensor imaging: And higher order models | |
JP4444413B2 (ja) | 4重フィールドエコーシーケンスを用いて水と脂肪を定量的にmr撮影する装置 | |
Osman et al. | Imaging longitudinal cardiac strain on short‐axis images using strain‐encoded MRI | |
US10768253B2 (en) | MR imaging with signal suppression of a spin series | |
RU2605516C2 (ru) | Магниторезонансная (мр) визуализация, использующая контрастность арт и дискретизацию с множеством времен эха | |
US8680861B1 (en) | System and method for multislice fast magnetic resonance imaging | |
US8030923B2 (en) | Method and system to perform phase correction for species separation in magnetic resonance imaging | |
Epstein et al. | Displacement‐encoded cardiac MRI using cosine and sine modulation to eliminate (CANSEL) artifact‐generating echoes | |
JP5848713B2 (ja) | 磁気共鳴イメージング装置及びコントラスト強調画像取得方法 | |
KR101084470B1 (ko) | 자기 공명 촬영 장치 | |
CN105556326A (zh) | 具有dixon类型的水/脂肪分离的mr成像 | |
JP2000135206A5 (ja) | 4重フィールドエコーシーケンスを用いて水と脂肪を定量的にmr撮影する装置 | |
CN102257399A (zh) | 具有cest对比度增强的mr成像 | |
US9523749B2 (en) | Magnetic resonance imaging of chemical species | |
CN104685368A (zh) | 耐金属mr成像参考扫描 | |
CN103472420B (zh) | 未知空间分布磁场下获取高分辨核磁共振异核谱图的方法 | |
CN105785298A (zh) | 一种高精度三维化学位移成像方法 | |
KR20160145103A (ko) | 낮은 고조파 펄스 시퀀스를 이용한 저-노이즈 자기 공명 이미징 | |
US9581667B2 (en) | Method and magnetic resonance system to implement a multi-echo measurement sequence | |
US20170192070A1 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
JP5004588B2 (ja) | 磁気共鳴イメージング装置 | |
US10114095B2 (en) | Multi-spectral MRI scan with magnetization recovery | |
US11585881B2 (en) | Magnetic resonance imaging of neuro-electro-magnetic oscillations | |
CN112244812B (zh) | 磁共振成像方法、磁共振成像系统和电子装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |