CN113970717B - 磁敏感加权成像方法、装置及磁共振成像系统 - Google Patents

磁敏感加权成像方法、装置及磁共振成像系统 Download PDF

Info

Publication number
CN113970717B
CN113970717B CN202010721447.3A CN202010721447A CN113970717B CN 113970717 B CN113970717 B CN 113970717B CN 202010721447 A CN202010721447 A CN 202010721447A CN 113970717 B CN113970717 B CN 113970717B
Authority
CN
China
Prior art keywords
echo
excitation
acquired
central
acquisition
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
Application number
CN202010721447.3A
Other languages
English (en)
Other versions
CN113970717A (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 CN202010721447.3A priority Critical patent/CN113970717B/zh
Priority to US17/383,563 priority patent/US11726158B2/en
Publication of CN113970717A publication Critical patent/CN113970717A/zh
Application granted granted Critical
Publication of CN113970717B publication Critical patent/CN113970717B/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/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/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5616—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI
    • 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/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5619—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences by temporal sharing of data, e.g. keyhole, block regional interpolation scheme for k-Space [BRISK]
    • 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/50—NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences
    • 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
    • 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
    • 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56536—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to magnetic susceptibility variations
    • 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56554—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by acquiring plural, differently encoded echo signals after one RF excitation, e.g. correction for readout gradients of alternating polarity in EPI
    • 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/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/482—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory
    • G01R33/4822—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory in three dimensions

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General 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)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明实施例公开了磁敏感加权成像方法、装置及磁共振系统。方法包括：在多次激发的平面回波成像中，针对每次激发的回波进行分层编码、相位编码和频率编码方向的流动补偿；确定每次激发后当采用线性重排序方式时将要采集的多个回波，针对每次激发，从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波，且对本次激发的回波采集方向与对上一次激发的回波采集方向相反；对采集到的回波进行磁敏感加权成像。本发明实施例降低了基于平面回波序列的磁敏感加权成像图像中的流动伪影。

Description

