CN111948590B - 磁共振成像方法及装置、电子设备、存储介质 - Google Patents

磁共振成像方法及装置、电子设备、存储介质 Download PDF

Info

Publication number
CN111948590B
CN111948590B CN202010669930.1A CN202010669930A CN111948590B CN 111948590 B CN111948590 B CN 111948590B CN 202010669930 A CN202010669930 A CN 202010669930A CN 111948590 B CN111948590 B CN 111948590B
Authority
CN
China
Prior art keywords
phase
image
narrow band
magnetic resonance
resonance imaging
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
CN202010669930.1A
Other languages
English (en)
Other versions
CN111948590A (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.)
Neusoft Medical Systems Co Ltd
Shanghai Neusoft Medical Technology Co Ltd
Original Assignee
Neusoft Medical Systems Co Ltd
Shanghai Neusoft Medical Technology 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 Neusoft Medical Systems Co Ltd, Shanghai Neusoft Medical Technology Co Ltd filed Critical Neusoft Medical Systems Co Ltd
Priority to CN202010669930.1A priority Critical patent/CN111948590B/zh
Publication of CN111948590A publication Critical patent/CN111948590A/zh
Application granted granted Critical
Publication of CN111948590B publication Critical patent/CN111948590B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/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
    • 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/5618—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 both RF and gradient refocusing, e.g. GRASE
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

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

本发明公开了磁共振成像方法及装置、电子设备、存储介质。方法包括：基于多窄带放射状填充轨迹将多个回波链填充为多个窄带，其中，多个回波链为利用包含多个shot的脉冲序列激发目标扫描对象得到，每个窄带中相邻的两条相位编码线来自于激发相位不同的两个shot；分别对各个窄带进行傅里叶变换，得到对应于各个窄带的初始图像；去除各幅初始图像中伪影所处的伪影区域，得到对应于各幅初始图像的中间图像，伪影区域的位置与shot的激发相位相关；根据各幅中间图像重建针对目标扫描对象的MR图像，能够消除因直流信号产生的伪影。

Description

