CN102697501B - 校正磁共振图像中的相位信息的方法和相应的磁共振设备 - Google Patents

校正磁共振图像中的相位信息的方法和相应的磁共振设备 Download PDF

Info

Publication number
CN102697501B
CN102697501B CN201210085167.3A CN201210085167A CN102697501B CN 102697501 B CN102697501 B CN 102697501B CN 201210085167 A CN201210085167 A CN 201210085167A CN 102697501 B CN102697501 B CN 102697501B
Authority
CN
China
Prior art keywords
phase
data
magnetic resonance
signal
space
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
CN201210085167.3A
Other languages
English (en)
Other versions
CN102697501A (zh
Inventor
J.菲尔佛
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 Healthineers AG
Original Assignee
Siemens AG
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 AG filed Critical Siemens AG
Publication of CN102697501A publication Critical patent/CN102697501A/zh
Application granted granted Critical
Publication of CN102697501B publication Critical patent/CN102697501B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56563Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image 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/5615Echo 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/5616Echo 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56509Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/567Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
    • G01R33/5676Gating or triggering based on an MR signal, e.g. involving one or more navigator echoes for motion monitoring and correction

Landscapes

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

Abstract

本发明涉及一种用于借助磁共振设备(5)校正检查对象的预定体积片段的MR图像中的相位信息的方法,包括以下步骤:施加基本磁场;记录预定的体积片段的MR数据;这样分析所述MR数据,使得对预定体积片段的每个像素计算相位信息;记录导航信号,该导航信号采集基本磁场的偶然改变,所述改变通过检查对象(O)的运动或通过磁共振设备(5)本身引起;并且在考虑该导航信号的条件下对相位相信进行校正。

Description

校正磁共振图像中的相位信息的方法和相应的磁共振设备
技术领域
本发明涉及一种用于校正在建立MR图像时的相位信息的方法,所述MR图像特别是利用平面回波方法拍摄的。此外,本发明还涉及一种相应构造的磁共振设备。
背景技术
例如显示血管的氧饱和度的功能MR成像(fMRI)是用于非介入式研究例如活动的大脑的普及的方法。血液的氧饱和度取决于大脑中局部活性改变,并且由于该改变的氧饱和度,T2*弛豫时间改变或者采集的MR数据改变。由此复数值的fMRI信号包含生理信息。
对于借助fMRI的分析的基础前提条件是,所测量的MR信号在空间上以及在时间上的稳定性。这些前提条件也适用于其他应用,例如在建立用来显示扩散或灌注图的MR图像的情况下。空间上的稳定性(即,假定,在检查对象的一个确定的位置上分别一个预定的磁场强度占主导)例如可以通过呼吸或心跳或通过检查对象的运动而受到干扰。时间上的稳定性(即,假定,在一个特定的时刻或者在预定的时间段期间(在一个特定的位置上)一个预定的磁场强度占主导)例如可以通过磁共振设备的相应漂移而受到干扰。
发明内容
本发明要解决的技术问题是,比按照现有技术更好地校正对稳定性的干扰。
在本发明的范围内,提供了一种用于借助磁共振设备校正检查对象的预定体积片段的MR图像中的相位信息的方法。在此,该方法包括以下步骤:
·施加基本磁场(B0)。
·通过特别地扫描与预定的体积片段相应的K空间,记录预定的体积片段的MR数据。
·这样分析所述MR数据,使得对预定体积片段的每个像素计算一个相位信息。由此每个像素例如除了振幅信息或除了振幅值还存在相应的相位值。
·记录导航信号,其中借助该导航信号采集基本磁场的(偶然)改变,所述改变例如通过检查对象或通过磁共振设备的相应的欠缺(漂移)本身引起。
·通过该导航信号这样校正相位信息,使得基本磁场的偶然改变的作用在相位信息中被考虑。
通过按照本发明这样校正相位信息(并且不是如现有技术中通常的那样仅校正大小信息),使得计算出或校正了基本磁场的偶然改变,对稳定性的干扰比按照现有技术更好地得到校正。
按照现有技术,例如在功能性MR成像中根据基本磁场的偶然改变校正振幅信息并且用于导出生理信息。按照本发明,此时还可以关于基本磁场的改变校正图像区域中的相位信息,使得作为生理信息有利地还提供MR数据的相应校正的相位信息或图像区域中整个复数多值的信息(相位信息和振幅信息)。
为了实现图像区域中的正确的相位信息,校正可以在K空间中(通过校正K空间信号)或在图像空间中(例如在K空间信号到MR图像数据的转换之后)进行。
校正可以实时进行,从而可以实时建立相应的MR图像(其表示校正了的相位信息)。
MR数据的记录在此特别地利用平面回波方法进行。在此,在平面回波方法中从单个选择性的HF激励中读出一个或多个回波信号。在此,平面回波信号也被理解为螺旋形或非笛卡尔(非直线)地扫描K空间的方法。
借助平面回波方法从HF激励出发扫描K空间中的多个行(或轨迹(例如在正弦形振荡的梯度场的情况下)),或者也扫描整个K空间,从而与每个HF脉冲仅扫描K空间的一行的其他方法相比有利地在时间上更快地采集K空间。由此,平面回波方法很好地适合于例如用于采集大脑中短期改变的功能性MR成像。
根据本发明实施方式,逐段地扫描与待扫描的预定体积片段相应的K空间。在此,在记录每个片段的MR数据时分开地考虑利用导航信号采集的基本磁场的偶然改变,以便在组合地重建片段的MR数据的特别是相位信息之前,根据采集的基本磁场改变校正各自的片段的MR数据。
通过在将K空间的MR数据组合并用于重建之前个别地校正K空间的每个片段的MR数据,可以与总体上校正K空间的MR数据的方案相比明显降低伪影。
作为导航信号,可以采用或者是内部的导航信号或者是外部的导航信号。在此,内部的导航信号被理解为在扫描K空间的范围内或者附加地(例如通过附加地采集特定的参考行)或者隐含地(例如通过使用来自于K空间中心的MR信号)所采集的MR信号。相应地,在正常扫描的范围内所采集的K空间信号,按照本发明同时也可以是内部的导航信号,从而为了采集导航信号无需采集附加的MR信号。特别地,内部的导航信号还可以用于其他校正(例如相位编码中的图像重影(“鬼影”))。外部的导航信号被理解为例如利用呼吸带(用来采集呼吸活动的设备)或者利用心脏监视器采集的信号。也可以利用磁场探针和相应的校准来采集外部的导航信号,利用该导航信号可以采集关于基本磁场的动态频率改变的信号。
按照本发明另一个实施方式,根据第一参考相位值、根据第二参考相位值、第一相位值和第二相位值,对相位信号进行校正。在此,为了采集该相位值,利用非相位编码的导航信号对横向磁化进行采集。换言之,利用导航信号测量自由感应衰减。在参考测量的范围内利用导航信号在HF激励脉冲之后的第一时间段测量第一参考相位值;而通过在参考测量时在K空间中心在HF激励脉冲之后的第二时间段测量从导航信号得到的回波信号的横向磁化,来根据导航信号确定第二参考相位值。如下地确定第一相位值:在正常测量的范围内在相应的HF激励之后的第一时间段之后,根据自由感应衰减,利用非相位编码的任意测量信号对横向磁化进行采集。如下地确定第二相位值:在正常测量的范围内在K空间的中心中在HF激励之后的所述第二时间段对从测量信号中得到的回波信号的横向磁化进行采集。
按照该实施方式,因为既采集在两个时刻的参考相位值又采集在相同的两个时刻的相位值,所以为了校正基本磁场的改变既考虑了偏置也考虑了该改变的每个时间单元的线性改变。
因为参考相位值的采集在上面描述的实施方式中意味着脉冲序列的改变,所以还存在这样的实施方式,在该实施方式中,通过在参考测量时在K空间的中心中在HF激励脉冲之后采集所产生的回波信号的横向磁化一个时间段来仅确定参考相位值并且通过在实际的测量中在K空间中心中在HF激励脉冲之后的该时间段之后采集所产生的回波信号的横向磁化来仅确定相位值。在该简化的实施方式中,仅根据参考相位值和相位值来校正相位信息。
按照该简化的实施方式,在校正基本磁场的改变的情况下仅考虑每个时间单位的线性改变(并且不考虑偏置)。
按照本发明,在建立MR图像时相位信息的校正出现在前面部分。然而,所采集的基本磁场改变也可以用于校正对预定体积片段的每个像素从MR数据中计算的振幅信息或大小信息。为此,在K空间中采集MR数据时或在图像空间中分析(重建)MR数据时可以考虑采集的基本磁场改变,以便相应地也校正振幅信息。
为了建立失真图(英语:“distortion map”)可以对K空间第一次和然后第二次进行扫描。其中在两种情况下采用平面回波方法。第二扫描中与第一扫描相比,平面回波方法的相位编码梯度具有小的附加梯度矩(即,该矩具有数值上比相位编码梯度的实际振幅小的振幅),通过该梯度矩进行在第一扫描和第二扫描的结果之间在空间中在与相位编码梯度相应的方向上的恒定移动。通过对每个像素计算在第一扫描的结果和第二扫描的结果之间的相位差,建立失真图。在此基本磁场的偶然改变的作用按照本发明在计算相位差的情况下被校正。由此,失真图对于预定体积片段的每个像素说明,要如何校正相应像素的相位值,以便校正采集的基本磁场改变。
为了建立场图(英语:“field map”)可以分别在第一回波时间之后和在第二回波时间之后(也就是双倍地)扫描K空间。在此,根据在第一回波时间采集的MR数据建立第一相位图,并且根据在第二回波时间采集的MR数据建立第二相位图。在此相位图对每个像素说明,对于相应的像素已经计算了哪些相位值。按照本发明,在此校正相位值,使得两个相位图的相位信息有利地考虑所采集的基本磁场改变。场图根据第一相位图和第二相位图的差而被建立并且对于每个像素说明,对于该像素基本磁场是多强。
通过采用本发明来建立场图,例如可以从场图中消除生理伪影。
在本发明的范围内,还提供一种用于采集MR图像数据的磁共振设备。在此磁共振设备包括基本场磁体、梯度场系统、HF天线和控制装置,用于控制梯度场系统和HF天线、接收由HF天线记录的测量信号、分析所述测量信号并由此采集MR图像。为此,磁共振设备施加基本磁场并且采集预定体积片段的MR数据。磁共振设备分析这些MR数据并且对于预定体积片段的每个像素计算一个相位信息。磁共振设备采集取决于基本磁场的偶然改变的导航信号并且借助该导航信号借助所述相位信息。
按照本发明的磁共振设备的优点在此基本上相应于按照本发明的方法的优点,所述方法在前面已经详细解释了,从而在此不再重复。
此外,本发明还描述了一种计算机程序产品,特别是计算机程序或软件,其可以被加载到磁共振设备的可编程控制器或计算单元的存储器中。当计算机程序产品在磁共振设备的控制器或控制装置中运行时,利用所述计算机程序产品可以执行按照本发明的方法的前述所有或不同实施方式。在此,计算机程序产品可能需要程序装置,例如数据库和辅助函数,用于实现该方法的相应实施方式。换言之,利用针对计算机程序产品的权利要求特别地要求保护一种用来执行按照本发明的方法的上述实施方式或执行这些实施方式的计算机程序或软件。在此,软件可以是尚需编译(翻译)和连接或仅需翻译的源代码(例如C++),或者是为了执行仅还需加载到相应的计算单元中的可执行软件代码。
最后本发明还公开了一种电子可读数据载体,例如DVD、磁带或USB棒,在其上存储了电子可读控制信息,特别是软件(参见以上)。当这些控制信息(软件)由数据载体读取并且存储到磁共振设备的控制器或计算单元中时,可以执行所述方法的所以按照本发明的实施方式。
本发明特别适合于功能性MR成像、产生扩散加权的MR图像或产生灌注加权的MR图像。当然,本发明不限于这些优选应用领域,因为本发明有利地可以在结果取决于正确的相位信息的任何应用中被采用。而且在功能性MR成像领域中的新的方法,例如,用来在处理通过大脑的特定的信息的情况下显示在特定的大脑区域之间的关系或者采用功能性连接的MR成像(英语:“connectivitymapping,连通映射”),得益于本发明,这一点通过按照本发明不仅校正振幅,而且校正相位来实现。特别是就稳定性和图像质量来说,MR成像方法如下地得益:按照本发明也关于相位信息来采集和校正基本磁场的偶然的动态的时间上和空间上的改变。
附图说明
以下借助附图结合按照本发明的实施方式详细解释本发明。其中,
图1示意性示出了按照本发明的MR设备。
图2对于在平面回波方法中的参考测量和正常测量示出了相位随时间的变化。
图3示出了在导航频率以及相位值和谱之间的关系。
图4示出了按照本发明的方法的流程图。
具体实施方式
图1示出了磁共振设备5(磁共振成像或核自旋断层造影设备)的示意图。在此,基本场磁体1产生在时间上恒定的强磁场,用于在对象O的检查区域中,例如人体的待检查的部位(例如头部)核自旋的极化或对齐,为了数据采集将置于卧榻23上的人体驶入磁共振设备5中。在典型的球形测量空间M中定义用于核自旋共振测量所需的基本磁场的高均匀性。为了支持均匀性要求并且特别是为了消除时间上不可变的影响,在合适的位置安装由铁磁材料构成的所谓的匀场片。通过匀场线圈2尽可能消除时间上可变的影响。
在基本场磁体1中采用由三个子线圈组成的圆柱形的梯度线圈系统3。由放大器给每个子线圈提供用于在笛卡尔坐标系的各个方向上产生线性梯度场的电流。在此,梯度场系统3的第一子线圈产生x方向上的梯度Gx,第二子线圈产生y方向上的梯度Gy,并且第三子线圈产生z方向上的梯度Gz。放大器包括数模转换器,该数模转换器由用于时间正确地产生梯度脉冲的序列控制装置18控制。
一个(或多个)高频天线4位于梯度线圈系统3内,所述高频天线4将高频功率放大器给出的高频脉冲转换为用于待检查的对象O或者对象O的待检查区域的核的激励以及核自旋的对齐的交变磁场。每个高频天线4由组件线圈的以环形优选线性或矩阵形布置的形式的一个或多个HF发送线圈和一个或多个HF接收线圈组成。各自的高频天线4的HF接收线圈也将从进动的核自旋出发的交变场、即通常由一个或多个高频脉冲和一个或多个梯度脉冲组成的脉冲序列引起的核自旋回波信号,转换为电压(测量信号),该电压经过放大器7被传输到高频系统22的高频接收通道8。高频系统22还包括发送信道9,在该发送信道9中产生用于激励核磁共振的高频脉冲。在此,将各个高频脉冲根据由设备计算机20预先给出的脉冲序列在序列控制装置18中数字地表示为复数的序列。该数列作为实部和虚部分别经过输入端12被传输到高频系统22中的数模转换器并且从该数模转换器被传输到发送信道9。在发送信道9中将脉冲序列加调制到高频载波信号上,其基频相应于测量空间中核自旋的谐振频率。
通过发送-接收转接器6进行发送运行和接收运行的切换。高频天线4的HF发送线圈将用于激励核自旋的高频脉冲入射到测量空间M,并且通过HF接收线圈扫描所得到的回波信号。相应获得的核共振信号在高频系统22的接收信道8′(第一解调器)中被相位敏感地解调到中频,并且在模拟-数字-转换器(ADC)中被数字化。该信号还被解调到频率0。在第二解调器8中在数字域中数字化之后进行到频率0的解调和到实部和虚部的分离。通过图像计算机17可以从这样获得的测量数据中重建MR图像。通过设备计算机20进行测量数据、图像数据和控制程序的管理。序列控制装置18根据利用控制程序的预定值来控制各个期望的脉冲序列的产生和k空间的相应扫描。在此,序列控制装置18特别地控制梯度的时间正确的接通、具有定义的相位振幅的高频脉冲的发送以及核共振信号的接收。由合成器19提供用于高频系统22和序列控制装置18的时间基准。通过包括键盘15、鼠标16和显示屏14的终端13选择用于产生MR图像的相应控制程序,所述MR图像例如存储在DVD 21中,以及显示所产生的MR图像。
图2对于参考测量42和用于获得MR图像的正常测量41示出了在K空间中采集的信号的相位Φ(t)关于时间t的关系,其中利用平面回波方法来采集MR数据。此外,以相同的时间比例示出了读出梯度Gro和相位编码梯度Gpe。按照假定的模型,相位值由于基本磁场的偶然改变而线性地关于时间改变,从而在不同的时刻需要两个测量点,以便确定(在时刻0的)相位值和频率偏移。如下地测量相位值Φ:利用非相位编码的导航信号(遵从自由感应衰减的信号)在HF激励之后在第一时刻测量横向磁化并且在K空间中心中在后面的第二时刻测量由相应的导航信号产生的回波信号。在此,在两种情况下扫描一个在整个待采集的层上综合的信号。如果在K空间中心采集回波信号,导航信号还可以在读出方向上编码。实践中,导航信号和回波信号对于参考测量根据信号最大值被采集并且被用于所有后面的测量。
根据参考测量42和实际测量41的在时刻T1的相位值ΦR,1和Φn,1和在时刻T2的相位值ΦR,2和Φn,2,可以借助以下等式(1)确定在时刻的0相位差ΔΦ。
ΔΦ=(T2·(Φn,1R,1)-T1·(Φn,2R,2))/(T2-T1)        (1)
关于相同的相位值和时刻(时间段)可以按照以下等式(2)确定在参考测量的频率(相位的斜率)和实际测量的频率(斜率)之间的频率改变(斜率改变)Δω。
Δω=((Φn,2R,2)-(Φn,1R,1))/(T2-T1)                (2)
根据相位差ΔΦ和频率改变Δω,现在可以通过乘以按照以下等式(3)的取决于时间t的系数F来校正在K空间中在实际测量期间所采集的信号。
F(t)=e-i(ΔΦ+Δω·t)                                (3)
按照简化的校正,在K空间中心中对于参考测量42和实际测量41在时刻T(相应于第二时间段T2)仅采集回波信号。从这样所采集的相位值ΦR,Φn和时间段T中,可以通过以下等式(4)确定在参考测量42的频率(斜率)和实际测量41的频率(斜率)之间的频率改变(斜率改变)Δω。
Δω=(ΦRn)/T                                      (4)
利用频率改变可以通过乘以按照以下等式(5)的取决于时间t的系数F来校正K空间信号。
F(t)=e-i·Δω·t                                     (5)
简化的校正的前提条件是,相位差ΔΦ在测量(K空间的扫描)期间基本上保持恒定。
从利用按照等式(3)或按照等式(5)的项所校正的K空间信号中,此时例如借助傅里叶变换对于预定体积片段的每个体素确定振幅值和相位值,所述预定体积片段与在其中采集K空间信号的K空间相应。换言之,通过两个前面描述的校正相应地对于预定体积片段的每个体素还对相位信息、也就是相位值进行校正。
图3示出本发明的结果。图3a)示出了导航频率31(即,利用导航信号所采集的频率改变)。导航频率31基本上相应于上面描述的频率改变Δω。
在图3b)中示出了在没有按照本发明的校正32以及利用按照本发明的校正33的情况下预定体积片段的体素的相位值关于时间的平均值。可以看出,没有校正的曲线32最大程度相应于图3a)中的所采集的频率改变31。
图3c)示出没有按照本发明的校正34以及利用按照本发明的校正35在K空间中所采集的测量信号的谱。通过按照本发明的校正消除了由于呼吸36、心跳37和由于磁共振设备的冷却头38所引起的伪影。
图4示出了按照本发明的方法的流程图。
在第一步骤S1中,施加基本磁场B0。
然后在接下来的步骤S2中,借助平面回波方法扫描与待扫描的预定体积片段相应的K空间,其中还采集导航信号。
根据通过导航信号采集的结果,在步骤S3中,直接校正K空间测量值或校正在图像空间中从中导出的MR数据。
从校正了的K空间测量值或校正了的MR数据出发,在步骤S4中,计算对于预定体积片段的每个像素的相位值(以及振幅值)。

