CN106646299A - 梯度回波三回波水脂分离方法及应用它的磁共振成像系统 - Google Patents

梯度回波三回波水脂分离方法及应用它的磁共振成像系统 Download PDF

Info

Publication number
CN106646299A
CN106646299A CN201611257940.4A CN201611257940A CN106646299A CN 106646299 A CN106646299 A CN 106646299A CN 201611257940 A CN201611257940 A CN 201611257940A CN 106646299 A CN106646299 A CN 106646299A
Authority
CN
China
Prior art keywords
echo
phase
water
data
corrected
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
CN201611257940.4A
Other languages
English (en)
Other versions
CN106646299B (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.)
Suzhou Lonwin Medical Systems Co Ltd
Original Assignee
Suzhou Lonwin Medical Systems Co 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 Suzhou Lonwin Medical Systems Co Ltd filed Critical Suzhou Lonwin Medical Systems Co Ltd
Priority to CN201611257940.4A priority Critical patent/CN106646299B/zh
Publication of CN106646299A publication Critical patent/CN106646299A/zh
Application granted granted Critical
Publication of CN106646299B publication Critical patent/CN106646299B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/4828—Resolving the MR signals of different chemical species, e.g. water-fat imaging
    • 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

Landscapes

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

Abstract

本发明公开了一种梯度回波（Gradient Echo,GRE）三回波水脂分离方法及应用它的磁共振成像系统。基于化学位移成像的梯度回波双回波技术成像序列，受到磁场不均匀性影响，导致数据包含额外的误差相位和T2*衰减等问题，从而造成水和脂肪的比例定量不精准，甚至水脂图像分离错误。本发明原始数据使用了GRAPPA（Generalized Auto‑calibrating Partially Parallel Acquisitions）技术提高采集速度；采用梯度回波三回波技术校正图像数据的T2*衰减；同时结合区域生长校正算法，消除额外相位误差带来的相位卷绕，避免了水脂分离错误。三者共同保证了水脂分离图像的高效性、准确性和稳定性。

Description