磁共振成像方法及装置、电子设备、存储介质
技术领域
本发明涉及医学成像技术领域，特别涉及磁共振成像方法及装置、电子设备、存储介质。
背景技术
磁共振成像(Magnetic Resonance Imaging，MRI)是现代医学影像中主要的成像方式之一，其基本工作原理是利用磁共振现象，采用激发目标扫描对象中的氢质子，利用梯度场进行位置编码，随后使用接收线圈接收带位置信息的信号，最终通过傅里叶变换重建出图像信息。
基于多窄带放射状填充轨迹(PROPELLER)的磁共振成像技术对运动(肌肉收缩、呼吸运动等)造成的伪影不敏感，目前使用的越来越频繁。但是PROPELLER采集不能解决实际成像过程中因设备噪音或外部干扰等原因引入的直流信号带来的伪影问题，该直流信号会出现在每个回波链的K空间中心位置，并在图像重建的过程中传播到整个MR图像，造成图像的放射状伪影或者致使成像较模糊。
发明内容
本发明提供磁共振成像方法及装置、电子设备、存储介质，以解决磁共振成像过程中，因设备噪音或外部干扰等原因引入的直流信号而带来的伪影问题。
具体地，本发明是通过如下技术方案实现的：
第一方面，提供一种磁共振成像方法，包括：
基于多窄带放射状填充轨迹将采集的多个回波链填充为多个窄带，每个窄带包含至少两条相位编码线，其中，所述多个回波链由利用包含多个射频激发shot的脉冲序列激发目标扫描对象得到，每个窄带中相邻的两条相位编码线来自于激发相位不同的两个shot；
分别对各个窄带进行傅里叶变换，得到对应于各个窄带的初始图像；
去除各幅初始图像中伪影所处的伪影区域，得到对应于各幅初始图像的中间图像，所述伪影区域的位置与所述shot的激发相位相关；
根据各幅中间图像重建针对所述目标扫描对象的MR图像。
可选地，去除各幅初始图像中伪影所处的伪影区域，包括：
根据所述shot的激发相位以及预定义的所述激发相位与成像视野的移动距离的对应关系确定所述初始化图像的成像视野的移动距离；
根据所述移动距离确定所述伪影区域，并裁剪所述初始图像中的伪影区域。
可选地，根据各幅中间图像重建针对所述目标扫描对象的MR图像，包括：
将对所述各幅中间图像进行傅里叶反变换得到的窄带填充于K空间；
根据所述K空间重建所述MR图像。
可选地，所述两个shot的激发相位的相位差为180°。
可选地，位于各个窄带中心的相位编码线对应的有效回波时间之差小于时间阈值。
可选地，
所述窄带的数量为偶数；
和/或，每个窄带包含的相位编码线的数量为偶数。
第二方面，提供一种磁共振成像装置，包括：
填充模块，用于基于多窄带放射状填充轨迹将采集的多个回波链填充为多个窄带，每个窄带包含至少两条相位编码线，其中，所述多个回波链由利用包含多个shot的脉冲序列激发目标扫描对象得到，每个窄带中相邻的两条相位编码线来自于激发相位不同的两个shot；
变换模块，用于分别对各个窄带进行傅里叶变换，得到对应于各个窄带的初始图像；
去除模块，用于去除各幅初始图像中伪影所处的伪影区域，得到对应于各幅初始图像的中间图像，所述伪影区域的位置与所述shot的激发相位相关；
重建模块，用于根据各幅中间图像重建针对所述目标扫描对象的MR图像。
可选地，所述去除模块具体用于：
根据所述shot的激发相位以及预定义的所述激发相位与成像视野的移动距离的对应关系确定所述初始化图像的成像视野的移动距离；
根据所述移动距离确定所述伪影区域，并裁剪所述初始图像中的伪影区域。
可选地，所述重建模块具体用于：
将对所述各幅中间图像进行傅里叶反变换得到的窄带填充于K空间；
根据所述K空间重建所述MR图像。
可选地，所述两个shot的激发相位的相位差为180°。
可选地，位于各个窄带中心的相位编码线对应的有效回波时间之差小于时间阈值。
可选地，
所述窄带的数量为偶数；
和/或，每个窄带包含的相位编码线的数量为偶数。
第三方面，提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述任一项所述的磁共振成像方法。
第四方面，提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任一项所述的磁共振成像方法的步骤。
本发明的实施例提供的技术方案可以包括以下有益效果：
本发明实施例，通过修改PROPELLER各个窄带与shot的循环采集顺序，可以消除直流信号伪影，提高MR图像的图像质量。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本发明。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本发明的实施例，并与说明书一起用于解释本发明的原理。
图1是本发明一示例性实施例示出的磁共振成像系统的结构示意图；
图2是本发明一示例性实施例示出的采用传统PROPELLER得到的一个K空间；
图3是本发明一示例性实施例示出的一种磁共振成像方法的流程图；
图4是本发明一示例性实施例示出的采用图3示出磁共振成像方法得到的K空间的示意图；
图5a是本发明一示例性实施例示出的将采用传统PROPELLER方式得到的单个窄带进行傅里叶变换后得到的初始图像的示意图；
图5b是本发明一示例性实施例示出的采用图3示出的方法得到的初始图像的示意图；
图5c是本发明一示例性实施例示出的进行一个fftshift操作的初始图像的示意图。
图6a是本发明一示例性实施例示出的基于传统的PROPELLER仿真得到的MR图像的示意图；
图6b是本发明一示例性实施例示出的采用图3示出的方法仿真得到的MR图的示意图；
图7是本发明一示例性实施例示出的一种磁共振成像装置的模块示意图；
图8是本发明一示例实施例示出的一种电子设备的结构示意图。
具体实施方式
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本发明的一些方面相一致的装置和方法的例子。
在本发明使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本发明。在本发明和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解，尽管在本发明可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本发明范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
图1是本发明一示例性实施例示出的磁共振成像(MRI)系统的结构示意图，如图1所示，MRI系包括MRI设备110、显示操作设备120、数据处理器130。MRI设备110包括磁体111、梯度线圈112、射频线圈113。
磁体111产生用于将目标扫描对象(可以是人或动物，或者人或动物的部位)中的原子核的磁偶极矩的方向调整为恒定方向的静磁场。
梯度线圈112包括用于沿着彼此相交成直角的X轴方向、Y轴方向和Z轴方向产生磁场梯度的X线圈、Y线圈和Z线圈。梯度线圈112可通过根据目标扫描对象的区域而不同地诱导共振频率来提供目标扫描对象的空间定位信息。具体的，用一个方向的梯度磁场作为层面选择梯度，确定扫描层面，然后用另外两个方向的梯度磁场来确定层面内的坐标位置，通过三个梯度磁场的不同组合，可以实现任意层面断层成像。
射频线圈113可向目标扫描对象发送RF(射频)信号，并采集从目标扫描对象发出的MR(磁共振)信号。具体的，射频线圈113产生具有与原子核的类型相应的RF信号，并且将RF信号施加到目标扫描对象上，使得目标扫描对象的原子核从低能态跃迁到高能态。当由射频线圈113产生的RF信号消失时，原子核从高能态跃迁到低能态，由此发射具有拉莫尔频率的电磁波(MR信号)。换句话讲，当停止向目标扫描对象施加RF信号时，目标扫描对象的原子核能级从高能级改变到低能级，因此发射具有拉莫尔频率的电磁波，此时射频线圈113可采集MR信号。
数据处理器130可对射频线圈113采集的MR信号进行处理，以重建成图像。显示操作设备120可显示数据处理器130重建的图像。
磁共振成像过程中，需要将采集的MR信号填充到K空间，也即将时间域的MR信号映射到频率域。K空间为MR信号的填充储存空间，对其进行傅里叶转换，即可重建出MR图像。
对于MR信号的采集，传统采用笛卡尔轨迹，它的磁共振成像方法虽然简单，但是逐行采样对运动(肌肉收缩、呼吸运动等)造成的伪影很敏感。为了解决伪影问题，基于多窄带放射状填充轨迹(PROPELLER)的磁共振成像技术发展越来越迅速，PROPELLER与快速自旋回波(turbo spin echo，TSE)序列结合应用得最为广泛，可以解决基于笛卡尔轨迹采集MR数据的运动伪影问题。
图2是本发明一示例性实施例示出的一个K空间，该K空间为将利用TSE序列对目标扫描对象进行射频激发，并基于PROPELLER采集的目标扫描对象的MR信号对K空间进行填充的结果。假设利用4次TSE序列激发目标扫描对象，该TSE序列包含一个偏转角为90°的激发脉冲以及偏转角为180°的6个重聚焦脉冲，参见图2，得到的K空间包含4个窄带(blade)分别为窄带a、窄带b、窄带c和窄带d。每次TSE序列激发称为一次shot或者一条回波链，一次shot产生6个回波(对应于6个重聚焦脉冲)，填充于一个窄带的6条相位编码线。
相关技术中，基于PROPELLER采集MR信号，通常将一次shot产生的包含多个回波信号(MR信号)的回波链填充于一个窄带，图2中相同类型的线条表示来自于相同shot的相位编码线。表1示出了图2中K空间的窄带信息与shot信息的对应关系。参见表1，窄带的编号与shot的编号一一对应，例如，编号为0的shot产生的6个回波对应于编号为a的窄带的6条相位编码线，该6条相位编码线的激发相位相同，均为0；编号为1的shot产生的6个回波对应于编号为b的窄带的6条相位编码线，该6条相位编码线的激发相位相同，均为π。表1中有效TE表示有效回波时间。
表1
Figure BDA0002581885260000071
虽然PROPELLER能够解决笛卡尔轨迹采集的MR信号的运动伪影问题，但是PROPELLER采集不能解决实际成像过程中因设备噪音或外部干扰等原因引入的直流信号带来的伪影问题，该直流信号会出现在每个回波链的K空间中心位置，并在网格化图像重建的过程中传播到整个MR图像，造成图像的放射状伪影或者致使成像较模糊。
基于上述情况，本发明提供了一种磁共振成像方法，以消除因直流信号导致的伪影。
图3是本发明一示例性实施例示出的一种磁共振成像方法的流程图，可以包括以下步骤：
步骤301、基于多窄带放射状填充轨迹将采集的多个回波链填充为多个窄带。
其中，多个回波链(MR信号)为利用包含多个shot的脉冲序列激发目标扫描对象得到，该多个shot的激发相位不完全相同。脉冲序列为适合多回波链的序列，可以但不限于采用FSE、TSE、FASE、EPI等实现。激发脉冲具有一个射频脉冲和多个重聚焦脉冲，射频脉冲的偏转角可以为90°，也可以是除了90°之外的偏转角，例如范围为0°至180°的任何大小的偏转角；同样的，重聚焦脉冲的偏转角可以为180°，也可以是除了180°之外的偏转角。
本实施例中，可以基于多窄带放射状填充轨迹将上述多个回波链填充于K空间，该K空间包含n个窄带，其中，n≥2。但是一次shot不仅仅采集一个窄带的相位编码线，而会采集多个窄带的相位编码线，并填充于K空间的对应位置，使得每个窄带中相邻的两条相位编码线来自于激发相位不同的两个shot，可以将上述K空间的填充方式称为相位循环模式。由于一次shot会采集多个窄带的相位编码线进行填充K空间，在设计脉冲序列时，最好使得K空间中窄带的数量为偶数，每个窄带包含的相位编码线的数量为偶数，使得采集到的MR信号能够完整的填充于K空间。
图4是一示例性实施例示出的K空间的示意图，图中相同类型的线条表示来自相同shot或回波链的相位编码线，从图中可以看出，每个窄带中，相邻相位编码线来自于不同的shot或回波链。
由于窄带中相邻的两条相位编码线来自于激发相位不同的两个shot，相对应的，相邻的两条相位编码线具有不同的激发相位，使得相邻两条相位编码线的激发相位具有一个相位差。基于傅里叶编码算子对相位编码线进行编码，相较于传统的PROPELLER(窄带中每条相位编码线具有相同的激发相位)，本实施例中相当于是将每条相位编码线均乘以exp(i*θ*n)，其中，i表示虚数单位，n表示第n条相位编码，θ表示激发相位。而若在频域将每条相位编码线均乘以exp(i*θ*n)，基于傅里叶变换的位移性(在频域乘以一个线性相位，相当于在时域内移动对应距离)，相当于将K空间转换为时域图像时，时域图像的FOV(成像视野)移动了与θ相对应的距离。而回波信号包含直流信号和目标扫描对象发出的质子信号，由于激发相位只会影响质子信号，而不会影响直流信号，时域图像的FOV移动了与θ相对应的距离，相当于成像于时域图像中的质子信号对应的感兴趣区域向时域图像的边缘移动了与θ相对应的距离，而直流信号产生的伪影区域的位置不变。
目前，为了便于图像重建，shot的激发相位一般采用0、π，若使用激发相位依次为0、π、0、π的4个shot对目标扫描对象进行射频激发，可以得到4个回波链，将该4个回波链填充于K空间可以形成4个窄带，每个窄带中相邻的两条相位编码线来自于激发相位不同的两个shot，也即相邻的两条相位编码线的激发相位不同。以窄带包含6条相位编码线为例，对于每个回波链，其包含的6条相位编码线的激发相位可以分别为0、π、0、π、0、π，相较于传统的PROPELLER(窄带中每条相位编码线具有相同的激发相位)，相当于在频域上将每条相位编码线均乘以exp(i*π*n)；将该回波链进行傅里叶变换后，相较于激发相同的6条相位编码线得到的时域图像，基于激发相位不同的6条相位编码线得到的时域图像的FOV移动了1/2FOV。而直流信号由于不受激发相位的影响，当时域图像移动了1/2FOV时，直流信号产生的伪影区域的位置不变。
表2示出了图4中K空间的窄带信息与shot信息的对应关系。不同于表1中每个窄带的相位编码的激发相位相同，本实施例中，参见表2，一个shot会采集多个窄带内的相位编码线，窄带的编号与shot的编号不存在一一对应关系，每个窄带中相邻的相位编码线的激发相位不同，激发相位的相位差为π。
需要说明的是，相邻的两条相位编码线不是必须来自于相邻两次shot，只需要保证相邻的两条相位编码线的激发相位不同即可。设计脉冲序列时，编号为奇数的shot一般具有相同的激发相位，编号为偶数的shot具有相同的激发相位，而奇数shot的激发相位与偶数shot的激发相位不同，从而对于相邻的两条相位编码线，可以设置为一条来自于编号为奇数的shot，一条来自于编号为偶数的shot。
表2
Figure BDA0002581885260000101
步骤302、分别对各个回波链进行傅里叶变换，得到对应于各个回波链的初始图像。
其中，傅立叶变换可以但不限于采用快速傅立叶变换FFT、二维FT、三维FT等等。该初始图像为低分辨率图像。
采集到的回波链包含目标扫描对象发出的质子信号和直流信号，参见对步骤301的解释部分，由于激发相位只会影响质子信号，而不会影响直流信号，若一次shot采集多个回波链的相位编码线进行填充K空间，相当于对质子信号乘以exp(i*θ*n)，在频域中相当于质子信号的位置移动了，而直流信号的位置不变。
参见图5a和图5b，图5a为将采用传统PROPELLER方式得到的回波链进行傅里叶变换后得到的初始图像，从图中可以看出直流信号对应的伪影区域位于图像中间位置，与质子信号对应的部分区域重叠，无法将伪影区域从图像中去除；图5b为将采用本实施例的方法得到的回波链进行傅里叶变换后得到的初始图像，该回波链中相邻两条相位编码线的激发相位的相位差为π，从图中可以看出质子信号对应的感兴趣区域发生了移动，也即感兴趣区域与伪影区域的相对位置发生了改变，两者不重合，可以将伪影区域重图中去除。
步骤303、去除各幅初始图像中伪影所在的伪影区域，得到对应于各幅初始图像的中间图像。
在一个实施例中，去除各幅初始图像中伪影所在的伪影区域时，可以采用人工手动的方式，图像重建过程中，可将各幅初始图像显示于显示操作设备，例如显示图5b，医护人员通过显示操作设备即可确定伪影区域，并将与裁剪伪影区域对应的指令发送给显示操作设备，即可实现对伪影区域的裁剪，得到中间图像，该中间图像同样为低分辨率图像。
在另一个实施例中，医护人员还可以手动把初始图像移动与θ对应的距离，例如在激发相位的相位差为π的情况下，移动1/2个FOV(或者说进行一个fftshift操作)，参见图5c，经过fftshift操作把感兴趣区域移到FOV中心位置，伪影区域移到图像边缘位置，然后对伪影区域进行裁剪。需要说明的是，本实施例中的图像移动是循环移动，比如往上移1/2个FOV，原来在图像最上面的不是移出图像，而且移到图像最下面。
在另一个实施例中，去除各幅初始图像中伪影所在的伪影区域时，也可以采用系统自动方式，首先根据相邻的两条相位编码线对应的两个shot的激发相位以及预定义的激发相位与成像视野的移动距离的对应关系确定FOV的移动距离。预定义的激发相位与成像视野的移动距离的对应关系例如可以是若两个shot的激发相位的相位差为π，对应的移动距离为1/2FOV；若两个shot的激发相位的相位差为1/2π，对应的移动距离为1/4FOV；若两个shot的激发相位的相位差为3/2π，对应的移动距离为-1/4FOV。假设当前两个shot的激发相位分别为0和π，两者相位差为π，可以确定FOV的移动距离的移动距离为1/2FOV，对初始图像进行一个fftshift操作，可以确定直流信号对应的伪影区域向图像边缘方向移动了1/2FOV，则可从初始图像中将1/2FOV对应的边缘区域裁剪掉。
由于伪影区域与感兴趣区域的相对位置与shot的激发相位相关，在另一个实施例中，可以先根据shot的激发相位以及预定义的激发相位与相对位置的对应关系确定伪影区域与感兴趣区域的相对位置，进而根据该相对位置确定FOV的移动距离。
步骤304、根据各幅中间图像重建针对目标扫描对象的MR图像。
在图像重建的过程中，可以先对各幅中间图像进行傅里叶反变换得到新的窄带，并将新的窄带重新填充于K空间，并根据该K空间重建MR图像。其中，对于MR图像的重建可以但不限于采用网格化算法重建图像。
图6a为一示例性实施例示出的基于传统的PROPELLER得到的MR图像，图6b为一示例性实施例示出采用上述实施例的方法得到的MR图，对比两幅图像可知，图6a的图像中心存在放射状伪影，而图6b中的放射状伪影被消除了，从而通过修改PROPELLER各个窄带与shot的循环采集顺序，可以消除直流信号伪影，提高MR图像的图像质量。
在另一个实施例中，由于TSE类图像的对比度由K空间中心的数据的采集时间(又称有效TE，effective TE)决定，所以每个窄带采集到的K空间中心的TE需要尽量保持一致，以保障得到T2加权图像的对比度最佳，在进行数据编码线填充时，需要确保位于各个窄带中心的相位编码线对应的有效回波时间之差小于时间阈值。
以表2为例，每个窄带包含6条相位编码线，可以将编码为0的相位编码线确定为中心编码线(K空间中心)，其决定T2图像的对比度。若每个中心编码线对应的有效TE均相同，可以得到与用户预期一致的T2加权对比度的MR图像。本实施例中，由于每个窄带的相位编码线来自于不同的shot，无法确保中心编码线对应的有效TE均相同，但是需要确保位于各个窄带中心的相位编码线对应的有效回波时间之差较小，例如可以安排各个窄带中心的相位编码线对应的有效回波时间一个是2*esp，一个是3*esp。具体的，填充K空间时，对于编号为奇数的窄带和编号为偶数的窄带的错开一个编码，例如对于第0个窄带，在第2个回波采集编码为0的相位编码线，此时对应的有效TE为3*esp；对于第1个窄带，在第1个回波采集编码为0的相位编码线，此时对应的有效TE为2*esp，这样可以使得在保证图像无运动伪影的基础上，有效回波时间之差尽量最小，使得T2加权受到的影响最小。
与前述磁共振成像方法实施例相对应，本发明还提供了磁共振成像装置的实施例。
图7是本发明一示例性实施例示出的一种磁共振成像装置的模块示意图，包括：
填充模块，用于基于多窄带放射状填充轨迹将采集的多个回波链填充为多个窄带，每个窄带包含至少两条相位编码线，其中，所述多个回波链由利用包含多个shot的脉冲序列激发目标扫描对象得到，每个窄带中相邻的两条相位编码线来自于激发相位不同的两个shot；
变换模块，用于分别对各个窄带进行傅里叶变换，得到对应于各个窄带的初始图像；
去除模块，用于去除各幅初始图像中伪影所处的伪影区域，得到对应于各幅初始图像的中间图像，所述伪影区域的位置与所述shot的激发相位相关；
重建模块，用于根据各幅中间图像重建针对所述目标扫描对象的MR图像。
可选地，所述去除模块具体用于：
根据所述shot的激发相位以及预定义的所述激发相位与成像视野的移动距离的对应关系确定所述初始化图像的成像视野的移动距离；
根据所述移动距离确定所述伪影区域，并裁剪所述初始图像中的伪影区域。
可选地，所述重建模块具体用于：
将对所述各幅中间图像进行傅里叶反变换得到的窄带填充于K空间；
根据所述K空间重建所述MR图像。
可选地，所述两个shot的激发相位的相位差为180°。
可选地，位于各个窄带中心的相位编码线对应的有效回波时间之差小于时间阈值。
可选地，
所述窄带的数量为偶数；
和/或，每个窄带包含的相位编码线的数量为偶数。
对于装置实施例而言，由于其基本对应于方法实施例，所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本发明方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
图8是本发明一示例实施例示出的一种电子设备的结构示意图，示出了适于用来实现本发明实施方式的示例性电子设备80的框图。图8显示的电子设备80仅仅是一个示例，不应对本发明实施例的功能和使用范围带来任何限制。
如图8所示，电子设备80可以以通用计算设备的形式表现，例如其可以为服务器设备。电子设备80的组件可以包括但不限于：上述至少一个处理器81、上述至少一个存储器82、连接不同系统组件(包括存储器82和处理器81)的总线83。
总线83包括数据总线、地址总线和控制总线。
存储器82可以包括易失性存储器，例如随机存取存储器(RAM)821和/或高速缓存存储器822，还可以进一步包括只读存储器(ROM)823。
存储器82还可以包括具有一组(至少一个)程序模块824的程序工具825(或实用工具)，这样的程序模块824包括但不限于：操作系统、一个或者多个应用程序、其它程序模块以及程序数据，这些示例中的每一个或某种组合中可能包括网络环境的实现。
处理器81通过运行存储在存储器82中的计算机程序，从而执行各种功能应用以及数据处理，例如上述任一实施例所提供的方法。
电子设备80也可以与一个或多个外部设备84(例如键盘、指向设备等)通信。这种通信可以通过输入/输出(I/O)接口85进行。并且，模型生成的电子设备80还可以通过网络适配器86与一个或者多个网络(例如局域网(LAN)，广域网(WAN)和/或公共网络，例如因特网)通信。如图所示，网络适配器86通过总线83与模型生成的电子设备80的其它模块通信。应当明白，尽管图中未示出，可以结合模型生成的电子设备80使用其它硬件和/或软件模块，包括但不限于：微代码、设备驱动器、冗余处理器、外部磁盘驱动阵列、RAID(磁盘阵列)系统、磁带驱动器以及数据备份存储系统等。
应当注意，尽管在上文详细描述中提及了电子设备的若干单元/模块或子单元/模块，但是这种划分仅仅是示例性的并非强制性的。实际上，根据本发明的实施方式，上文描述的两个或更多单元/模块的特征和功能可以在一个单元/模块中具体化。反之，上文描述的一个单元/模块的特征和功能可以进一步划分为由多个单元/模块来具体化。
本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任一实施例所提供的方法。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明保护的范围之内。

