CN106353706B - 磁共振成像序列相位校正方法及磁共振设备 - Google Patents

磁共振成像序列相位校正方法及磁共振设备 Download PDF

Info

Publication number
CN106353706B
CN106353706B CN201610879030.3A CN201610879030A CN106353706B CN 106353706 B CN106353706 B CN 106353706B CN 201610879030 A CN201610879030 A CN 201610879030A CN 106353706 B CN106353706 B CN 106353706B
Authority
CN
China
Prior art keywords
signal
echo
displacement
phase
echo signal
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
CN201610879030.3A
Other languages
English (en)
Other versions
CN106353706A (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 Neusoft Medical Technology Co Ltd
Original Assignee
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 Shanghai Neusoft Medical Technology Co Ltd filed Critical Shanghai Neusoft Medical Technology Co Ltd
Priority to CN201610879030.3A priority Critical patent/CN106353706B/zh
Publication of CN106353706A publication Critical patent/CN106353706A/zh
Priority to US15/725,065 priority patent/US10905353B2/en
Application granted granted Critical
Publication of CN106353706B publication Critical patent/CN106353706B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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/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
    • 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/5617—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 RF refocusing, e.g. RARE
    • 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/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • 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]
    • 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

Landscapes

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

Abstract

本申请提供一种磁共振成像序列相位校正方法及磁共振设备，所述方法包括：通过磁共振预扫描，获取成像序列在多通道条件下每个通道的回波信号，所述回波信号包括第一回波信号和第二回波信号；对所述多个通道的第二回波信号进行多次平移，确定使所述多个通道的成像相位一致性最高的位移量，作为校正位移量；根据所述校正位移量，确定所述多个通道的回波信号的一阶相位校正值和零阶相位校正值。本申请磁共振成像序列相位校正方法，不利用信号形状直接计算一阶相位校正值，而是通过平移信号曲线确定多通道的一阶相位校正值和零阶相位校正值，能避免先计算一阶相位校正值，再计算零阶相位校正值的过程容易受各种因素的影响而导致计算错误的缺陷。

Description

