CN113534030B - 磁共振成像的方法、装置、医疗设备及存储介质 - Google Patents

磁共振成像的方法、装置、医疗设备及存储介质 Download PDF

Info

Publication number
CN113534030B
CN113534030B CN202010285038.3A CN202010285038A CN113534030B CN 113534030 B CN113534030 B CN 113534030B CN 202010285038 A CN202010285038 A CN 202010285038A CN 113534030 B CN113534030 B CN 113534030B
Authority
CN
China
Prior art keywords
gradient field
data set
space data
group
echo signals
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
CN202010285038.3A
Other languages
English (en)
Other versions
CN113534030A (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.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare 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 Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN202010285038.3A priority Critical patent/CN113534030B/zh
Publication of CN113534030A publication Critical patent/CN113534030A/zh
Application granted granted Critical
Publication of CN113534030B publication Critical patent/CN113534030B/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/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
    • 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
    • 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

Landscapes

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

Abstract

本发明实施例公开了一种磁共振成像的方法、装置、医疗设备及存储介质，该方法包括：基于第一成像序列获取第一组回波信号，其中，第一成像序列包括第一射频脉冲和第一极性频率编码梯度场；基于第二成像序列获取第二组回波信号，其中，第二成像序列包括第二射频脉冲和第二极性频率编码梯度场，第一极性频率编码梯度场和第二极性频率编码梯度场的极性相反；将第一组回波信号和第二组回波信号分别填充至K空间，以获取第一K空间数据集和第二K空间数据集；将第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集；根据平均K空间数据集生成重建图像。本发明实施例的技术方案，有效去除了磁共振图像的伪影。

Description

磁共振成像的方法、装置、医疗设备及存储介质
技术领域
本发明实施例涉及磁共振技术领域，尤其涉及一种磁共振成像的方法、装置、医疗设备及存储介质。
背景技术
磁共振成像(MRI，Magnetic Resonance Imaging)是断层成像的一种，它利用磁共振现象从人体中获得电磁信号,从而重建出人体信息。
梯度自旋回波序列(GARSE，gradient and spin echo)的一次激发中既包括自旋回波和梯度回波，具备自旋回波的磁敏感特点的同时，缩短了扫描时间，同时序列选择具有较大的灵活性，广泛应用于磁共振成像中。
基于GRASE采集的ASL(Arterial Spin Labeling，动脉自旋标记)序列，由于GRASE序列正负梯度交替进行数据采集，因此，会出现由于系统性能不完美带来的伪影。现有技术往往通过PE(相位方向)参考线进行伪影去除，具体为：由参考线算出校准项，在图像重建阶段，利用校准参数校准后续的采集数据。然而，由于PE参考线与实际采集的数据结果不一致，或者实际的相位不一致，不能完全的去除由于系统变化带来的图像伪影问题。
发明内容
本发明实施例提供了一种磁共振成像的方法、装置、医疗设备及存储介质，以去除基于梯度自旋回波序列进行数据采集所产生的图像伪影。
第一方面，本发明实施例提供了一种磁共振成像的方法，该磁共振成像的方法包括：
基于第一成像序列激发检测对象的目标区域，以获取第一组回波信号，其中，所述第一成像序列包括第一射频脉冲和第一极性频率编码梯度场；
基于第二成像序列激发所述检测对象的目标区域，以获取第二组回波信号，其中，所述第二成像序列包括第二射频脉冲和第二极性频率编码梯度场，所述第一极性频率编码梯度场和所述第二极性频率编码梯度场的极性相反；
将所述第一组回波信号和第二组回波信号分别填充至K空间，以获取第一K空间数据集和第二K空间数据集；
将所述第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集；
根据所述平均K空间数据集生成重建图像。
第二方面，本发明实施例还提供了一种磁共振成像的方法，该磁共振成像的方法包括：
获取第一组回波信号，并将所述第一组回波信号填充至K空间，获取第一K空间数据集，其中，所述第一组回波信号为基于第一射频脉冲激发目标区域且在第一极性频率编码梯度场下编码获得；
获取第二组回波信号，并将所述第二组回波信号填充至K空间，获取第二K空间数据集，其中，所述第二组回波信号为基于第二射频脉冲激发目标区域且在第二极性频率编码梯度场下编码获得，所述第一极性频率编码梯度场和所述第二极性频率编码梯度场的极性相反；
将所述第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集；
根据所述平均K空间数据集生成重建图像。
第三方面，本发明实施例还提供了一种磁共振成像的装置，该磁共振成像的装置包括：
第一成像序列激发模块，用于基于第一成像序列激发检测对象的目标区域，以获取第一组回波信号，其中，所述第一成像序列包括第一射频脉冲和第一极性频率编码梯度场；
第二成像序列激发模块，用于基于第二成像序列激发所述检测对象的目标区域，以获取第二组回波信号，其中，所述第二成像序列包括第二射频脉冲和第二极性频率编码梯度场，所述第一极性频率编码梯度场和所述第二极性频率编码梯度场的极性相反；
K空间填充模块，用于将所述第一组回波信号和第二组回波信号分别填充至K空间，以获取第一K空间数据集和第二K空间数据集；
数据取平均模块，用于将所述第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集；
图像重建模块，用于根据所述平均K空间数据集生成重建图像。
第四方面，本发明实施例还提供了一种磁共振成像的装置，该磁共振成像的装置包括：
第一组回波获取模块，用于获取第一组回波信号，并将所述第一组回波信号填充至K空间，获取第一K空间数据集，其中，所述第一组回波信号为基于第一射频脉冲激发目标区域且在第一极性频率编码梯度场下编码获得；
第二组回波获取模块，用于获取第二组回波信号，并将所述第二组回波信号填充至K空间，获取第二K空间数据集，其中，所述第二组回波信号为基于第二射频脉冲激发目标区域且在第二极性频率编码梯度场下编码获得，所述第一极性频率编码梯度场和所述第二极性频率编码梯度场的极性相反；
数据平均处理模块，用于将所述第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集；
图像生成模块，用于根据所述平均K空间数据集生成重建图像。
第五方面，本发明实施例还提供了一种医疗设备，该医疗设备包括：
一个或多个处理器；
存储器，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如本发明任意实施例所提供的磁共振成像的方法。
第六方面，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现本发明任意实施例所提供的磁共振成像的方法。
本发明实施例的技术方案，通过第一成像序列以及与第一成像序列具有极性序列相反的频率编码梯度磁场的第二成像序列采集目标区域的两组回波信号，将该两组回波信号分别填充至K空间，得到两组K空间数据集，通过将两组K空间数据集进行取平均，利用取平均后的K空间数据进行图像重建，消除了由于采用不同梯度磁场交替采集数据带来的伪影，提高了成像质量，且算法复杂度低，易于实现。
附图说明
图1A是本发明实施例一中的一种磁共振成像的方法的流程图；
图1B是本发明实施例一中的第一成像序列的示意图；
图1C是本发明实施例一中的第二成像序列的示意图；
图1D是本发明实施例一中的动脉自旋标记序列的示意图；
图2是本发明实施例二中的一种磁共振成像的方法的流程图；
图3是本发明实施例三中的一种磁共振成像的装置的结构示意图；
图4是本发明实施例四中的一种磁共振成像的装置的结构示意图；
图5是本发明实施例五中的一种医疗设备的结构示意图；
图6是本发明实施例六中的一种磁共振成像系统的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
实施例一
图1A是本发明实施例一提供的一种磁共振成像的方法的流程图，本实施例可适用于基于梯度自旋回波序列进行磁共振成像的情况，该方法可以由磁共振成像的装置来执行，如图1A所示，该磁共振成像的方法，包括：
步骤110、基于第一成像序列激发检测对象的目标区域，以获取第一组回波信号。
其中，第一成像序列可以是梯度自旋回波成像序列。所述第一成像序列包括第一射频脉冲和第一极性频率编码梯度场。
其中，目标区域即待检测的区域，可以是脑部、心脏、乳腺、骨关节等部位。梯度场指的是磁场强度沿某一给定的方向变化的磁场，第一极性频率编码梯度场用于以梯度场的形式对第一射频脉冲进行空间编码，以生成第一组回波信号。第一射频脉冲可以高斯脉冲、矩形脉冲或者其他形式的脉冲。
具体的，磁共振成像中的成像序列(第一成像序列或后续的第二成像序列)通常用于确定激励所需的射频脉冲的参数以及各个梯度场的参数。第一成像序列中可以包括第一射频脉冲RF1、层面选择梯度场GSS、相位编码梯度场GPE和第一极性频率编码梯度场GRO1，第一射频脉冲RF1由射频线圈发射至受检对象的目标区域，从而激发目标区域的核自旋以产生磁共振信号，并通过在层面选择梯度GSS方向、相位编码梯度GPE方向和频率编码梯度GRO方向对磁共振信号进行空间编码，从而生成空间编码后的磁共振信号，该空间编码后的磁共振信号经过接收线圈，转换为MR(磁共振，Magnetic Resonance)模拟信号，再经过模数转换得到相应的数字信号，该数字信号便为上述第一组回波信号，接着，便可以按照预设方式将第一组回波信号填充至K空间，以得到第一K空间数据集。第二组回波信号的生成过程与第一组回波信号的过程相同，仅在空间编码时所采用的频率编码梯度场的极性不同。
步骤120、基于第二成像序列激发所述检测对象的目标区域，以获取第二组回波信号。
其中，第二成像序列可以是梯度自旋回波成像序列。所述第二成像序列包括第二射频脉冲和第二极性频率编码梯度场，所述第一极性频率编码梯度场和所述第二极性频率编码梯度场的极性相反。
具体的，第一极性频率编码梯度场的极性可以是正向，而第二极性频率编码梯度磁场的极性则为反向，或者，第一极性频率编码梯度场的极性为反向，而第二极性频率编码梯度磁场的极性为正向。其中，正向可以表示为梯度从左向右增大，相应的，反向则表示梯度从右向左增大。当然也可以用其他方向表示正向。
可选的，第一射频脉冲和第二射频脉冲的脉冲参数可以相同，均包括90°激发脉冲和随后施加的一个或多个180°聚集脉冲。
具体的，第一射频脉冲或第二射频脉冲的一个脉冲序列为一个90°激发脉冲和随后施加的预设数量的180°聚集脉冲。
进一步地，第一成像序列与第二成像序列的不同点可以仅为时间顺序不同以及频率编码梯度场的极性不同，其他参数均相同。
可选的，所述第一成像序列和第二成像序列的相位编码梯度场包括主相位编码梯度场和尖波形/尖峰(blip)梯度场，其中，所述尖波形梯度场的梯度场的波形为尖波形。
具体的，第一成像序列或第二成像序列中的相位编码梯度场相同，包括主相位编码梯度场和尖波形梯度场，该主相位编码梯度场作为主梯度场用于对磁共振信号进行相位编码，以进行空间识别，而尖波形梯度场的主要作用则是增加K空间中相邻的两线数据的间隔，以便于进行后续的数据处理。
可选的，还可以在所述第一成像序列前施加动脉自旋标记序列，以基于施加自旋标记序列的第一成像序列激发检测对象的目标区域，以获取动脉自旋标记的第一组回波信号，和/或，在所述第二成像序列前施加动脉自旋标记序列，以基于施加自旋标记序列的第二成像序列激发检测对象的目标区域，以获取动脉自旋标记的第二组回波信号。
其中，动脉自旋标记(Arterial Spin Labeling，ASL)是一种利用动脉血液中的水分子作为内源性对比剂实现血流成像的磁共振技术。
示例性的，请参见图1B是本发明实施例一中的第一成像序列示意图，第一成像序列为第一梯度自旋回波成像序列，采用GRASE1表示。其中，RF表示由射频线圈发射的第一射频脉冲；梯度线圈分别形成沿层面选择梯度GSS方向、相位编码梯度GPE方向和频率编码梯度GRO方向的梯度场；Echo表示在信号采集窗内采集得到的回波信号。该第一射频脉冲包括90°激发脉冲和随后施加的一个或者更多个180°聚集脉冲；第一射频脉冲的对应时序位置施加有沿层面选择梯度GSS方向的选层梯度，且在沿层面选择梯度GSS方向上分别施加有第一极性反转梯度GSS1、第二极性反转梯度GSS2、第三极性反转梯度GSS3和第四极性反转梯度GSS4。本实施例图中的上、下方向的箭头对应的梯度则为极性反转梯度。在第一极性反转梯度GSS1和第二极性反转梯度GSS2之间的时序间隙，沿相位编码梯度GPE方向上施加尖波形梯度场；同时，在第一极性反转梯度GSS1和第二极性反转梯度GSS2之间的时序间隙，沿频率编码梯度GRO方向施加连续切换的第一极性频率编码梯度场GRO1，该第一极性频率编码梯度场GRO1的起始极性为正向。90°激发脉冲的中心与回波信号中心的时序间隔即为有效回波时间(Effective TE)。进一步地，在第三极性反转梯度GSS3和第四极性反转梯度GSS4之间的时序间隙，沿频率编码梯度GRO方向施加连续切换的第一极性频率编码梯度场GRO1，该第一极性频率编码梯度场GRO1的起始极性为正向；在第三极性反转梯度GSS3和第四极性反转梯度GSS4之间的时序间隙，沿相位编码梯度GPE方向上施加尖波形梯度场。在此实施例中，第一组回波中心的回波为自旋回波，位于自旋回波两侧的为梯度回波，且梯度回波信号的强度随着距离自旋回波的距离增大而衰减。
示例性的，请参见图1C是本发明实施例一中的第二成像序列的示意图，其中，第二成像序列为第二梯度自旋回波序列，采用GRASE2表示。该第二射频脉冲包括90°激发脉冲和随后施加的一个或者更多个180°聚集脉冲；第二射频脉冲的对应时序位置施加有沿层面选择梯度GSS方向的选层梯度，且在沿层面选择梯度GSS方向上分别施加有第五极性反转梯度GSS5、第六极性反转梯度GSS6、第七极性反转梯度GSS7和第八极性反转梯度GSS8。在第五极性反转梯度GSS5和第六极性反转梯度GSS6之间的时序间隙，沿相位编码梯度GPE方向上施加尖波形梯度场；同时，在第五极性反转梯度GSS5和第六极性反转梯度GSS6之间的时序间隙，沿频率编码梯度GRO方向施加连续切换的第二极性频率编码梯度场GRO1，该第二极性频率编码梯度场GRO1的起始极性为负向。90°激发脉冲额中心与回波信号中心的时序间隔即为有效回波时间(Effective TE)。进一步地，在第七极性反转梯度GSS7和第八极性反转梯度GSS8之间的时序间隙，沿频率编码梯度GRO方向施加连续切换的第二极性频率编码梯度场GRO2，该第二极性频率编码梯度场GRO2的起始极性为负向；在第七极性反转梯度GSS7和第八极性反转梯度GSS8之间的时序间隙，沿相位编码梯度GPE方向上施加尖波形梯度场。
本实施例中，第一成像序列包括如图1B所示的第一梯度自旋回波序列和在该第一梯度自旋回波序列前施加的动脉自旋标记(ASL)序列。第二成像序列包括如图1C所示的第二梯度自旋回波序列和在该第二梯度自旋回波序列前施加的动脉自旋标记(ASL)。通过向目标区域施加第一成像序列获取动脉自旋标记的第一组回波信号，通过向目标区域施加第二成像序列获取动脉自旋标记的第二组回波信号。
在一个实施例中，如图1D是本发明实施例一提供的动脉自旋标记序列的示意图，如图1D所示，其包括多个翻转角较小的射频脉冲和沿层面选择梯度GSS方向的连续反转选层梯度。当然，第一成像序列、第二成像序列中所包含的ASL序列可包括沿相位编码梯度GPE方向和/或沿频率编码梯度GRO方向的连续反转选层梯度。
步骤130、将所述第一组回波信号和第二组回波信号分别填充至K空间，以获取第一K空间数据集和第二K空间数据集。
具体的，K空间又称为数字点阵。具体的填充过程为：根据第一组回波信号和第二组回波信号的空间编码(频率编码和相位编码)方式，进行K空间填充，可以是循序对称方式、中央优先填充方式等填充方式。
步骤140、将所述第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集。
具体的，可以将第一K空间数据集和第二K空间数据集均包括相同线数的K空间数据，将编码位置相同的第一K空间数据集和第二K空间数据集的每线数据进行相加并取平均，从而得到与第一K空间数据集或第二K空间数据集线数相同的平均K空间数据集。
步骤150、根据所述平均K空间数据集生成重建图像。
具体的，在得到平均K空间数据集后，对该平均K空间数据集进行傅里叶变换，便可以解析出磁共振信号的频率、相位和信号强度，将各个不同频率、相位以及信号强度的磁共振信号分配至相应的像素中，变得到重建图像。
本发明实施例的技术方案，通过第一成像序列以及与第一成像序列具有极性序列相反的频率编码梯度磁场的第二成像序列采集目标区域的两组回波信号，将该两组回波信号分别填充至K空间，得到两组K空间数据集，通过将两组K空间数据集进行取平均，利用取平均后的K空间数据进行图像重建，消除了由于采用不同梯度磁场交替采集数据带来的伪影，提高了成像质量，且算法复杂度低，易于实现。
实施例二
图2是本发明实施例二提供的一种磁共振成像的方法的流程图，实施例可适用于基于梯度自旋回波序列进行磁共振成像的情况，该方法可以由磁共振成像的装置来执行，如图2所示，该磁共振成像的方法，包括：
步骤210、获取第一组回波信号，并将所述第一组回波信号填充至K空间，获取第一K空间数据集。
其中，所述第一组回波信号为基于第一射频脉冲激发目标区域且在第一极性频率编码梯度场下编码获得。
可选的，第一组回波信号为动脉自旋标记的第一组回波信号，所述动脉自旋标记的第一组回波信号为基于施加动脉自旋标记序列的第一射频脉冲激发目标区域且在第一极性频率编码梯度下编码获得的。
步骤220、获取第二组回波信号，并将所述第二组回波信号填充至K空间，获取第二K空间数据集。
其中，所述第二组回波信号为基于第二射频脉冲激发目标区域且在第二极性频率编码梯度场下编码获得，所述第一极性频率编码梯度场和所述第二极性频率编码梯度场的极性相反。
可选的，第二组回波信号为动脉自旋标记的第二组回波信号，所述动脉自旋标记的第二组回波信号为基于施加动脉自旋标记序列的第二射频脉冲激发目标区域且在第二极性频率编码梯度下编码获得的。
可选的，所述第一组回波信号和第二组回波信号的相位编码梯度场包括主相位编码梯度场和尖波形梯度场，其中，所述尖波形梯度场的梯度场的波形为尖波形。
可选的，所述第一射频脉冲和第二射频脉冲的脉冲参数相同，均包括90°激发脉冲和180°聚集脉冲。
步骤230、将所述第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集。
步骤240、根据所述平均K空间数据集生成重建图像。
本发明实施例的技术方案，通过具有极性序列相反的频率编码梯度磁场采集的目标区域的两组回波信号，第一组回波信号和第二组回波信号，将该两组回波信号分别填充至K空间，得到两组K空间数据集，通过将两组K空间数据集进行取平均，利用取平均后的K空间数据进行图像重建，消除了由于采用不同梯度磁场交替采集数据带来的伪影，提高了成像质量，且算法复杂度低，易于实现。
实施例三
图3是本发明实施例三提供的一种磁共振成像的装置的结构示意图，如图3所示，该磁共振成像的装置包括：第一成像序列激发模块310、第二成像序列激发模块320、K空间填充模块330、数据取平均模块340和图像重建模块350。
其中，第一成像序列激发模块310，用于基于第一成像序列激发检测对象的目标区域，以获取第一组回波信号，其中，所述第一成像序列包括第一射频脉冲和第一极性频率编码梯度场；第二成像序列激发模块320，用于基于第二成像序列激发所述检测对象的目标区域，以获取第二组回波信号，其中，所述第二成像序列包括第二射频脉冲和第二极性频率编码梯度场，所述第一极性频率编码梯度场和所述第二极性频率编码梯度场的极性相反；K空间填充模块330，用于将所述第一组回波信号和第二组回波信号分别填充至K空间，以获取第一K空间数据集和第二K空间数据集；数据取平均模块340，用于将所述第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集；图像重建模块350，用于根据所述平均K空间数据集生成重建图像。
本发明实施例的技术方案，通过第一成像序列以及与第一成像序列具有极性序列相反的频率编码梯度磁场的第二成像序列采集目标区域的两组回波信号，将该两组回波信号分别填充至K空间，得到两组K空间数据集，通过将两组K空间数据集进行取平均，利用取平均后的K空间数据进行图像重建，消除了由于采用不同梯度磁场交替采集数据带来的伪影，提高了成像质量，且算法复杂度低，易于实现。
可选的，所述第一成像序列和第二成像序列的相位编码梯度场包括主相位编码梯度场和尖波形梯度场，其中，所述尖波形梯度场的梯度场的波形为尖波形。
可选的，第一射频脉冲和第二射频脉冲的脉冲参数相同，均包括90°激发脉冲和180°聚集脉冲。
可选的，该磁共振成像的装置，还包括：
动脉自旋标记施加模块，用于在所述第一成像序列前施加动脉自旋标记序列，以基于施加自旋标记序列的第一成像序列激发检测对象的目标区域，以获取动脉自旋标记的第一组回波信号，和/或，在所述第二成像序列前施加动脉自旋标记序列，以基于施加自旋标记序列的第二成像序列激发检测对象的目标区域，以获取动脉自旋标记的第二组回波信号。
本发明实施例所提供的磁共振成像的装置可执行本发明任意实施例所提供的磁共振成像的方法，具备执行方法相应的功能模块和有益效果。
实施例四
图4是本发明实施例四提供的一种磁共振成像的装置的结构示意图，如图4所示，该磁共振成像的装置包括：第一组回波获取模块410、第二组回波获取模块420、数据平均处理模块430和图像生成模块440。
其中，第一组回波获取模块410，用于获取第一组回波信号，并将所述第一组回波信号填充至K空间，获取第一K空间数据集，其中，所述第一组回波信号为基于第一射频脉冲激发目标区域且在第一极性频率编码梯度场下编码获得；第二组回波获取模块420，用于获取第二组回波信号，并将所述第二组回波信号填充至K空间，获取第二K空间数据集，其中，所述第二组回波信号为基于第二射频脉冲激发目标区域且在第二极性频率编码梯度场下编码获得，所述第一极性频率编码梯度场和所述第二极性频率编码梯度场的极性相反；数据平均处理模块430，用于将所述第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集；图像重建模块440，用于根据所述平均K空间数据集生成重建图像。
本发明实施例的技术方案，通过具有极性序列相反的频率编码梯度磁场采集的目标区域的两组回波信号，第一组回波信号和第二组回波信号，将该两组回波信号分别填充至K空间，得到两组K空间数据集，通过将两组K空间数据集进行取平均，利用取平均后的K空间数据进行图像重建，消除了由于采用不同梯度磁场交替采集数据带来的伪影，提高了成像质量，且算法复杂度低，易于实现。
可选的，所述第一组回波信号和第二组回波信号的相位编码梯度场包括主相位编码梯度场和尖波形梯度场，其中，所述尖波形梯度场的梯度场的波形为尖波形。
可选的，所述第一射频脉冲和第二射频脉冲的脉冲参数相同，均包括90°激发脉冲和180°聚集脉冲。
可选的，第一组回波信号为动脉自旋标记的第一组回波信号，所述动脉自旋标记的第一组回波信号为基于施加动脉自旋标记序列的第一射频脉冲激发目标区域且在第一极性频率编码梯度下编码获得的。
可选的，第二组回波信号为动脉自旋标记的第二组回波信号，所述动脉自旋标记的第二组回波信号为基于施加动脉自旋标记序列的第二射频脉冲激发目标区域且在第二极性频率编码梯度下编码获得的。
本发明实施例所提供的磁共振成像的装置可执行本发明任意实施例所提供的磁共振成像的方法，具备执行方法相应的功能模块和有益效果。
实施例五
图5是本发明实施例五提供的一种医疗设备的结构示意图，如图5所示，该设备包括处理器510、存储器520、输入装置530和输出装置540；设备处理器510的数量可以是一个或多个，图5中以一个处理器510为例；设备中的处理器510、存储器520、输入装置530和输出装置540可以通过总线或其他方式连接，图5中以通过总线连接为例。
存储器520作为一种计算机可读存储介质，可用于存储软件程序、计算机可执行程序以及模块，如本发明实施例中的业务数据的转换方法对应的程序指令/模块(例如，磁共振成像装置中的第一成像序列激发模块310、第二成像序列激发模块320、K空间填充模块330、数据取平均模块340和图像重建模块350，或者第一组回波获取模块410、第二组回波获取模块420、数据平均处理模块430和图像生成模块440)。处理器510通过运行存储在存储器520中的软件程序、指令以及模块，从而执行设备的各种功能应用以及数据处理，即实现上述的磁共振成像的方法。
存储器520可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序；存储数据区可存储根据终端的使用所创建的数据等。此外，存储器520可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中，存储器520可进一步包括相对于处理器510远程设置的存储器，这些远程存储器可以通过网络连接至设备/终端/服务器。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
输入装置530可用于接收输入的数字或字符信息，以及产生与设备的用户设置以及功能控制有关的键信号输入。输出装置540可包括显示屏等显示设备。
实施例六
图6是本发明实施例六提供的一种磁共振成像系统的结构示意图，如图6所示，该磁共振成像系统包括：扫描成像装置610和显示装置620，其中，扫描成像装置610包括：磁体单元611和射频单元612和控制模块613。
其中，磁体单元611主要包括产生主磁场B0的主磁体和产生梯度场的梯度组件，梯度场用于生成相应空间编码信号的梯度场，以对磁共振信号进行空间定位。射频单元612主要包括射频发射线圈和射频接收线圈，射频发射线圈用于向受检者或人体发射射频脉冲信号，射频接收线圈用于接收从人体采集的磁共振信号。控制模块613用于执行本发明任意实施例所提供的磁共振成像的方法，以生成重建图像。显示装置620用于接收所述重建图像并进行显示。
具体的，显示装置620可以是显示屏。主磁体可以是永磁体或超导磁体，梯度组件主要包含梯度电流放大器和梯度线圈。射频线圈可分为体线圈和局部线圈，体线圈或局部线圈的种类可以是鸟笼形线圈、螺线管形线圈、马鞍形线圈、亥姆霍兹线圈、相控阵列线圈、回路线圈等。局部线圈可以设置为相控阵列线圈，该相控阵列线圈可设置为4通道模式、8通道模式或16通道模式。
实施例七
本发明实施例七还提供一种包含计算机可执行指令的存储介质，所述计算机可执行指令在由计算机处理器执行时用于执行一种磁共振成像的方法，该方法包括：
基于第一成像序列激发检测对象的目标区域，以获取第一组回波信号，其中，所述第一成像序列包括第一射频脉冲和第一极性频率编码梯度场；
基于第二成像序列激发所述检测对象的目标区域，以获取第二组回波信号，其中，所述第二成像序列包括第二射频脉冲和第二极性频率编码梯度场，所述第一极性频率编码梯度场和所述第二极性频率编码梯度场的极性相反；
将所述第一组回波信号和第二组回波信号分别填充至K空间，以获取第一K空间数据集和第二K空间数据集；
将所述第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集；
根据所述平均K空间数据集生成重建图像。
该磁共振成像的方法，还可以包括：
获取第一组回波信号，并将所述第一组回波信号填充至K空间，获取第一K空间数据集，其中，所述第一组回波信号为基于第一射频脉冲激发目标区域且在第一极性频率编码梯度场下编码获得；
获取第二组回波信号，并将所述第二组回波信号填充至K空间，获取第二K空间数据集，其中，所述第二组回波信号为基于第二射频脉冲激发目标区域且在第二极性频率编码梯度场下编码获得，所述第一极性频率编码梯度场和所述第二极性频率编码梯度场的极性相反；
将所述第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集；
根据所述平均K空间数据集生成重建图像。
当然，本发明实施例所提供的一种包含计算机可执行指令的存储介质,其计算机可执行指令不限于如上所述的方法操作，还可以执行本发明任意实施例所提供的磁共振成像的方法所执行的步骤的相关操作。
通过以上关于实施方式的描述，所属领域的技术人员可以清楚地了解到，本发明可借助软件及必需的通用硬件来实现，当然也可以通过硬件实现，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如计算机的软盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(RandomAccess Memory,RAM)、闪存(FLASH)、硬盘或光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法、步骤。
值得注意的是，上述磁共振成像的装置的实施例中，所包括的各个单元、子单元和模块只是按照功能逻辑进行划分的，但并不局限于上述的划分，只要能够实现相应的功能即可；另外，各功能单元的具体名称也只是为了便于相互区分，并不用于限制本发明的保护范围。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。

Claims (10)

1.一种磁共振成像的方法，其包括：
基于第一成像序列激发检测对象的目标区域，以获取第一组回波信号，其中，所述第一成像序列为GRASE序列且包括第一射频脉冲和第一极性频率编码梯度场；
基于第二成像序列激发所述检测对象的目标区域，以获取动第二组回波信号，其中，所述第二成像序列为GRASE序列且包括第二射频脉冲和第二极性频率编码梯度场，所述第一极性频率编码梯度场和所述第二极性频率编码梯度场分别沿频率编码梯度方向连续切换，且第一极性频率编码梯度场和所述第二极性频率编码梯度场的起始极性相反；
将所述第一组回波信号和第二组回波信号分别填充至K空间，以获取第一K空间数据集和第二K空间数据集；
将所述第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集；
根据所述平均K空间数据集生成重建图像。
2.根据权利要求1所述的方法，其特征在于，所述第一成像序列和第二成像序列的相位编码梯度场包括主相位编码梯度场和尖波形梯度场，其中，所述尖波形梯度场的梯度场的波形为尖波形。
3.根据权利要求1所述的方法，其特征在于，所述第一射频脉冲和第二射频脉冲的脉冲参数相同，均包括90°激发脉冲和180°聚集脉冲。
4.根据权利要求1所述的方法，其特征在于，还包括：
在所述第一成像序列前施加动脉自旋标记序列，以基于施加自旋标记序列的第一成像序列激发检测对象的目标区域，以获取动脉自旋标记的第一组回波信号，和/或，
在所述第二成像序列前施加动脉自旋标记序列，以基于施加自旋标记序列的第二成像序列激发检测对象的目标区域，以获取动脉自旋标记的第二组回波信号。
5.一种磁共振成像的方法，其包括：
获取第一组回波信号，并将所述第一组回波信号填充至K空间，获取第一K空间数据集，其中，所述第一组回波信号为基于GRASE序列的第一射频脉冲激发目标区域且在第一极性频率编码梯度场下编码获得；
获取第二组回波信号，并将所述第二组回波信号填充至K空间，获取第二K空间数据集，其中，所述第二组回波信号为基于GRASE序列的第二射频脉冲激发目标区域且在第二极性频率编码梯度场下编码获得，所述第一极性频率编码梯度场和所述第二极性频率编码梯度场分别沿频率编码梯度方向连续切换，且第一极性频率编码梯度场和所述第二极性频率编码梯度场的起始极性相反；
将所述第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集；
根据所述平均K空间数据集生成重建图像。
6.根据权利要求5所述的方法，其特征在于，所述第一组回波信号和第二组回波信号的相位编码梯度场包括主相位编码梯度场和尖波形梯度场，其中，所述尖波形梯度场的梯度场的波形为尖波形。
7.一种磁共振成像装置，其特征在于，包括：
第一成像序列激发模块，用于基于第一成像序列激发检测对象的目标区域，以获取第一组回波信号，其中，所述第一成像序列为GRASE序列且包括第一射频脉冲和第一极性频率编码梯度场；
第二成像序列激发模块，用于基于第二成像序列激发所述检测对象的目标区域，以获取第二组回波信号，其中，所述第二成像序列为GRASE序列且包括第二射频脉冲和第二极性频率编码梯度场，所述第一极性频率编码梯度场和所述第二极性频率编码梯度场分别沿频率编码梯度方向连续切换，且第一极性频率编码梯度场和所述第二极性频率编码梯度场的起始极性相反；
K空间填充模块，用于将所述第一组回波信号和第二组回波信号分别填充至K空间，以获取第一K空间数据集和第二K空间数据集；
数据取平均模块，用于将所述第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集；
图像重建模块，用于根据所述平均K空间数据集生成重建图像。
8.一种磁共振成像装置，其特征在于，包括：
第一组回波获取模块，用于获取第一组回波信号，并将所述第一组回波信号填充至K空间，获取第一K空间数据集，其中，所述第一组回波信号为基于GRASE序列的第一射频脉冲激发目标区域且在第一极性频率编码梯度场下编码获得；
第二组回波获取模块，用于第二组回波信号，并将所述第二组回波信号填充至K空间，获取第二K空间数据集，其中，所述第二组回波信号为基于GRASE序列的第二射频脉冲激发目标区域且在第二极性频率编码梯度场下编码获得，所述第一极性频率编码梯度场和所述第二极性频率编码梯度场分别沿频率编码梯度方向连续切换，且第一极性频率编码梯度场和所述第二极性频率编码梯度场的起始极性相反；
数据平均处理模块，用于将所述第一K空间数据集和第二K空间数据集相同编码位置的数据相加并取平均，以获取平均K空间数据集；
图像生成模块，用于根据所述平均K空间数据集生成重建图像。
9.一种医疗设备，其特征在于，包括：
一个或多个处理器；
存储器，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如权利要求1-4或5-6中任一所述的磁共振成像的方法。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现如权利要求1-4或5-6中任一项所述的磁共振成像的方法。
CN202010285038.3A 2020-04-13 2020-04-13 磁共振成像的方法、装置、医疗设备及存储介质 Active CN113534030B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010285038.3A CN113534030B (zh) 2020-04-13 2020-04-13 磁共振成像的方法、装置、医疗设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010285038.3A CN113534030B (zh) 2020-04-13 2020-04-13 磁共振成像的方法、装置、医疗设备及存储介质

Publications (2)

Publication Number Publication Date
CN113534030A CN113534030A (zh) 2021-10-22
CN113534030B true CN113534030B (zh) 2023-03-21

Family

ID=78087908

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010285038.3A Active CN113534030B (zh) 2020-04-13 2020-04-13 磁共振成像的方法、装置、医疗设备及存储介质

Country Status (1)

Country Link
CN (1) CN113534030B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1769921A (zh) * 2004-11-02 2006-05-10 西门子公司 在磁共振断层造影中避免皱褶伪影的优化方法
CN103389481A (zh) * 2012-05-11 2013-11-13 上海联影医疗科技有限公司 磁共振频率相位双编码采样方法及其图像重建方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0633480A1 (en) * 1993-07-09 1995-01-11 Koninklijke Philips Electronics N.V. Artefact suppression in GRASE MR imaging
JPH07323021A (ja) * 1994-05-31 1995-12-12 Shimadzu Corp Ｍｒイメージング装置
US6876199B2 (en) * 2003-05-30 2005-04-05 General Electric Company Method and system for accelerated imaging using parallel MRI
CN101109791B (zh) * 2006-07-19 2010-09-29 西门子（中国）有限公司 回波平面成像序列的纠正方法
AU2015215941A1 (en) * 2011-03-23 2015-09-10 Millikelvin Technologies Llc Improved Techniques, Systems and Machine Readable Programs for Magnetic Resonance
CN104597420B (zh) * 2015-02-02 2018-01-16 清华大学 基于多次激发的磁共振扩散成像方法
CN106646302B (zh) * 2015-10-29 2019-11-19 上海联影医疗科技有限公司 一种磁共振成像方法
CN109310361B (zh) * 2016-04-07 2022-04-26 上海联影医疗科技股份有限公司 用于快速磁共振成像的系统和方法
CN106324537B (zh) * 2016-08-12 2018-11-27 江亲瑜 一种超快速分段式单次激发水脂分离方法
CN107037386B (zh) * 2016-11-01 2019-08-23 上海联影医疗科技有限公司 一种平面回波成像方法以及系统
CN110215209B (zh) * 2019-04-29 2023-03-28 上海联影医疗科技股份有限公司 一种磁共振成像方法和磁共振成像系统
CN113534031B (zh) * 2020-04-21 2023-05-05 上海联影医疗科技股份有限公司 图像域数据生成方法、计算机设备和可读存储介质

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1769921A (zh) * 2004-11-02 2006-05-10 西门子公司 在磁共振断层造影中避免皱褶伪影的优化方法
CN103389481A (zh) * 2012-05-11 2013-11-13 上海联影医疗科技有限公司 磁共振频率相位双编码采样方法及其图像重建方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
康立丽,林意群.MRI常见伪影及解决方法.1999,(02),全文. *
张中伟 ; .磁共振成像中的伪影(一):运动相关伪影.2018,(01),全文. *
李仰康,刘国瑞,Michael Amann.R成像新技术――螺旋磁共振成像.2004,(12),全文. *
胡熙宁.MRI快速成像技术原理及发展趋势.2001,(12),全文. *
马晓晖 ; 李石玲 ; 张敏 ; 赵建 ; 郭智萍 ; .关节软骨和纤维软骨的MR扫描序列技术分析.2008,(12),全文. *

Also Published As

Publication number Publication date
CN113534030A (zh) 2021-10-22

Similar Documents

Publication Publication Date Title
US6515476B1 (en) Magnetic field inhomogeneity measurement method and apparatus, phase correction method and apparatus, and magnetic resonance imaging apparatus
JP6037652B2 (ja) 拡散強調磁気共鳴データの生成方法、磁気共鳴システムおよびコンピュータ読み取り可能な記憶媒体
US6885763B2 (en) Image processing method and apparatus, recording medium, and imaging apparatus
CN102621510B (zh) 用于抑制mr成像中伪影的系统
US9746539B2 (en) MR imaging with suppresion of flow artifacts
CN103477238A (zh) 采用来自先验采集的约束的压缩感测mr图像重建
CN105572613B (zh) 磁共振化学交换饱和转移成像方法和系统
CN107076819B (zh) 具有对流伪影的抑制的Dixon MR成像
EP3201643B1 (en) Magnetic resonance imaging with enhanced bone visualization
CN111329479B (zh) 用于提供mr图像的方法、磁共振成像设备和介质
US9880250B2 (en) Method and magnetic resonance apparatus for image reconstruction with multiple virtual coils
US7418286B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
US6804383B2 (en) Image processing method and apparatus, recording medium, and imaging apparatus
US6498485B2 (en) Method and apparatus for preventing artifacts due to magnetic field outside of field of view
KR100610991B1 (ko) 자기 공명 촬상 장치
CN106175765B (zh) 磁共振成像系统及方法
CN106137198B (zh) 一种磁共振成像方法及装置
CN113534030B (zh) 磁共振成像的方法、装置、医疗设备及存储介质
US6407549B1 (en) Magnetic resonance signal acquiring method and apparatus, recording medium and magnetic resonance imaging apparatus
EP3102103B1 (en) Magnetic resonance imaging method with t2* mapping based on echo planar imaging
CN103502832A (zh) 使用具有流动敏化梯度的受激回波脉冲序列的黑血mri
CN113466768B (zh) 磁共振成像方法及磁共振成像系统
US11119173B2 (en) Dynamic imaging based on echo planar imaging sequence
CN112244812A (zh) 磁共振成像方法、磁共振成像系统和电子装置
CN112826491A (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