磁敏感加权成像方法、装置及磁共振成像系统
技术领域
本发明涉及MRI(Magnetic Resonance Imaging，磁共振成像)技术领域，特别是一种SWI(Susceptibility-Weighted Imaging，磁敏感加权成像)方法、装置及MRI系统。
背景技术
MRI通过对静磁场中的人体施加某种特定频率的射频(Radio Frequency，RF)脉冲，使人体中的氢质子受到激发而发生MR(Magnetic Resonance，磁共振)现象，停止脉冲后，质子在弛豫过程中产生MR信号，通过对MR信号的接收、空间编码和图像重建等处理过程，即产生MR影像。
成像过程中，在RF激发之后，首先通过Z方向梯度完成选层，即选择成像的一个层；接着打开Y方向的梯度磁场，这样不同Y位置的磁矩的进动速度也会不同，再关闭Y梯度，这样各个位置的磁矩速度又恢复成一样，但是由于之前进动速度不同造成的不同Y位置的相位偏移不同，该过程称为相位编码；接下来，打开X方向的梯度，不同X位置的磁矩的速度不同，该过程称为频率编码。至此，通过相位编码和频率编码，一个2D图像的各个位置可以被确定下来，接收线圈采集到的信号是k空间信号，通过傅里叶变换即可得到图像；再通过Z方向梯度完成一轮新的选层，再重复上述过程，即可得到3D的MRI图像。
通过90°RF和180°RF的组合可以构成不同的RF。常用的MRI RF有：SE(Spin-Echosequence，自旋回波序列)、TSE(Turbo Spin-Echo sequence，快速自旋回波序列)、GRE(Gradient-Echo sequence，梯度回波序列)和EPI(Echo Planar Imaging，回波平面成像)。其中，GRE和EPI都属于梯度回波。梯度回波的产生是在一次RF激发后，利用读出梯度场方向(即X方向)正反向切换产生一个梯度回波。EPI与一般梯度回波不同的是在一次RF激发后，利用读出梯度场的连续正反向切换，每次切换产生一个梯度回波，因此有回波链的存在，与TSE类似。
在磁共振成像中，静止组织在经历了大小相同、方向相反的选层梯度(即Z方向梯度)后，获得的相位累积为零；而对于沿着读出梯度场方向(即X方向)移动的运动组织如流动的血液、脑脊液等的相位累积则不为零，接下来的相位编码梯度(即Y方向梯度)施加时，已经有相位的运动组织会被错误编码，出现在其他位置，成为流动伪影。
FC(Flow Compensation，流动补偿)也叫GMN(Gradient Moment Nulling，梯度力矩消除)，是利用特殊设计的梯度场来减少或消除流动伪影的技术。FC技术的梯度组合模式有很多种。通过多次不同面积的正、反向梯度场的变换，各种速度流体的相位偏移最终都能接近于零，从而消除运动伪影。一般情况下FC是通过消除一阶梯度矩来补偿匀速流动带来的影响，文中的GMN也指的是一阶矩补偿。
SWI(Susceptibility-Weighted Imaging，磁敏感加权成像)是一种磁共振对比增强成像技术，以T2*加权梯度回波序列作为序列基础，利用采集到的MR信号的相位信息，为具有磁化率差异的组织提供增强的图像对比度，已广泛应用于各种临床环境下的含铁组织、静脉血管的评价。
目前的一种方案是通过全流量补偿T2*加权3D-GRE序列获取SWI图像，该方法对每一个回波都做全流动补偿，即对每一个回波都在选层方向(即Z方向)、相位编码方向(即Y方向)和频率编码方向(即X方向，也称读出梯度场方向)做流动补偿，因此通常需要较长的扫描时间。
为了提高扫描效率，目前提出了一种3D-多次激发平面回波成像(iEPI，interleaved Echo Planar Imaging)技术作为3D-GRE的快速替代方案。其通过使用较短的EPI回波链长度，典型的EPI相关伪影(失真和模糊)是有限的，并且图像将在信噪比和效率上都获得增益，同时保持与传统3D-GRE相比的幅度和相位图像的相似对比度。然而，与GRE相比，EPI的流动补偿设计更为复杂。因为与GRE一次激发采集一个回波所不同，EPI序列是一次激发采集多个回波，这样在做流动补偿时很难做到三个方向的全流补偿。除非是采用回扫法(Flyback)对频率编码方向上的读出梯度极性相同的回波都进行补偿，但这又会牺牲掉扫描效率。
因此，在将GMN应用于EPI回波链时，为了充分利用GMN的扫描效率，目前提出了一种基于线性重排序的3D-iEPI方案，该方案在选层梯度的流动补偿实现与3D-GRE相同的情况下，对所有方向即分层编码、相位编码和频率编码方向仅补偿每次激发的中心回波，然后采用线性重排序方式采集流动补偿后的回波，来实现更快的SWI，提高信噪比效率。但是对该方案实现的SWI图像进行分析发现，图像中仍存在明显的动脉流动污染。
发明内容
有鉴于此，本发明实施例一方面提出了一种SWI方法，另一方面提出了一种SWI装置及MRI系统，以降低流动效应对基于平面回波序列的SWI图像的影响。
一种磁敏感加权成像方法，包括：
在多次激发的平面回波成像中，针对每次激发的回波进行分层编码、相位编码和频率编码方向的流动补偿；
确定每次激发后当采用线性重排序方式时将要采集的多个回波，针对每次激发，从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波，且对本次激发的回波采集方向与对上一次激发的回波采集方向相反；
对采集到的回波进行磁敏感加权成像。
所述对本次激发的回波采集方向与对上一次激发的回波采集方向相反包括：
当对上一次激发的回波采集方向为：从中心回波开始向k空间正向采集到相当于以线性重排序方式采集时的最后一个回波时，对本次激发的回波采集方向为：从中心回波开始向k空间负向采集到相当于以线性重排序方式采集时的第一个回波；或者，
当对上一次激发的回波采集方向为：从中心回波开始向k空间负向采集到相当于以线性重排序方式采集时的第一个回波时，对本次激发的回波采集方向为：从中心回波开始向k空间正向采集到相当于以线性重排序方式采集时的最后一个回波。
所述从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波包括：
当本次激发的选层编码梯度脉冲和相位编码梯度脉冲发射完毕时，即从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波。
所述从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波之后、对采集到的回波进行磁敏感加权成像之前，进一步包括：
将第m次和第m+1次采集的两个中心回波叠加后取平均，将得到的平均数据作为第m次和第m+1次采集的中心回波，1≤m<M，M为本次平面回波成像的激发总次数。
所述针对每次激发的回波进行分层编码、相位编码和频率编码方向的流动补偿包括：
针对每次激发的中心回波进行分层编码和相位编码方向的流动补偿，所述中心回波为每次激发后当采用线性重排序方式时将要采集的多个回波的中心回波；
针对每次激发的第奇数个回波进行频率编码方向的流动补偿。
所述针对每次激发的中心回波进行分层编码和相位编码方向的流动补偿包括：
按照如下公式进行分层编码和相位编码方向的流动补偿：
M1,par＝M0,parΔtpar
M1,phase＝M0,phaseprephaseΔtp
其中，M1,par为分层编码方向在中心回波处的一阶矩，M1,phase为相位编码方向在中心回波处的一阶矩；M0,par为分层编码方向上分层编码梯度的零阶矩，M0,phaseprephase为相位编码方向上预相位梯度的零阶矩，M0,par、M0,phaseprephase在不同的激发中是不同的；Δtpar为分层编码方向上分层编码梯度中心到中心回波的时间，Δtp为相位编码方向上预相位梯度中心到中心回波的时间。
一种磁敏感加权成像装置，包括：
流动补偿模块(51)：在多次激发的平面回波成像中，针对每次激发的回波进行分层编码、相位编码和频率编码方向的流动补偿；
回波采集模块(52)：确定每次激发后当采用线性重排序方式时将要采集的多个回波，针对每次激发，从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波，且对本次激发的回波采集方向与对上一次激发的回波采集方向相反；
加权成像模块(53)：用于对采集到的回波进行磁敏感加权成像。
所述回波采集模块(52)对本次激发的回波采集方向与对上一次激发的回波采集方向相反包括：
当对上一次激发的回波采集方向为：从中心回波开始向k空间正向采集到相当于以线性重排序方式采集时的最后一个回波时，对本次激发的回波采集方向为：从中心回波开始向k空间负向采集到相当于以线性重排序方式采集时的第一个回波；或者，
当对上一次激发的回波采集方向为：从中心回波开始向k空间负向采集到相当于以线性重排序方式采集时的第一个回波时，对本次激发的回波采集方向为：从中心回波开始向k空间正向采集到相当于以线性重排序方式采集时的最后一个回波。
所述回波采集模块(52)从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波包括：
当本次激发的选层编码梯度脉冲和相位编码梯度脉冲发射完毕时，即从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波。
所述流动补偿模块(51)针对每次激发的回波进行分层编码、相位编码和频率编码方向的流动补偿包括：
针对每次激发的中心回波进行分层编码和相位编码方向的流动补偿，所述中心回波为每次激发后当采用线性重排序方式时将要采集的多个回波的中心回波；
针对每次激发的第奇数个回波进行频率编码方向的流动补偿。
一种磁共振成像系统，包括如上任一项所述的磁敏感加权成像装置。
本发明实施例中，通过将采用线性重排序方式时一次激发采集的多个回波拆分成采用中心重排序方式时的2次相邻激发来采集，从而通过两次相邻激发采集得到和采用线性重排序方式一次激发采集得到的一样的完整回波，且因为每次激发后从中心回波开始采集，从而降低了流动效应对SWI图像的影响，且由于相邻两次激发过程都采集了一次中心回波，提高了SWI图像的信噪比。
附图说明
下面将通过参照附图详细描述本发明的优选实施例，使本领域的普通技术人员更清楚本发明的上述及其它特征和优点，附图中：
图1为本发明一实施例提供的SWI方法的流程图；
图2为本发明一应用示例中在相邻两次激发过程(第m次和第m+1次)中采用的RF和梯度脉冲序列的示意图；
图3为对采用图2所示的RF和梯度脉冲序列进行激发后，回波采集过程在K空间的轨迹图；
图4给出了用于3DGRE、基于线性重排序和基于中心重排序的3DiEPI的SWI处理的mIP图像的比较示意图；
图5为本发明一实施例提供的SWI装置的结构示意图。
其中，附图标记如下：
标号 含义
101～103 步骤
51 流动补偿模块
52 回波采集模块
53 加权成像模块
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，以下举实施例对本发明进一步详细说明。
图1为本发明一实施例提供的SWI方法的流程图，其具体步骤如下：
步骤101：在多次激发的平面回波成像中，针对每次激发的回波进行分层编码、相位编码和频率编码方向的流动补偿。
一可选实施例中，针对每次激发的中心回波进行分层编码和相位编码方向的流动补偿，中心回波为每次激发后当采用线性重排序方式时将要采集的多个回波的中心回波；针对每次激发的第奇数个(即第1、3、5、7、...个)回波进行频率编码方向的流动补偿。
其中，可按照如下公式进行分层编码和相位编码方向的流动补偿：
M1,par＝M0,parΔtpar
M1,phase＝M0,phaseprephaseΔtp
其中，M1,par为分层编码方向在中心回波处的一阶矩，M1,phase为相位编码方向在中心回波处的一阶矩；M0,par为分层编码方向上分层编码梯度的零阶矩，M0,phaseprephase为相位编码方向上预相位梯度的零阶矩，M0,par、M0,phaseprephase在不同的激发中是不同的；Δtpar为分层编码方向上分层编码梯度中心到中心回波的时间，Δtp为相位编码方向上预相位梯度中心到中心回波的时间。本发明实施例中，回波链中的第一个回波即为当采用线性重排序方式采集时的中心回波。
线性重排序方式的回波采集属于已有成熟技术，其具体过程为：每次激发后向k空间的正向或者负向采集多个回波得到部分k空间数据，多次激发采集得到的K空间数据在Y方向上以交错的方式结合起来，形成完整的K空间数据。例如：每次激发后向k空间的正向采集n个回波(n>1)，激发次数为m，则m次激发共得到m*n个回波，该m*n个回波在k空间的Y方向正向上的排列顺序为：
第1次激发采集的第1个回波、第2次激发采集的第1个回波、第3次激发采集的第1个回波、...、第m-1次激发采集的第1个回波、第m次激发采集的第1个回波；第1次激发采集的第2个回波、第2次激发采集的第2个回波、第3次激发采集的第2个回波、...、第m-1次激发采集的第2个回波、第m次激发采集的第2个回波；......；第1次激发采集的第n-1个回波、第2次激发采集的第n-1个回波、第3次激发采集的第n-1个回波、...、第m-1次激发采集的第n-1个回波、第m次激发采集的第n-1个回波；第1次激发采集的第n个回波、第2次激发采集的第n个回波、第3次激发采集的第n个回波、...、第m-1次激发采集的第n个回波、第m次激发采集的第n个回波。
本发明方案是对线性重排序方式的改进，本发明方案中的中心回波指的是：假设每次激发后以线性重排序方式采集多个回波时，该多个回波的中心回波。例如：每次激发后以线性重排序方式采集9个回波，则第5个回波即为中心回波。
步骤102：确定每次激发后当采用线性重排序方式时将要采集的多个回波，针对每次激发，从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波，且对本次激发的回波采集方向与对上一次激发的回波采集方向相反。
其中，对本次激发的回波采集方向与对上一次激发的回波采集方向相反可以是如下两种之一：
一、当对上一次激发的回波采集方向为：从中心回波开始向k空间正向采集到相当于以线性重排序方式采集时的最后一个回波时，对本次激发的回波采集方向为：从中心回波开始向k空间负向采集到相当于以线性重排序方式采集时的第一个回波；
二、当对上一次激发的回波采集方向为：从中心回波开始向k空间负向采集到相当于以线性重排序方式采集时的第一个回波时，对本次激发的回波采集方向为：从中心回波开始向k空间正向采集到相当于以线性重排序方式采集时的最后一个回波。
步骤103：对采集到的回波进行磁敏感加权成像。
由于相邻两次激发的回波采集过程都采集了中心回波，因此，在进行磁敏感加权成像之前，可将第m次和第m+1次采集的两个中心回波叠加后取平均，将得到的平均数据作为第m次和第m+1次采集的中心回波，且1≤m<M，M为本次平面回波成像的激发总次数。
上述实施例中，通过将采用线性重排序方式时一次激发采集的多个回波拆分成采用中心重排序方式时的2次相邻激发来采集，从而通过两次相邻激发采集得到和采用线性重排序方式一次激发采集得到的一样的完整回波，且因为每次激发后从中心回波开始采集，从而降低了流动效应对SWI图像的影响，且由于相邻两次激发过程都采集了一次中心回波，提高了SWI图像的信噪比。
为了进一步降低流动效应对SWI图像的影响，一较佳实施例中，在每次激发的选层编码梯度脉冲和相位编码梯度脉冲发射完毕时，即进行从中心回波开始向k空间正向或者负向采集每一回波。也就是说，回波采集过程在选层编码梯度脉冲和相位编码梯度脉冲发射完毕后立即开始。
本发明实施例中，可限制每次激发在采用线性重排序方式时采集的回波总数(即EPI因子)为奇数，这样中心回波就是中间那个回波，例如：确定每次激发后当采用线性重排序方式时将要采集的回波数为9，则中心回波就是第5个回波。此外，需要说明的是，本发明实施例中采用中心重排序的方式采集回波时，需要保证相邻两次激发采集的对称性，即相邻两次激发需要从中心回波开始分别向k空间的正向和负向对称采集相同数目的回波。
图2为本发明一应用示例中在相邻两次激发过程(第m次和第m+1次)中采用的RF和梯度脉冲序列的示意图。可见，针对每次激发，在分层编码PA、相位编码PE和频率编码FE方向均施加有流动补偿梯度Sg、Pg、Fg。本示例中一次激发当采用线性重排序方式采集时采集的回波数为9，则其中心回波EC，即第5个回波为本发明实施例中在每次激发后采集的第一个回波。
图3为对采用图2所示的RF和梯度脉冲序列进行激发后，回波采集过程在K空间的轨迹图。如图3所示，一次激发后若采用线性重排序方式要采集9个回波，则本发明实施例中采用中心重排序方法后，将拆分成2次激发，即在第m次激发中，从中心回波a5开始，依次采集a6、a7、a8和a9，在第m+1次激发中，从中心回波b5开始，依次采集b4、b3、b2、b1。
为了评估本发明实施例提供的基于中心重排序的3D-iEPI方案降低流动效应影响的效果，在一个例子中，采用一台配置有20通道的头颈线圈的商用的1.5T扫描仪，分别采用3D-GRE、基于线性重排序的3D-iEPI和本发明实施例提供的基于中心重排序的3D-iEPI方案对一健康志愿者进行试验，具体成像参数见表1。
表1
其中，在基于线性重排序的3D-iEPI和本发明方案提供的基于中心重排序的3D-iEPI方案中调节TR和翻转角度以匹配3D-GRE中的背景对比度。采集回波数据后，对3D-GRE、基于线性重排序的3D-iEPI数据和基于中心重排序的3D-iEPI数据进行标准化SWI处理。
图4给出了用于3DGRE、基于线性重排序和基于中心重排序的3DiEPI的SWI处理的mIP图像的比较示意图。其中，a对应的是有流动补偿的3D-GRE方案，b对应的无流动补偿的3D-GRE方案，c对应的是有流动补偿和线性重排序的3D-iEPI方案，d对应的是有流动补偿和中心重排序的3D-iEPI方案。
从图4中的分图a、b、c、d可以看出：
一、与未采用流动补偿的图b相比，采用流动补偿的图a、c和d的图像质量提高了，血管描绘清晰，且图a、c、d的图像质量相似。
二、将图d与图a和图c相比，可以看出：采用中心重排序的图d能够进一步减轻血流引起的动脉污染，更好地刻画血管。
三、此外，在采用3D-iEPI(图c、图d)的图像中显示了更多的小血管，这可能是由于高频信号的TE增加。
综上，本发明实施例提供的基于中心重排序的3D-iEPI方案与传统的3D-GRE方案相比，可以使扫描时间减少约一半，图像质量和对比度与传统的3D-GRE方案相当。同时，与基于线性重排序的3D-iEPI方案相比，本发明实施例提供的基于中心重排序的3D-iEPI方案，其图像质量能够进一步减轻血流引起的动脉污染，更好地刻画血管。
图5为本发明一实施例提供的磁敏感加权成像装置的结构示意图，该装置主要包括：
流动补偿模块51：在多次激发的平面回波成像中，针对每次激发的回波进行分层编码、相位编码和频率编码方向的流动补偿。
回波采集模块52：确定每次激发后当采用线性重排序方式时将要采集的多个回波，针对每次激发，从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波，且对本次激发的回波采集方向与对上一次激发的回波采集方向相反。
加权成像模块53：用于对采集到的回波进行磁敏感加权成像。
一可选实施例中，回波采集模块52对本次激发的回波采集方向与对上一次激发的回波采集方向相反包括：当对上一次激发的回波采集方向为：从中心回波开始向k空间正向采集到以线性重排序方式采集时的最后一个回波时，对本次激发的回波采集方向为：从中心回波开始向k空间负向采集到以线性重排序方式采集时的第一个回波；或者，当对上一次激发的回波采集方向为：从中心回波开始向k空间负向采集到以线性重排序方式采集时的第一个回波时，对本次激发的回波采集方向为：从中心回波开始向k空间正向采集到以线性重排序方式采集时的最后一个回波。
一可选实施例中，回波采集模块52从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波包括：当本次激发的选层编码梯度脉冲和相位编码梯度脉冲发射完毕时，即从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波。
一可选实施例中，流动补偿模块51针对每次激发的回波进行分层编码、相位编码和频率编码方向的流动补偿包括：针对每次激发的中心回波进行分层编码和相位编码方向的流动补偿，所述中心回波为每次激发后当采用线性重排序方式时将要采集的多个回波的中心回波；针对每次激发的第奇数个回波进行频率编码方向的流动补偿。
本发明实施例中提出的一种磁共振成像系统，可包括上述的快速磁敏感加权成像装置。
综上，本发明实施例的有益技术效果如下：
1)与已有的SWI序列相比，由于采用了高效的回波采集方案，大大缩短了回波采集时间。
2)与采用回扫的全流动补偿EPI序列相比，由于回波间隔较短，本发明实施例可以保持较高的扫描效率和较小的失真。
3)与具有线性重排序的流动补偿EPI序列相比，本发明实施例可以进一步减少SWI图像中来自动脉流的污染信号。
4)用相反的相位编码方向在连续两次激发中分别向k空间正向和负向采样获得部分回波，这提供了用标准的反向梯度方法进行失真校正的可行性，这将允许对回波序列长度进行更灵活的选择。
需要说明的是，虽然本发明方案中采用中心重排序的部分流动补偿3D-iEPI能够更好地抑制动脉信号，但中心重排序方案与线性重排序方案相比，由于使用相同的EPI因子，将采用线性重排序方式时在一次激发中采集的回波拆分成本发明实施例中采用中心重排序方式时需要在2次激发中采集，每次激发获得的回波较少，因此必然导致较低的信噪比效率。然而，通过中心重排序采集和对中心回波进行两次采样，可以减轻流动效应并且提高信噪比。此外，通过连续两次采集部分回波，可以将运动影响降到最低，并有利于后续的数据拼接。最终，从实验结果来看，由于使用的回波序列较短，采用本发明方案得到的SWI图像中的失真和模糊伪影不明显。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (11)

1.一种磁敏感加权成像方法，其特征在于，包括：
在多次激发的平面回波成像中，针对每次激发的回波进行分层编码、相位编码和频率编码方向的流动补偿；
确定每次激发后当采用线性重排序方式时将要采集的多个回波，针对每次激发，从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波，且对本次激发的回波采集方向与对上一次激发的回波采集方向相反；
对采集到的回波进行磁敏感加权成像。
2.根据权利要求1所述的方法，其特征在于，所述对本次激发的回波采集方向与对上一次激发的回波采集方向相反包括：
当对上一次激发的回波采集方向为：从中心回波开始向k空间正向采集到相当于以线性重排序方式采集时的最后一个回波时，对本次激发的回波采集方向为：从中心回波开始向k空间负向采集到相当于以线性重排序方式采集时的第一个回波；或者，
当对上一次激发的回波采集方向为：从中心回波开始向k空间负向采集到相当于以线性重排序方式采集时的第一个回波时，对本次激发的回波采集方向为：从中心回波开始向k空间正向采集到相当于以线性重排序方式采集时的最后一个回波。
3.根据权利要求1所述的方法，其特征在于，所述从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波包括：
当本次激发的选层编码梯度脉冲和相位编码梯度脉冲发射完毕时，即从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波。
4.根据权利要求1所述的方法，其特征在于，所述从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波之后、对采集到的回波进行磁敏感加权成像之前，进一步包括：
将第m次和第m+1次采集的两个中心回波叠加后取平均，将得到的平均数据作为第m次和第m+1次采集的中心回波，1≤m<M，M为本次平面回波成像的激发总次数。
5.根据权利要求1至4任一所述的方法，其特征在于，所述针对每次激发的回波进行分层编码、相位编码和频率编码方向的流动补偿包括：
针对每次激发的中心回波进行分层编码和相位编码方向的流动补偿，所述中心回波为每次激发后当采用线性重排序方式时将要采集的多个回波的中心回波；
针对每次激发的第奇数个回波进行频率编码方向的流动补偿。
6.根据权利要求5所述的方法，其特征在于，所述针对每次激发的中心回波进行分层编码和相位编码方向的流动补偿包括：
按照如下公式进行分层编码和相位编码方向的流动补偿：
M1,par＝M0,parΔtpar
M1,phase＝M0,phaseprephaseΔtp
其中，M1,par为分层编码方向在中心回波处的一阶矩，M1,phase为相位编码方向在中心回波处的一阶矩；M0,par为分层编码方向上分层编码梯度的零阶矩，M0,phaseprephase为相位编码方向上预相位梯度的零阶矩，M0,par、M0,phaseprephase在不同的激发中是不同的；Δtpar为分层编码方向上分层编码梯度中心到中心回波的时间，Δtp为相位编码方向上预相位梯度中心到中心回波的时间。
7.一种磁敏感加权成像装置，其特征在于，包括：
流动补偿模块(51)：在多次激发的平面回波成像中，针对每次激发的回波进行分层编码、相位编码和频率编码方向的流动补偿；
回波采集模块(52)：确定每次激发后当采用线性重排序方式时将要采集的多个回波，针对每次激发，从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波，且对本次激发的回波采集方向与对上一次激发的回波采集方向相反；
加权成像模块(53)：用于对采集到的回波进行磁敏感加权成像。
8.根据权利要求7所述的装置，其特征在于，所述回波采集模块(52)对本次激发的回波采集方向与对上一次激发的回波采集方向相反包括：
当对上一次激发的回波采集方向为：从中心回波开始向k空间正向采集到相当于以线性重排序方式采集时的最后一个回波时，对本次激发的回波采集方向为：从中心回波开始向k空间负向采集到相当于以线性重排序方式采集时的第一个回波；或者，
当对上一次激发的回波采集方向为：从中心回波开始向k空间负向采集到相当于以线性重排序方式采集时的第一个回波时，对本次激发的回波采集方向为：从中心回波开始向k空间正向采集到相当于以线性重排序方式采集时的最后一个回波。
9.根据权利要求7所述的装置，其特征在于，所述回波采集模块(52)从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波包括：
当本次激发的选层编码梯度脉冲和相位编码梯度脉冲发射完毕时，即从所述多个回波的中心回波开始向k空间正向或者负向采集每一回波。
10.根据权利要求7所述的装置，其特征在于，所述流动补偿模块(51)针对每次激发的回波进行分层编码、相位编码和频率编码方向的流动补偿包括：
针对每次激发的中心回波进行分层编码和相位编码方向的流动补偿，所述中心回波为每次激发后当采用线性重排序方式时将要采集的多个回波的中心回波；
针对每次激发的第奇数个回波进行频率编码方向的流动补偿。
11.一种磁共振成像系统，包括如权利要求7至10中任一项所述的磁敏感加权成像装置。
CN202010721447.3A 2020-07-24 2020-07-24 磁敏感加权成像方法、装置及磁共振成像系统 Active CN113970717B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010721447.3A CN113970717B (zh) 2020-07-24 2020-07-24 磁敏感加权成像方法、装置及磁共振成像系统
US17/383,563 US11726158B2 (en) 2020-07-24 2021-07-23 3D flow compensated interleaved EPI with a centric reordering scheme for fast high-resolution susceptibility-weighted imaging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010721447.3A CN113970717B (zh) 2020-07-24 2020-07-24 磁敏感加权成像方法、装置及磁共振成像系统