磁共振成像序列相位校正方法及磁共振设备
技术领域
本申请涉及磁共振技术领域，尤其涉及磁共振成像序列相位校正方法及磁共振设备。
背景技术
在磁共振成像序列中，通过多次采集不同编码梯度下磁共振信号来重建图像，对于大部分序列来说，重建图像过程中通常都会使用FFT(Fast Fourier Transformation，快速傅氏变换)算法，该算法要求多次采集的信号有相同的初相位，以避免或降低信噪比或是伪影的产生。对于TSE(快速自旋回波)序列而言，当系统采集的信号为自旋回波信号(SE)和受激回波信号(STE)的叠加时，由于两种信号的产生机理不一致，容易导致两信号不同相，出现信号强度相互抵消，导致图像变暗等情况。而对于平面回波序列(EPI)，系统采集信号时使用交替的正负梯度，正负梯度的时间幅度等的不一致也能导致采集信号的相位不同步，出现N/2伪影。
针对这些情况，往往会在正式序列扫描之前进行一个预扫描，利用预扫描数据计算出相位，再进行相位校正，从而消除相位不一致带来的影响。一般情况，由于高阶相位差相对较小，且补偿方法相对复杂，而一阶和零阶相位差出现偏差的情况较为常见，因此校正过程一般只考虑两信号的一阶和零阶相位差，其中一阶相位差表现为信号时间上的延迟，零阶相位差表现为两信号的初始相差，实际相位校正过程就是校正零阶和一阶相位差。传统的校正方法一般是通过求信号最大值附近取一段范围的重心对应的时刻，作为信号完全聚相时刻，那么两信号重心时刻的差值即为一阶相位差在时间对应的时间差。零级相差为两个信号完全聚相时刻对应的相位值的差。但是，该校正方法对一阶相位差的准确性的依赖性非常高，当一阶相位差计算错误时，零阶相位校正值也必然是错误的，而由于一阶相差计算与整个扫描区域的信号均有关系，受磁场均匀性等影响，在扫描复杂组织时有时会导致计算错误，特别是使用多通道接收线圈进行成像时。
发明内容
本申请提供一种磁共振成像序列相位校正方法，包括：
通过磁共振预扫描，获取成像序列在多通道条件下每个通道的回波信号，所述回波信号包括第一回波信号和第二回波信号；
对所述多个通道的第二回波信号的信号曲线进行多次平移，确定使所述多个通道的成像相位一致性最高的位移量，作为校正位移量；
根据所述校正位移量，确定所述多个通道的回波信号的一阶相位校正值和零阶相位校正值。
本申请磁共振成像序列相位校正方法，不利用信号形状直接计算一阶相位校正值，而是通过平移信号曲线确定多通道的一阶相位校正值和零阶相位校正值，能够提高多通道成像相位校正时计算一阶和零阶相位的正确性与鲁棒性，从而提高磁共振成像序列的图像质量，避免了先计算一阶相位校正值，再计算零阶相位校正值的过程容易受各种因素的影响而导致计算错误的缺陷。
本申请磁共振成像序列相位校正方法的进一步改进在于，对所述多个通道的第二回波信号的信号曲线进行多次平移，包括：提供多个备选的位移量，对所述多个通道的第二回波信号的信号曲线根据所述多个备选的位移量进行相应次数的平移；其中，同一次平移过程中所述多个通道的第二回波信号的信号曲线平移的位移量相同。
本申请磁共振成像序列相位校正方法的进一步改进在于，在每次平移过程中，均包括数据处理步骤以供确定所述校正位移量，所述数据处理步骤包括：
确定在每次平移过程对应的位移量条件下所述多个通道各自的第一回波信号和第二回波信号之间的相位差；
确定在每次平移过程对应的位移量条件下所述多个通道的第一回波信号和第二回波信号的相位差的标准差，组成标准差集合。
本申请磁共振成像序列相位校正方法的进一步改进在于，确定使所述多个通道的成像相位一致性最高的位移量，作为校正位移量，包括：将所述标准差集合中标准差最小值对应的位移量以及该标准差最小值对应的位移量的前后两个位移量与及它们各自对应的标注差的值进行曲线拟合，将拟合曲线中的最小值点对应的位移量作为所述校正位移量。
本申请磁共振成像序列相位校正方法的进一步改进在于，根据所述校正位移量，确定所述多个通道的回波信号的一阶相位校正值和零阶相位校正值，包括：
将所述校正位移量作为所述多个通道的回波信号的一阶相位校正值；
将在所述校正位移量条件下所述多个通道的第一回波信号和第二回波信号的相位差的平均值作为所述多个通道的回波信号的零阶相位校正值。
本申请磁共振成像序列相位校正方法的进一步改进在于，采用二次曲线拟合的方法将所述标准差集合中标准差最小值对应的位移量以及该标准差最小值对应的位移量的前后两个位移量与及它们各自对应的标注差的值进行曲线拟合。
本申请磁共振成像序列相位校正方法的进一步改进在于，确定在每次平移过程对应的位移量条件下所述多个通道各自的第一回波信号和第二回波信号之间的相位差，包括：
在每次平移过程对应的位移量条件下，根据预设的阈值条件，在所述多个通道各自的第一回波信号和第二回波信号中选取位于同一时间段内的一段数据区域，使所述数据区域内的数据信息符合所述预设的阈值条件；
通过所述数据区域内的数据信息，确定对应的位移量条件下对应的通道的第一回波信号和第二回波信号之间的相位差。
本申请磁共振成像序列相位校正方法的进一步改进在于，在所述数据区域内进一步选取一段数据窗，在所述数据窗内选取多个数据作为所述数据信息。
本申请磁共振成像序列相位校正方法的进一步改进在于，将所述阈值条件预设为所述第一回波信号或所述第二回波信号的最大幅值的10％～20％范围内。
本申请磁共振成像序列相位校正方法的进一步改进在于，采用加权求和的方法确定在每次平移过程对应的位移量条件下所述多个通道各自的第一回波信号和第二回波信号之间的相位差。
本申请磁共振成像序列相位校正方法的进一步改进在于，所述数据处理步骤还包括：设定信噪比阈值，仅选取信噪比低于所述信噪比阈值的通道进行所述数据处理步骤。
本申请磁共振成像序列相位校正方法的进一步改进在于，所述成像序列为快速自旋回波序列，所述第一回波信号为自旋回波信号，所述第二回波信号为受激回波信号。
本申请还提供一种磁共振设备，包括：
回波信号获取模块，其构造为用于通过预扫描获取成像序列在多通道条件下每个通道的回波信号，所述回波信号包括第一回波信号和第二回波信号；
校正位移量确定模块，其构造为用于对所述回波信号获取模块获取的多个通道的第二回波信号的信号曲线进行多次平移，确定使所述多个通道的成像相位一致性最高的位移量，作为校正位移量；
相位校正值确定模块，其构造为用于根据所述校正位移量确定模块确定的校正位移量，确定所述多个通道的回波信号的一阶相位校正值和零阶相位校正值。
本申请的磁共振设备，不利用信号形状直接计算一阶相位校正值，而是通过校正位移量确定模块采用平移信号曲线的方法，进而确定多通道的一阶相位校正值和零阶相位校正值，能够提高多通道成像相位校正时计算一阶和零阶相位的正确性与鲁棒性，从而提高磁共振成像序列的图像质量，避免了先计算一阶相位校正值，再计算零阶相位校正值的过程容易受各种因素的影响而导致计算错误的缺陷。
本申请磁共振设备的进一步改进在于，所述校正位移量确定模块进一步构造为：用于根据提供的多个备选的位移量，对所述多个通道的第二回波信号的信号曲线根据所述多个备选的位移量进行相应次数的平移；其中，同一次平移过程中所述多个通道的第二回波信号的信号曲线平移的位移量相同。
本申请磁共振设备的进一步改进在于，所述校正位移量确定模块包括：数据处理模块，其构造为用于执行以下数据处理步骤：
在每次平移过程中，均进行数据处理步骤以供确定所述校正位移量，所述数据处理步骤包括：
确定在每次平移过程对应的位移量条件下所述多个通道各自的第一回波信号和第二回波信号之间的相位差；
确定在每次平移过程对应的位移量条件下所述多个通道的第一回波信号和第二回波信号的相位差的标准差，组成标准差集合。
本申请磁共振设备的进一步改进在于，所述数据处理模块进一步构造为：用于将所述标准差集合中标准差最小值对应的位移量以及该标准差最小值对应的位移量的前后两个位移量与及它们各自对应的标注差的值进行曲线拟合，将拟合曲线中的最小值点对应的位移量作为所述校正位移量。
本申请磁共振设备的进一步改进在于，所述相位校正值确定模块进一步构造为：用于将所述校正位移量作为所述多个通道的回波信号的一阶相位校正值，以及将在所述校正位移量条件下所述多个通道的第一回波信号和第二回波信号的相位差的平均值作为所述多个通道的回波信号的零阶相位校正值。
本申请磁共振设备的进一步改进在于，所述数据处理模块进一步构造为：用于在每次平移过程对应的位移量条件下，根据预设的阈值条件，在所述多个通道各自的第一回波信号和第二回波信号中选取位于同一时间段内的一段数据区域，使所述数据区域内的数据信息符合所述预设的阈值条件，以及通过所述数据区域内的数据信息，确定对应的位移量条件下对应的通道的第一回波信号和第二回波信号之间的相位差。
本申请磁共振设备的进一步改进在于，所述数据处理模块进一步构造为：用于在所述数据区域内进一步选取一段数据窗，在所述数据窗内选取多个数据作为所述数据信息。
本申请磁共振设备的进一步改进在于，所述数据处理模块还包括：信噪比阈值设定模块，其构造为用于仅选取信噪比低于所述信噪比阈值的通道进行所述数据处理步骤。
附图说明
图1是本申请示出的一种磁共振成像序列相位校正方法的流程示意图。
图2是图1所示的磁共振成像序列相位校正方法中的数据处理步骤的流程示意图。
图3是图2中步骤S201的详细流程图。
图4是图1中步骤S103的详细流程图。
图5是本申请示出的一种磁共振成像序列相位校正方法的回波信号的波形图。
图6是本申请示出的一种磁共振成像序列相位校正方法在回波信号的波形图上选取数据区域的示意图。
图7和图8是采用现有技术的相位校正方法和采用本申请示出的一种磁共振成像序列相位校正方法的结果对比图。
图9是本申请示出的一种磁共振设备的结构框图。
图10是本申请示出的一种磁共振设备的数据处理模块的结构框图。
具体实施方式
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。
在本申请使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本申请。在本申请和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解，尽管在本申请可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本申请范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
下面结合附图，对本申请磁共振成像序列相位校正方法及磁共振设备进行详细介绍。在不冲突的情况下，下述的实施例及实施方式中的特征可以相互组合。
参见图1所示，图1是本申请示出的一种磁共振成像序列相位校正方法的流程示意图。本申请的磁共振成像序列相位校正方法，包括：
步骤S101：通过磁共振预扫描，获取成像序列在多通道条件下每个通道的回波信号，所述回波信号包括第一回波信号和第二回波信号。
步骤S102：对所述多个通道的第二回波信号进行多次平移，确定使所述多个通道的成像相位一致性最高的位移量，作为校正位移量。
步骤S103：根据所述校正位移量，确定所述多个通道的回波信号的一阶相位校正值和零阶相位校正值。
本申请磁共振成像序列相位校正方法，不利用信号形状直接计算一阶相位校正值，而是通过平移信号曲线确定多通道的一阶相位校正值和零阶相位校正值，能够提高多通道成像相位校正时计算一阶和零阶相位的正确性与鲁棒性，从而提高磁共振成像序列的图像质量，避免了先计算一阶相位校正值，再计算零阶相位校正值的过程容易受各种因素的影响而导致计算错误的缺陷。
其中，在一方面，在上述步骤S102中，对所述多个通道的第二回波信号的信号曲线进行多次平移，包括：提供多个备选的位移量，对所述多个通道的第二回波信号的信号曲线根据所述多个备选的位移量进行相应次数的平移；其中，同一次平移过程中所述多个通道的第二回波信号的信号曲线平移的位移量相同。也就是说，同一个通道的第二回波信号的信号曲线每次平移的位移量是不同的，但是同一次平移过程中多个通道的第二回波信号的信号曲线平移的位移量是相同的。
进一步地，参见图2所示，图2是图1所示的磁共振成像序列相位校正方法中的数据处理步骤的流程示意图。在每次平移过程中，均包括数据处理步骤以供确定所述校正位移量，所述数据处理步骤包括：
步骤S201：确定在每次平移过程对应的位移量条件下所述多个通道各自的第一回波信号和第二回波信号之间的相位差。
步骤S202：确定在每次平移过程对应的位移量条件下所述多个通道的第一回波信号和第二回波信号的相位差的标准差，组成标准差集合。
优选地，所述数据处理步骤还可以包括步骤：设定信噪比阈值，仅选取信噪比低于所述信噪比阈值的通道进行所述数据处理步骤。这样，可以筛选掉一部分无效的通道的回波信号，以减少数据处理步骤的计算量，提高效率。
参见图3所示，图3是图2中步骤S201的详细流程图。在本申请的一个实施方式中，上述步骤S201，进一步包括：
步骤S2011：在每次平移过程对应的位移量条件下，根据预设的阈值条件，在所述多个通道各自的第一回波信号和第二回波信号中选取位于同一时间段内的一段数据区域，使所述数据区域内的数据信息符合所述预设的阈值条件。
在本申请的一个实施方式中，可以在所述数据区域内进一步选取一段数据窗，在所述数据窗内选取多个数据作为所述数据信息。优选地，将所述阈值条件预设为所述第一回波信号或所述第二回波信号的最大幅值的10％～20％范围内。
步骤S2012：通过所述数据区域内的数据信息，确定对应的位移量条件下对应的通道的第一回波信号和第二回波信号之间的相位差。
在本申请的一个实施方式中，可以采用加权求和的方法确定在每次平移过程对应的位移量条件下所述多个通道各自的第一回波信号和第二回波信号之间的相位差。但是，需要说明的是，本申请并不仅限于采用加权求和的方法，也可以采用其他的求和方法，只要能够确定在每次平移过程对应的位移量条件下所述多个通道各自的第一回波信号和第二回波信号之间的相位差即可。
在另一方面，在上述步骤S102中，确定使所述多个通道的成像相位一致性最高的位移量，作为校正位移量，包括：将所述标准差集合中标准差最小值对应的位移量以及该标准差最小值对应的位移量的前后两个位移量与及它们各自对应的标注差的值进行曲线拟合，将拟合曲线中的最小值点对应的位移量作为所述校正位移量。
在本申请的一个实施方式中，可以采用二次曲线拟合的方法将所述标准差集合中标准差最小值对应的位移量以及该标准差最小值对应的位移量的前后两个位移量与及它们各自对应的标注差的值进行曲线拟合。但是，需要说明的是，本申请并不仅限于采用二次曲线拟合的方法，也可以采用其他的曲线拟合方法，只要能实现将所述标准差集合中标准差最小值对应的位移量以及该标准差最小值对应的位移量的前后两个位移量与及它们各自对应的标注差的值进行曲线拟合即可。
参见图4所示，图4是图1中步骤S103的详细流程图，在本申请的一个实施方式中，上述步骤S103，进一步包括：
步骤S1031：将所述校正位移量作为所述多个通道的回波信号的一阶相位校正值。
步骤S1032：将在所述校正位移量条件下所述多个通道的第一回波信号和第二回波信号的相位差的平均值作为所述多个通道的回波信号的零阶相位校正值。
一般，基于k空间的数据处理方法，由于所有信号在回波中心时汇聚，信号信噪比较大，且在一阶相位差正确的时候，两信号相位差为一个常数，与成像体的分布等无关，因此，此时不同通道的相位差应该具有高度的一致性。本申请磁共振成像序列相位校正方法，正是利用了上述特点，不利用信号形状直接计算一阶相位校正值，而是通过平移信号曲线确定多通道的一阶相位校正值和零阶相位校正值，类似于在不同的一阶相位校正值下，计算零阶相位值，因为不同线圈通道的零阶相位应该保持一致，所以通过平移信号曲线的位置，找到接近正确的一阶相位校正值的值，最后通过拟合方法找出最可能的校正值。因此，本申请磁共振成像序列相位校正方法，能够提高多通道成像相位校正时计算一阶和零阶相位的正确性与鲁棒性，从而提高磁共振成像序列的图像质量，避免了先计算一阶相位校正值，再计算零阶相位校正值的过程容易受各种因素的影响而导致计算错误的缺陷。
以下通过具体实施例，对本申请磁共振成像序列相位校正方法进行更为详细的介绍。在下述的实施例中，磁共振的成像序列以快速自旋回波序列、通道的数量为7个(及通道1～7)、平移的次数为5次、5次平移的位移量分别为位移量1～5为例。则相应地，所述第一回波信号为自旋回波信号，所述第二回波信号为受激回波信号。
1)关闭相位编码梯度进行磁共振预扫描，获取快速自旋回波序列不加相位编码梯度的自旋回波信号和受激回波信号两组信号。参见图5所示，假设图5所示为通道1在位移量1的条件下的自旋回波信号和受激回波信号的波形图。其中，singal1为自旋回波信号的波形图，singal1为受激回波信号的波形图。
2)根据预设的阈值条件，在所述多个通道各自的第一回波信号和第二回波信号中选取位于同一时间段内的一段数据区域，使所述数据区域内的数据信息符合所述阈值条件。参见图6所示，假设将阈值条件设定为受激回波信号的最大幅值的10％，那么根据该阈值条件，选取出如图6中虚线范围内所表示的数据区域。将自旋回波信号在该数据区域内的数据表示为s1(i)，将受激回波信号在该数据区域内的数据表示为s2(i)，i表示数据段的位置。
3)在所述数据区域内取一定长度的数据窗，对数据窗内两信号加权求和，确定通道1在位移量1的条件下的自旋回波信号和受激回波信号之间的相位差。再对所有通道进行相同的计算，统计得到所有通道在位移量1的条件下各自的自旋回波信号和受激回波信号之间的相位差。假设数据窗的长度为5，数据中点索引为0，那么该相位差的计算公式如下：
s＝s1(-2)*s2*(-2)+s1(-1)*s2*(-1)…+s1(2)*s2*(2)
然后，对所有通道的受激回波信号的信号曲线分别在位移量2～5的条件下进行平移，计算得到所有通道在位移量2～5的条件下各自的自旋回波信号和受激回波信号之间的相位差。
最后再统计出所有通道在位移量1～5的条件下的相位值的标准差，组成标准差集合。统计结果见下表1所示。
表1
4)将所述标准差集合中标准差最小值对应的位移量以及该标准差最小值对应的位移量的前后两个位移量及它们各自对应的标注差的值进行曲线拟合，将拟合曲线中的最小值点对应的位移量作为校正位移量。如表1所示，在标准差集合中找出的最小值为2，对应的位移量为位移量2。那么，位移量2的前后两个位移量就是位移量1和位移量3，将位移量1～3与位移量1～3各自对应的标注差的值(即2.6273、2以及3.01696)进行曲线拟合，将拟合曲线中的最小值点对应的位移量作为校正位移量。
5)将所述校正位移量作为所述多个通道的回波信号的一阶相位校正值。将在所述校正位移量条件下所述多个通道的自旋回波信号和受激回波信号的相位差的平均值作为所述多个通道的回波信号的零阶相位校正值。
参见图7和图8所示，图7和图8是采用现有技术的相位校正方法和采用本申请示出的一种磁共振成像序列相位校正方法的结果对比图。图7所示为某次实测数据的结果展示图(图像数据仅为展示改善，扫描参数是随意设置，未考虑图像质量)，极坐标图中A处所指实心点表示采用本申请示出的一种磁共振成像序列相位校正方法的校正值(有多个点重合)，B处所指星点表示采用现有技术的相位校正方法得到的各通道的校正值(极轴长度为一阶相位校正值的长度，相位为零阶相位校正值)。
对应的，图8所示的两幅腰椎图是在同样扫描条件下，采用现有技术的相位校正方法和采用本申请示出的一种磁共振成像序列相位校正方法得到的腰椎对比图，图8中左边所示那幅图是采用现有技术的相位校正方法得到的腰椎图，图8中右边所示那幅图是采用本申请示出的一种磁共振成像序列相位校正方法得到的腰椎图。可见，当现有技术的相位校正方法进行相位校正时，由于相位差未达到完全一致而引起图像变暗，而采用本申请示出的一种磁共振成像序列相位校正方法可以得到正确校正值，因此得到的图像亮度相对较好。
与前述本申请的磁共振成像序列相位校正方法的实施例相对应，本申请还提供了一种磁共振设备的实施例。
参见图9所示，图9是本申请示出的一种磁共振设备的结构框图。本申请的磁共振设备，包括：
回波信号获取模块10，其构造为用于通过预扫描获取成像序列在多通道条件下每个通道的回波信号，所述回波信号包括第一回波信号和第二回波信号；
校正位移量确定模块20，其构造为用于对回波信号获取模块10获取的多个通道的第二回波信号的信号曲线进行多次平移，确定使所述多个通道的成像相位一致性最高的位移量，作为校正位移量；
相位校正值确定模块30，其构造为用于根据校正位移量确定模块20确定的校正位移量，确定所述多个通道的回波信号的一阶相位校正值和零阶相位校正值。
本申请的磁共振设备，不利用信号形状直接计算一阶相位校正值，而是通过校正位移量确定模块采用平移信号曲线的方法，进而确定多通道的一阶相位校正值和零阶相位校正值，能够提高多通道成像相位校正时计算一阶和零阶相位的正确性与鲁棒性，从而提高磁共振成像序列的图像质量，避免了先计算一阶相位校正值，再计算零阶相位校正值的过程容易受各种因素的影响而导致计算错误的缺陷。
其中，校正位移量确定模块20进一步构造为：用于根据提供的多个备选的位移量，对所述多个通道的第二回波信号的信号曲线根据所述多个备选的位移量进行相应次数的平移；其中，同一次平移过程中所述多个通道的第二回波信号的信号曲线平移的位移量相同。
相位校正值确定模块30进一步构造为：用于将所述校正位移量作为所述多个通道的回波信号的一阶相位校正值，以及将在所述校正位移量条件下所述多个通道的第一回波信号和第二回波信号的相位差的平均值作为所述多个通道的回波信号的零阶相位校正值。
进一步地，校正位移量确定模块20包括：数据处理模块210，其构造为用于执行以下数据处理步骤：
在每次平移过程中，均进行数据处理步骤以供确定所述校正位移量，所述数据处理步骤包括：
确定在每次平移过程对应的位移量条件下所述多个通道各自的第一回波信号和第二回波信号之间的相位差；
确定在每次平移过程对应的位移量条件下所述多个通道的第一回波信号和第二回波信号的相位差的标准差，组成标准差集合。
在本申请的一个实施方式中，数据处理模块210进一步构造为：用于将所述标准差集合中标准差最小值对应的位移量以及该标准差最小值对应的位移量的前后两个位移量与及它们各自对应的标注差的值进行曲线拟合，将拟合曲线中的最小值点对应的位移量作为所述校正位移量。
在本申请的一个实施方式中，数据处理模块210进一步构造为：用于在每次平移过程对应的位移量条件下，根据预设的阈值条件，在所述多个通道各自的第一回波信号和第二回波信号中选取位于同一时间段内的一段数据区域，使所述数据区域内的数据信息符合所述预设的阈值条件，以及通过所述数据区域内的数据信息，确定对应的位移量条件下对应的通道的第一回波信号和第二回波信号之间的相位差。
在本申请的一个实施方式中，数据处理模块210进一步构造为：用于在所述数据区域内进一步选取一段数据窗，在所述数据窗内选取多个数据作为所述数据信息。
参见图10所示，图10是本申请示出的一种磁共振设备的数据处理模块的结构框图。在本申请的一个实施方式中，数据处理模块210还包括：信噪比阈值设定模块220，其构造为用于仅选取信噪比低于所述信噪比阈值的通道进行所述数据处理步骤。这样，可以筛选掉一部分无效的通道的回波信号，以减少数据处理步骤的计算量，提高效率。
上述装置中各个单元的功能和作用的实现过程具体详见上述方法中对应步骤的实现过程，在此不再赘述。
对于装置实施例而言，由于其基本对应于方法实施例，所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本申请方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
以上所述仅是本申请的较佳实施例而已，并非对本申请做任何形式上的限制，虽然本申请已以较佳实施例揭露如上，然而并非用以限定本申请，任何熟悉本专业的技术人员，在不脱离本申请技术方案的范围内，当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例，但凡是未脱离本申请技术方案的内容，依据本申请的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本申请技术方案的范围内。