Claims (12)

1.一种磁共振成像方法，其特征在于，包括：
基于多窄带放射状填充轨迹将采集的多个回波链填充为多个窄带，每个窄带包含至少两条相位编码线，其中，所述多个回波链由利用包含多个射频激发shot的脉冲序列激发目标扫描对象得到，每个窄带中相邻的两条相位编码线来自于激发相位不同的两个shot；
分别对各个窄带进行傅里叶变换，得到对应于各个窄带的初始图像；
去除各幅初始图像中伪影所处的伪影区域，得到对应于各幅初始图像的中间图像，所述伪影区域的位置与所述shot的激发相位相关；
根据各幅中间图像重建针对所述目标扫描对象的MR图像；
其中，所述去除各幅初始图像中伪影所处的伪影区域，包括：
根据所述shot的激发相位以及预定义的激发相位与成像视野的移动距离的对应关系确定所述初始图像的成像视野的移动距离；
根据所述移动距离确定所述伪影区域，并裁剪所述初始图像中的伪影区域。
2.根据权利要求1所述的磁共振成像方法，其特征在于，根据各幅中间图像重建针对所述目标扫描对象的MR图像，包括：
将对所述各幅中间图像进行傅里叶反变换得到的窄带填充于K空间；
根据所述K空间重建所述MR图像。
3.根据权利要求1所述的磁共振成像方法，其特征在于，所述两个shot的激发相位的相位差为180°。
4.根据权利要求1所述的磁共振成像方法，其特征在于，位于各个窄带中心的相位编码线对应的有效回波时间之差小于时间阈值。
5.根据权利要求1所述的磁共振成像方法，其特征在于，
所述窄带的数量为偶数；
和/或，每个窄带包含的相位编码线的数量为偶数。
6.一种磁共振成像装置，其特征在于，包括：
填充模块，用于基于多窄带放射状填充轨迹将采集的多个回波链填充为多个窄带，每个窄带包含至少两条相位编码线，其中，所述多个回波链由利用包含多个shot的脉冲序列激发目标扫描对象得到，每个窄带中相邻的两条相位编码线来自于激发相位不同的两个shot；
变换模块，用于分别对各个窄带进行傅里叶变换，得到对应于各个窄带的初始图像；
去除模块，用于去除各幅初始图像中伪影所处的伪影区域，得到对应于各幅初始图像的中间图像，所述伪影区域的位置与所述shot的激发相位相关；其中，所述去除模块具体用于：根据所述shot的激发相位以及预定义的激发相位与成像视野的移动距离的对应关系确定所述初始图像的成像视野的移动距离；根据所述移动距离确定所述伪影区域，并裁剪所述初始图像中的伪影区域；
重建模块，用于根据各幅中间图像重建针对所述目标扫描对象的MR图像。
7.根据权利要求6所述的磁共振成像装置，其特征在于，所述重建模块具体用于：
将对所述各幅中间图像进行傅里叶反变换得到的窄带填充于K空间；
根据所述K空间重建所述MR图像。
8.根据权利要求6所述的磁共振成像装置，其特征在于，所述两个shot的激发相位的相位差为180°。
9.根据权利要求6所述的磁共振成像装置，其特征在于，位于各个窄带中心的相位编码线对应的有效回波时间之差小于时间阈值。
10.根据权利要求6所述的磁共振成像装置，其特征在于，
所述窄带的数量为偶数；
和/或，每个窄带包含的相位编码线的数量为偶数。
11.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5中任一项所述的磁共振成像方法。
12.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至5中任一项所述的磁共振成像方法的步骤。
CN202010669930.1A 2020-07-13 2020-07-13 磁共振成像方法及装置、电子设备、存储介质 Active CN111948590B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010669930.1A CN111948590B (zh) 2020-07-13 2020-07-13 磁共振成像方法及装置、电子设备、存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010669930.1A CN111948590B (zh) 2020-07-13 2020-07-13 磁共振成像方法及装置、电子设备、存储介质