Publications (2)

Publication Number Publication Date
CN113970717A CN113970717A (zh) 2022-01-25
CN113970717B true CN113970717B (zh) 2023-11-24

Family

ID=79585808

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010721447.3A Active CN113970717B (zh) 2020-07-24 2020-07-24 磁敏感加权成像方法、装置及磁共振成像系统

Country Status (2)

Country Link
US (1) US11726158B2 (zh)
CN (1) CN113970717B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20010074354A (ko) * 2000-01-25 2001-08-04 이민화 Ｍｒａ의 영상구현시 발생하는 아티펙트제거방법
US6271665B1 (en) * 1998-09-02 2001-08-07 Stuart S. Berr Extra slice spin tagging (EST) magnetic resonance imaging for determining perfusion
WO2010030389A2 (en) * 2008-09-15 2010-03-18 University Of Pittsburgh -Of The Common Wealth System Of Higher Education An echo-specific k-space reordering approach to compatible dual-echo arteriovenography
DE102009014461A1 (de) * 2009-03-23 2010-09-30 Siemens Aktiengesellschaft Verfahren, Magnetresonanzgerät und Computerprogramm zur Erstellung von Bildern mittels paralleler Akquistionstechnik
CN102920455A (zh) * 2011-01-27 2013-02-13 美国西门子医疗解决公司 用于在磁共振成像中独立地操作脂肪和水分量的方法
CN105232046A (zh) * 2015-11-26 2016-01-13 厦门大学 一种基于重叠回波的单扫描定量磁共振t2成像方法
CN105738846A (zh) * 2014-12-12 2016-07-06 西门子(深圳)磁共振有限公司 K空间数据采集方法及其磁共振成像方法
KR20170024243A (ko) * 2015-08-25 2017-03-07 한국과학기술원 고속 자기공명영상을 위한 그룹화된 위상부호화 경사자장 배치 방법
CN108567422A (zh) * 2018-03-02 2018-09-25 上海东软医疗科技有限公司 一种血管成像方法和装置

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7046004B2 (en) * 2004-05-07 2006-05-16 University Of Basel Generic eddy-current compensation in magnetic resonance imaging
US8120358B2 (en) * 2006-04-13 2012-02-21 The Regents Of The University Of California Magnetic resonance imaging with high spatial and temporal resolution
WO2010112039A1 (en) * 2009-03-31 2010-10-07 Max-Planck-Gesellschaft zur Förderung der Wissenschaften e. V. Magnetic resonance imaging with improved imaging contrast
US10222439B2 (en) * 2014-04-24 2019-03-05 Dignity Health System and method for spiral multislab magnetic resonance imaging
US11163028B2 (en) * 2017-04-05 2021-11-02 Beth Israel Deaconess Medical Center, Inc. System and method for improved spin-echo-based magnetic resonance imaging

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6271665B1 (en) * 1998-09-02 2001-08-07 Stuart S. Berr Extra slice spin tagging (EST) magnetic resonance imaging for determining perfusion
KR20010074354A (ko) * 2000-01-25 2001-08-04 이민화 Ｍｒａ의 영상구현시 발생하는 아티펙트제거방법
WO2010030389A2 (en) * 2008-09-15 2010-03-18 University Of Pittsburgh -Of The Common Wealth System Of Higher Education An echo-specific k-space reordering approach to compatible dual-echo arteriovenography
DE102009014461A1 (de) * 2009-03-23 2010-09-30 Siemens Aktiengesellschaft Verfahren, Magnetresonanzgerät und Computerprogramm zur Erstellung von Bildern mittels paralleler Akquistionstechnik
CN102920455A (zh) * 2011-01-27 2013-02-13 美国西门子医疗解决公司 用于在磁共振成像中独立地操作脂肪和水分量的方法
CN105738846A (zh) * 2014-12-12 2016-07-06 西门子(深圳)磁共振有限公司 K空间数据采集方法及其磁共振成像方法
KR20170024243A (ko) * 2015-08-25 2017-03-07 한국과학기술원 고속 자기공명영상을 위한 그룹화된 위상부호화 경사자장 배치 방법
CN105232046A (zh) * 2015-11-26 2016-01-13 厦门大学 一种基于重叠回波的单扫描定量磁共振t2成像方法
CN108567422A (zh) * 2018-03-02 2018-09-25 上海东软医疗科技有限公司 一种血管成像方法和装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
3D Flow Compensated Interleaved EPI for aFast High-Resolution Susceptibility-WeightedI maging at 1.5T;Wei Liu 等;《ISMRM 2019》;第1-4 页 *
Reducing oblique flow effects in interleaved EPI with a centric reordering technique;Gabriele Beck 等;《Magnetic Resonance in Medicine》(第45期);623–629 *

Also Published As

Publication number Publication date
US11726158B2 (en) 2023-08-15
CN113970717A (zh) 2022-01-25
US20220026515A1 (en) 2022-01-27

Similar Documents

Publication Publication Date Title
EP2414861B1 (en) Magnetic resonance imaging with improved imaging contrast
JP6188597B2 (ja) 磁気共鳴システムの作動方法および制御装置
US9488711B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP4832786B2 (ja) 磁気共鳴断層撮影の拡散強調画像のためのマルチショット方法
US8760163B2 (en) Diffusion-weighted magnetic resonance imaging using 3D mosaic segmentation and 3D navigator phase correction
US7098662B2 (en) Method and system of MR imaging with variable readout gradient filtering
US7102352B2 (en) Method and system for image artifact reduction using nearest-neighbor phase correction for echo planar imaging
JPH09276243A (ja) 磁気共鳴診断装置
US5459400A (en) Method to enhance the sensitivity of MRI for magnetic susceptibility effects
CN105699923A (zh) 一种无创动态测量组织r2，r2*和r2`参数图像的磁共振成像方法
US5869964A (en) Magnetic resonance imaging apparatus in which gradient echo signals are acquired at a time distant from the center of a gradient echo
JPH0871054A (ja) 磁気共鳴映像装置
CN113970717B (zh) 磁敏感加权成像方法、装置及磁共振成像系统
WO2023035930A1 (zh) 一种减少图像失真的磁共振平面回波成像方法
US12111377B2 (en) Method and device for acquiring and reconstructing a sequence of diffusion-weighted magnetic resonance images covering a volume
JP2004089275A (ja) 磁気共鳴イメージング装置における位相補正方法
CN115144801A (zh) 飞跃时间成像方法、装置及磁共振成像系统
Zhang et al. Hybrid-space reconstruction with add-on distortion correction for simultaneous multi-slab diffusion MRI
CN117388781A (zh) 纵向弛豫时间加权和磁敏感加权成像图像获取方法及装置
Newbould et al. PERMEATE: High temporal resolution multi-echo/multi-slice dynamic susceptibility contrast perfusion imaging using GRAPPA EPI
JP2023552783A (ja) Ｄｉｘｏｎ式水／脂肪分離ＭＲ撮像
Dietrich General Advantages of Parallel Imaging
CN118707417A (zh) 回波平面成像扫描方法、装置及磁共振成像系统
WO2021259494A1 (en) Method and device for acquiring and reconstructing a sequence of diffusion-weighted magnetic resonance images covering a volume
JP2591405B2 (ja) Ｍｒイメージング装置のバイポーラグラジェントの調整法

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