梯度回波三回波水脂分离方法及应用它的磁共振成像系统
技术领域
本发明属于磁共振成像领域，具体涉及梯度回波三回波水脂分离方法及应用它的磁共振成像系统。
背景技术
磁共振成像（Magnetic resonance imaging，MRI）技术已经成为现代医学诊断中一种常用的技术手段。在MRI中，由于人体组织水和脂肪中的氢质子所处的分子环境不同，使得其共振频率存在一定差异，称之为化学位移。Dixon最早利用该原理获得了MRI中的水图和脂肪图像。假设人体组织只含有水和脂肪两种成分，即人体组织在MRI中只含有两种共振频率，当人体受到射频脉冲激励之后，两种组织的弛豫时间不同，在不同回波时间采集信号，水和脂肪也会表现出信号强度的差异。Dixon方法分别采集水和脂肪的同相位（InPhase，简称IP）和反相位（Out Phase，简称OP）两种回波信号，通过两个回波的数据后处理，得到纯水和纯脂肪的图像。
早期Dixon方法受磁场不均匀性影响较大，成像过程容易受到人体呼吸等运动影响，计算方法复杂且容错性较差。实际应用中存在着系统噪声、涡流、主磁场不均匀性和人体组织之间磁导率变化等都会导致磁场不均匀，如何克服硬件和人为误差带来的影响是Dixon算法的重点。随着水脂分离技术的发展，三点法、不对称三点法和多点法相继提出：其基本原理都是依据三个及三个以上多个回波，计算出磁场不均匀性导致的附加相位，校正各个回波相位之后再求得水和脂肪图像。
三点法的不足在于，附加相位的计算和校正无法避免需要对相位进行去卷绕，现有技术对图像中每个点的附加相位去卷绕速度很慢，且本身在数学上无解，可能存在奇点等导致算法错误。而不对称三点法和其他多点法，不仅一定程度需要更多扫描时间，每点后处理计算也相对更加复杂，增加数据处理时间。
因此，需要提供一种快速准确稳定的水脂分离算法，在高效的扫描时间内，获得正确的水脂图像。
发明内容
本发明目的是：提供梯度回波三回波水脂分离方法及应用它的磁共振成像系统，避免使用相位解卷绕的同时，有效克服磁场不均匀和人为运动等因素带来的数据误差。
本发明的技术方案是：提供一种梯度回波三回波水脂分离方法，应用于磁共振成像系统中，其特征在于，包括以下步骤：
步骤1：采用三维梯度回波三回波序列对磁共振成像区域进行磁共振成像扫描，并对三个回波进行数据采集，同时在采集过程中利用GRAPPA技术加速数据采集；
步骤2：利用GRAPPA技术拟合还原三个回波数据各自K空间内欠采的相位编码数据，然后做初步的数据处理，求得完整的三个回波的图像域数据；
步骤3：将三个回波的图像域数据进行T2*衰减的拟合估计，并对第二个回波进行校正，以消除T2*的衰减影响；
步骤4：利用第一回波数据和校正之后的第二个回波数据，利用区域生长法对前两个回波进行逐个像素点的相位校正，并最终计算获得水和脂肪图像。
进一步的，所述步骤1中对磁共振成像区域使用三维梯度回波序列进行磁共振成像扫描，一次激发后采集三个回波，其中OP的回波时间TE为2.2ms、IP1的回波时间TE为4.4ms、IP2的回波时间TE为8.8ms。
进一步的，所述步骤2的具体步骤如下：
A：将采集到的三个回波的数据分为三组分别进行GRAPPA后处理，拟合补充完整的K空间数据，三个回波的K空间数据分别为SIOP、SIIP1和SIIP1；
B：对三组K空间数据各个通道的信号进行滤波去噪处理，傅里叶变化求得图像域数据之后，进行多通道合成进一步增加信噪比，同时减少后处理运算量。
进一步的，所述步骤3的具体步骤如下：
通过三个回波的图像域数据，拟合每一个像素点的T2*值，并根据公式：IIP1-corrected=IIP1×e(△TE/ T2*)，校正IIP1，其中△TE=2.2ms。
进一步的，所述步骤4具体步骤如下：
a：利用区域生长算法获得正相位IIP1-corrected和反相位IOP图的正确相位差值图Pcorrected，即：以正相位图为参考，IIP1-corrected的相位为0，IIP=abs（IIP1-corrected）；同时需要依据正确相位差值图Pcorrected得到反相位图像IOP=abs(IOP)×ei×Pcorrected；
b：根据如下公式最终计算得到水和脂肪图像：Iwater=(IIP＋IOP)/2，Ifat=(IIP－IOP)/2。
进一步的，所述步骤a的具体步骤如下：
（1）首先计算正相位图IIP1-corrected和反相位图IOP各自的相位图，分别为PIP和POP；
（2）根据PIP和POP每一个像素与相邻像素的相位差异大小，获得正反相位的相位变化差异图QIP和QOP，每个像素点其值越小，说明相位变化越小，该点处的磁场均匀性越好；
（3）根据QIP和QOP综合考量找到磁场均匀性最好的区域作为区域生长的起始位置，即种子点Seed；
（4）以Seed为起始，优先考量磁场均匀性好的区域，每次校正某个像素的相位时，先参考周围相邻已经校正过的像素点相位，以此保证磁场均匀性好的区域的正确相位能够被最优先利用，最终获得整个扫描区域的校正相位图Pcorrected。
本发明为了解决上述问题，还设计了一种磁共振成像系统，其特征在于，应用了上文所述的梯度回波三回波水脂分离方法来获得水脂分离图像。
本发明的优点是：
利用GRAPPA技术提高了水脂分离成像序列的实际扫描速度，减少潜在人为运动带来的伪影；三回波的T2*衰减校正进一步提高了水和脂肪信号比例的准确性；区域生长法对正反相图像数据的相位校正，回避了相位解卷绕这个无解的数学问题，节省大量时间，提高了水脂分离后处理的准确性和稳定性。
附图说明
下面结合附图及实施例对本发明作进一步描述：
图1为本发明的流程示意图；
图2为本发明中扫描序列及回波演变示意简图；
图3 为本发明区域生长过程示意图；
图4 为本发明校正之后的IOP相位图；
图5 为本发明最终计算得到的水和脂肪图像。
具体实施方式
下面结合附图对本发明进行详细描述。
本发明提供一种磁共振多通道成像水脂分离方法，其特征在于，包括以下步骤：
（1）采用GRAPPA加速技术的GRE三回波序列扫描被试。参见图2所示，扫描中每个相位编码过程序列及回波演变。所述GRAPPA技术可参见论文“Griswold M A, Jakob P M,Heidemann R M, et al. Generalized autocalibrating partially parallelacquisitions (GRAPPA)[J]. Magnetic Resonance in Medicine, 2002, 47(6):1202–1210.”；
（2）利用GRAPPA技术拟合还原三个回波各自K空间内欠采的相位编码数据行，然后做初步的数据处理，包括滤波、傅里叶变换和多通道合并；
（3）将完整的三组回波数据进行T2*衰减的拟合估计，并对第二个回波进行校正，消除T2*衰减影响；
（4）使用校正后的两个回波数据，利用区域生长法对前两个回波进行逐个像素点的相位校正。图3所示为区域生长各个像素点相位校正的先后次序，暗色区域为先，高亮区域靠后，图4所示为校正之后IOP的相位图Pcorrected；
（5）根据公式：Iwater=(IIP＋IOP)/2，Ifat=(IIP－IOP)/2，计算水和脂肪图像，图5为最终图像示例。
当然上述实施例只为说明本发明的技术构思及特点，其目的在于让熟悉此项技术的人能够了解本发明的内容并据以实施，并不能以此限制本发明的保护范围。凡根据本发明主要技术方案的精神实质所做的修饰，都应涵盖在本发明的保护范围之内。