Publications (2)

Publication Number Publication Date
CN111948590A CN111948590A (zh) 2020-11-17
CN111948590B true CN111948590B (zh) 2023-06-16

Family

ID=73340942

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010669930.1A Active CN111948590B (zh) 2020-07-13 2020-07-13 磁共振成像方法及装置、电子设备、存储介质

Country Status (1)

Country Link
CN (1) CN111948590B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103505210A (zh) * 2012-06-28 2014-01-15 西门子(深圳)磁共振有限公司 一种实现水脂分离的磁共振成像方法和装置
CN103584864A (zh) * 2012-08-15 2014-02-19 深圳迈瑞生物医疗电子股份有限公司 一种磁共振成像方法和装置

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7023207B1 (en) * 2005-02-16 2006-04-04 General Electric Company Method and system of MR imaging with reduced radial ripple artifacts
DE102009014498B4 (de) * 2009-03-23 2011-03-24 Siemens Aktiengesellschaft Verfahren, Magnetresonanzgerät und Computerprogramm zur Erstellung von Bildern mittels paralleler Akquisitionstechnik
DE102009014461B4 (de) * 2009-03-23 2011-06-22 Siemens Aktiengesellschaft, 80333 Verfahren, Magnetresonanzgerät und Computerprogramm zur Erstellung von Bildern mittels paralleler Akquistionstechnik
JP6050041B2 (ja) * 2012-07-06 2016-12-21 株式会社日立製作所 磁気共鳴イメージング装置及びｆｓｅ撮像方法
US10928473B2 (en) * 2018-09-17 2021-02-23 General Electric Company Methods and systems for reduced shading and blurring in magnetic resonance imaging
EP3650880B1 (en) * 2018-11-07 2023-03-22 Siemens Healthcare GmbH Optimized blade magnetic resonance imaging
CN109907759B (zh) * 2019-04-01 2023-05-05 上海联影医疗科技股份有限公司 磁共振成像方法和系统
CN110542872B (zh) * 2019-09-30 2022-05-13 上海东软医疗科技有限公司 磁共振成像方法及设备、相位信息的获取方法及装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103505210A (zh) * 2012-06-28 2014-01-15 西门子(深圳)磁共振有限公司 一种实现水脂分离的磁共振成像方法和装置
CN103584864A (zh) * 2012-08-15 2014-02-19 深圳迈瑞生物医疗电子股份有限公司 一种磁共振成像方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
New Algorithm for Extracting Motion Information from PROPELLER Data and Head Motion Correction in T1-Weighted MRI;Yanqiu Feng;2005 IEEE Engineering in Medicine and Biology 27th Annual Conference;第1378-1381页 *

Also Published As

Publication number Publication date
CN111948590A (zh) 2020-11-17

Similar Documents

Publication Publication Date Title
JP6998218B2 (ja) 動き検出を用いるｍｒ撮像
Gallichan et al. Retrospective correction of involuntary microscopic head movement using highly accelerated fat image navigators (3D FatNavs) at 7T
US9396562B2 (en) MRI reconstruction with incoherent sampling and redundant haar wavelets
US7840045B2 (en) Method and system for parallel reconstruction in the K-space domain for application in imaging systems
US7358730B2 (en) Diffusion tensor imaging using highly constrained image reconstruction method
CN104204838B (zh) 具有Dixon类型水/脂肪分离和关于主磁场的不均匀性的先验知识的MRI
EP2798364B1 (en) Mr imaging with suppression of flow artefacts
CN105934683B (zh) 具有对主磁场变化的估计的dixon型水/脂肪分离的mri
JP2014503290A (ja) 校正スキャン、コイル感度マップ及びナビゲータを使用する、剛体動き補償のためのパラレルｍｒｉ方法
EP2893363A1 (en) Propeller with dixon water fat separation
CN106796274B (zh) 具有伪迹抑制的propeller-mr成像
US9500734B2 (en) Magnetic resonance phase contrast angiography with rotating coding gradients
Chen et al. Self‐calibrating wave‐encoded variable‐density single‐shot fast spin echo imaging
WO2014154544A1 (en) Real-time motion correction for mri using fat navigators
CN109564268A (zh) 在多梯度回波mri中对场波动的回顾性校正
JP2001309903A (ja) Ｍｒ投影画像の自動修正
JP2020522344A (ja) 並列マルチスライスｍｒ撮像
CN110542872B (zh) 磁共振成像方法及设备、相位信息的获取方法及装置
Honal et al. Compensation of breathing motion artifacts for MRI with continuously moving table
CN114167334A (zh) 磁共振图像的重建方法、装置及电子设备
US10416265B2 (en) Method and system for generating MR images of a moving object in its environment
CN111948590B (zh) 磁共振成像方法及装置、电子设备、存储介质
Hakky et al. Application of basic physics principles to clinical neuroradiology: differentiating artifacts from true pathology on MRI
CN113466768B (zh) 磁共振成像方法及磁共振成像系统
CN111292388A (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
CB03 Change of inventor or designer information

Inventor after: Ning Jia

Inventor after: Zhang Xingxing

Inventor after: Gao Mingda

Inventor before: Ning Jia

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant