CN100414316C - 能降低声学噪声的磁共振成像方法 - Google Patents

能降低声学噪声的磁共振成像方法 Download PDF

Info

Publication number
CN100414316C
CN100414316C CNB028232585A CN02823258A CN100414316C CN 100414316 C CN100414316 C CN 100414316C CN B028232585 A CNB028232585 A CN B028232585A CN 02823258 A CN02823258 A CN 02823258A CN 100414316 C CN100414316 C CN 100414316C
Authority
CN
China
Prior art keywords
pulse
signal
magnetic resonance
magnetic
sweep time
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
CNB028232585A
Other languages
English (en)
Other versions
CN1589410A (zh
Inventor
C·L·G·哈姆
P·R·哈维伊
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
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 Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN1589410A publication Critical patent/CN1589410A/zh
Application granted granted Critical
Publication of CN100414316C publication Critical patent/CN100414316C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

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/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • G01R33/385Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using gradient magnetic field coils
    • G01R33/3854Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using gradient magnetic field coils means for active and/or passive vibration damping or acoustical noise suppression in gradient magnet coil systems
    • 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/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/34Constructional details, e.g. resonators, specially adapted to MR
    • G01R33/341Constructional details, e.g. resonators, specially adapted to MR comprising surface coils
    • G01R33/3415Constructional details, e.g. resonators, specially adapted to MR comprising surface coils comprising arrays of sub-coils, i.e. phased-array coils with flexible receiver channels
    • 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/5611Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE

Landscapes

  • Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

一种磁共振成像方法包括应用包括一或多个脉冲的脉冲序列。该脉冲序列具有在k-空间内基于全采样速率按预定的全“观察场”的固有扫描时间和具有磁场梯度脉冲的参考时间脉冲的形状。利用具有空间灵敏度剖面图的接收器天线系统来采集一系列磁共振信号。运用欠采样信号采集在k-空间内以预定的降低了的采样速率来采集欠采样磁共振信号,其采样速率相对于全采样速率降低一个缩减系数。在实际信号扫描时间中运用脉冲序列。该实际信号扫描时间大于固有信号扫描时间乘以缩减系数。欠采样允许磁共振信号有较小的采集速率和磁梯度脉冲较小的转换速率和振幅,重调焦脉冲的较低的峰值RF-场。这样就能实现降低声学噪声和特定吸收率。

Description

能降低声学噪声的磁共振成像方法
本发明涉及到一种磁共振成像方法,该方法包括运用含有一或多个脉冲的脉冲序列.
在美国专利US4680545中公开了这样一种磁共振成像方法.
磁共振成像方法一般要涉及包括射频(RF)脉冲和磁梯度脉冲的脉冲序列.RF-脉冲例如用作励磁脉冲和重调焦脉冲.磁梯度脉冲被特别用于磁共振信号的空间编码.RF-脉冲和磁梯度脉冲都可以用来操纵磁共振信号的相位.
磁梯度脉冲的导通和关断会使用来施加磁梯度的梯度线圈产生声学振动.如果不采取措施,梯度线圈的声学振动就会造成诸如噪声、喀哒声和嗡嗡声等恼人的声响.
公知的磁共振成像方法采用降低转换速率并具有正弦形状脉冲边沿的磁梯度脉冲.尽管公知的磁共振成像方法能降低梯度脉冲造成的声学噪声,但降低转换速率和正弦形状的脉冲边沿会延长扫描时间,也就是延长了为获得用来重组磁共振图像所需的所有磁共振信号所需要的时间.因此,公知的磁共振成像方法不适合按高速率产生连续的磁共振图像.
本发明的目的是提供一种不会产生恼人的高音噪声并且能缩短扫描时间的磁共振成像方法.
用来实现本发明这一目的的一种磁共振成像方法包括利用具有空间灵敏度剖面图的接收器天线系统来获得一系列磁共振信号,其中:
-应用欠采样信号采集,以便在k-空间内用预定的降低了的采样速率来采集欠采样磁共振信号,其采样速率相对于全采样速率降低一个缩减系数,
-在实际信号扫描时间中应用脉冲序列,该实际信号扫描时间大于固有信号扫描时间乘以该缩减系数.
为了进行比较,引入了脉冲的参考时间脉冲形状.磁梯度脉冲的参考时间脉冲形状可以采用常用的时间脉冲形状(例如是梯形时间脉冲形状),但参考时间脉冲形状也可以采用任何预定的时间脉冲形状.对于RF-脉冲,特别是RF-励磁脉冲和重调焦脉冲,参考脉冲形状可以用其时间脉冲持续时间、峰值RF-场振幅、和所讨论的该RF-脉冲发生自旋翻转处的翻转角来代表.应该注意到翻转角取决于时间脉冲持续时间和峰值RF-振幅.
此外,对于脉冲序列中各个别的磁梯度脉冲和/或RF-脉冲,参考时间脉冲形状还可以包含不同或相同的脉冲形状.
按照本发明,磁共振信号是用k-空间内的采样速率获得的,它比预定的全‘视场(field of view)’所需的速率要低.降低采样速率的结果是,这些欠采样磁共振信号仅仅具有局部空间编码,因为各个单独的欠采样磁共振信号包括预定的全‘观察场’内不同空间位置的影响.根据接收器天线系统的空间灵敏度剖面图由欠采样磁共振信号计算出有关分散的空间位置的信号影响.换句话说,接收器天线系统的空间灵敏度剖面图提供了欠采样磁共振信号中没有的附加空间编码.由欠采样磁共振信号配合使用接收器天线系统的空间灵敏度剖面图来重组磁共振图像.
欠采样允许获得的欠采样磁共振信号比全采样所需的信号要少.按照本发明的磁共振成像方法,与为了获得同样数量的磁共振信号而在k-空间中降低采样速率相比,脉冲序列的扫描时间缩短较少或是根本不缩短.因此,本发明能够按较低的时间速率获得磁共振信号,因而脉冲序列中磁梯度脉冲的时间变化比参考时间脉冲形状的时间变化更慢.为了执行公平的比较,为脉冲序列所规定的固有信号扫描时间仅仅不同于为磁梯度脉冲和RF-脉冲的时间参考脉冲形状实际采用的脉冲序列.也就是说,固有信号扫描时间是为和实际脉冲序列具有相同的磁梯度脉冲和RF-脉冲顺序的脉冲序列而计算的,但是其中的磁梯度脉冲具有参考时间脉冲的形状.因而较少产生声学噪声,而扫描时间却等于或是仍然小于如果是对磁梯度脉冲采用参考脉冲形状时所具有的固有扫描时间.特别是,声学噪声比实际信号扫描时间缩短到小于k-空间中所允许的采样速率时更小.特别是,本发明能够针对不同环境在缩短实际信号扫描时间和降低声学噪声之间获得-种最佳折衷.例如,若是采用在1-3ms内上升到10mTm-1的梯度斜率,与参考脉冲形状在0.2ms内的21mTm-1的梯度斜率相比,声学噪声能降低20-30dB.
此外,欠采样还允许使用降低的梯度值,由于信号扫描时间相对较长可使用较低的接收器带宽.这样能提高磁共振信号的信噪比.
欠采样程度是由缩减系数给定的,它被定义为实际采样速率或k-空间内的密度相对于全采样速率的比例.因此,如果例如实际采样仅仅是全采样所需磁共振信号的一半,则缩减系数就是1/2.缩减系数通常是间隔[0,1]中的一个有理数.
除了降低声学噪声之外,还能降低脉冲序列在病人组织中积聚的电磁能量.电磁能量的积聚通常是用‘特定吸收率’表示的.特定吸收率(或SAR)被定义为每单位质量吸收的射频功率.具体说,特定吸收率因降低了磁梯度脉冲的转换速率(slew rate)而降低.较低的RF-脉冲的峰值RF-场也会降低特定吸收率.由于扫描时间缩短的程度比磁共振信号欠采样所允许的程度要低,所以RF-脉冲的持续时间可以增加,从而在采用低峰值RF-脉冲的同时获得理想的翻转角度,而不需要增加相对于参考脉冲序列的扫描时间.按照本发明的磁共振成像方法获得的特定吸收率的下降不会无意中加热病人的局部体温.因此,由于病人体温升高所带来的生理上的不适甚至有害影响被减轻乃至避免.降低特定吸收率配合着采用比较短扫描时间的可能性特别有利于采用更强的静磁场,例如是磁场强度达到1.5T以上,甚至达到3T的磁场.另外,降低磁梯度脉冲的转换速率和/或梯度值能够降低对受检病人的外围神经刺激.
以下要参照在权利要求书中限定的实施例具体描述本发明的各个方面.
声学噪声和/或特定吸收率的具体而有效的降低是在实际信号扫描时间等于或大约等于固有信号扫描时间时获得的.在本实施例中,所获得的磁共振信号数量的减少被完全用于降低信号采集速率,从而减少磁梯度脉冲的时间变化和/或减少重调焦RF-脉冲的数量和峰值RF-场.然而,实际信号扫描时间仍然不超过固有信号扫描时间.
具体说,采用下述磁梯度脉冲能够有效地降低噪声,这种脉冲的时间变化比参考时间变化要慢.尽管用来施加磁梯度脉冲所需的时间延长了,但获得磁共振信号所需的时间却不会增加,因为欠采样所获得的磁共振信号比较少.
例如,通过采用其速率比参考时间脉冲形状中的转换速率要慢的磁梯度脉冲的转换速率可实现噪声降低.由于转换速率是声学噪声的一个主要来源,降低转换速率可导致明显降低声学噪声.
为单个磁共振图像获得的磁共振(MR)信号的数量因对MR-信号采用欠采样而减少了.这种欠采样包括在k-空间内减少采样点的数量,这可以通过各种途径来实现.已经有许多包含欠采样的磁共振成像方法.这些磁共振成像方法一般属于‘平行成像方法’,因为欠采样的磁共振信号同时关系到k-空间内的若干条线,实际上是平行地获得k-空间内两条以上的线.MR信号是通过属于若干接收器天线的信号通道获取的,最好是表面线圈等接收器线圈.通过若干条信号通道能够平行地获得信号,从而进一步缩短信号扫描时间.
由于欠采样,采样的数据包含成像物体内从若干位置来的信息.利用有关信号通道的灵敏度剖面图来重组从欠采样MR-信号获得的磁共振图像.特别是,灵敏度剖面图是接收器天线例如是接收器线圈的空间灵敏度剖面图.接收器天线最好是采用表面线圈.重组的磁共振图像可以认为是由与各个波长上的亮度/对比度变化有关的大量空间谐波成分组成的.磁共振图像的分辨率是由最小波长也就是最高波数(k-值)确定的.最大波长即最小波数就是磁共振图像的观察场(FOV).分辨率是由观察场和采样数量确定的.
可以这样来获得欠采样,用各相应的接收器天线获得MR信号,使它们在k-空间内的分辨率比起磁共振图像所需的分辨率更加粗糙.采样的最小波数也就是在k-空间内的最小步幅被增大,同时采样的最大波数则保持不变.因此,在采用欠采样时图象分辨率保持不变,而与此同时最小k-空间步幅增加,也就是减小了FOV.降低k-空间内的采样密度就能实现欠采样,例如是在扫描k-空间时跳过一些线,使得k-空间内对线的扫描间隔大于磁共振图像的分辨率所需的间隔.在缩减观察场的同时维持最大k-值,使得采样点数量相应减少,就能实现欠采样.由于缩减了观察场,采样的数据中包含成像物体内的若干个位置的信息.
特别是,若是从各个接收器线圈的欠采样MR信号来重组接收器线圈图像,这种接收器图像会包含了因缩减观察场造成的人为混淆现象.根据接收器线圈图像和灵敏度剖面图来查清接收器线圈图像中个别位置的影响,并且重组成磁共振图像。这个MR成象法是按缩写而称作SENSE法,并是已知的.在国际专利申请WO99/54746-A1号中具体描述了这一SENSE方法.
也可以将欠采样的MR-信号组合成组合MR-信号,这些信号提供对应于全观察场的k-空间的采样.特别是,按照所谓SMASH-方法,欠采样MR-信号近似于按照灵敏度剖面图合成的低阶球形谐波.参阅国际申请WO98/21600号可以了解这种SMASH-方法.
还可以在空间上执行欠采样.在这种情况下,MR-信号的空间分辨率要小于磁共振图像的分辨率,并且根据灵敏度剖面图形成对应着磁共振图像的MR-信号.空间欠采样具体是这样获得的,由例如来自各个接收器线圈的各独立的信号通道的MR-信号形成来自物体中若干部分的组合影响.这些部位例如是同时受激励的切片(slice).各个信号通道中的MR-信号通常形成来自各个部位例如是切片的线性的组合影响.这种线性组合包括有关信号通道即接收器线圈的灵敏度剖面图.因此,各相应的信号通道的MR信号和各相应部位(切片)的MR信号是通过一个灵敏度矩阵而相关联的,这个矩阵代表物体若干部位的影响在各相应的信号通道中因灵敏度剖面图而产生的权重.这一灵敏度矩阵的反演可导出属于物体各部位的MR-信号.特别是,这样就能导出来自各个切片的MR-信号,并且重组成这些切片的磁共振图像.
在按照本发明的磁共振成像方法的脉冲序列中,磁梯度脉冲的脉冲形状有利地可以调节.因此可以实现缩短信号扫描时间,而代价是声学噪声电平和/或特定吸收率较高.在某些检查环境中,较短的信号扫描时间比降低声学噪声和/或特定吸收率更有价值.具体例子是:涉及病人屏住呼吸的扫描的磁共振成像方法;以在病人大脑内灌注成像为目标的磁共振成像方法;以及对受检病人在时间上的生理过程所引起的广泛的磁共振成像方法.可调节的实际信号扫描时间为本发明的磁共振成像方法提供了比常规磁共振成像方法更大的灵活性.
以下面所述的实施例并参照附图来解释本发明的上述及其他方面,其中
图1表示采用本发明的一种磁共振成像系统的示意图;以及
图2到5表示实际信号扫描时间对固有扫描时间对照的简单示例.
图1表示采用本发明的一种磁共振成像系统的示意图.磁共振成像系统包括一组主线圈10,用来产生稳定均匀的磁场.主线圈例如是这样构成的,使它们包围一个隧道形检查空间.将受检病人推入这一隧道形检查空间.磁共振成像系统还包括许多梯度线圈11,12,由它们产生主要以各个方向上时间梯度的形式而呈现为空间上变化的磁场以叠加在均匀的磁场上.梯度线圈11,12连接到一个可控供电单元21.利用供电单元21提供一个电流来激励梯度线圈11,12.通过控制供电单元来控制梯度的强度,方向和持续时间.磁共振成像系统还包括分别用来产生RF激励脉冲和用来拾取磁共振信号的发送和接收线圈13,16.发送线圈13最好是构成一个用来包围受检物体(的一部分)的体线圈13.体线圈在磁共振成像系统中通常是这样布置的,在他或她被送入磁共振成像系统时,用体线圈13包围受检的病人30.体线圈13作为发射RF激励脉冲和RF重调焦脉冲的发射天线.体线圈13最好能提供空间上均匀分布的发射RF脉冲(RFS).同一线圈或天线通常被交替用作发射天线和接收天线.体线圈最好是采用所谓的协同线圈.发射和接收线圈通常采用线圈的形状,但是能够作为RF电磁信号的发射和接收天线的其他形状也是可以的.发射和接收线圈13被连接到一个电子发射和接收电路15.
应该注意到最好是采用独立的接收线圈16,具体说接收线圈可以采用表面线圈16.这种表面线圈在比较小的范围内具有高灵敏度.表面线圈的空间灵敏度剖面图要相对于体线圈的均匀灵敏度剖面图来校准.
发射线圈例如是表面线圈被连接到解调器24,利用解调器24对接收的磁共振信号(MS)解调.解调后的磁共振信号(DMS)被提供给一个重组单元.接收线圈被连接到一个前置放大器23.前置放大器23放大由接收线圈16接收的RF磁共振信号(MS)并将放大的RF磁共振信号提供给解调器24.解调器24把放大的RF磁共振信号解调.解调的磁共振信号包含有关被成像物体组织内局部自旋密度的实际信息.此外,发射和接收电路15被连接到调制器22.调制器22及发射和接收电路15激励发射线圈13发射RF激励和重调焦脉冲.重组单元从解调的磁共振信号(DMS)导出一或多个图像信号,图像信号代表受检物体成像部分的图像信息.实际的重组单元25是用数字图像处理单元25构成的,它被编程以便从解调的磁共振信号导出代表被成像物体该部分图像信息的图像信号.信号输出到重组监视器26,以便监视器能显示出磁共振图像.也可以将来自重组单元25的信号存储在缓冲器单元27中等待进一步处理.
按照本发明的磁共振成像系统还设有一个控制单元20,例如是一台包括(微)处理器的计算机.控制单元20控制RF激励的执行并且施加时间的梯度场.具体说,按照本发明,控制单元被用来调节磁场梯度脉冲的脉冲形状,例如读梯度脉冲和相位编码梯度脉冲.控制单元用来从用户接收设置脉冲序列的指令,以便能够降低声学噪声.控制单元按照这一指令设置脉冲序列,以使对k-空间采用欠采样扫描,并且可以例如对磁共振图像的重组采用SENSE技术.此外,还能调节脉冲形状,使得脉冲序列的信号扫描时间大致等于固有信号扫描时间.为此,可以将按照本发明的计算机程序加载到例如控制单元20和重组单元25中.
图2到5表示实际信号扫描时间与固有扫描时间对照的简单例子.在图2到5中用图形表示的脉冲序列是TSE序列.图2表示本例中的参考TSE脉冲序列.图2中用曲线表示的‘Gread’代表读方向上的读梯度脉冲,同时还表示了磁共振信号.曲线‘Gslice’表示切片磁选择梯度脉冲.Genc表示磁场相位编码梯度脉冲.曲线RF表示RF-脉冲,例如是射频激励和重调焦脉冲.脉冲序列从一个激励RF-脉冲(RFer,α1)开始,例如是激励在受检物体中磁自旋的90°脉冲.RF-激励脉冲(RFer)通过一个和RF-激励脉冲(RFer)同时施加的磁切片选择梯度脉冲(Gs1r)起着有选择地切片的作用.接着,施加连续的RF重调焦脉冲(Rf1r,Rf2r,…Rf4r)(β1)以便对横向磁化分量重调焦,以产生自旋-回波磁共振信号.由于这些RF重调焦脉冲产生了连续的自旋回波磁共振信号(MR1r,MR2r,MR3r,MR4r).为了获得空间编码,如图2中曲线‘Genc’所示施加相位编码磁梯度脉冲(PE1r,PE2r,…PE8r)和读梯度脉冲(RG1r,RG2r,RG3r,RG4r和RG5r).读梯度脉冲在kx-方向上实现在k-空间中的扫描.连续相位编码脉冲(PE1r,PE2r,…PE8r)提供kx-方向上相应的位移.读梯度脉冲和相位编码梯度脉冲的各种脉冲形状的时间形状具有参考时间脉冲的形状.尤其是,参考时间脉冲形状具有比较高的转换速率,以允许相对较短的固有扫描时间(TR0).
图2表示按照固定重复时间等于TR的典型的四ETL TSE序列的一种梯度和RF-波形的示意图.因此总的信号扫描时间是Tscan=N×TR,其中的N是获得一个完整图像矩阵所需的交错次数.对于256编码的图像矩阵,N=64.参考脉冲序列的这一信号扫描时间Tscan是这一ETLTSE序列的固有信号扫描时间.本例中的梯度转换速率被设在最大.
SAR ∝ Σ T scan α 1 ( B 1 max , t ) 2 + Σ T scan β 1 ( B 1 max , t ) 2 T scan - - - [ 1 ]
典型的β1=2*α1,其中β1是一个180度RF脉冲,它是脉冲形状,持续时间t的函数,而峰值RF场振幅等于某一最大有效值BImaxμT.在本例中,对于全编码的扫描,α1脉冲的总数是64,而β1脉冲的总数是4×64=256.对于固定的TR,声学噪声与梯度值,转换速率和梯度斜率的总数成正比.假设所有图都是按相同比例绘制的,则在图2中,梯度值代表TSE序列的参数所需的初始幅值.梯度上升时间是时间表示的最大转换速率.一般表示为:
Acoustic_Noise_Level ∝G·S·Nslopes    [2]
其中G=梯度幅值,S=梯度转换速率,而Nslopes=梯度斜率数.降低G,S,和Nslopes之一或是同时降低三者就能明显降低声学噪声.
图3中表明为‘Gread’的曲线表示在采用欠采样的一例实际脉冲序列中读梯度脉冲的脉冲形状.曲线Galice代表切片磁选择梯度脉冲.曲线GENC代表磁相位编码梯度脉冲.曲线RF代表RF-脉冲,例如是射频激励和重调焦脉冲.图3中所示的脉冲序列相对于图2中具有固有信号扫描时间和参考脉冲形状的脉冲序列而言具有较少的梯度脉冲,它们有较低的幅值和转换速率.用来对横向磁化分量重调焦的RF重调焦脉冲(Rf1a,…Rf2a)数量较少,大约是图2中所示参考TSE脉冲序列中RF重调焦脉冲数量的一半.图3中的读梯度脉冲(RG1a,RG2a,RG3a)和磁切片选择脉冲(Gs1a)相对于图2中所示的参考脉冲形状有小得多的转换速率.此外,产生的自旋回波磁共振信号(MR1a,MR2a)比较少.相位编码磁梯度脉冲(PE1a,…PE4a)也具有低转换速率和振幅.重复时间TR和图2中的参考序列一样.由于采集的磁共振信号的数量较少,在本例中仅有一半,实际信号扫描时间等于固有信号扫描时间(Tscan).
图3类似于图2,不同的只是磁读梯度脉冲和磁切片选择梯度脉冲的梯度转换速率都被降低.这样的效果是降低声学噪声并降低SAR.由于梯度转换速率被降低了,有必要进一步加大梯度波形的间距.由于重复时间TR是固定的并有一个有效的空余时间,因此这时有可能不影响总的信号扫描时间,仍然维持如图2中的Tscan.在本例中,SENSE被用来完成对全矩阵的采集.降低转换速率的进一步优点是具有降低磁(梯度)场的时间导数dB/dt的效果,它在采用强梯度的快速成像序列中有可能造成外围精神刺激(PNS).
图4中用‘Gread’表示的曲线在采用欠采样的一个实际脉冲序列的例子中表明读梯度脉冲的脉冲形状.曲线Gslice代表切片磁选择梯度脉冲.曲线GENC代表磁相位编码梯度脉冲.曲线RF代表RF-脉冲,例如是射频激励和重调焦脉冲.图4中所示的脉冲序列相对于图2中具有固有信号扫描时间和参考脉冲形状的脉冲序列而言具有较少的梯度脉冲,这些梯度脉冲有较低的幅值和转换速率.图4中读梯度脉冲(RG1a,RG2b,RG3b)和磁切片选择脉冲(Gs1a)的转换速率相对于图2中所示的参考脉冲形状要小得多.另外,读梯度脉冲(RG2b,RG3b)的振幅被降低并且时间被延长.实际信号扫描时间(Tscan)等于固有信号扫描时间(Tscan).
图4表示的序列与图3类似,不同的只是‘Gread’梯度波形的振幅被降低并且时间被延长.TR中可利用的多余时间被用于容纳所增加的回波间隔.降低和延长‘Gread’波形的净效果有三方面.
1.降低梯度振幅能进一步降低声学噪声.
2.降低梯度振幅能进一步降低PNS的可能性.
3.‘Gread’波形的延长通过采用在多数情况下能获得更高的SNR的较低的采集带宽而得到补偿.
再次利用SENSE技术以便从实际采集用来重组磁共振图像的64×2=128编码步幅而恢复出全256矩阵.
SAR与图3中相同(低于采用参考序列的扫描),声学噪声和PNS可能性甚至比图3更低,SNR由于对各个回波采样所用的低带宽而得到改善.同样,总扫描时间Ts can和图像分辨率/特性基本上与图2中使用参考脉冲序列的扫描时的相同.
图5中用‘Gread’表示的曲线在采用欠采样的一个实际脉冲序列例子中表明读梯度脉冲的脉冲形状.曲线Gslice代表切片磁选择梯度脉冲.曲线GENC代表磁相位编码梯度脉冲.曲线RF代表RF-脉冲,例如是射频激励和重调焦脉冲.图5中所示的脉冲序列相对于图2的具有固有信号扫描时间和参考脉冲形状的脉冲序列而言,具有较少的梯度脉冲,这些脉冲具有较低的幅值和转换速率.图5中读梯度脉冲(RG1a,RG2b,RG3b)和磁切片选择脉冲(Gs1b)的转换速率相对于图2中所示参考脉冲形状的转换速率要小得多.另外,RF-激励脉冲(Rfeb)和重调焦RF-脉冲(Rf1b,RF2b)的峰值RF-场与参考脉冲序列相比被降低.同样,与图4的例子一样,读梯度脉冲(RG2b,RG3b)的振幅被降低并且时间被延长.实际信号扫描时间大约等于固有信号扫描时间.
图5类似于图4.在这种情况下所采用的RF脉冲使用较低的B1(用α2和β2表示).
使用低B1进一步降低序列中的SAR.另外还需要缩小RF脉冲带宽,这将导致RF脉冲持续时间延长.缩小带宽的结果是‘Gslice’梯度的振幅也会降低.这种振幅降低导致声学噪声和PNS可能性进一步降低.
由于RF脉冲的持续时间被延长,它们的间隔更大.利用TR内留下的剩余时间有可能使RF脉冲进一步分离,并且额外延长和降低‘Gread’梯度波形.这种从属的效果能够进一步同时降低声学噪声电平和进一步提高SNR.
按照前述所有的实施例,总采集时间保持不变,并且利用SENSE方法产生全256矩阵的数据.
本实施例中的PNS可能性比前面的实施例更低.

Claims (6)

1. 一种磁共振成像方法包括
-应用包括一或多个脉冲的脉冲序列,
-该脉冲序列具有对预定的全‘观察场’的在k-空间内基于全采样速率的固有扫描时间并且具有磁梯度脉冲的参考时间脉冲的形状,
-利用具有空间灵敏度剖面图的接收器天线系统采集一系列磁共振信号,其中,
-应用欠采样的信号采集,以便在k-空间内用预定的降低了的采样速率来采集欠采样的磁共振信号,该采样速率相对于全采样速率降低一个缩减系数,
-在实际信号扫描时间中应用该脉冲序列,该实际信号扫描时间大于固有信号扫描时间乘以缩减系数。
2. 按照权利要求1的磁共振成像方法,其中实际信号扫描时间基本上等于固有信号扫描时间。
3. 按照权利要求1的磁共振成像方法,其中所应用的脉冲的脉冲形状所具有的时间上的变化比参考时间脉冲形状所涉及的时间上的变化要慢。
4. 按照权利要求1的磁共振成像方法,其中脉冲序列包括磁梯度脉冲,并且以一个小于参考脉冲形状所涉及的转换速率的转换速率来应用该磁梯度脉冲。
5. 按照权利要求4的磁共振成像方法,其中脉冲序列包括射频(RF)脉冲,并且以一个小于参考脉冲形状所涉及的RF-磁场分量(B1)的RF-磁场分量(B1)来应用该射频(RF)脉冲。
6. 按照权利要求1的磁共振成像方法,其中该脉冲的脉冲形状是可以调节的,以使得在固有信号扫描时间乘以缩减系数与固有信号扫描时间之间的范围内适配脉冲序列的实际信号扫描时间。
CNB028232585A 2001-11-26 2002-11-20 能降低声学噪声的磁共振成像方法 Expired - Fee Related CN100414316C (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP01204540 2001-11-26
EP01204540.7 2001-11-26

Publications (2)

Publication Number Publication Date
CN1589410A CN1589410A (zh) 2005-03-02
CN100414316C true CN100414316C (zh) 2008-08-27

Family

ID=8181301

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB028232585A Expired - Fee Related CN100414316C (zh) 2001-11-26 2002-11-20 能降低声学噪声的磁共振成像方法

Country Status (6)

Country Link
US (1) US7061240B2 (zh)
EP (1) EP1454156A1 (zh)
JP (1) JP4360912B2 (zh)
CN (1) CN100414316C (zh)
AU (1) AU2002348830A1 (zh)
WO (1) WO2003046597A1 (zh)

Families Citing this family (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1506421A1 (en) * 2002-05-13 2005-02-16 Koninklijke Philips Electronics N.V. Inherent limitation of the reduction factor in parallel imaging as a function of field strength
US6989673B2 (en) * 2003-11-26 2006-01-24 General Electric Company Method and apparatus to reduce RF power deposition during MR data acquisition
US7463920B2 (en) * 2005-03-25 2008-12-09 Siemens Medical Solutions Usa, Inc. Rapid MR visualization of interventional devices
WO2008134891A1 (en) * 2007-05-03 2008-11-13 National Research Counsil Of Canada Rf based spatially selective excitation in mri
CN101308202B (zh) * 2007-05-17 2011-04-06 西门子公司 并行采集图像重建的方法和装置
WO2009111027A2 (en) * 2008-03-04 2009-09-11 Oni Medical Systems, Inc. Optimization of rf transmit and gradient magnetic field imaging using radio frequency and gradient coils
WO2010014686A1 (en) 2008-07-30 2010-02-04 Aleva Neurotherapeutics, S.A. Apparatus and method for optimized stimulation of a neurological target
EP3231476B1 (en) 2008-11-12 2019-06-19 Ecole Polytechnique Fédérale de Lausanne Microfabricated neurostimulation device
EP2414857B1 (en) * 2009-04-01 2019-12-04 Koninklijke Philips N.V. Sar hotspot reduction by temporal averaging in parallel transmission mri
US8774937B2 (en) 2009-12-01 2014-07-08 Ecole Polytechnique Federale De Lausanne Microfabricated surface neurostimulation device and methods of making and using the same
AU2011234422B2 (en) 2010-04-01 2015-11-05 Ecole Polytechnique Federale De Lausanne (Epfl) Device for interacting with neurological tissue and methods of making and using the same
WO2012143821A1 (en) * 2011-04-21 2012-10-26 Koninklijke Philips Electronics N.V. Magnetic resonance imaging of object in motion
DE102012203512B4 (de) 2012-02-03 2014-02-13 Siemens Aktiengesellschaft Ermittlung einer MR-Messsequenz mittels eines Gradienten-Optimierungsverfahrens
KR101967243B1 (ko) * 2012-12-28 2019-04-09 삼성전자주식회사 고속 자기 공명 영상 방법 및 장치
WO2014111782A1 (en) * 2013-01-15 2014-07-24 Koninklijke Philips N.V. Arterial spin labeling with reduced acoustic noise emissions and method of operation thereof
US9557397B2 (en) * 2013-11-04 2017-01-31 Aspect Imaging Ltd. Method for manipulating the MRI's protocol of pulse-sequences
JP6419730B2 (ja) * 2014-01-27 2018-11-07 株式会社日立製作所 磁気共鳴イメージング装置及び騒音低減方法
WO2015173787A1 (en) 2014-05-16 2015-11-19 Aleva Neurotherapeutics Sa Device for interacting with neurological tissue and methods of making and using the same
US11311718B2 (en) 2014-05-16 2022-04-26 Aleva Neurotherapeutics Sa Device for interacting with neurological tissue and methods of making and using the same
US9474894B2 (en) 2014-08-27 2016-10-25 Aleva Neurotherapeutics Deep brain stimulation lead
US9403011B2 (en) 2014-08-27 2016-08-02 Aleva Neurotherapeutics Leadless neurostimulator
DE102015223658B4 (de) * 2015-11-30 2017-08-17 Siemens Healthcare Gmbh Verfahren zum Erfassen von Magnetresonanz-Signalen eines Untersuchungsobjekts
EP3411111A1 (en) 2016-02-02 2018-12-12 Aleva Neurotherapeutics SA Treatment of autoimmune diseases with deep brain stimulation
CN105738843B (zh) * 2016-02-19 2018-10-09 上海联影医疗科技有限公司 Pns水平检测及其有效梯度刺激时间计算方法
JP6944816B2 (ja) * 2017-06-07 2021-10-06 株式会社日立製作所 磁気共鳴イメージング装置
US10702692B2 (en) 2018-03-02 2020-07-07 Aleva Neurotherapeutics Neurostimulation device
CN109917315B (zh) * 2019-04-30 2021-09-28 上海联影医疗科技股份有限公司 磁共振成像扫描方法、装置、计算机设备和存储介质
CN112858975B (zh) * 2019-11-26 2023-05-30 上海联影医疗科技股份有限公司 用于磁共振检查的噪声控制方法、及扫描序列确定方法
US11340323B2 (en) * 2020-01-06 2022-05-24 General Electric Company Low acoustic noise magnetic resonance image acquisition
CN112068054A (zh) 2020-08-19 2020-12-11 上海东软医疗科技有限公司 数据监控方法、装置、设备及磁共振系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1120673A (zh) * 1994-05-31 1996-04-17 株式会社岛津制作所 磁共振成像方法和装置
CN1263601A (zh) * 1998-04-17 2000-08-16 皇家菲利浦电子有限公司 磁共振成象方法和装置
US6198282B1 (en) * 1999-10-07 2001-03-06 General Electric Company Optimized MRI gradient system for providing minimum-duration gradient pulses
US20010010464A1 (en) * 2000-01-21 2001-08-02 Hiromitsu Takamori Magnetic resonance imaging apparatus

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4680545A (en) * 1985-01-04 1987-07-14 General Electric Company Method for reduction of acoustical noise generated by magnetic field gradient pulses
FR2620827B1 (fr) * 1987-09-21 1990-01-26 Crouzet Sa Procede de mesure de la frequence d'un signal periodique et frequence metre pour la mise en oeuvre du procede
DE3811066A1 (de) * 1988-03-31 1989-10-12 Philips Patentverwaltung Kernresonanz-spektrometer
US5910728A (en) 1996-11-12 1999-06-08 Beth Israel Deaconess Medical Center Simultaneous acquisition of spatial harmonics (SMASH): ultra-fast imaging with radiofrequency coil arrays
US7003343B2 (en) * 2001-03-30 2006-02-21 Cornell Research Foundation, Inc. Method and apparatus for anatomically tailored k-space sampling and recessed elliptical view ordering for bolus-enhanced 3D MR angiography
JP2005508238A (ja) * 2001-11-07 2005-03-31 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 高速ダイナミック画像を形成する磁気共鳴方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1120673A (zh) * 1994-05-31 1996-04-17 株式会社岛津制作所 磁共振成像方法和装置
CN1263601A (zh) * 1998-04-17 2000-08-16 皇家菲利浦电子有限公司 磁共振成象方法和装置
US6198282B1 (en) * 1999-10-07 2001-03-06 General Electric Company Optimized MRI gradient system for providing minimum-duration gradient pulses
US20010010464A1 (en) * 2000-01-21 2001-08-02 Hiromitsu Takamori Magnetic resonance imaging apparatus

Also Published As

Publication number Publication date
US7061240B2 (en) 2006-06-13
EP1454156A1 (en) 2004-09-08
CN1589410A (zh) 2005-03-02
AU2002348830A1 (en) 2003-06-10
WO2003046597A1 (en) 2003-06-05
JP2005510320A (ja) 2005-04-21
JP4360912B2 (ja) 2009-11-11
US20040245987A1 (en) 2004-12-09

Similar Documents

Publication Publication Date Title
CN100414316C (zh) 能降低声学噪声的磁共振成像方法
US6897655B2 (en) Moving table MRI with frequency-encoding in the z-direction
CN102576061B (zh) 使用多通道rf激励进行mr成像
DE4436801A1 (de) MRI-Vorrichtung mit dynamischer Empfängerverstärkung
GB2230866A (en) Computed nmr image with enhanced signal to noise ratio
US6448770B1 (en) Gain selection for magnetic resonance imaging and spectroscopy
EP0686855B1 (en) MR imaging apparatus
US8022702B2 (en) MR imaging with an RF pulse producing reduced magnetization transfer
EP3191862B1 (en) Zero echo time mr imaging
KR101310706B1 (ko) 선택적 회질 영상을 획득할 수 있는 자기공명영상 장치 및 이를 이용한 자기공명영상 획득방법
US6169398B1 (en) Magnetic resonance imaging apparatus and imaging method
US5168227A (en) High resolution imaging using short te and tr pulse sequences with asymmetric nmr echo acquisition
US7047062B2 (en) Magnetic resonance imaging with nested gradient pulses
US6255820B1 (en) Variable bandwidth MRI data collection
DE19814677A1 (de) Korrektur einer durch Maxwell-Terme verursachten Verschlechterung eines Axial-Bild-Signals
US7075299B1 (en) Method and apparatus to correct amplitude modulation in multi-echo magnetic resonance imaging
EP0753759A1 (en) Diffusion sensitizing imaging method, dynamic imaging method, and MRI apparatus
Weishaupt et al. Factors affecting the signal-to-noise ratio
JP2000175882A (ja) Mrイメージング装置
US6753683B2 (en) System and method for generating spectra
JP2667706B2 (ja) 磁気共鳴イメージングシステム
CN116523744A (zh) 一种超高场磁共振系统中超高分辨率功能成像方法
CN113544526A (zh) 使用波编码的并行mr成像
CN116635733A (zh) Dixon型水/脂肪分离MR成像
JPH10234708A (ja) Mrイメージング装置

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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20080827

Termination date: 20121120