Claims (14)

1.一种用于借助磁共振设备(5)校正检查对象(O)的预定体积片段的MR图像中的相位信息的方法,包括以下步骤:
施加基本磁场,
记录预定的体积片段的MR数据,
这样分析所述MR数据,使得对所述预定体积片段的每个像素计算相位信息,
记录导航信号,该导航信号采集了基本磁场的偶然改变,所述改变通过检查对象的运动或通过所述磁共振设备(5)本身引起,并且
在考虑该导航信号的条件下对相位信息进行校正。
2.根据权利要求1所述的方法,其特征在于,
利用平面回波方法来进行对所述MR数据的记录,
在该平面回波方法中,从单个选择性的HF激励中读出一个或多个回波信号。
3.根据权利要求2所述的方法,其特征在于,
逐段地扫描与所述预定体积片段相应的K空间,
在采集每个片段的MR数据时分开地考虑采集的基本磁场改变,以便在分析组合的片段的MR数据以重建相位信息之前,相应校正各自的片段的MR数据。
4.根据权利要求1的方法,其特征在于,借助于呼吸带或者借助于心脏监视器来采集所述导航信号。
5.根据权利要求1所述的方法,其特征在于,
如下地确定第一参考相位值(ΦR,1):在参考测量(42)中根据在HF激励脉冲之后的第一时间段(T1)的自由感应衰减,利用非相位编码的导航信号对横向磁化进行采集,
如下地确定第二参考相位值(ΦR,2):在该参考测量(42)中在K空间中心在HF激励脉冲之后的第二时间段(T2),对从导航信号得到的回波信号的横向磁化进行采集,
如下地确定第一相位值(Φn,1):在测量(41)中在HF激励之后的第一时间段(T1)之后根据自由感应衰减,利用非相位编码的测量信号对横向磁化进行采集,
如下地确定第二相位值(Φn,2):在该测量(41)中在K空间的中心中在HF激励之后第二时间段(T2)之后,对从测量信号中得到的回波信号的横向磁化进行采集,并且
根据第一参考相位值(ΦR,1)、根据第二参考相位值(ΦR,2)、第一相位值(Φn,1)和第二相位值(Φn,2),在扫描K空间时对MR数据进行校正。
6.根据权利要求5所述的方法,其特征在于,
通过将MR数据乘以取决于时间的系数F,校正MR数据,
按照以下等式确定所述系数F,
F(t)=e-i(ΔΦ+Δω·t)
按照以下等式计算ΔΦ和Δω,
ΔΦ=(T2·(Φn,1R,1)-T1·(Φn,2R,2))/(T2-T1)
Δω=((Φn,2R,2)-(Φn,1R,1))/(T2-T1)
其中,ΦR,1是第一参考相位值,ΦR,2是第二参考相位值,Φn,1是第一相位值,Φn,2是第二相位值,T1是第一时间段并且T2是第二时间段。
7.根据权利要求1至6中任一项所述的方法,其特征在于,
如下地确定参考相位值(ΦR,2):在参考测量中在K空间中心中在HF激励之后的一个时间段(T2)对所得到的回波信号的横向磁化进行采集,
如下地确定相位值(Φn,2):在测量中在K空间中心中在HF激励之后所述时间段(T2)之后的该时间段(T2)之后对所产生的回波信号的横向磁化进行采集,并且
根据所述参考相位值和相位值,在扫描K空间时对MR数据进行校正。
8.根据权利要求7所述的方法,其特征在于,
如下地校正MR数据:将MR数据乘以取决于时间的系数F,
按照以下等式来确定所述系数F,
F(t)=e-i·Δω·t
按照以下等式来计算Δω,
Δω=(ΦRn)/T
其中,ΦR是参考相位值,Φn是相位值并且T是时间段。
9.根据权利要求1至6中任一项所述的方法,其特征在于,
这样分析MR数据,使得对所述预定体积片段的每个像素计算振幅信息,并且
在K空间中采集MR数据时和/或在分析MR数据时考虑所采集的基本磁场改变,以便也相应校正振幅信息。
10.根据权利要求1至6中任一项所述的方法,其特征在于,
将所述方法用于功能性MR成像、用于功能性连接的MR成像、用于产生扩散加权的MR图像或用于产生灌注加权的MR图像。
11.根据权利要求1至6中任一项所述的方法,其特征在于,
分别利用平面回波方法对K空间进行第一次并且然后第二次扫描,
在第二扫描中与第一扫描相比,平面回波方法的相位编码梯度具有小的附加梯度矩,由此实现了在第一扫描和第二扫描的结果之间在K空间中沿着与相位编码梯度相应的方向的恒定移动,
通过第一扫描的结果和第二扫描的结果的逐像素的相位差建立失真图,所述失真图对于预定体积片段的每个像素说明,如何校正对于相应的像素计算的相位值,以便考虑基本磁场的改变。
12.根据权利要求1至6中任一项所述的方法,其特征在于,
如下地双倍扫描K空间:在第一和第二回波时间之后采集MR数据,
从在第一回波时间所采集的MR数据中建立第一相位图,并且从在第二回波时间所采集的MR数据中建立第二相位图,其中相位图对每个像素给出相应的相位值,并且从第一相位图和第二相位图的差中建立场图,利用所述场图能够对于每个像素确定基本磁场的强度。
13.一种用于校正在检查对象(O)的预定的体积片段的MR图像中的相位信息的磁共振设备,
其中,所述磁共振设备(5)包括基本场磁体(1)、梯度场系统(3)、至少一个HF天线(4)和控制装置(10),所述控制装置被用于控制所述梯度场系统(3)和所述至少一个HF天线(4),用于接收由该至少一个HF天线(4)所记录的测量信号并且用于分析所述测量信号,
所述磁共振设备(5)被构造为,所述磁共振设备(5)施加基本磁场,所述磁共振设备(5)记录预定的体积片段的MR数据,所述磁共振设备(5)这样分析所述MR数据,使得对所述预定体积片段的每个像素计算相位信息,所述磁共振设备(5)记录导航信号,该导航信号采集基本磁场的偶然改变,所述改变通过检查对象(O)的运动或通过磁共振设备(5)本身引起,并且所述磁共振设备(5)在考虑该导航信号的条件下对相位信息进行校正。
14.根据权利要求13所述的磁共振设备,其特征在于,
所述磁共振设备(5)被构造为用于执行按照权利要求1至12中任一项所述的方法。
CN201210085167.3A 2011-03-28 2012-03-28 校正磁共振图像中的相位信息的方法和相应的磁共振设备 Active CN102697501B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102011006230.0 2011-03-28
DE102011006230A DE102011006230B4 (de) 2011-03-28 2011-03-28 Pixelweise Korrektur einer Phaseninformation in MR-Bildern unter Verwendung eines Navigatorsignals

