CN1576875B - 在自旋回波图像中避免外围干扰的方法及磁共振设备 - Google Patents

在自旋回波图像中避免外围干扰的方法及磁共振设备 Download PDF

Info

Publication number
CN1576875B
CN1576875B CN200410054959XA CN200410054959A CN1576875B CN 1576875 B CN1576875 B CN 1576875B CN 200410054959X A CN200410054959X A CN 200410054959XA CN 200410054959 A CN200410054959 A CN 200410054959A CN 1576875 B CN1576875 B CN 1576875B
Authority
CN
China
Prior art keywords
layer
gradient
frequency
pulse
amplitude
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.)
Expired - Fee Related
Application number
CN200410054959XA
Other languages
English (en)
Other versions
CN1576875A (zh
Inventor
彼得·休贝丝
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens 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 CN1576875A publication Critical patent/CN1576875A/zh
Application granted granted Critical
Publication of CN1576875B publication Critical patent/CN1576875B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, 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
    • 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

Landscapes

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

Abstract

本发明涉及一种用于在采用自旋回波序列的磁共振断层造影中避免外围噪声信号的方法,该自旋回波序列具有高频激励脉冲、高频重聚焦脉冲、层选择梯度脉冲、相位编码梯度脉冲和读出梯度脉冲,其中,这样区分高频激励脉冲的平均频率(ω1)和带宽(Δω1)与高频重聚焦脉冲的平均频率(ω2)和带宽(Δω2),和在高频激励脉冲期间接通的层选择梯度的振幅(g1)与高频重聚焦脉冲期间接通的层选择梯度的振幅(g2):将HF激励脉冲的激励层和HF重聚焦脉冲的重聚焦层在MRT磁场的均匀区域(FOV)内重叠,而使HF激励脉冲的激励层和HF重聚焦脉冲的重聚焦层在MRT磁场的非均匀区域内空间上分离,由此避免在非均匀区域内出现回波信号。

Description

在自旋回波图像中避免外围干扰的方法及磁共振设备
技术领域
本发明总的涉及在医疗中用于检查患者的核自旋断层造影(同义词:磁共振断层造影,MRT)。在此,本发明尤其涉及一种用于在自旋回波图像中避免外围干扰信号、如双重伪影的方法。
背景技术
MRT基于核自旋共振的物理现象,作为成像方法成功地应用在医疗和生物物理领域中已有超过15年的历史。在这种检查方法中,将对象置于一个强的稳定磁场中。由此使对象体内先前无规则旋转的原子核的自旋定向。高频波可以将这些“定向”的核自旋激励为特定的振荡。该振荡在MRT中产生实际的测量信号,并借助合适的接收线圈接收。采用由梯度线圈产生的非均匀磁场,可以在所有3个方向对各感兴趣区(也称为视场FOV)中的测量对象进行空间编码,这一般称为“位置编码”。
在MRT中数据的记录是在所谓的k域(同义词:频域)中进行的。所谓图像域中的MRT图像借助付立叶变换与k域的MRT数据关联。在k域中展开的对象的位置编码借助所有3个方向的梯度进行。在此,要区分层选择(确定对象中要拍摄的层,通常为Z轴)、频率编码(确定该层的方向,通常为x轴)和相位编码(确定在该层内的第二维,通常为y轴)。
因此,首先通过层选择梯度Gs或Gz选择性地在例如z方向激励一个断层。借助两个业已提到的正交梯度场GR或GP,通过组合相位和频率编码来对该层中的位置信息进行编码,所述正交梯度场对于在z方向上激励的断层的例子来说,是通过同样已经提到的在x和y方向上的梯度线圈产生的。
在图2a和2b中示出在MRT实验中记录数据的可能形式。所采用的序列是自旋回波序列。在这种序列中,通过90°激励脉冲(具有确定的振幅和带宽)将自旋的磁化翻转到x-y平面上。随着时间的变化,出现构成x-y平面Mxy上的横向磁化的磁化分量的相位差。在一定时间(例如1/2TE,TE是回波时间)之后在x-y平面上这样入射一个180°脉冲(同样具有确定的振幅和带宽),使得有相位差的磁化分量被反射,而不改变各磁化分量的精确方向和精确速度。在另一个持续时间1/2TE之后,所述磁化分量又指向相同的方向,即重新产生横向磁化,称为“再相位化”,该横向磁化可以相应地通过读出来采集。横向磁化的完全再生称为自旋回波。
为了测量待检查对象的整个断层,针对相位编码梯度GP和GY的不同值将成像序列重复N次,其中,在每个序列以等距离的时距Δt通过以Δt为时钟的ADC(模拟数字转换器)N次时,将核共振信号(自旋回波信号)的频率在具有读出梯度GR和Gx的情况下进行扫描、数字化和存储。通过这种方式,根据图2b获得逐行产生的数据矩阵(k域中的矩阵或k矩阵),该数据矩阵具有N×N个数据点(具有N×N个点的对称矩阵只是一个示例,也可以产生非对称矩阵)。可以通过付立叶变换,从该数据组中以N×N个点的分辨率直接再现有关断层的MR图像。
对于自旋回波序列,为了具有可用于诊断的图像质量,对k矩阵(在拍摄多个断层时为多个k矩阵)的扫描需要若干分钟的测量时间,而这对于很多临床应用是个问题。例如患者无法在所需的时间段内保持不动。在检查胸腔和骨盆部位时,骨骼的运动一般是无法避免的(心脏和呼吸运动,蠕动)。一种用于加速自旋回波序列的途径在1986年作为Turbo自旋回波序列(TES序列)或以缩写词RARE(驰豫增强的快速采集)公开(J.Hennig等人,Magn.Reson.Med.3,823-833,1986)。在这种成像方法中(与上述传统自旋回波方法相比更为快速),在一个90°激励脉冲之后产生多个多重回波,其中对这些回波中的每一个都单独进行相位编码。在图3a中示出在分别产生7个回波的情况下的相应序列图。在每个回波之前和之后,必须根据待选择的付立叶行接通相位编码梯度。通过这种方式,在唯一的HF激励脉冲(90°)之后,对k矩阵进行逐行扫描,如图3b所示。在该例中,所需的总测量时间缩短了七分之一。在图3a中示出了信号的理想变化。实际上,后面的回波通过横向磁化的T2衰减而具有逐渐减小的振幅。
一个更快的成像序列表示RARE和半付立叶技术的组合,其在1994年作为所谓的HASTE序列(半付立叶采集的单点快速自旋回波)公开(B.Kiefer等人,J.Magn.Reson.Imaging,4(P),86,1994)。HASTE采用与RARE相同的基本技术,但只扫描k矩阵的一半。k矩阵的另一半借助半付立叶算法,以计算方式重建。在此,利用了k矩阵的数据点关于k矩阵的中点镜像对称这个事实。因此,只要测量一半k矩阵的数据点、并通过在中点上的反射(和复数共轭)以计算方式补充原始矩阵就足够了。通过这种方式,可以将测量时间缩短一半。但是拍摄时间的减少是与信噪比(S/R)恶化了倍紧密相连的。
在自旋回波序列(SE序列)中一般存在这样的问题,在通过高频脉冲进行高频激励期间,共振条件不仅在FOV(视场,以基本磁场的均匀性和梯度场的线性为特征),还在FOV的非均匀边界区域中给出。由于在MRT设备的边界区域中基本磁场实际的非均匀性和梯度场的非线性,损害了每个空间点与一个确定磁场强度的唯一可反转对应的原则。这意味着,按照伪影的方式给实际测量场的图像叠加了一幅来自非均匀区域的、一般具有干扰的图像。该不期望的伪影称为“双重伪影”,并尤其是在自旋回波序列中由于自旋重聚焦而以显现出的方式出现。基本磁场在z方向上的延伸越短,“双重伪影”的可能性就越大。因此,在将来趋向于更短磁铁的MRT系统中,这个问题会加剧,并且不能再用当前用于抑制伪影的措施解决。
目前用于减少伪影类型的策略一方面在于硬件措施,另一方面在于脉冲序列的修改。
在HF系统中的硬件措施在于,在给定磁铁和梯度设计时,用有用空间之外的双重场来确定空间位置。然后,将HF线圈的实现置于这样一个限制下,即在该临界空间位置上使其灵敏度足够最小化,使得无法形成明显的伪影。但是HF场分布不是可任意形成的。由此,作为不期望的负效应,还在有用空间内出现损害图像质量的HF场非均匀性。
设计未来MR设备的重点是更短的磁铁、更大的空间和尽可能大的患者进入时间(例如由于介入)。在一个这样的磁场几何中,损害了必需的、空间与场的可反转的唯一对应关系,并加剧了由此产生的伪影问题,使得目前的硬件措施无用。因此,对于具有大直径的短磁铁,无法预见有用的HF线圈设计。
脉冲序列修改(建立序列的新途径)通常是唯一的实用解决办法,如果直接避免触及原理或技术限制边界的弊病(这里是非单调磁场变化)或所需要的很大花费从而使该产品的经济性受到质疑的弊病的话。
一个根据US 6486668的脉冲序列修改的可能形式是,通过施加或入射附加的所谓准备脉冲来起到抑制伪影的作用。这种措施的缺点在于,明显减小了时间效率,同时通过准备脉冲而按照寄生自旋回波信号分量的形式产生了另外的图像质量问题。
脉冲序列修改的另一种形式在US 2002/0101237中实现,其中通过HF激励脉冲和HF重聚焦脉冲的选择梯度的极性变换来抑制伪影。与传统的用SE序列激励断层的方法不同,在层激励期间通过(90°)HF脉冲接通的层选择梯度与(180°)重聚焦脉冲期间接通的层选择梯度相比,前者的符号或极性是相反的。这使得在位置空间中,在不同的非重叠区域内出现(90°)HF脉冲以及(180°)重聚焦脉冲的错误选择(Fehlselektion)。通过这种方式无法形成干扰回波信号。但是这种方法的缺点在于,它提高了HF脉冲和梯度脉冲要在时间上精确同步的技术要求,以及对系统补偿(Shimmung)的要求。此外,只有在信号损失明显时才能同时显示不同的化学成分(例如脂肪和水)。
发明内容
因此本发明要解决的技术问题在于,提供一种新的成像方法或在自旋回波序列的框架内进行脉冲序列修改,使得外围干扰信号(例如双重伪影)被抑制,同时具有最广泛的应用和出现最少的缺陷性的并因此不期望的负效应。
根据本发明的用于在采用自旋回波序列的磁共振断层造影中避免外围噪声信号的方法,该自旋回波序列具有高频激励脉冲、高频重聚焦脉冲、层选择梯度脉冲、相位编码梯度脉冲和读出梯度脉冲,其中,这样区分高频激励脉冲的平均频率和带宽与高频重聚焦脉冲的平均频率和带宽,以及在高频激励脉冲期间接通的层选择梯度的振幅与高频重聚焦脉冲期间接通的层选择梯度的振幅,即将HF激励脉冲的激励断层和HF重聚焦脉冲的重聚焦层在MRT磁场的均匀区域(FOV)内重叠,而使HF激励脉冲的激励断层和HF重聚焦脉冲的重聚焦断层在MRT磁场的非均匀区域内空间上分离,并由此避免在非均匀区域内出现回波信号。
在本发明方法的第一实施方式中,提高HF激励脉冲的平均频率、带宽和振幅以及在HF激励脉冲期间接通的层选择梯度的振幅,同时与HF重聚焦脉冲期间接通的层选择梯度的持续时间相比,相应减少在HF激励脉冲期间接通的层选择梯度的持续时间。
在本发明方法的第二实施方式中,保持HF脉冲以及层选择梯度脉冲的总持续时间,并与在HF重聚焦脉冲期间接通的层选择梯度的振幅相比,提高HF激励脉冲的平均频率、带宽和振幅以及在HF激励脉冲期间接通的层选择梯度的振幅。
优选的,将HF激励脉冲的不变总持续时间用于改善被激励层的层特性,其中,补充其它的HF激励脉冲的旁瓣最大值。
在横向断层选择(GS=Gz)以及径向或冠状断层选择(GS≠Gz)的情况下,按照差异因子的形式计算层选择梯度振幅的最小差异。
在HF激励脉冲期间接通的层选择梯度的振幅g1与在HF重聚焦脉冲期间接通的层选择梯度的振幅g2相比较的差别可以根据算式g2=(1+ε)g1来公式化。在横向层选择的情况下,根据以下不等式计算差异因子ε:
ϵ > 1 ( R - r ) / d - 1 / 2
其中R表示FOV的半径,r表示在FOV中受激励的层的层位置,d表示在FOV中受激励横向层的层厚。
在径向或冠状层选择的情况下,根据以下不等式计算差异因子ε:
ϵ > 1 - r / d - 1 / 2
其中r表示在FOV中受激励层的层位置,d表示在FOV中受激励的径向或冠状层的层厚。
此外还提供了一种磁共振断层造影设备,包括具有附属梯度线圈的梯度放大器、输入-显示终端、序列控制器和设备计算机以及模拟数字转换器(ADC),该磁共振断层造影设备适用于实施按照上述方法特征的方法。
附图说明
下面借助基于附图所示的实施例详细解释本发明的其它优点、特征和特性。
图1示意性示出核自旋断层造影设备,
图2a示意性示出公知的自旋回波序列的梯度脉冲流函数的时间变化过程,
图2b示意性示出通过根据图2a的自旋回波序列对k矩阵的时间扫描,
图3a示意性示出公知快速自旋回波序列(TSE或RARE序列)的梯度脉冲流函数的时间变化过程,其中产生7个自旋回波,
图3b示意性示出在按照图3a的TSE或RARE序列中对k矩阵的时间扫描,
图4a示意性示出HF脉冲的时间变化过程以及本发明第一修改SE序列的梯度脉冲流函数的时间变化过程,
图4b示意性示出HF脉冲的时间变化过程以及本发明第二修改SE序列的梯度脉冲流函数的时间变化过程,
图5a示意性示出在具有相同选择梯度强度(左侧)以及较少但足够的不同选择梯度强度(右侧)的横向层选择时在z方向上的磁场情况,
图5b示意性示出在按照图5a的出现双重伪影的情况(左侧图)和按照图5a的按照本发明抑制双重伪影的情况(右侧图)进行横向层定向时的多层测量的叠加的总和图像,
图6a示意性示出在具有相同选择梯度强度(左侧)以及较少但足够的不同选择梯度强度(右侧)下的径向或冠状层选择时在z方向上的磁场情况,
图6b示意性示出在按照图6a的出现双重伪影的情况(左侧图)和按照图6a的按照本发明抑制双重伪影的情况(右侧图)进行径向或冠状层定向时的多层测量的叠加的总和图像,
图7a示意性示出在两个根据本发明的较少但足够的不同层选择梯度(左侧是分离的磁场和梯度场,右侧是所有的场)时的横向层选择的简化场模型,
图7b示意性示出在横向层选择中,对于不同的层厚取决于到FOV中心的距离的参数ε的曲线簇,
图8a示意性示出在径向或冠状层选择(左侧)时以及在根据本发明的考虑较少但足够的不同层选择梯度(右侧)下的FOV外寄生层选择的基本磁场磁铁的简化场模型,
图8b示意性示出在径向或冠状层选择中,对于不同的层厚取决于到FOV中心的距离的参数ε的曲线簇。
具体实施方式
图1示出用于产生根据本发明的梯度脉冲的核自旋断层造影设备的示意图。在此,核自旋断层造影设备的结构对应于传统断层造影设备的结构。基本磁场磁铁1产生时间上稳定的强磁场,用于极化或定向对象检查区域内的核自旋,该对象例如是人体的待检查部分。在圆形测量空间域M中定义核自旋共振测量所需的基本磁场的高均匀性,人体的待检查部分被送入该测量区域。为了支持均匀性要求,特别是为了消除不随时间变化的影响,在合适的位置上设置由铁磁材料制成的所谓填隙片。随时间变化的影响可以通过由补偿电源15控制的补偿线圈2消除。
在基本磁场磁铁1中采用圆柱形的梯度线圈系统3,其由3个部分绕组构成。每个部分绕组都由一个放大器14提供电流,从而在笛卡尔坐标系的3个方向上分别产生一个线性梯度场。梯度场系统3的第一部分绕组产生x方向上的梯度Gx,第二绕组产生y方向上的梯度Gy,第三绕组产生z方向上的梯度Gz。每个放大器14包括数字模拟转换器,其由序列控制器18控制,以及时产生梯度脉冲。
在梯度场系统3内设置了高频天线4,该高频天线4将高频功率放大器30输出的高频脉冲转换为交变磁场,以激励待检查对象或对象的待检查区域中的原子核以及使核自旋定向。高频天线4还将由确定的核自旋发出的交变场、即一般由一个或多个高频脉冲和一个或多个梯度脉冲组成的脉冲序列引起的核自旋回波信号,转换为电压,该电压通过放大器7输入高频系统22的高频接收信道8。高频系统22还包括一个发送信道9,其中产生用于激励磁核共振的高频脉冲。在此,根据设备计算机20预先给定的脉冲序列,在序列控制器18中将各高频脉冲数字化地表示为复数序列。该数列作为实部和虚部,分别通过输入端12输入到高频系统22中的数字模拟转换器,并由该数字模拟转换器输入到发送信道9。在发送信道9中,将该脉冲序列调制为高频载波信号,其基本频率对应于测量空间内的核自旋的共振频率。
通过发送-接收转接器6将发送运行切换到接收运行。高频天线4将用于激励核自旋的高频脉冲入射到测量空间M中,并对产生的回波信号进行扫描。在高频系统22的接收信道8中对相应获得的核共振信号进行相敏解调,并通过模拟数字转换器分别转换为测量信号的实部和虚部。通过图像计算机17,从获得的测量数据中再现图像。测量数据、图像数据和控制程序的管理都通过设备计算机20进行。根据预先给定的控制程序,序列控制器18检查各期望的脉冲序列的产生以及相应的k域扫描。特别的,序列控制器18在此控制梯度的及时接通、具有限定相位和振幅的高频脉冲的发送以及核共振信号的接收。高频系统22和序列控制器18的时间基准由合成器19提供。通过终端21选择用于产生核自旋图像的相应控制程序,并显示所产生的核自旋图像,该终端21包括键盘以及一个或多个显示屏。
本发明在于在自旋回波序列下产生新的HF激励机制。在此存在随着距离增大的处理要求,因为所有成像过程的大多数都建立在该序列类型的基础上,同时双重伪影由于自旋重聚焦而以显现的方式出现。本发明涉及这样的措施,其在只进行边缘的序列修改时保证最大有效地抑制伪影。在序列控制器18或设备计算机20中产生根据本发明的序列变化。同样在设备计算机20中进行图像再现。
本发明的思想基本上如下所述:为了产生自旋回波需要至少两个层选择过程。与通常的实践相反,这些层选择过程是不相同的,而是考虑选择梯度强度和HF脉冲带宽而不同地进行,其预先给定的条件是空间分离FOV(也称为“有用立体区域”)外的错误选择的立体区域。通过适当选择平均频率ω1、ω2以及HF激励脉冲和HF重聚焦脉冲在同一地点(同一层)的带宽Δω1、Δω2来实现在FOV内的选择,该选择产生或使得从该断层中获得期望的回波信号成为可能,如目前已知的。通过空间分离FOV外的错误选择的立体区域(寄生激励)(根据ω1、ω2、Δω1、Δω2、以及HF激励脉冲期间的层选择梯度GS1和HF重聚焦脉冲期间的层选择梯度GS2的不同振幅g1和g2进行选择,图5a和6a),无法在非均匀区域内形成回波信号,并抑制了图像中的双重伪影。优选的,在序列开始时(90度脉冲)这样改变HF自旋激励,使得较少地(20%至30%)提高选择梯度强度和HF脉冲带宽。
在此,尽可能少地选择差异,以保证广泛的应用和最小化不期望的负面效应。这种要考虑的负面效应主要是化学移动和SAR(特殊吸收率)方面:
由于3.4ppm的化学脂肪-水移动,对于脂肪和水来说层位置始终不同。对于两个化学成分中的每一个,激励特性和所有重聚焦特性只有当选择梯度强度不变时才是配合的。严格地说,这在下面介绍的序列设计中不再适用,即由于较少相对移动的信号特性而将不同化学成分的信号损失限制为很小的比例并由此使得该信号损失是无关紧要的。
在层厚或投影厚度增加时,完全分离有用立体区域外的错误选择的立体区域要求相应增加选择梯度强度和HF脉冲带宽的差别。通过一方面(上限)可用的HF功率和SAR方面,另一方面(下限)通过降低序列在例如回波时间和重复时间上的时间性能而限制HF带宽的有用范围。在更实用的例如100%的差别下,可以利用根据本发明的序列修改将层厚或投影厚度调整到最大约100mm而不会出现明显的伪影。
下面借助两个可能的实施例(图4a和4b)解释本发明:
两幅图4a和4b的序列图的上半部示出传统的SE序列。在最上一行示出发送模式TX下的平均频率为ω、带宽为Δω的初始α激励脉冲(HF激励脉冲,α一般为90°),同时正层选择梯度GS1接通。在HF激励脉冲之后,将层选择梯度反向。在此期间一方面接通读出梯度GR,另一方面通过相位编码梯度GP进行相位编码。接着,通过180°重聚焦脉冲,同时与通过层选择梯度GS2进行的第二次层激励一起进行重聚焦。对于传统的SE序列,两个层选择梯度GS1和GS2具有相同的振幅。为了能在后面的MRT系统的接收模式RX中,在回波时间TE之后通过接通读出梯度GR获得MR回波信号,需要将180度重聚焦脉冲的平均频率和带宽选择为分别等于HF激励脉冲的平均频率ω和带宽Δω。后面还将借助其它附图详细解释这种情况。
根据本发明的序列是修改的自旋回波序列(SE序列)。该修改在于这样改变HF激励脉冲(α脉冲)、HF重聚焦脉冲(180度脉冲)和层选择梯度(GS)的参数,使得在FOV内与以往一样激励一个断层,并由该层获得自旋回波,但在FOV外从空间上分离两个HF脉冲的寄生激励,并由此不产生自旋回波信号。
在本发明的第一可能实施方式(图4a)中,参数的改变例如是这样进行的,即改变两个HF脉冲的振幅、平均频率以及带宽。同时相应地改变层选择梯度GS1和GS2的持续时间和振幅。根据图4a将α脉冲在时间上进行可见的压缩,加大其振幅。该α脉冲具有平均频率ω1和带宽Δω1。180度脉冲具有平均频率ω2和带宽Δω2。缩短第一层选择梯度GS1,其振幅g1相对于第二层选择梯度GS2的振幅g2增大。
在本发明的第二可能实施方式(图4b)中,参数的改变例如是这样进行的,即保持α脉冲和第一层选择梯度GS1的总持续时间不变。脉冲振幅g1相对于GS2的振幅g2增大。压缩α脉冲,但是最初的脉冲持续时间保持不变,其中给激励脉冲的类似正弦函数又补充了几个旁瓣最大值(Nebenmaxima,英语:Side-Lobes),并具有改善的层特性的积极旁效应。这对SAR值只有很小的影响,因为在此未经改变的重聚焦脉冲占据主导地位。相应地改变平均频率ω1和ω2以及带宽Δω1和Δω2,并相应地相互调谐。
图5a和6a(分别是左侧的图像)显示了对于传统自旋回波序列在FOV内或FOV外的期望的或寄生的层激励。由于人的体形和通常患者躺下的位置,场双重效应主要出现在身体的长轴上(Z方向)。这产生了情况差异,即是利用Z梯度(横向断层)还是不用Z梯度(径向和冠状断层)进行断层选择的差异。图5a示出横向层选择(GS=Gz),图6a示出径向或冠状层选择(GS≠Gz)。
在横向层激励(图5a)中,层选择梯度GS(=GZ)在FOV区域内首先线性变化。接着的剧烈的场下降既是由于基本磁场的非均匀性,又是由于梯度场的非线性,并引起寄生激励。
在径向或冠状层激励(图6a)中,在z方向上不存在层选择梯度。但是,由于基本磁场在FOV外的剧烈下降,也产生了寄生激励。
总之,这产生了这种现象,即在利用传统2D自旋回波进行测量时,对于横向层定向和径向或冠状层定向都会在图像中的不同位置出现双重伪影。该事实在图5b和6b中示出,其中利用传统2D自旋回波进行的多层测量的总和图像重叠显示。图5b(左侧)示出横向层定向时的双重伪影,图6b(左侧)示出径向或冠状层定向时的双重伪影。如果在信号产生和信号指示过程中还考虑不同的场情况,则可以看到,根据各自的层位置而在图像中出现不同的伪影。后者影响图像中伪影的位置,并且不应在后续变化中加深,因为本发明只考虑层选择过程。
如上所述(图4a和4b的描述),利用不同的选择梯度强度进行层选择,也就是至少利用第一梯度强度g1和第二梯度强度g2:
g2=(1+ε)g2                                    (1)
通常的自旋回波序列采用ε=0。ε=-2的情况对应于上面提到的并在现有技术中已经采用的序列修改(选择梯度的极性转换)。在本发明中,优选的序列特性在选择ε>0时表现出来。
在激励和重聚焦特定断层时HF脉冲的平均频率ω和带宽Δω当然必须对同一位置区域起作用,该位置区域通过期望的层位置r和层厚d预先给出,即
ω1=γg1r                    (2a)
ω2=γg2r=(1+ε)ω        (2b)
Δω1=γg1d                  (3a)
Δω2=γg2d=(1+ε)Δω1     (3b)
图5a就激励和重聚焦(右侧图像)将传统横向断层激励(左侧图像)的情况与不同选择梯度强度时的横向断层选择的情况进行了比较。右侧的图像示出在没有形成双重伪影下较少但足够的不同选择梯度强度(g2=(1+ε)g2)。
图5b基于图5a中2D自旋回波测量例子的配置而示出这种序列修改的效果。在中度不同的层选择(ε=0.2至0.3)中,已经非常有效地抑制了这种情况下的伪影(图5b右侧)。
图6a关于激励和重聚焦(α脉冲或重聚焦脉冲;右侧图像)将传统径向或冠状层激励(左侧图像)的情况与不同HF脉冲时的径向或冠状层选择的情况进行比较。右侧图像示出在没有形成双重伪影下较少但足够的不同HF脉冲(ω2=(1+ε)ω1)。
图6b基于图6a中2D自旋回波测量例子的配置而示出这种序列修改的效果。在稍微不同的层选择(ε=0.2至0.3)中,已经有效地抑制了这种情况下的伪影(图6b右侧)。
下面,借助空间场分布的简单数学模型进一步解释按照目标对具有决定意义的参数ε的选择。
图7a示出横向断层选择的场模型。在图7a左侧图像的情况下(实际情况的理想化),Z梯度的场与Z坐标成正比地一直上升到FOV的边界,然后保持不变,而基本磁场首先为恒定的,并在FOV外(>R)以不变的场梯度-G下降。图7a中的右侧图像示出在选择梯度强度g1和g2不同时根据本发明的总和场的变化。在两种情况下都选择期望的层(位置r、厚度d),但有用立体区域外的寄生区域是不同的(具有厚度D1和D2的位置ρ1和ρ2)。
图7a中的理想化模型这样简化了实际情况,即在半径为R的FOV内假定场分布是理想的(均匀基本磁场和线性梯度场)。在>R的外部区域,梯度场不会继续上升,磁场下降应当以恒定的场梯度-G进行。
对于利用梯度强度g1和g2进行的两个选择过程(图7a,右侧图像),下式成立:
ω1=γg1r=γg1R-γGρ1            (4a)
ω2=γg2r=γg2R-γGρ2            (4b)
Δω1=γg1d=γGD1                 (5a)
Δω2=γg2d=γGD2                 (5b)
在位置ρ1和ρ2上进行的寄生选择的相对位移为:
ρ21=(R-r)(g2-g1)/G             (6)
目标是完全抑制双重伪影,也就是将具有厚度D1和D2的错误选择的立体区域完全分离:
ρ21>(D1+D2)/2                  (7)
从中,通过将方程(1)到(6)代入不等式(7)得出参数ε的条件:
ϵ > 1 ( R - r ) / d - 1 / 2 - - - ( 8 )
其中r的值满足
r<R-d/2                            (9)
可以观察到,公式(8)仅取决于层和FOV的边界之间以层厚d为单位测量的距离(R-r)。该距离越大,对ε的要求越小,这可以借助图7b中的双曲线函数族示出。图7b示出在横向层选择时的ε值变化:如果参数ε超过所显示的最小值,则对应于图7a的场模型完全抑制双重伪影。利用通常层厚的2D自旋回波方法满足ε=0.2到0.3时的该条件。
FOV边界的极点对应于双重性在该点的消失(理想断层和错误选择在此相互转变(uebergehen))。
对g1、g2和G的不依赖性也是值得注意的,即在外部区域中的梯度强度和HF脉冲带宽的绝对大小以及磁场下降是不起作用的,可以任意构造。
为完整起见还针对径向或冠状断层定向示出了具有上述断言的该模型。图8a示出涉及该层定向的场模型。在图8a的左侧图中,基本磁场首先是恒定的,并在FOV外(>R)以不变的场梯度-G下降。在图8a的右侧图中,在FOV外选择不同的寄生区域(具有厚度D1和D2的位置ρ1和ρ2),象在期望的层外示出的那样(该层由于X或Y相关而未在此示出)。
由于所指出的措施只涉及层选择过程,并且在此Z梯度没有参与,因此不用考虑有关的假设。又假定磁场的变化在有用立体区域内恒定,而在>R的外部区域内以恒定的场梯度-G描述。
对于具有梯度强度g1和g2的两个选择过程下式成立:
ω1=γg1r=-γGρ1                    (10a)
ω2=γg2r=-γGρ2                    (10b)
Δω1=γg1d=γGD1                    (11a)
Δω2=γg2d=γGD2                    (11b)
在位置ρ1和ρ2上进行的寄生选择的相对位移为:
ρ21=-r(g2-g1)/G                   (12)
目标还是完全抑制双重伪影,也就是将具有厚度D1和D2的错误选择的立体区域完全分离:
ρ21>(D1+D2)/2                     (13)
从中,通过将方程(10a)到(12)代入不等式(13)得出参数ε的条件:
ϵ > 1 - r / d - 1 / 2 - - - ( 14 )
其中r的值满足
r<-d/2                                (15)
在这种情况下,参数ε取决于断层和磁场中心之间以层厚d为单元测量的距离r。该距离越大,对ε的要求越小。在图8b中示出径向或冠状层选择时的双曲线函数:如果参数ε超过所示出的最小值,则完全抑制对应于图8a的场模型的双重伪影。利用通常层厚的2D自旋回波方法满足ε=0.2到0.3时的条件。磁场中心的极点对应于双重性在该点的消失(理想断层和错误选择在此相互转变)。
下面指出该序列方法抑制伪影的边界:
在ε图(图7b和8b)中的极点不只是这种简单数学模型的现象。如已经示出的,这些极点描述了一种实际现象,即从双重伪影到几何符号的转变。该转变是在期望层和寄生选择的立体区域直接与磁场换向点相互紧密相邻时产生的。在重叠或完全一致时当然不可能实现分离,这就在形成数学模型时产生了极点。
因此,在该边界情况下不能完全避免伪影,为此这种现象的图像与典型的双重伪影有很大的偏差,并更近似于一种图像失真。该伪影结构与实际对象的对应是可以识别的,本身并不包含错误解释的风险。
由于简化了所采用的场模型而未能示出的另一方面在于FOV外的实际空间场分布。在横向断层选择时假定的、Z梯度在>R的外部区域中的恒定场变化至少对于>>R的大距离来说是不会遇到的。更合理的是,该梯度场自身又单独(即不涉及磁场)下降,并具有场双重效应。
但是这种情况并没有负面影响。同时假定的在>R的外部区域中的磁场下降也必须排除。
因此,仅在磁铁同时根本不具有场偏差的地方才有出现伪影的风险。取代目前的一维场变化,在此必须观察和考虑磁场的整个空间场分布。只有这种巧合与在该点具有足够高的HF灵敏度的发送和接收线圈一起,才会包含不能用所介绍的序列修改完全抑制双重伪影的风险。
如果需要,在设计HF线圈时必须考虑该方面,与传统自旋回波序列时的限制相反,该方面表示很少且可以满足的附加条件。

Claims (5)

1.一种用于在采用自旋回波序列的磁共振断层造影中避免外围噪声信号的方法,该自旋回波序列具有高频激励脉冲、高频重聚焦脉冲、层选择梯度脉冲、相位编码梯度脉冲和读出梯度脉冲,
其特征在于,这样区分高频激励脉冲的平均频率(ω1)和带宽(Δω1)与高频重聚焦脉冲的平均频率(ω2)和带宽(Δω2),以及区分在高频激励脉冲期间接通的层选择梯度的振幅(g1)与在高频重聚焦脉冲期间接通的层选择梯度的振幅(g2),即将高频激励脉冲的激励层和高频重聚焦脉冲的重聚焦层在磁共振断层造影磁场的均匀区域(FOV)内重叠,而使高频激励脉冲的激励层和高频重聚焦脉冲的重聚焦层在磁共振断层造影磁场的非均匀区域内空间上分离,并由此避免在非均匀区域内出现回波信号,以及保持所述高频脉冲以及层选择梯度脉冲的总持续时间,并相对于在高频重聚焦脉冲期间接通的层选择梯度(GS2)的振幅(g2),提高高频激励脉冲的平均频率(ω1)、带宽(Δω1)和振幅以及在高频激励脉冲期间接通的层选择梯度(GS1)的振幅(g1)。
2.根据权利要求1所述的方法,其特征在于,引入其它的高频激励脉冲的旁瓣最大值。
3.根据权利要求1或2所述的方法,其特征在于,在高频激励脉冲期间接通的层选择梯度(GS1)的振幅(g1)与在高频重聚焦脉冲期间接通的层选择梯度(GS2)的振幅(g2)相比较的差别可以根据关系式g2=(1+ε)g1来公式化,并且在横向层选择的情况下,根据以下不等式计算差异因子ε:
ϵ > 1 ( R - r ) / d - 1 / 2
其中R表示FOV的半径,r表示在FOV中受到激励的层的层位置,d表示在FOV中受到激励的横向断层的层厚。
4.根据权利要求1或2所述的方法,其特征在于,在高频激励脉冲期间接通的层选择梯度(GS1)的振幅(g1)与在高频重聚焦脉冲期间接通的层选择梯度(GS2)的振幅(g2)相比较的差别可以根据关系式g2=(1+ε)g1来公式化,并且在径向或冠状断层选择的情况下,根据以下不等式计算差异因子ε:
ϵ > 1 - r / d - 1 / 2
其中r表示在FOV中受到激励的层的层位置,d表示在FOV中受到激励的径向或冠状层的层厚。
5.一种磁共振断层造影设备,包括具有附属梯度线圈(3)的梯度放大器、输入-显示终端(21)、序列控制器(18)和设备计算机(20)以及模拟数字转换器,该磁共振断层造影设备适用于实施按照权利要求1所述的方法。
CN200410054959XA 2003-07-24 2004-07-26 在自旋回波图像中避免外围干扰的方法及磁共振设备 Expired - Fee Related CN1576875B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE10333795A DE10333795B4 (de) 2003-07-24 2003-07-24 Verfahren und Vorrichtung zur Vermeidung von peripheren Störsignalen in Spin-Echo-Bildern bei nicht monotonem Magnetfeldverlauf in der Magnetresonanz-Tomographie-Bildgebung
DE10333795.4 2003-07-24

Publications (2)

Publication Number Publication Date
CN1576875A CN1576875A (zh) 2005-02-09
CN1576875B true CN1576875B (zh) 2010-06-16

Family

ID=34071888

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200410054959XA Expired - Fee Related CN1576875B (zh) 2003-07-24 2004-07-26 在自旋回波图像中避免外围干扰的方法及磁共振设备

Country Status (5)

Country Link
US (1) US7091720B2 (zh)
JP (1) JP4546179B2 (zh)
KR (1) KR20050012172A (zh)
CN (1) CN1576875B (zh)
DE (1) DE10333795B4 (zh)

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5121219B2 (ja) * 2006-12-07 2013-01-16 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置および磁気共鳴イメージング方法
WO2008134891A1 (en) * 2007-05-03 2008-11-13 National Research Counsil Of Canada Rf based spatially selective excitation in mri
WO2010097094A1 (en) * 2009-02-25 2010-09-02 MAX-PLANCK-Gesellschaft zur Förderung der Wissenschaften e.V. Magnetic resonance imaging with chemical species selectivity
US8970217B1 (en) 2010-04-14 2015-03-03 Hypres, Inc. System and method for noise reduction in magnetic resonance imaging
DE102010062290A1 (de) * 2010-12-01 2012-06-06 Siemens Aktiengesellschaft Pulssequenz, Magnetresonanzanlage und Verfahren zum Erzeugen einer Pulssequenz sowie Computerprogrammprodukt und elektronisch lesbarer Datenträger
EP2500742A1 (en) * 2011-03-17 2012-09-19 Koninklijke Philips Electronics N.V. Restriction of the imaging region for MRI in an inhomogeneous magnetic field
DE102011082010B4 (de) 2011-09-01 2013-03-14 Siemens Aktiengesellschaft Rephasierung von Spinsystemen in einer ersten und zweiten Schicht in Schicht-Multiplexing-Messsequenzen zur Magnetresonanzbildgebung
EP2799849B1 (en) 2011-12-21 2020-12-02 Japan Science And Technology Agency Nmr imaging device and nmr imaging method
DE102012204434B3 (de) 2012-03-20 2013-07-11 Siemens Aktiengesellschaft Mehrschicht-MRI-Anregung mit simultaner Refokussierung aller angeregten Schichten
DE102012205587B4 (de) * 2012-04-04 2013-12-24 Siemens Aktiengesellschaft Schichtspezifische Phasenkorrektur bei Schicht-Multiplexing
DE102013201616B3 (de) 2013-01-31 2014-07-17 Siemens Aktiengesellschaft TSE-basierte, gegen lokale B0-Feldvariationen unempfindliche MR-Mulitschicht-Anregung
DE102013219754B4 (de) * 2013-09-30 2018-11-22 Siemens Healthcare Gmbh Optimierung der Geräuschentwicklung einer 3D-Gradientenchosequenz in einer Magnetresonanzanlage
CN105433944B (zh) * 2014-07-31 2018-07-03 西门子公司 用于获取对象的磁共振数据的方法及装置
JP6407620B2 (ja) * 2014-08-11 2018-10-17 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
US10132903B2 (en) * 2014-11-26 2018-11-20 Toshiba Medical Systems Corporation Magnetic resonance imaging method, magnetic resonance imaging apparatus and magnetic resonance imaging system
CN107536609B (zh) * 2016-06-28 2020-06-19 上海联影医疗科技有限公司 磁共振成像设备及磁共振成像的扫描方法
DE102016212632A1 (de) 2016-07-12 2018-01-18 Siemens Healthcare Gmbh Reduzierung von Artefakten in der Magnetresonanztechnik
WO2018175807A1 (en) * 2017-03-22 2018-09-27 Viewray Technologies, Inc. Reduction of artifacts in magnetic resonance imaging by creating inhomogeneity in the magnetic field at gradient null position of an mri system
CN108652624A (zh) * 2017-03-31 2018-10-16 西门子(深圳)磁共振有限公司 一种前瞻性心电触发方法和装置
US11454690B2 (en) * 2017-12-08 2022-09-27 Rensselaer Polytechnic Institute Synergized pulsing-imaging network (SPIN)
DE102018201810A1 (de) * 2018-02-06 2019-08-08 Siemens Healthcare Gmbh Multischicht-Bildgebung mit Binomialpulsen
WO2021247857A1 (en) * 2020-06-03 2021-12-09 The Brigham And Women's Hospital, Inc. System and methods for ultra-fast multi-dimensional diffusion-relaxation mri using time-division multiplexing sequences

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5386190A (en) * 1991-07-31 1995-01-31 Takeuchi; Hiroyuki Selective excitation method of nuclear spins in magnetic resonance imaging and MRI apparatus using the same
US5530354A (en) * 1994-07-29 1996-06-25 Medical Advances, Inc. Non-monotonic gradient coil system for magnetic resonance imaging
CN1286960A (zh) * 1999-06-24 2001-03-14 通用电器横河医疗系统株式会社 梯度磁场施加方法和设备以及磁共振成像设备
US6486668B1 (en) * 2000-11-20 2002-11-26 Ge Medical Systems Global Technology Company, Llc Preparatory pulse sequence for suppression of artifacts in MR images

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH025926A (ja) * 1988-06-24 1990-01-10 Hitachi Medical Corp Mri装置
US4952876A (en) * 1988-11-23 1990-08-28 General Electric Company Variable bandwidth multiecho NMR imaging
US4959611A (en) * 1989-06-29 1990-09-25 Picker International, Inc. Out-of-slice artifact reduction technique for magnetic resonance imagers
JP3033851B2 (ja) * 1991-03-25 2000-04-17 株式会社日立メディコ 磁気共鳴イメージング装置
DE19616388A1 (de) * 1996-04-24 1997-11-06 Siemens Ag Pulssequenz zur Selektion mehrerer Schichten in einem Kernspintomographiegerät
GB2348005A (en) * 1999-03-13 2000-09-20 Marconi Electronic Syst Ltd NMR apparatus with additional RF receive coil to reduce aliasing
DE19928110A1 (de) * 1999-06-19 2000-12-21 Philips Corp Intellectual Pty MR-Verfahren

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5386190A (en) * 1991-07-31 1995-01-31 Takeuchi; Hiroyuki Selective excitation method of nuclear spins in magnetic resonance imaging and MRI apparatus using the same
US5530354A (en) * 1994-07-29 1996-06-25 Medical Advances, Inc. Non-monotonic gradient coil system for magnetic resonance imaging
CN1286960A (zh) * 1999-06-24 2001-03-14 通用电器横河医疗系统株式会社 梯度磁场施加方法和设备以及磁共振成像设备
US6486668B1 (en) * 2000-11-20 2002-11-26 Ge Medical Systems Global Technology Company, Llc Preparatory pulse sequence for suppression of artifacts in MR images

Also Published As

Publication number Publication date
US7091720B2 (en) 2006-08-15
DE10333795A1 (de) 2005-03-03
CN1576875A (zh) 2005-02-09
US20050017719A1 (en) 2005-01-27
JP4546179B2 (ja) 2010-09-15
JP2005040612A (ja) 2005-02-17
DE10333795B4 (de) 2008-01-31
KR20050012172A (ko) 2005-01-31

Similar Documents

Publication Publication Date Title
CN1576875B (zh) 在自旋回波图像中避免外围干扰的方法及磁共振设备
US7495437B2 (en) Method for MRT imaging on the basis of conventional PPA reconstruction methods
Adriany et al. Transmit and receive transmission line arrays for 7 Tesla parallel imaging
Staroswiecki et al. In vivo sodium imaging of human patellar cartilage with a 3D cones sequence at 3 T and 7 T
CN101034153B (zh) 基于梯度回波序列的磁共振成像方法及设备
KR101625557B1 (ko) 자기 공명 시스템을 동작시키는 방법 및 제어 디바이스
US6445184B1 (en) Multiple gradient echo type projection reconstruction sequence for MRI especially for diffusion weighted MRI
CN104101853B (zh) 用于确定共振频率偏差的方法和磁共振设备
CN107430175B (zh) 具有多个独立的发射接收信道的磁共振体积线圈及其操作方法
KR101625715B1 (ko) 자기 공명 시스템을 제어하는 방법 및 제어 장치
CN1573352A (zh) 用于在磁共振断层造影成像中避免双重含义伪影的方法
US7254435B2 (en) Method and magnetic resonance apparatus for calibrating coil sensitivities
US7609058B2 (en) Method and apparatus for generating a magnetic resonance data file
CN105445684B (zh) 重建mr图像的方法、磁共振设备和数据载体
CN103083020A (zh) 磁共振成像设备及其控制方法
US7429862B2 (en) PPA MRT method and apparatus with scaled contribution of k-space data to the reconstructed image
JP2003010148A (ja) Mrイメージング方法およびmri装置
JPWO2004060156A1 (ja) 磁気共鳴イメージング装置
US5544653A (en) Magnetic resonance imaging method and apparatus
CN110215209B (zh) 一种磁共振成像方法和磁共振成像系统
US6906516B2 (en) Artifact reduction in SSFP MRI using weighted sum of combined signals
CN101201392B (zh) 产生磁共振数据文件的方法和设备
US9772390B2 (en) Magnetic resonance imaging device and method for generating image using same
US6369569B1 (en) Magnetic resonance tomography apparatus and operating method for displaying tissue contrasts with a short measuring time
Somai et al. A multi spin echo pulse sequence with optimized excitation pulses and a 3D cone readout for hyperpolarized 13C imaging

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20100616

Termination date: 20180726