Claims (7)

1.一种梯度回波三回波水脂分离方法，应用于磁共振成像系统中，其特征在于，包括以下步骤：
步骤1：采用三维梯度回波三回波序列对磁共振成像区域进行磁共振成像扫描，并对三个回波进行数据采集，同时在采集过程中利用GRAPPA技术加速数据采集；
步骤2：利用GRAPPA技术拟合还原三个回波数据各自K空间内欠采的相位编码数据，然后做初步的数据处理，求得完整的三个回波的图像域数据；
步骤3：将三个回波的图像域数据进行T2*衰减的拟合估计，并对第二个回波进行校正，以消除T2*的衰减影响；
步骤4：利用第一回波数据和校正之后的第二个回波数据，利用区域生长法对前两个回波进行逐个像素点的相位校正，并最终计算获得水和脂肪图像。
2.根据权利要求1所述的梯度回波三回波水脂分离方法，其特征在于，所述步骤1中对磁共振成像区域使用三维梯度回波序列进行磁共振成像扫描，一次激发后采集三个回波，其中OP的回波时间TE为2.2ms、IP1的回波时间TE为4.4ms、IP2的回波时间TE为8.8ms。
3.根据权利要求1所述的梯度回波三回波水脂分离方法，其特征在于，所述步骤2的具体步骤如下：
A：将采集到的三个回波的数据分为三组分别进行GRAPPA后处理，拟合补充完整的K空间数据，三个回波的K空间数据分别为SIOP、SIIP1和SIIP1；
B：对三组K空间数据各个通道的信号进行滤波去噪处理，傅里叶变化求得图像域数据之后，进行多通道合成进一步增加信噪比，同时减少后处理运算量。
4.根据权利要求1所述的梯度回波三回波水脂分离方法，其特征在于，所述步骤3的具体步骤如下：
通过三个回波的图像域数据，拟合每一个像素点的T2*值，并根据公式：IIP1-corrected=IIP1×e(△TE/ T2*)，校正IIP1，其中△TE=2.2ms。
5.根据权利要求1所述的梯度回波三回波水脂分离方法，其特征在于，所述步骤4具体步骤如下：
a：利用区域生长算法获得正相位IIP1-corrected和反相位IOP图的正确相位差值图Pcorrected，即：以正相位图为参考，IIP1-corrected的相位为0，IIP=abs（IIP1-corrected）；同时需要依据正确相位差值图Pcorrected得到反相位图像IOP=abs(IOP)×ei×Pcorrected；
b：根据如下公式最终计算得到水和脂肪图像：Iwater=(IIP＋IOP)/2，Ifat=(IIP－IOP)/2。
6.根据权利要求5所述的梯度回波三回波水脂分离方法，其特征在于，所述步骤a的具体步骤如下：
（1）首先计算正相位图IIP1-corrected和反相位图IOP各自的相位图，分别为PIP和POP；
（2）根据PIP和POP每一个像素与相邻像素的相位差异大小，获得正反相位的相位变化差异图QIP和QOP，每个像素点其值越小，说明相位变化越小，该点处的磁场均匀性越好；
（3）根据QIP和QOP综合考量找到磁场均匀性最好的区域作为区域生长的起始位置，即种子点Seed；
（4）以Seed为起始，优先考量磁场均匀性好的区域，每次校正某个像素的相位时，先参考周围相邻已经校正过的像素点相位，以此保证磁场均匀性好的区域的正确相位能够被最优先利用，最终获得整个扫描区域的校正相位图Pcorrected。
7.一种磁共振成像系统，其特征在于，应用了权利要求1-6中任一种的梯度回波三回波水脂分离方法来获得水脂分离图像。
CN201611257940.4A 2016-12-30 2016-12-30 梯度回波三回波水脂分离方法及应用它的磁共振成像系统 Active CN106646299B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611257940.4A CN106646299B (zh) 2016-12-30 2016-12-30 梯度回波三回波水脂分离方法及应用它的磁共振成像系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611257940.4A CN106646299B (zh) 2016-12-30 2016-12-30 梯度回波三回波水脂分离方法及应用它的磁共振成像系统