Publications (2)

Publication Number Publication Date
CN102697501A CN102697501A (zh) 2012-10-03
CN102697501B true CN102697501B (zh) 2015-10-28

Family

ID=46844685

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210085167.3A Active CN102697501B (zh) 2011-03-28 2012-03-28 校正磁共振图像中的相位信息的方法和相应的磁共振设备

Country Status (5)

Country Link
US (1) US9329254B2 (zh)
JP (1) JP2012205897A (zh)
KR (1) KR101657331B1 (zh)
CN (1) CN102697501B (zh)
DE (1) DE102011006230B4 (zh)

Families Citing this family (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2807494A1 (en) * 2012-03-22 2014-12-03 TDK Corporation Movable coil scanner systems and methods
DE102012222413B4 (de) 2012-12-06 2023-10-26 Siemens Healthcare Gmbh Verfahren zum Erzeugen eines HF-Anregungspulses zur Anregung eines beliebig geformten Volumens, Verfahren zum gezielten Anregen von Spins innerhalb eines Gefäßes und Verfahren zur Erstellung von MR-Angiographiebildern sowie entsprechende Magnetresonanzanlage
DE102012223789B4 (de) * 2012-12-19 2014-07-17 Siemens Aktiengesellschaft Ermittlung einer Phasendifferenz-Karte
KR101593480B1 (ko) 2013-01-04 2016-02-15 연세대학교 원주산학협력단 영상강도 최소화를 통한 자기공명 확산텐서 영상정합 및 왜곡보정 방법 및 그 시스템
DE102013201670B3 (de) 2013-02-01 2014-07-03 Siemens Aktiengesellschaft Verfahren zum Erfassen von MR-Daten und zur Bestimmung eines B1-Magnetfelds sowie entsprechend ausgestaltete Magnetresonanzanlage
DE102013209295B4 (de) * 2013-05-21 2016-11-17 Siemens Healthcare Gmbh Korrektur von MR-Bilddatensätzen unter Nutzung einer Ähnlichkeit zeitlich aufeinanderfolgender Datensätze
DE102014209364B3 (de) * 2014-05-16 2015-10-08 Siemens Aktiengesellschaft Bestimmung von komplexen Sensitivitätsfaktoren von HF-Empfangsspulen
US10162037B2 (en) 2015-09-29 2018-12-25 Siemens Healthcare Gmbh Navigator-based data correction for simultaneous multislice MR imaging
DE102015222835B4 (de) * 2015-11-19 2019-06-13 Siemens Healthcare Gmbh Magnetresonanzbildgebungsverfahren mit simultaner Bildaufnahme von mehreren Teilvolumen mit einer synchronen Bildaufnahme von Navigatoren
DE102015224054B4 (de) 2015-12-02 2017-11-23 Julius-Maximilians-Universität Würzburg Modifizierte TrueFISP-Sequenz zur parallelen MR-Daten-Erfassung
US10466381B2 (en) * 2015-12-28 2019-11-05 Baker Hughes, A Ge Company, Llc NMR logging in formation with micro-porosity by using first echoes from multiple measurements
DE102016200603B4 (de) 2016-01-19 2018-02-01 Siemens Healthcare Gmbh Mehrschicht gradientenecho magnetresonanz-bildgebung
DE102016217223A1 (de) * 2016-09-09 2018-03-15 Siemens Healthcare Gmbh Überprüfung einer zeitlichen Änderung eines Magnetfeldes in einer Magnetresonanzvorrichtung
CN106443533B (zh) * 2016-09-21 2019-08-09 清华大学 基于多次激发的导航磁共振扩散成像方法及装置
CN107037386B (zh) * 2016-11-01 2019-08-23 上海联影医疗科技有限公司 一种平面回波成像方法以及系统
EP3388855A1 (de) * 2017-04-12 2018-10-17 Siemens Healthcare GmbH Vorrichtung und verfahren zur rückgewinnung des zeitlichen bezuges in freilaufenden mr-empfangsketten
CN109901088B (zh) * 2017-12-11 2023-08-22 通用电气公司 用于磁共振成像的运动追踪方法、计算机程序、存储设备
US10816624B2 (en) * 2018-03-26 2020-10-27 Siemens Healthcare Gmbh Method and device for correcting a B0 inhomogeneity by a radio frequency signal
EP3627172B1 (de) 2018-09-18 2022-02-09 Siemens Healthcare GmbH Verfahren und gerät für das mrt-schicht-multiplexing
DE102018216774A1 (de) 2018-09-28 2020-04-02 Siemens Healthcare Gmbh Korrekturverfahren für Schicht-Multiplexing-EPI-Verfahren
CN109521383B (zh) 2018-10-17 2019-08-30 浙江大学 一种基于频率稳定模块的磁共振cest成像序列及装置
CA3117901A1 (en) * 2018-10-26 2020-04-30 Aspect Imaging Ltd. Systems and methods for mri motion correction during mri image acquisition
US11163029B2 (en) * 2019-08-14 2021-11-02 GE Precision Healthcare LLC MRI system with improved navigator
CN110473271B (zh) 2019-08-20 2022-12-06 上海联影医疗科技股份有限公司 一种图像数据处理方法、系统、装置及存储介质
EP4081820A4 (en) * 2019-12-26 2024-01-24 Brigham & Womens Hospital Inc SYSTEM AND METHOD FOR CONTROLLING PHYSIOLOGICAL NOISE IN IMAGE PRODUCTION USING MAGNETIC RESONANCE IMAGING
CN113466765A (zh) * 2020-03-31 2021-10-01 通用电气精准医疗有限责任公司 磁共振扫描方法及系统、计算机可读存储介质
CN113625209B (zh) * 2020-05-09 2024-02-27 上海联影医疗科技股份有限公司 磁共振系统频率漂移量的确定方法、装置和计算机设备
DE102020209382A1 (de) * 2020-07-24 2022-01-27 Siemens Healthcare Gmbh Verfahren zur Aufnahme von Messdaten mittels einer Magnetresonanzanlage mit einer Korrektur der verwendeten k-Raumtrajektorien
DE102020209911A1 (de) 2020-08-05 2022-02-10 Siemens Healthcare Gmbh Korrektur von durch Fluktuationen eines Grundmagnetfelds verursachten Einflüssen auf eine Magnetresonanztomographie eines Untersuchungsobjekts
WO2023034044A1 (en) * 2021-08-30 2023-03-09 Children's Medical Center Corporation Dynamic distortion correction for mri using fid navigators
DE102022202094A1 (de) 2022-03-01 2023-09-07 Siemens Healthcare Gmbh Bildrekonstruktion aus Magnetresonanzmessdaten mit einer trainierten Funktion
DE102022207892B4 (de) 2022-07-29 2024-03-07 Siemens Healthcare Gmbh Verfahren zur schichtspezifischen Korrektur von mittels einer echo-planaren simultanen-Mehrschicht-Technik simultan für mindestens zwei Schichten aufgenommenen Messdaten
DE102022207891A1 (de) 2022-07-29 2024-02-01 Siemens Healthcare Gmbh Verfahren zur schichtspezifischen Korrektur von mittels einer echo-planaren simultanen-Mehrschicht-Technik simultan für mindestens zwei Schichten aufgenommenen Messdaten

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1216242A (zh) * 1997-10-17 1999-05-12 日立医药株式会社 磁共振成像方法及实施该方法的装置
CN1378816A (zh) * 2001-04-05 2002-11-13 Ge医疗系统环球技术有限公司 相位校正方法和磁共振成像系统
US6586935B1 (en) * 2000-03-31 2003-07-01 Ge Medical Technology Services, Inc. Magnetic resonance image artifact correction using navigator echo information
US6853191B1 (en) * 2003-12-10 2005-02-08 The Board Of Trustees Of The Leland Stanford Junior University Method of removing dynamic nonlinear phase errors from MRI data
CN1625366A (zh) * 2002-02-01 2005-06-08 株式会社日立医药 核磁共振成像方法及装置

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10117752C1 (de) 2001-04-09 2003-02-27 Siemens Ag Kernspintomographievorrichtung mit einer Einrichtung zur Bewegungskorrektur
DE10330926B4 (de) 2003-07-08 2008-11-27 Siemens Ag Verfahren zur absoluten Korrektur von B0-Feld-Abweichungen in der Magnetresonanz-Tomographie-Bildgebung
US7408345B2 (en) * 2006-02-06 2008-08-05 The Board Of Trustees Of The Leland Stanford Junior University Generalized MRI reconstruction with correction for multiple image distortion
US8406849B2 (en) * 2006-03-31 2013-03-26 University Of Utah Research Foundation Systems and methods for magnetic resonance imaging
US8483457B2 (en) * 2010-07-07 2013-07-09 General Electric Company System and method of image artifact reduction using self-navigated real-time phase correction in echo planar imaging
JP5835989B2 (ja) * 2010-08-11 2015-12-24 株式会社東芝 磁気共鳴イメージング装置および磁気共鳴イメージング方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1216242A (zh) * 1997-10-17 1999-05-12 日立医药株式会社 磁共振成像方法及实施该方法的装置
US6586935B1 (en) * 2000-03-31 2003-07-01 Ge Medical Technology Services, Inc. Magnetic resonance image artifact correction using navigator echo information
CN1378816A (zh) * 2001-04-05 2002-11-13 Ge医疗系统环球技术有限公司 相位校正方法和磁共振成像系统
CN1625366A (zh) * 2002-02-01 2005-06-08 株式会社日立医药 核磁共振成像方法及装置
US6853191B1 (en) * 2003-12-10 2005-02-08 The Board Of Trustees Of The Leland Stanford Junior University Method of removing dynamic nonlinear phase errors from MRI data

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Correction of physiologically induced global off-resonance effects in dynamic echo-planar and spiral functional imaging;Josef Pfeuffer et al.;《Magnetic Resonance in Medicine》;20020228;第47卷(第2期);344-353 *

Also Published As

Publication number Publication date
DE102011006230A1 (de) 2012-10-04
KR20120110056A (ko) 2012-10-09
US20120249138A1 (en) 2012-10-04
CN102697501A (zh) 2012-10-03
JP2012205897A (ja) 2012-10-25
KR101657331B1 (ko) 2016-09-13
US9329254B2 (en) 2016-05-03
DE102011006230B4 (de) 2013-01-24

Similar Documents

Publication Publication Date Title
CN102697501B (zh) 校正磁共振图像中的相位信息的方法和相应的磁共振设备
US9687172B2 (en) System for motion corrected MR diffusion imaging
JP4229487B2 (ja) マクスウェル項誤差を補正する方法
CN103238082B (zh) 使用多点Dixon技术和低分辨率校准的MR成像
US8093894B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method for improving uniformity in sensitivity map
US9983284B2 (en) MRI with dixon-type water/fat separation and prior knowledge about inhomogeneity of the main magnetic field
CN102342833B (zh) 用于建立磁共振图像的方法和相应的磁共振设备
US10203394B2 (en) Metal resistant MR imaging
CN104068859A (zh) 用于确定多个磁共振图像的方法和磁共振设备
CN104101853A (zh) 用于确定共振频率偏差的方法和磁共振设备
US11585884B2 (en) Continual trajectory correction in magnetic resonance imaging
JP5828763B2 (ja) 磁気共鳴イメージング装置および磁気共鳴イメージング方法
CN103018694A (zh) 校正失真的方法、磁共振设备、程序产品以及数据载体
US9846215B2 (en) MRI embodiments for controlling an arrangement order of multiple echoes in a k-space
CN106419917A (zh) 磁共振成像中用体积导航做前瞻性运动校正的方法和装置
EP1102082B1 (en) Method and apparatus for reducing image artifacts caused by magnet vibration in an MR imaging system
EP2635914B1 (en) Method of characterizing an rf transmit chain
CN109983358A (zh) Propeller mr成像
US11402454B2 (en) Correction of distorted diffusion-weighted magnetic resonance image data
US20220057467A1 (en) Epi mr imaging with distortion correction
CN102565736B (zh) 确定层关于相对其运动的部位的位置的方法和磁共振设备
WO2012143540A1 (en) Determining positions of a magnetic field probe in a magnetic resonance measurement
JP2024506575A (ja) Dixonタイプ水/脂肪分離MR画像
US20190128979A1 (en) Magnetic resonance imaging apparatus and pulse setting method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20220127

Address after: Erlangen

Patentee after: Siemens Healthineers AG

Address before: Munich, Germany

Patentee before: SIEMENS AG

TR01 Transfer of patent right