Claims (20)

1.一种磁共振成像序列相位校正方法，其特征在于，包括：
通过磁共振预扫描，获取成像序列在多通道条件下每个通道的回波信号，所述回波信号包括第一回波信号和第二回波信号；
提供多个备选的位移量，对多个通道的第二回波信号的信号曲线根据所述多个备选的位移量进行相应次数的平移，确定使所述多个通道的成像相位一致性最高的位移量，作为校正位移量；
根据所述校正位移量，确定所述多个通道的回波信号的一阶相位校正值和零阶相位校正值。
2.根据权利要求1所述的磁共振成像序列相位校正方法，其特征在于，同一次平移过程中所述多个通道的第二回波信号的信号曲线平移的位移量相同。
3.根据权利要求2所述的磁共振成像序列相位校正方法，其特征在于，在每次平移过程中，均包括数据处理步骤以供确定所述校正位移量，所述数据处理步骤包括：
确定在每次平移过程对应的位移量条件下所述多个通道各自的第一回波信号和第二回波信号之间的相位差；
确定在每次平移过程对应的位移量条件下所述多个通道的第一回波信号和第二回波信号的相位差的标准差，组成标准差集合。
4.根据权利要求3所述的磁共振成像序列相位校正方法，其特征在于，确定使所述多个通道的成像相位一致性最高的位移量，作为校正位移量，包括：将所述标准差集合中标准差最小值对应的位移量以及该标准差最小值对应的位移量的前后两个位移量与及它们各自对应的标注差的值进行曲线拟合，将拟合曲线中的最小值点对应的位移量作为所述校正位移量。
5.根据权利要求4所述的磁共振成像序列相位校正方法，其特征在于，根据所述校正位移量，确定所述多个通道的回波信号的一阶相位校正值和零阶相位校正值，包括：
将所述校正位移量作为所述多个通道的回波信号的一阶相位校正值；
将在所述校正位移量条件下所述多个通道的第一回波信号和第二回波信号的相位差的平均值作为所述多个通道的回波信号的零阶相位校正值。
6.根据权利要求4所述的磁共振成像序列相位校正方法，其特征在于，采用二次曲线拟合的方法将所述标准差集合中标准差最小值对应的位移量以及该标准差最小值对应的位移量的前后两个位移量与及它们各自对应的标注差的值进行曲线拟合。
7.根据权利要求3所述的磁共振成像序列相位校正方法，其特征在于，确定在每次平移过程对应的位移量条件下所述多个通道各自的第一回波信号和第二回波信号之间的相位差，包括：
在每次平移过程对应的位移量条件下，根据预设的阈值条件，在所述多个通道各自的第一回波信号和第二回波信号中选取位于同一时间段内的一段数据区域，使所述数据区域内的数据信息符合所述预设的阈值条件；
通过所述数据区域内的数据信息，确定对应的位移量条件下对应的通道的第一回波信号和第二回波信号之间的相位差。
8.根据权利要求7所述的磁共振成像序列相位校正方法，其特征在于，在所述数据区域内进一步选取一段数据窗，在所述数据窗内选取多个数据作为所述数据信息。
9.根据权利要求7或8所述的磁共振成像序列相位校正方法，其特征在于，将所述阈值条件预设为所述第一回波信号或所述第二回波信号的最大幅值的10％～20％范围内。
10.根据权利要求3所述的磁共振成像序列相位校正方法，其特征在于，采用加权求和的方法确定在每次平移过程对应的位移量条件下所述多个通道各自的第一回波信号和第二回波信号之间的相位差。
11.根据权利要求3所述的磁共振成像序列相位校正方法，其特征在于，所述数据处理步骤还包括：设定信噪比阈值，仅选取信噪比低于所述信噪比阈值的通道进行所述数据处理步骤。
12.根据权利要求1所述的磁共振成像序列相位校正方法，其特征在于，所述成像序列为快速自旋回波序列，所述第一回波信号为自旋回波信号，所述第二回波信号为受激回波信号。
13.一种磁共振设备，其特征在于，包括：
回波信号获取模块，其构造为用于通过预扫描获取成像序列在多通道条件下每个通道的回波信号，所述回波信号包括第一回波信号和第二回波信号；
校正位移量确定模块，其构造为用于根据提供的多个备选的位移量，对所述回波信号获取模块获取的多个通道的第二回波信号曲线进行相应次数的平移，确定使所述多个通道的成像相位一致性最高的位移量，作为校正位移量；
相位校正值确定模块，其构造为用于根据所述校正位移量确定模块确定的校正位移量，确定所述多个通道的回波信号的一阶相位校正值和零阶相位校正值。
14.根据权利要求13所述的磁共振设备，其特征在于，同一次平移过程中所述多个通道的第二回波信号的信号曲线平移的位移量相同。
15.根据权利要求14所述的磁共振设备，其特征在于，所述校正位移量确定模块包括：数据处理模块，其构造为用于执行以下数据处理步骤：
在每次平移过程中，均进行数据处理步骤以供确定所述校正位移量，所述数据处理步骤包括：
确定在每次平移过程对应的位移量条件下所述多个通道各自的第一回波信号和第二回波信号之间的相位差；
确定在每次平移过程对应的位移量条件下所述多个通道的第一回波信号和第二回波信号的相位差的标准差，组成标准差集合。
16.根据权利要求15所述的磁共振设备，其特征在于，所述数据处理模块进一步构造为：用于将所述标准差集合中标准差最小值对应的位移量以及该标准差最小值对应的位移量的前后两个位移量与及它们各自对应的标注差的值进行曲线拟合，将拟合曲线中的最小值点对应的位移量作为所述校正位移量。
17.根据权利要求16所述的磁共振设备，其特征在于，所述相位校正值确定模块进一步构造为：用于将所述校正位移量作为所述多个通道的回波信号的一阶相位校正值，以及将在所述校正位移量条件下所述多个通道的第一回波信号和第二回波信号的相位差的平均值作为所述多个通道的回波信号的零阶相位校正值。
18.根据权利要求15所述的磁共振设备，其特征在于，所述数据处理模块进一步构造为：用于在每次平移过程对应的位移量条件下，根据预设的阈值条件，在所述多个通道各自的第一回波信号和第二回波信号中选取位于同一时间段内的一段数据区域，使所述数据区域内的数据信息符合所述预设的阈值条件，以及通过所述数据区域内的数据信息，确定对应的位移量条件下对应的通道的第一回波信号和第二回波信号之间的相位差。
19.根据权利要求18所述的磁共振设备，其特征在于，所述数据处理模块进一步构造为：用于在所述数据区域内进一步选取一段数据窗，在所述数据窗内选取多个数据作为所述数据信息。
20.根据权利要求15所述的磁共振设备，其特征在于，所述数据处理模块还包括：信噪比阈值设定模块，其构造为用于仅选取信噪比低于所述信噪比阈值的通道进行所述数据处理步骤。
CN201610879030.3A 2016-10-08 2016-10-08 磁共振成像序列相位校正方法及磁共振设备 Active CN106353706B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201610879030.3A CN106353706B (zh) 2016-10-08 2016-10-08 磁共振成像序列相位校正方法及磁共振设备
US15/725,065 US10905353B2 (en) 2016-10-08 2017-10-04 Magnetic resonance imaging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610879030.3A CN106353706B (zh) 2016-10-08 2016-10-08 磁共振成像序列相位校正方法及磁共振设备

Publications (2)

Publication Number Publication Date
CN106353706A CN106353706A (zh) 2017-01-25
CN106353706B true CN106353706B (zh) 2019-07-23

Family

ID=57867205

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610879030.3A Active CN106353706B (zh) 2016-10-08 2016-10-08 磁共振成像序列相位校正方法及磁共振设备

Country Status (2)

Country Link
US (1) US10905353B2 (zh)
CN (1) CN106353706B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107479015B (zh) * 2017-08-21 2020-07-24 上海联影医疗科技有限公司 磁共振校准扫描序列配置和图像获取的方法及系统
CN109581254B (zh) * 2017-09-29 2021-07-30 西门子(深圳)磁共振有限公司 相位偏差获取方法及系统、相位校准方法及系统
US11249163B2 (en) 2017-11-27 2022-02-15 Shenzhen Institutes Of Advanced Technology Eddy-current correction method and apparatus, mobile terminal and readable storage medium
CN109709503B (zh) 2019-02-13 2021-03-23 上海联影医疗科技股份有限公司 一种磁共振系统控制方法、磁共振系统、计算机设备
CN110118950B (zh) * 2019-06-19 2021-01-01 华东师范大学 一种腹部定量磁化率成像中双极读出梯度的相位校正方法
CN110501664B (zh) * 2019-08-23 2021-07-16 苏州朗润医疗系统有限公司 磁共振波谱相位校正方法及应用该方法的磁共振成像系统
EP3798662A1 (en) * 2019-09-30 2021-03-31 Siemens Healthcare GmbH Trained image processing for dwi and/or tse with focus on body applications
CN112932492B (zh) * 2021-01-29 2022-08-05 上海跃磁生物科技有限公司 心磁图仪多个扫描数据同步方法、系统、介质及装置
CN115607133B (zh) * 2022-09-26 2023-11-28 中科微影(泰州)医疗科技有限公司 一种磁共振图像相位解卷绕方法
CN115951281B (zh) * 2022-12-08 2023-10-20 无锡鸣石峻致医疗科技有限公司 磁共振梯度延迟和梯度切换率同时检测方法、装置和介质

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07323021A (ja) * 1994-05-31 1995-12-12 Shimadzu Corp Ｍｒイメージング装置
JP4330247B2 (ja) * 2000-04-14 2009-09-16 株式会社日立メディコ 核磁気共鳴イメージング装置
JP3875479B2 (ja) 2000-10-20 2007-01-31 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴撮影装置
JP2002143115A (ja) 2000-10-30 2002-05-21 Ge Medical Systems Global Technology Co Llc Ｍｒイメージング方法、位相エラー測定方法およびｍｒｉ装置
JP3878429B2 (ja) 2001-04-05 2007-02-07 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｍｒｉ装置
JP5624273B2 (ja) * 2007-11-22 2014-11-12 株式会社東芝 磁気共鳴イメージング装置
EP2063282B1 (en) * 2007-11-22 2013-01-16 Kabushiki Kaisha Toshiba MR imaging of flowing matter using steady state free precession
DE102008007048B4 (de) * 2008-01-31 2010-06-17 Siemens Aktiengesellschaft Dynamische Verzeichnungskorrektur bei EPI-Messungen in der medizinischen Magnet-Resonanz-Bildgebung
JP5686660B2 (ja) * 2010-05-27 2015-03-18 株式会社東芝 磁気共鳴イメージング装置及び方法
JP5611882B2 (ja) * 2010-05-31 2014-10-22 株式会社東芝 磁気共鳴イメージング装置
CN105247382B (zh) * 2013-03-29 2018-07-24 皇家飞利浦有限公司 具有重影伪影降低的磁共振成像系统和其操作的方法
CN104569882B (zh) * 2014-12-29 2017-08-22 苏州朗润医疗系统有限公司 磁共振成像系统快速自旋回波的相位校正方法及其应用
JP6591208B2 (ja) * 2015-06-09 2019-10-16 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置

Also Published As

Publication number Publication date
US10905353B2 (en) 2021-02-02
US20180098712A1 (en) 2018-04-12
CN106353706A (zh) 2017-01-25

Similar Documents

Publication Publication Date Title
CN106353706B (zh) 磁共振成像序列相位校正方法及磁共振设备
Yu et al. Field map estimation with a region growing scheme for iterative 3‐point water‐fat decomposition
Funai et al. Regularized field map estimation in MRI
Yu et al. Combination of complex‐based and magnitude‐based multiecho water‐fat separation for accurate quantification of fat‐fraction
Jhooti et al. Phase ordering with automatic window selection (PAWS): a novel motion‐resistant technique for 3D coronary imaging
EP3869225A1 (en) System and method for estimating time of arrival (toa)
US7112964B2 (en) Eddy current measurement and correction in magnetic resonance imaging systems with a static phantom
Liu et al. Sliding interleaved kY (SLINKY) acquisition: a novel 3D MRA technique with suppressed slab boundary artifact
Griswold et al. Autocalibrated coil sensitivity estimation for parallel imaging
CN107907846B (zh) 涡流校正方法、装置、移动终端及可读存储介质
US20200300947A1 (en) System and method for reducing artifacts in echo planar magnetic resonance imaging
US10989782B2 (en) Method and magnetic resonance imaging device for providing an MR-image, computer program, and computer-readable storage medium
Benkhedah et al. Evaluation of adaptive combination of 30‐channel head receive coil array data in 23 N a MR imaging
CN104520728A (zh) 具有重影稳定性校正的多次激发磁共振成像
Zun et al. Pseudocontinuous arterial spin labeling with prospective motion correction (PCASL‐PROMO)
WO2012145754A1 (en) Correcting constant, linear and oblique phase errors in mr images acquired using propeller echo planar imaging
Robison et al. Three‐dimensional ultrashort echo‐time imaging using a FLORET trajectory
Moran et al. High‐resolution 3D radial bSSFP with IDEAL
Vegh et al. Selective channel combination of MRI signal phase
CN109544652B (zh) 基于深度生成对抗神经网络的核磁共振多加权成像方法
CN109729032A (zh) 一种频偏估计值的校正方法、装置及计算机可读存储介质
CN115629347B (zh) 一种磁共振成像系统中获得梯度轨迹的方法、装置和介质
CN108333545B (zh) 基于高通滤波的磁共振图像重建方法
Slawig et al. Multifrequency reconstruction for frequency‐modulated b SSFP
Schneider et al. B0‐informed variable density trajectory design for enhanced correction of off‐resonance effects in parallel transmission

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