Publications (2)

Publication Number Publication Date
CN106646299A true CN106646299A (zh) 2017-05-10
CN106646299B CN106646299B (zh) 2019-09-10

Family

ID=58837750

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611257940.4A Active CN106646299B (zh) 2016-12-30 2016-12-30 梯度回波三回波水脂分离方法及应用它的磁共振成像系统

Country Status (1)

Country Link
CN (1) CN106646299B (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 苏州朗润医疗系统有限公司 一种梯度回波多回波水脂分离方法及应用该方法的磁共振成像系统
CN112669406A (zh) * 2020-12-31 2021-04-16 苏州朗润医疗系统有限公司 一种基于投影估计的相控阵线圈磁共振图像非均匀性校正方法
CN112763955A (zh) * 2020-12-31 2021-05-07 苏州朗润医疗系统有限公司 一种基于图割算法的磁共振图像水脂分离方法
CN113419202A (zh) * 2020-12-29 2021-09-21 苏州朗润医疗系统有限公司 一种采用Radial 3DTOF获取颈动脉磁共振血管图像的方法及磁共振成像系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103513202A (zh) * 2012-06-16 2014-01-15 上海联影医疗科技有限公司 一种磁共振成像中的dixon水脂分离方法
CN104918546A (zh) * 2013-02-12 2015-09-16 株式会社日立医疗器械 磁共振成像装置及其处理方法
US20160187447A1 (en) * 2014-12-29 2016-06-30 General Electric Company Method and apparatus for separating chemical species in magnetic resonance imaging

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103513202A (zh) * 2012-06-16 2014-01-15 上海联影医疗科技有限公司 一种磁共振成像中的dixon水脂分离方法
CN104918546A (zh) * 2013-02-12 2015-09-16 株式会社日立医疗器械 磁共振成像装置及其处理方法
US20160187447A1 (en) * 2014-12-29 2016-06-30 General Electric Company Method and apparatus for separating chemical species in magnetic resonance imaging

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ESIN OZTURK 等: "Partially Parallel MR Spectroscopic Imaging of Gliomas at 3T", 《PROCEEDINGS OF THE 28TH IEEE EMBS ANNUAL INTERNATIONAL CONFERENCE》 *
MARK A. GRISWOLD 等: "Generalized Autocalibrating Partially Parallel Acquisitions (GRAPPA)", 《MAGNETIC RESONANCE IN MEDICINE》 *
许林 等: "基于各向异性扩散的GRAPPA重建算法", 《计算机工程》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108720834A (zh) * 2018-02-06 2018-11-02 苏州朗润医疗系统有限公司 一种梯度回波多回波水脂分离方法及应用该方法的磁共振成像系统
CN108720834B (zh) * 2018-02-06 2021-06-01 苏州朗润医疗系统有限公司 一种梯度回波多回波水脂分离方法及应用该方法的磁共振成像系统
CN113419202A (zh) * 2020-12-29 2021-09-21 苏州朗润医疗系统有限公司 一种采用Radial 3DTOF获取颈动脉磁共振血管图像的方法及磁共振成像系统
CN113419202B (zh) * 2020-12-29 2022-06-21 苏州朗润医疗系统有限公司 一种采用Radial 3DTOF获取颈动脉磁共振血管图像的方法及磁共振成像系统
CN112669406A (zh) * 2020-12-31 2021-04-16 苏州朗润医疗系统有限公司 一种基于投影估计的相控阵线圈磁共振图像非均匀性校正方法
CN112763955A (zh) * 2020-12-31 2021-05-07 苏州朗润医疗系统有限公司 一种基于图割算法的磁共振图像水脂分离方法

Also Published As

Publication number Publication date
CN106646299B (zh) 2019-09-10

Similar Documents

Publication Publication Date Title
CN102870000B (zh) 使用单个和多个通道接收器线圈的用于同时多切片磁共振成象的方法
US7847546B2 (en) Magnetic resonance imaging apparatus
US8023705B2 (en) Method for reconstructing image from echo planar imaging sequence
US8154294B2 (en) Magnetic resonance imaging apparatus and method for correcting error due to gradient magnetic field
JP5502682B2 (ja) 磁気共鳴イメージング装置
US9050018B2 (en) Means and methods for providing high resolution MRI
CN106646299A (zh) 梯度回波三回波水脂分离方法及应用它的磁共振成像系统
CN104068859A (zh) 用于确定多个磁共振图像的方法和磁共振设备
CN107997763A (zh) 一种水脂分离磁共振成像方法
US20110274331A1 (en) Magnetic resonance imaging method for achieving water-fat separation
US10254370B2 (en) Modified pulse sequence for magnetic resonance imaging using MRI-visible markers for motion correction
US20080218169A1 (en) Methods for fat signal suppression in magnetic resonance imaging
CN102973271B (zh) 磁共振动态成像方法及系统
CN108720834A (zh) 一种梯度回波多回波水脂分离方法及应用该方法的磁共振成像系统
JP2012090957A (ja) 磁気共鳴イメージング装置
CN103747728B (zh) 磁共振成像装置以及磁共振成像方法
CN107076819A (zh) 具有对流伪影的抑制的Dixon MR成像
JP4047553B2 (ja) 磁気共鳴イメージング装置
JP2017529960A (ja) アーチファクト抑制を有するプロペラｍｒイメージング
CN103890602A (zh) 利用定制信号激励模块进行快速mri采集（rate）的方法与系统
US10175328B2 (en) System and method for reconstructing ghost-free images from data acquired using simultaneous multislice magnetic resonance imaging
US7626388B2 (en) Magnetic resonance imager
CN106137198B (zh) 一种磁共振成像方法及装置
CN102841329A (zh) 磁共振信号处理方法及装置
CN104918546A (zh) 磁共振成像装置及其处理方法

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