CN104739409A - 磁共振成像方法和装置 - Google Patents

磁共振成像方法和装置 Download PDF

Info

Publication number
CN104739409A
CN104739409A CN201310753743.1A CN201310753743A CN104739409A CN 104739409 A CN104739409 A CN 104739409A CN 201310753743 A CN201310753743 A CN 201310753743A CN 104739409 A CN104739409 A CN 104739409A
Authority
CN
China
Prior art keywords
image
fat
mrow
msub
space data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201310753743.1A
Other languages
English (en)
Other versions
CN104739409B (zh
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.)
Siemens Shenzhen Magnetic Resonance Ltd
Original Assignee
Siemens Shenzhen Magnetic Resonance Ltd
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 Siemens Shenzhen Magnetic Resonance Ltd filed Critical Siemens Shenzhen Magnetic Resonance Ltd
Priority to CN201310753743.1A priority Critical patent/CN104739409B/zh
Priority to US14/587,069 priority patent/US10042024B2/en
Publication of CN104739409A publication Critical patent/CN104739409A/zh
Application granted granted Critical
Publication of CN104739409B publication Critical patent/CN104739409B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/4828Resolving the MR signals of different chemical species, e.g. water-fat imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7221Determining signal validity, reliability or quality

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Veterinary Medicine (AREA)
  • Signal Processing (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Psychiatry (AREA)
  • Physiology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)

Abstract

本发明是一种磁共振成像方法和装置。其中,所述磁共振成像方法包括如下步骤:进行扫描序列获得一正相位图像和一反相位图像;修正所述正相位图像和所述反相位图像的磁场误差;根据磁共振图像信号的波谱的多个脂肪峰值利用所述正相位图像和所述反相位图像重建水图像和脂图像。利用磁共振图像信号的波谱的多个脂肪峰值重建水图像和脂图像,可以减少化学位移产生的伪影。

Description

磁共振成像方法和装置
技术领域
本发明涉及磁共振成像方法和装置,特别涉及DIXON技术的磁共振成像方法和装置。
背景技术
磁共振成像是随着计算机技术、电子电路技术、超导体技术的发展而迅速发展起来的一种生物磁学核自旋成像技术。在磁共振成像中,人体组织被置于静磁场B0中,随后用频率与氢原子核的进动频率相同的射频脉冲激发人体组织内的氢原子核,引起氢原子核共振,并吸收能量;在停止射频脉冲后,氢原子核按特定频率发出射电信号,并将吸收的能量释放出来,由体外的接受器收录,经计算机处理后获得图像。
在磁共振成像中,为了获得较好的成像质量,往往需要对例如脂肪信号、水信号、硅胶信号(乳房植入物)等特定谱成分的信号进行抑制。例如,在用于腹部和胸部等的磁共振成像检查中,普遍需要抑制脂肪信号,以突出感兴趣组织或病灶的显示。至今为止,已经提出了许多脂肪抑制技术,例如CHESS(CHEmical Shift Suppression,化学位移选择)技术,FatSat(fat saturation,脂肪饱和)技术,SPAIR(Spectral PresaturationAttenuated Inversion Recovery,频谱预饱和衰减反转恢复)技术,STIR(short inversion timeinversion recovery,短反转时间反转恢复)技术,DIXON技术等。
其中,在CHESS技术失效的情况下,DIXON技术是一种广为接受的具有自防故障的方法。与CHESS技术比较,DIXON技术得到的水图像具有较多的脂肪残余信号;而脂肪残余信号在T2加权图像(脂肪组织和非脂肪组织之间的对比度很低)中是个棘手的技术问题,T2是横向弛豫时间。
现有技术的DIXON技术具有以下技术特征:
s ( t ) = ( ρ w + ρ f Σ p = 1 P e j 2 πf p t ) · e j 2 πψt - - - [ 1 ]
具体而言,式[1]是DIXON技术的信号模型:其中,磁共振图像信号s(t)包括脂信号ρf和水信号ρw;同时,在磁共振图像信号的波谱中包括P个脂信号峰值;由于化学位移,脂信号的各个峰值经其子身频率调制,该频率反映了由于本地磁场的不均匀性导致的相位偏移。
发明内容
有鉴于此,本发明提出了一种磁共振成像方法,其特征在于,包括如下步骤:进行扫描序列获得一正相位图像和一反相位图像;修正所述正相位图像和所述反相位图像的磁场误差;根据磁共振图像信号的波谱的多个脂肪峰值利用经修正的所述正相位图像和经修正的所述反相位图像重建水图像和脂图像。
优选地,所述扫描序列是一多回波扫描序列。
优选地,所述多回波扫描序列包括:一多回波TSE序列、一多回波mDIXON序列或一多回波vibe DIXON序列。
优选地,所述根据磁共振图像信号的波谱的多个脂肪峰值利用经修正的所述正相位图像和经修正的所述反相位图像重建水图像和脂图像的步骤包括,利用一第一通式和一第二通式根据磁共振图像信号的波谱的多个脂肪峰值利用经修正的所述正相位图像和经修正的所述反相位图像重建水图像和脂图像,其中,
所述第一通式是:
W F = 1 1 1 c - 1 S in S op
其中,Sin是经修正的所述同相位图像,Sop是经修正的所述反相位图像,W是所述水图像,F是所述脂图像,
所述第二通式是:
c = Σ n = 1 N α n e - R n * t op + iΔ f n t op
其中,αn是各个脂肪峰值所占所述多个脂肪峰值的信号强度的比例,△fn是频率偏移,是衰变率,top是反相位图像的回波时间,N是所述多个脂肪峰值的数量。
优选地,所述根据磁共振图像信号的波谱的多个脂肪峰值利用经修正的所述正相位图像和经修正的所述反相位图像重建水图像和脂图像的步骤包括:
将经修正的所述正相位图像和经修正的所述反相位图像进行傅里叶变换得到一正相位图像的K空间数据和一反相位图像的K空间数据;
根据磁共振图像信号的波谱的多个脂肪峰值利用所述正相位图像的K空间数据和所述反相位图像的K空间数据得出一水图像的K空间数据和一脂图像的K空间数据;
对所述水图像的K空间数据和所述脂图像的K空间数据进行傅里叶逆比变换得出所述水图像和所述脂图像。
优选地,所述根据磁共振图像信号的波谱的多个脂肪峰值利用所述正相位图像的K空间数据和所述反相位图像的K空间数据得出一水图像的K空间数据和一脂图像的K空间数据的步骤包括,利用一第三通式和一第四通式,根据磁共振图像信号的波谱的多个脂肪峰值利用所述正相位图像的K空间数据和所述反相位图像的K空间数据得出一水图像的K空间数据和一脂图像的K空间数据,其中,所述第三通式是:
W KSPACE ( t ) F KSPACE ( t ) = 1 1 1 c ( t ) - 1 KS in ( t ) KS op ( t )
其中,KSin(t)是所述正相位图像的K空间数据,KSop(t)是所述反相位图像的K空间数据,WKSPACE(t)是所述水图像的K空间数据,FKSPACE(t)是所述脂图像的K空间数据;
其中,所述第四通式是:
c ( t ) = Σ n = 1 N α n e - R n * t op + iΔ f n t op
其中,αn是各个脂肪峰值所占所述多个脂肪峰值的信号强度的比例,△fn是频率偏移,是衰变率,top是反相位图像的回波时间,N是所述多个脂肪峰值的数量。
优选地,所述扫描序列是一基本DIXON序列。
优选地,所述根据磁共振图像信号的波谱的多个脂肪峰值利用经修正的所述正相位图像和经修正的所述反相位图像重建水图像和脂图像的步骤包括,利用一第五通式得出一未经修正的水图像和一未经修正的脂图像,
其中,所述第五通式是:
W BDM F BDM = 0.5 0.5 0.5 - 0.5 S in S op
其中,WBDM是所述未经修正的水图像,FBDM是所述未经修正的脂图像,
利用一第六通式对所述未经修正的水图像和所述未经修正的脂图像FBDM进行修正从而得出一最终的水图像和一最终的脂图像,
其中,所述第六通式是:
W BDMF F BDMF = K W BDM F BDM
其中,WBDMF是所述最终的水图像,FBDMF是所述最终的脂图像,
K = 1 1 1 c - 1 1 1 1 - 1
其中, c = Σ n = 1 N α n e - R n * t op + iΔ f n t op ,
其中,αn是各个脂肪峰值所占所述多个脂肪峰值的信号强度的比例,△fn是频率偏移,是衰变率,top是反相位图像的回波时间,N是所述多个脂肪峰值的数量。
本发明还提供一种磁共振成像装置,其特征在于,包括如下模块:
一扫描模块,用于进行扫描序列获得一正相位图像和一反相位图像;
一修正模块,用于修正所述正相位图像和所述反相位图像的磁场误差;
一重建模块,用于根据磁共振图像信号的波谱的多个脂肪峰值利用经修正的所述正相位图像和所述反相位图像重建水图像和脂图像。
优选地,所述扫描序列是一多回波扫描序列。
优选地,所述多回波扫描序列包括:一多回波TSE序列、一多回波mDIXON序列或一多回波vibe DIXON序列。
优选地,所述重建模块,用于利用一第一通式和一第二通式根据磁共振图像信号的波谱的多个脂肪峰值利用经修正的所述正相位图像和所述反相位图像重建水图像和脂图像,其中,
所述第一通式是:
W F = 1 1 1 c - 1 S in S op
其中,Sin是经修正的所述同相位图像,Sop是经修正的所述反相位图像,W是所述水图像,F是所述脂图像,
所述第二通式是:
c = Σ n = 1 N α n e - R n * t op + iΔ f n t op
其中,αn是各个脂肪峰值所占所述多个脂肪峰值的信号强度的比例,△fn是频率偏移,是衰变率,top是反相位图像的回波时间,N是所述多个脂肪峰值的数量。
优选地,所述重建模块,用于:
将经修正的所述正相位图像和经修正的所述反相位图像进行傅里叶变换得到一正相位图像的K空间数据和一反相位图像的K空间数据;
根据磁共振图像信号的波谱的多个脂肪峰值利用所述正相位图像的K空间数据和所述反相位图像的K空间数据得出一水图像的K空间数据和一脂图像的K空间数据;
对所述水图像的K空间数据和所述脂图像的K空间数据进行傅里叶逆比变换得出所述水图像和所述脂图像。
优选地,所述重建模块,用于利用一第三通式和一第四通式,根据磁共振图像信号的波谱的多个脂肪峰值利用所述正相位图像的K空间数据和所述反相位图像的K空间数据得出一水图像的K空间数据和一脂图像的K空间数据,其中,所述第三通式是:
W KSPACE ( t ) F KSPACE ( t ) = 1 1 1 c ( t ) - 1 KS in ( t ) KS op ( t )
其中,KSin(t)是所述正相位图像的K空间数据,KSop(t)是所述反相位图像的K空间数据,WKSPACE(t)是所述水图像的K空间数据,FKSPACE(t)是所述脂图像的K空间数据;
其中,所述第四通式是:
c ( t ) = Σ n = 1 N α n e - R n * t op + iΔ f n t op
其中,αn是各个脂肪峰值所占所述多个脂肪峰值的信号强度的比例,△fn是频率偏移,是衰变率,top是反相位图像的回波时间,N是所述多个脂肪峰值的数量。
优选地,所述扫描序列是一基本DIXON序列。
优选地,所述重建模块,用于利用一第五通式得出一未经修正的水图像和一未经修正的脂图像,
其中,所述第五通式是:
W BDM F BDM = 0.5 0.5 0.5 - 0.5 S in S op
其中,WBDM是所述未经修正的水图像,FBDM是所述未经修正的脂图像,
利用一第六通式对所述未经修正的水图像和所述未经修正的脂图像FBDM进行修正从而得出一最终的水图像和一最终的脂图像,
其中,所述第六通式是:
W BDMF F BDMF = K W BDM F BDM
其中,WBDMF是所述最终的水图像,FBDMF是所述最终的脂图像,
K = 1 1 1 c - 1 1 1 1 - 1
其中, c = Σ n = 1 N α n e - R n * t op + iΔ f n t op ,
其中,αn是各个脂肪峰值所占所述多个脂肪峰值的信号强度的比例,△fn是频率偏移,是衰变率,top是反相位图像的回波时间,N是所述多个脂肪峰值的数量。
从上述方案中可以看出,利用磁共振图像信号的波谱的多个脂肪峰值重建水图像和脂图像,可以减少化学位移产生的伪影。以多回波TSE序列为例,在正相位回波信号中,由于补偿了CPMG自旋回波磁振脉冲序列的频率偏移和衰变效果,所以修正仅影响反相位图像,因此通过伪逆运算可以重建出真正的水图像和脂图像。
附图说明
下面将通过参照附图详细描述本发明的优选实施例,使本领域的普通技术人员更清楚本发明的上述及其它特征和优点,附图中:
图1为脂肪组织的磁共振波谱示意图。
图2是根据本发明的第一具体实施例的磁共振成像方法的示意图。
图3是根据本发明的第二具体实施例的磁共振成像方法的示意图。
图4是根据本发明的第三具体实施例的磁共振成像方法的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,以下举实施例对本发明进一步详细说明。
根据本发明的磁共振成像方法的实质在于:对现有技术的DIXON技术的信号模型进行了校正,并将该信号模型用于DIXON技术,从而得到更好的水脂对比度。
DIXON技术具有以下技术特征:
s ( t ) = ( ρ w + ρ f Σ p = 1 P e j 2 πf p t ) · e j 2 πψt - - - [ 1 ]
具体而言,式[1]是DIXON技术的信号模型:磁共振图像信号s(t)包括脂信号ρf和水信号ρw;同时,在磁共振图像信号的波谱中包括P个脂信号峰值;由于化学位移,脂信号的各个峰值经其子身频率调制,该频率反映了由于本地磁场的不均匀性导致的相位偏移。
图1是脂肪组织的磁共振图像信号的波谱示意图,其中脂肪组织的磁共振波谱包括十个脂峰值A、B、C、D、E、F、G、H、I和J,以及一个水峰值water,其中,不同脂肪峰就是脂肪里不同种类的物质的磁共振信号,由图1可见,经验值exp和拟合值fit相差无几。现有技术的DIXON技术默认水频谱和脂频谱都只具有一个峰值,及最大脂肪峰值B,也就是在式[1]中P=1。通过不同的回波时间t得到两个回波信号,现有技术的DIXON技术即可利用式[1]计算出脂信号ρf和水信号ρw从而重建脂图像和水图像。由此可见,现有技术的DIXON技术忽略了最大峰值B之外的其他峰值A、C、D、E、F、G、H、I和J,从而导致脂信号在现有技术的DIXON技术中被低估,因此在水图像中可以观测到大量残留的脂信号:在T2加权像中,脂肪组织和非脂肪组织之间的低对比度问题尤其严重。
第一具体实施例
图2是根据本发明的第一具体实施例的磁共振成像方法的示意图。如图2所示,根据本发明的第一具体实施例的磁共振成像方法包括如下步骤:
步骤S101,进行多回波扫描序列获得一正相位图像和一反相位图像。
其中,多回波扫描序列包括:多回波TSE序列、多回波mDIXON序列或多回波vibeDIXON序列等等。
步骤S102:修正所述正相位图像和所述反相位图像的磁场误差。
根据本发明的第一具体实施例的磁共振成像方法利用中国专利申请200510009873.0所公开的方法对上述多回波信号的磁场误差进行修正;除此之外,还可以通过其他多种方式对上述多回波信号的磁场误差进行修正,例如,两点DIXON法,三点DIXON法等等。具体而言:
(1)取得一幅同相位图像和两幅反相位图像。
假设分别表示i通道的同相位图像和两反相位图像,其中i=1…n,n为总通道数。
(2)求各通道的线圈灵敏度分布(Profile)。
本步骤的求解过程可以分为以下几步:
(a)求各通道同相位图像模值平方和,然后开方:
SOS ( x , y ) = Σ i = 1 n | S 0 i ( x , y ) | 2
(b)取同相位图像中一个信噪比较好的通道信号对模值进行归一化:
S 0 ′ ( x , y ) = S 0 I ( x , y ) | S 0 I ( x , y ) |
在本发明的另外一个实施例中,(b)步骤的运算也可以为:取各通道信号的一种线性组合,对模值进行归一化:
S 0 ′ ( x , y ) = Σ i = 1 n S 0 i ( x , y ) · W i | Σ i = 1 n S 0 i ( x , y ) · W i | , 式中Wi表示权重。
(c)计算各通道的线圈灵敏度分布(Profile):
P i ( x , y ) = S 0 i ( x , y ) · [ S 0 ′ ( x , y ) ] · SOS ( x , y ) , 式中[]·表示求共轭。
由步骤(a)~(c)实现的求各通道的线圈灵敏度分布的过程还可以由其它优化算法代替,比如修正空间匹配滤波器法,即利用各通道图像每个象素点及其邻域点信号,计算信号及噪声相关矩阵,通过求信号相关矩阵的特征值及特征向量,从而得出各通道线圈的灵敏度分布。
(3)合成各通道图像。
S 0 ( x , y ) = Σ i = 1 n [ P i ( x , y ) ] · · S 0 i ( x , y )
S 1 ( x , y ) = Σ i = 1 n [ P i ( x , y ) ] · · S 1 i ( x , y )
S 2 ( x , y ) = Σ i = 1 n [ P i ( x , y ) ] · · S 2 i ( x , y )
本步骤把各通道的图像合成得到最优信噪比合成图像,同时保存了图像的相位。
(4)求两幅反相位图像的相位差。
式中angle[]表示取相角算子。
(5)使用相位反缠绕算法对进行相位反缠绕,得到修正后的相位
检测同相位图像中的一些特征区域以作为修正相位的判据。
在本步骤中,可利用大部分图像中脂肪信号比较高的特点,检测图像的脂肪,或者用边缘检测方法,检测皮下脂肪,利用这些脂肪点做为判据,修正相位
(6)修正反相位图像的相位。
步骤S103:根据磁共振图像信号的波谱的多个脂肪峰值利用所述正相位图像和所述反相位图像重建水图像和脂图像。
具体而言,各个脂肪峰值都具有已知的共振频率和T2数值;而且各个脂肪峰值所占比例也可以认为是固定的。因此,根据步骤S102,分别表示i通道的同相位图像和两反相位图像,如果选取作为经修正的同相位图像Sin,选取之中任一作为经修正的反相位图像Sop,那么通过式[2]重建水图像W和脂图像F。
W F = 1 1 1 c - 1 S in S op - - - [ 2 ]
其中,通过式[3]得出式[2]中的c。
c = Σ n = 1 N α n e - R n * t op + iΔ f n t op - - - [ 3 ]
其中,αn是各个脂肪峰值所占所述多个脂肪峰值的信号强度的比例,△fn是频率偏移,是衰变率,top是反相位图像的TE时间,N是所述多个脂肪峰值的数量,在根据本发明的第一具体实施例中N是10。在Composition of adipose tissue and marrow fat in humans by1H NMR at7Tesla J.Lipid Res.2008.49:2055-2062中对上述数据有相应公开。
化学位移伪影是水图像和脂图像之间的对位偏差导致的伪影,这种伪影源自水图像和脂图像具有不同的共振频率(由于在读出阶段脂信号的相位会发生变化,所以会产生相位位移)。在根据本发明的第一具体实施例的磁共振成像方法中,利用磁共振图像信号的波谱的多个脂肪峰值重建水图像和脂图像,可以减少化学位移产生的伪影。以多回波TSE序列为例,在正相位回波信号中,由于补偿了CPMG自旋回波磁振脉冲序列的频率偏移和衰变效果,所以修正仅影响反相位图像,因此通过伪逆运算可以重建出真正的水图像和脂图像。
第二具体实施例
图3是根据本发明的第二具体实施例的磁共振成像方法的示意图。如图3所示,根据本发明的第二具体实施例的磁共振成像方法包括如下步骤:
步骤S201,进行多回波扫描序列获得一正相位图像和一反相位图像。
其中,多回波扫描序列包括:两点DIXON序列、多回波TSE序列、多回波mDIXON序列,多回波vibe DIXON序列等等。
步骤S202:修正所述正相位图像和所述反相位图像的磁场误差。
步骤S203:对经修正的所述正相位图像和经修正的所述反相位图像进行傅里叶变换得到一正相位图像的K空间数据和一反相位图像的K空间数据。
根据步骤S102,分别表示i通道的同相位图像和两反相位图像,如果选取作为经修正的同相位图像Sin,选取之中任一作为经修正的反相位图像Sop,那么通过傅里叶变换得到一正相位图像的K空间数据KSin(t)和一反相位图像的K空间数据KSop(t)。
步骤S204:根据磁共振图像信号的波谱的多个脂肪峰值利用所述正相位图像的K空间数据KSin(t)和所述反相位图像的K空间数据KSop(t)得出一水图像的K空间数据WKSPACE(t)和一脂图像的K空间数据FKSPACE(t)。
具体而言,通过式[4]得出一水图像的K空间数据WKSPACE(t)和一脂图像的K空间数据FKSPACE(t)。
W KSPACE ( t ) F KSPACE ( t ) = 1 1 1 c ( t ) - 1 KS in ( t ) KS op ( t ) - - - [ 4 ]
其中,并且通过式[5]得出式[4]中的c(t)。
c ( t ) = Σ n = 1 N α n e - R n * t op + iΔ f n t op - - - [ 5 ]
其中,αn是各个脂肪峰值所占所述多个脂肪峰值的信号强度的比例,△fn是频率偏移,是衰变率,top是反相位图的TE时间,N是所述多个脂肪峰值的数量,在根据本发明的第二具体实施例中N是10。
化学位移伪影是水图像和脂图像之间的对位偏差导致的伪影,这种伪影源自水图像和脂图像具有不同的共振频率(由于在读出阶段脂信号的相位会发生变化,所以会产生相位位移)。在根据本发明的第二具体实施例的磁共振成像方法中,利用磁共振图像信号的波谱的多个脂肪峰值重建水图像和脂图像,可以减少化学位移产生的伪影。
步骤S205:对所述水图像的K空间数据和所述脂图像的K空间数据进行傅里叶逆比变换得出所述水图像和所述脂图像。
第三具体实施例
图4是根据本发明的第三具体实施例的磁共振成像方法的示意图。如图4所示,根据本发明的第三具体实施例的磁共振成像方法包括如下步骤:
步骤S301:进行DIXON序列得到一正相位图像和一反相位图像。
步骤S302:修正所述正相位图像和所述反相位图像的磁场误差。
步骤S303:利用式[6]得出一水图像WBDM和一脂图像FBDM,所述水图像WBDM和所述脂图像FBDM未经修正。
W BDM F BDM = 0.5 0.5 0.5 - 0.5 S in S op - - - [ 6 ]
其中Sin是经修正的同相位图像,Sop是经修正的反相位图像,WBDM是所述水图像,FBDM是所述脂图像。
步骤S304:利用式[7]对所述水图像WBDM和所述脂图像FBDM进行修正从而得出最终的经修正的水图像WBDMF和经修正的脂图像FBDMF
W BDMF F BDMF = K W BDM F BDM - - - [ 7 ]
其中,利用式[8]得出K。
K = 1 1 1 c - 1 1 1 1 - 1 - - - [ 8 ]
其中,利用式[3]得出 c = Σ n = 1 N α n e - R n * t op + iΔ f n t op .
化学位移伪影是水图像和脂图像之间的对位偏差导致的伪影,这种伪影源自水图像和脂图像具有不同的共振频率(由于在读出阶段脂信号的相位会发生变化,所以会产生相位位移)。在根据本发明的第三具体实施例的磁共振成像方法中,利用磁共振图像信号的波谱的多个脂肪峰值重建水图像和脂图像,可以减少化学位移产生的伪影。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (16)

1.一种磁共振成像方法,其特征在于,包括如下步骤:
进行扫描序列获得一正相位图像和一反相位图像;
修正所述正相位图像和所述反相位图像的磁场误差;
根据磁共振图像信号的波谱的多个脂肪峰值利用经修正的所述正相位图像和经修正的所述反相位图像重建水图像和脂图像。
2.如权利要求1所述的磁共振成像方法,其特征在于,所述扫描序列是一多回波扫描序列。
3.如权利要求2所述的磁共振成像方法,其特征在于,所述多回波扫描序列包括:一多回波TSE序列、一多回波mDIXON序列或一多回波vibe DIXON序列。
4.如权利要求2所述的磁共振成像方法,其特征在于,所述根据磁共振图像信号的波谱的多个脂肪峰值利用经修正的所述正相位图像和经修正的所述反相位图像重建水图像和脂图像的步骤包括,利用一第一通式和一第二通式根据磁共振图像信号的波谱的多个脂肪峰值利用经修正的所述正相位图像和经修正的所述反相位图像重建水图像和脂图像,其中,
所述第一通式是:
W F = 1 1 1 c - 1 S in S op
其中,Sin是经修正的所述同相位图像,Sop是经修正的所述反相位图像,W是所述水图像,F是所述脂图像,
所述第二通式是:
c = Σ n = 1 N α n e - R n * t op + iΔ f n t op
其中,αn是各个脂肪峰值所占所述多个脂肪峰值的信号强度的比例,△fn是频率偏移,是衰变率,top是反相位图像的回波时间,N是所述多个脂肪峰值的数量。
5.如权利要求2所述的磁共振成像方法,其特征在于,所述根据磁共振图像信号的波谱的多个脂肪峰值利用经修正的所述正相位图像和经修正的所述反相位图像重建水图像和脂图像的步骤包括:
将经修正的所述正相位图像和经修正的所述反相位图像进行傅里叶变换得到一正相位图像的K空间数据和一反相位图像的K空间数据;
根据磁共振图像信号的波谱的多个脂肪峰值利用所述正相位图像的K空间数据和所述反相位图像的K空间数据得出一水图像的K空间数据和一脂图像的K空间数据;
对所述水图像的K空间数据和所述脂图像的K空间数据进行傅里叶逆比变换得出所述水图像和所述脂图像。
6.如权利要求5所述的磁共振成像方法,其特征在于,所述根据磁共振图像信号的波谱的多个脂肪峰值利用所述正相位图像的K空间数据和所述反相位图像的K空间数据得出一水图像的K空间数据和一脂图像的K空间数据的步骤包括,利用一第三通式和一第四通式,根据磁共振图像信号的波谱的多个脂肪峰值利用所述正相位图像的K空间数据和所述反相位图像的K空间数据得出一水图像的K空间数据和一脂图像的K空间数据,其中,所述第三通式是:
W KSPACE ( t ) F KSPACE ( t ) = 1 1 1 c ( t ) - 1 KS in ( t ) KS op ( t )
其中,KSin(t)是所述正相位图像的K空间数据,KSop(t)是所述反相位图像的K空间数据,WKSPACE(t)是所述水图像的K空间数据,FKSPACE(t)是所述脂图像的K空间数据;
其中,所述第四通式是:
c ( t ) = Σ n = 1 N α n e - R n * t op + iΔ f n t op
其中,αn是各个脂肪峰值所占所述多个脂肪峰值的信号强度的比例,△fn是频率偏移,是衰变率,top是反相位图像的回波时间,N是所述多个脂肪峰值的数量。
7.如权利要求1所述的磁共振成像方法,其特征在于,所述扫描序列是一基本DIXON序列。
8.如权利要求7所述的磁共振成像方法,其特征在于,所述根据磁共振图像信号的波谱的多个脂肪峰值利用经修正的所述正相位图像和经修正的所述反相位图像重建水图像和脂图像的步骤包括,利用一第五通式得出一未经修正的水图像和一未经修正的脂图像,
其中,所述第五通式是:
W BDM F BDM = 0.5 0.5 0.5 - 0.5 S in S op
其中,WBDM是所述未经修正的水图像,FBDM是所述未经修正的脂图像,
利用一第六通式对所述未经修正的水图像和所述未经修正的脂图像FBDM进行修正从而得出一最终的水图像和一最终的脂图像,
其中,所述第六通式是:
W BDMF F BDMF = K W BDM F BDM
其中,WBDMF是所述最终的水图像,FBDMF是所述最终的脂图像,
K = 1 1 1 c - 1 1 1 1 - 1
其中, c = Σ n = 1 N α n e - R n * t op + iΔ f n t op ,
其中,αn是各个脂肪峰值所占所述多个脂肪峰值的信号强度的比例,△fn是频率偏移,是衰变率,top是反相位图像的回波时间,N是所述多个脂肪峰值的数量。
9.一种磁共振成像装置,其特征在于,包括如下模块:
一扫描模块,用于进行扫描序列获得一正相位图像和一反相位图像;
一修正模块,用于修正所述正相位图像和所述反相位图像的磁场误差;
一重建模块,用于根据磁共振图像信号的波谱的多个脂肪峰值利用经修正的所述正相位图像和所述反相位图像重建水图像和脂图像。
10.如权利要求9所述的磁共振成像装置,其特征在于,所述扫描序列是一多回波扫描序列。
11.如权利要求10所述的磁共振成像装置,其特征在于,所述多回波扫描序列包括:一多回波TSE序列、一多回波mDIXON序列或一多回波vibe DIXON序列。
12.如权利要求10所述的磁共振成像装置,其特征在于,所述重建模块,用于利用一第一通式和一第二通式根据磁共振图像信号的波谱的多个脂肪峰值利用经修正的所述正相位图像和所述反相位图像重建水图像和脂图像,其中,
所述第一通式是:
W F = 1 1 1 c - 1 S in S op
其中,Sin是经修正的所述同相位图像,Sop是经修正的所述反相位图像,W是所述水图像,F是所述脂图像,
所述第二通式是:
c = Σ n = 1 N α n e - R n * t op + iΔ f n t op
其中,αn是各个脂肪峰值所占所述多个脂肪峰值的信号强度的比例,△fn是频率偏移,是衰变率,top是反相位图像的回波时间,N是所述多个脂肪峰值的数量。
13.如权利要求10所述的磁共振成像装置,其特征在于,所述重建模块,用于:
将经修正的所述正相位图像和经修正的所述反相位图像进行傅里叶变换得到一正相位图像的K空间数据和一反相位图像的K空间数据;
根据磁共振图像信号的波谱的多个脂肪峰值利用所述正相位图像的K空间数据和所述反相位图像的K空间数据得出一水图像的K空间数据和一脂图像的K空间数据;
对所述水图像的K空间数据和所述脂图像的K空间数据进行傅里叶逆比变换得出所述水图像和所述脂图像。
14.如权利要求13所述的磁共振成像装置,其特征在于,所述重建模块,用于利用一第三通式和一第四通式,根据磁共振图像信号的波谱的多个脂肪峰值利用所述正相位图像的K空间数据和所述反相位图像的K空间数据得出一水图像的K空间数据和一脂图像的K空间数据,其中,所述第三通式是:
W KSPACE ( t ) F KSPACE ( t ) = 1 1 1 c ( t ) - 1 KS in ( t ) KS op ( t )
其中,KSin(t)是所述正相位图像的K空间数据,KSop(t)是所述反相位图像的K空间数据,WKSPACE(t)是所述水图像的K空间数据,FKSPACE(t)是所述脂图像的K空间数据;
其中,所述第四通式是:
c ( t ) = Σ n = 1 N α n e - R n * t op + iΔ f n t op
其中,αn是各个脂肪峰值所占所述多个脂肪峰值的信号强度的比例,△fn是频率偏移,是衰变率,top是反相位图像的回波时间,N是所述多个脂肪峰值的数量。
15.如权利要求9所述的磁共振成像装置,其特征在于,所述扫描序列是一基本DIXON序列。
16.如权利要求15所述的磁共振成像装置,其特征在于,所述重建模块,用于利用一第五通式得出一未经修正的水图像和一未经修正的脂图像,
其中,所述第五通式是:
W BDM F BDM = 0.5 0.5 0.5 - 0.5 S in S op
其中,WBDM是所述未经修正的水图像,FBDM是所述未经修正的脂图像,
利用一第六通式对所述未经修正的水图像和所述未经修正的脂图像FBDM进行修正从而得出一最终的水图像和一最终的脂图像,
其中,所述第六通式是:
W BDMF F BDMF = K W BDM F BDM
其中,WBDMF是所述最终的水图像,FBDMF是所述最终的脂图像,
K = 1 1 1 c - 1 1 1 1 - 1
其中, c = Σ n = 1 N α n e - R n * t op + iΔ f n t op ,
其中,αn是各个脂肪峰值所占所述多个脂肪峰值的信号强度的比例,△fn是频率偏移,是衰变率,top是反相位图像的回波时间,N是所述多个脂肪峰值的数量。
CN201310753743.1A 2013-12-31 2013-12-31 磁共振成像方法和装置 Active CN104739409B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201310753743.1A CN104739409B (zh) 2013-12-31 2013-12-31 磁共振成像方法和装置
US14/587,069 US10042024B2 (en) 2013-12-31 2014-12-31 Magnetic resonance imaging method and apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310753743.1A CN104739409B (zh) 2013-12-31 2013-12-31 磁共振成像方法和装置

Publications (2)

Publication Number Publication Date
CN104739409A true CN104739409A (zh) 2015-07-01
CN104739409B CN104739409B (zh) 2018-02-13

Family

ID=53481411

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310753743.1A Active CN104739409B (zh) 2013-12-31 2013-12-31 磁共振成像方法和装置

Country Status (2)

Country Link
US (1) US10042024B2 (zh)
CN (1) CN104739409B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108720834A (zh) * 2018-02-06 2018-11-02 苏州朗润医疗系统有限公司 一种梯度回波多回波水脂分离方法及应用该方法的磁共振成像系统
WO2019119506A1 (zh) * 2017-12-20 2019-06-27 深圳先进技术研究院 磁共振温度成像方法与装置
CN110471014A (zh) * 2018-05-10 2019-11-19 西门子医疗有限公司 磁共振成像中的条纹伪影减少
CN110501664A (zh) * 2019-08-23 2019-11-26 苏州朗润医疗系统有限公司 磁共振波谱相位校正方法及应用该方法的磁共振成像系统

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US12007458B2 (en) * 2019-07-08 2024-06-11 Canon Medical Systems Corporation Magnetic resonance imaging apparatus and magnetic resonance imaging method
CN114155234A (zh) * 2021-12-08 2022-03-08 沈阳东软智能医疗科技研究院有限公司 病灶肺段位置的识别方法、装置、存储介质及电子设备

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050033153A1 (en) * 2003-04-25 2005-02-10 Case Westem Reserve University Dixon Techniques in spiral trajectories with off-resonance correction
US20050122105A1 (en) * 2003-11-17 2005-06-09 Toshiba America Mri, Inc. Water-fat separation for fast spin echo imaging in an inhomogeneous field with progressive encoding
CN1827038A (zh) * 2005-02-28 2006-09-06 西门子(中国)有限公司 磁共振多通道成像水脂分离重建算法
US20080048659A1 (en) * 2006-08-28 2008-02-28 Scott Brian Reeder Generalized Method for MRI Chemical Species Separation Using Arbitrary k-space Trajectories
US20090261823A1 (en) * 2008-04-17 2009-10-22 Huanzhou Yu Method for water-fat separation and t2* estimation in an mri system employing multiple fat spectral peaks and fat spectrum self-calibration
CN102232831A (zh) * 2010-04-30 2011-11-09 西门子(深圳)磁共振有限公司 一种实现水脂分离的磁共振成像方法
CN102841329A (zh) * 2011-06-24 2012-12-26 西门子(深圳)磁共振有限公司 磁共振信号处理方法及装置
CN102859387A (zh) * 2010-02-22 2013-01-02 皇家飞利浦电子股份有限公司 采用谱模型的化学物类的磁共振成像
CN103257333A (zh) * 2012-02-17 2013-08-21 西门子(深圳)磁共振有限公司 一种磁共振成像中的水脂分离成像方法及装置
CN103608693A (zh) * 2011-04-21 2014-02-26 皇家飞利浦有限公司 具有用于脂肪抑制的化学位移编码的对比度增强磁共振血管造影

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6091243A (en) * 1997-11-13 2000-07-18 The University Of British Columbia Water-fat imaging with direct phase encoding (DPE)
WO2005055136A2 (en) * 2003-11-26 2005-06-16 Board Of Regents, The University Of Texas System Method and apparatus for phase-sensitive magnetic resonance imaging
US20070098298A1 (en) * 2005-11-02 2007-05-03 The University Of British Columbia Imaging methods, apparatus, systems, media and signals
CN101172036B (zh) * 2006-11-02 2011-08-03 西门子公司 相位校正的方法
US8761464B2 (en) * 2007-07-13 2014-06-24 Boad of Regents, The University of Texas System Methods of efficient and improved phase-sensitive MRI
JP5675362B2 (ja) * 2007-11-07 2015-02-25 ゼネラル・エレクトリック・カンパニイ 磁気共鳴撮像における脂肪/水の分離を強化するためのシステム、方法及びマシン読み取り可能プログラム
US8373415B2 (en) * 2010-04-15 2013-02-12 Wisconsin Alumni Research Foundation Method for separating magnetic resonance imaging signals using spectral distinction of species
US9449386B2 (en) * 2012-02-28 2016-09-20 Board Of Regents, The University Of Texas System Method and apparatus for extended phase correction in phase sensitive magnetic resonance imaging
DE102012223789B4 (de) * 2012-12-19 2014-07-17 Siemens Aktiengesellschaft Ermittlung einer Phasendifferenz-Karte
US9519040B2 (en) * 2014-04-04 2016-12-13 Board Of Regents, The University Of Texas System Methods and apparatus for phase correction in magnetic resonance imaging
WO2015190508A1 (ja) * 2014-06-12 2015-12-17 株式会社 日立メディコ 磁気共鳴イメージング装置及び水脂肪分離画像作成方法
DE102014226052A1 (de) * 2014-12-16 2016-06-16 Siemens Healthcare Gmbh Ermittlung einer phasenbeschreibenden Karte mit Hilfe eines min-cut/max-flow-Verfahrens

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050033153A1 (en) * 2003-04-25 2005-02-10 Case Westem Reserve University Dixon Techniques in spiral trajectories with off-resonance correction
US20050122105A1 (en) * 2003-11-17 2005-06-09 Toshiba America Mri, Inc. Water-fat separation for fast spin echo imaging in an inhomogeneous field with progressive encoding
CN1827038A (zh) * 2005-02-28 2006-09-06 西门子(中国)有限公司 磁共振多通道成像水脂分离重建算法
US20080048659A1 (en) * 2006-08-28 2008-02-28 Scott Brian Reeder Generalized Method for MRI Chemical Species Separation Using Arbitrary k-space Trajectories
US20090261823A1 (en) * 2008-04-17 2009-10-22 Huanzhou Yu Method for water-fat separation and t2* estimation in an mri system employing multiple fat spectral peaks and fat spectrum self-calibration
CN102859387A (zh) * 2010-02-22 2013-01-02 皇家飞利浦电子股份有限公司 采用谱模型的化学物类的磁共振成像
CN102232831A (zh) * 2010-04-30 2011-11-09 西门子(深圳)磁共振有限公司 一种实现水脂分离的磁共振成像方法
CN103608693A (zh) * 2011-04-21 2014-02-26 皇家飞利浦有限公司 具有用于脂肪抑制的化学位移编码的对比度增强磁共振血管造影
CN102841329A (zh) * 2011-06-24 2012-12-26 西门子(深圳)磁共振有限公司 磁共振信号处理方法及装置
CN103257333A (zh) * 2012-02-17 2013-08-21 西门子(深圳)磁共振有限公司 一种磁共振成像中的水脂分离成像方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HUANZHOU YU,ET AL.,: "Multiecho Water-Fat Separation and Simultaneous R2* Estimation With Multifrequency Fat Spectrum Modeling", 《MAGNETIC RESONANCE IN MEDICINE》 *
JIMIN REN,ET AL.,: "Composition of adipose tissue and marrow fat in humans by 1H NMR at 7 Tesla", 《JOURNAL OF LIPID RESEARCH》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019119506A1 (zh) * 2017-12-20 2019-06-27 深圳先进技术研究院 磁共振温度成像方法与装置
CN109938704A (zh) * 2017-12-20 2019-06-28 深圳先进技术研究院 磁共振温度成像方法与装置
CN109938704B (zh) * 2017-12-20 2020-05-19 深圳先进技术研究院 磁共振温度成像方法与装置
CN108720834A (zh) * 2018-02-06 2018-11-02 苏州朗润医疗系统有限公司 一种梯度回波多回波水脂分离方法及应用该方法的磁共振成像系统
CN108720834B (zh) * 2018-02-06 2021-06-01 苏州朗润医疗系统有限公司 一种梯度回波多回波水脂分离方法及应用该方法的磁共振成像系统
CN110471014A (zh) * 2018-05-10 2019-11-19 西门子医疗有限公司 磁共振成像中的条纹伪影减少
US11435419B2 (en) 2018-05-10 2022-09-06 Siemens Healthcare Gmbh Streak artifact reduction in magnetic resonance imaging
CN110471014B (zh) * 2018-05-10 2023-01-24 西门子医疗有限公司 磁共振成像中的条纹伪影减少
CN110501664A (zh) * 2019-08-23 2019-11-26 苏州朗润医疗系统有限公司 磁共振波谱相位校正方法及应用该方法的磁共振成像系统

Also Published As

Publication number Publication date
US20150185304A1 (en) 2015-07-02
CN104739409B (zh) 2018-02-13
US10042024B2 (en) 2018-08-07

Similar Documents

Publication Publication Date Title
Zaitsev et al. Point spread function mapping with parallel imaging techniques and high acceleration factors: Fast, robust, and flexible method for echo‐planar imaging distortion correction
CN104739409B (zh) 磁共振成像方法和装置
US8143890B2 (en) Spectral resolution enhancement of magnetic resonance spectroscopic imaging
US10317493B1 (en) System and method for multislice fast magnetic resonance imaging
CN103238082B (zh) 使用多点Dixon技术和低分辨率校准的MR成像
JP6243522B2 (ja) 正則化された検出再構成を使用するマルチエコーディクソン水−脂肪分離及びb0歪補正による並列mri
US9983284B2 (en) MRI with dixon-type water/fat separation and prior knowledge about inhomogeneity of the main magnetic field
US20110260726A1 (en) Phase labeling using sensitivity encoding: data acquisition and image reconstruction for geometric distortion correction in epi
US9753109B2 (en) Magnetic resonance imaging of chemical species with a spectral model
US8890526B2 (en) Method and apparatus for making distinction in a magnetic resonance imaging water-fat image
US10262385B2 (en) Method and apparatus for acquiring magnetic resonance data and entering the data into k-space
US20080211499A1 (en) Low Power Decoupling for Multi-Nuclear Spectroscopy
US9523749B2 (en) Magnetic resonance imaging of chemical species
CN104685368A (zh) 耐金属mr成像参考扫描
US20160161580A1 (en) Magnetic resonance imaging device and water-fat separation method
US8217656B2 (en) Motion compensation in magnetic resonance imaging
JP6814325B2 (ja) ディクソン型水/脂肪分離mrイメージング
US9535148B2 (en) Dynamic contrast enhanced magnetic resonance imaging with high spatial-temporal resolution
WO2009135167A1 (en) Phase labeling using sensitivity encoding: data acquisition and image reconstruction for geometric distortion correction in epi.
JP7128972B2 (ja) 歪み補正を伴うepi mr画像
CN111164444B (zh) 具有经改进的脂肪位移校正的Dixon型水/脂肪分离MR成像
RU2773486C2 (ru) Магнитно-резонансная визуализация с разделением вода/жир методом диксона
Bhattacharya Pushing the limits of spectroscopic imaging using novel low-rank based reconstruction algorithm
Granlund et al. Multi-frequency off-resonance correction and water/fat separation for spiral breast imaging
Pawar Compressive parallel MRI for accelerated imaging

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant