CN101879063B - 用于扩散加权地拍摄磁共振信号的方法和磁共振设备 - Google Patents
用于扩散加权地拍摄磁共振信号的方法和磁共振设备 Download PDFInfo
- Publication number
- CN101879063B CN101879063B CN201010169936.9A CN201010169936A CN101879063B CN 101879063 B CN101879063 B CN 101879063B CN 201010169936 A CN201010169936 A CN 201010169936A CN 101879063 B CN101879063 B CN 101879063B
- Authority
- CN
- China
- Prior art keywords
- gradient
- momentum
- diffuse
- sequence
- magnetic resonance
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
- G01R33/56341—Diffusion imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
-
- 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/5616—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 gradient refocusing, e.g. EPI
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Vascular Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明涉及一种用于利用拍摄序列(50)扩散加权地拍摄磁共振信号的方法,该拍摄序列包括具有多个扩散编码梯度(56-59)的扩散模块(51)和具有用于拍摄磁共振信号的读出梯度(60)的读出模块(52),其中,所述拍摄序列被配置为用于拍摄相应于预先确定的信号相干路径的磁共振信号,其中所述方法包括利用拍摄序列(50)拍摄磁共振信号,其中,在拍摄序列(50)期间以预先确定的梯度动量(M0-M3)接通扩散编码梯度(56-59)。在此,这样设置所述扩散编码梯度(56-59)的梯度动量(M0-M3),使得相应于与预先确定的相干路径不同的相干路径的磁共振信号被减小,其中,基于阈值进行梯度动量(M0-M3)的设置,以实现预先确定的减小。
Description
技术领域
本发明涉及一种利用包括多个扩散编码梯度的拍摄序列扩散加权地拍摄MR信号的方法,以及一种为此的磁共振设备。
背景技术
扩散加权的磁共振(MR)图像在临床例行程序中提供重要的诊断信息,例如在中风诊断和肿瘤诊断中。在扩散加权的成像(DWI)中使用具有高振幅和长的脉冲持续时间结合合适的读出或图像拍摄模块的扩散-编码梯度。相对于扩散的灵敏度从如下得出:水分子沿着施加的场梯度的扩散减小磁共振(MR)信号。由此,在检查的对象的具有小的扩散的区域发生更小的信号损失,并且相应地该区域被以更高的信号强度成像。在此扩散加权的强度与扩散编码梯度的强度相关。在此扩散加权通过所谓的b值来表征,该b值是梯度参数的函数,例如梯度强度、持续时间或者在梯度之间的间隔的函数。为了避免运动伪影,例如可以使用“单次激发”读出模块,利用其在一个拍摄序列期间拍摄一个完整的原始数据组。例如可以使用平面回波成像序列(EPI)。
在此,用于扩散编码的常用方法使用由Stejskal和Tanner描述的单极自旋回波拍摄序列(参见“Spin Diffusion Measurements:Spin Echoes in the Presence of aTime-Dependent Field Gradient”in The Journal of Chemical Physics 42,1965)。在该方法中在一个自旋回波序列中的180°重聚焦脉冲的每一侧定位两个强对称梯度。该对称扩散编码梯度的目的是,通过去相位的输送加速由于扩散导致的信号损失。在此去相位通常与在其中接通梯度的时间(梯度脉冲长度)的平方、与两个梯度脉冲的时间间隔和与接通的梯度场的强度的平方成比例。
在此,对于所拍摄的图像的质量来说,主要是信噪比(SNR)和几何失真是决定性的。此外使用的梯度系统限制了最大可接通的磁场梯度的强度。对于信噪比重要的成像参数是回波时间TE。几何失真基于梯度基本场振幅B0的空间变化,其中EPI相对该变化是特别敏感的。静态失真由于基本场非线性和检查对象的不同区域的磁化系数而引起。诸如涡流效应的动态失真特别受到梯度脉冲的时间顺序影响。场梯度的每个接通和断开可以引起这样的涡流,该涡流以不同的时间常数衰减。在缓慢衰减的情况下,场分量可以保持直到读出,从而会造成拍摄的MR数据的干扰和失真。特别是在组合以不同的扩散梯度方向和扩散梯度振幅拍摄的图像时值得期望的是,保持动态失真尽可能小,以便减少在产生的数据(例如各向异性图、扩散图、张量数据等)中的错误。
由Stejskal和Tanner描述的单极自旋回波拍摄序列具有强的几何失真(特别是高分量的剩余涡流场),以及梯度信号的高的负担。其它开发的拍摄序列为了避免强的失真使用双极的双自旋回波方案(Doppel-Spin-Echo-Schemata),其能够隐含地补偿涡流场。这样的方案例如在T.G.Reese等人的:“MagneticResonance in Medicine 49:177(2003)”中描述。利用这样的方案能够减小单极方案的缺点,然而代价是更长的回波时间。例如附加要使用的高频(HP)脉冲需要5ms。由于检查的组织的弛豫特征,该更长的回波时间会不利地影响图像质量。此外,通过附加的HF脉冲插入其它不期望的信号相干路径(Signal-),其又会导致其它不期望的自旋回波信号或激励的回波信号,以及导致自由感应衰减信号。该附加的相干路径导致由拍摄的MR数据重建的图像数据中的伪影。已知用于抑制这样的不期望的相干路径的扰相梯度。然而这样的附加的梯度又延长回波时间(TE),并且由此不利地影响信噪比。由此还需要更麻烦的拍摄序列。
值得期望的是,有效地抑制通过附加的信号相干路径导致的图像伪影。在此不应该延长回波时间,以实现良好的信噪比。同样值得期望的是,在最佳利用梯度系统的情况下保持失真小。
发明内容
相应地,本发明要解决的技术问题是,提供一种用于扩散加权地拍摄MR信号的改进的方法,以及一种磁共振设备,利用其能够进行改进的扩散加权的拍摄。
按照本发明的第一方面,提供一种用于利用拍摄序列扩散加权地拍摄MR信号的方法,该拍摄序列包括具有多个扩散编码梯度的扩散模块和具有用于拍摄MR信号的读出梯度的读出模块。拍摄序列被配置为用于拍摄相应于预先确定的信号相干路径的MR信号。该方法包括利用拍摄序列拍摄MR信号,其中在拍摄序列期间以预先确定的梯度动量接通扩散编码梯度。这样设置扩散编码梯度的梯度动量,使得相应于与预先确定的相干路径不同的相干路径的MR信号被减小。在此基于阈值进行梯度动量的设置,以实现预先确定的减小。
扩散模块还可以被称为扩散序列,并且例如除了多个扩散编码梯度还包括用于激励或重聚相位的HF脉冲。借助扩散编码梯度的梯度动量,设置期望的扩散加权和扩散方向。读出模块还可以被称为读出序列并且例如包括以相位编码梯度或频率编码梯度形式的读出梯度。相干路径可以描述由于例如在拍摄序列期间被入射的HF脉冲,不同信号的形成。每个HF脉冲可以产生自由的感应衰减信号,该自由的感应衰减信号在存在其它重聚相位梯度的情况下又会产生可以用这样的相干路径来描述的回波信号或激励的回波信号。
通过设置用于抑制不期望的相干路径的扩散编码梯度,可以避免在由MR信号重建的图像数据中的伪影。因为该不期望的信号分量的抑制不是通过附加的梯度,而是通过设置扩散编码梯度实现的,所以不会延长拍摄序列的持续时间,从而利用该方法可以实现改善的信噪比。通过使用阈值还可以确保,实现不期望的相干路径的预先确定的减少。
在本发明的一种实施方式中,可以根据使用的读出模块来确定所述阈值。由此在使用“单次激发”方案的读出模块中也能实现对不期望的信号分量的有效抑制。
例如,可以基于读出模块的读出梯度的产生的梯度动量(Mreadout)来确定所述阈值。在使用更强的读出梯度的情况下例如还可以提供所述阈值,使得能够确保有效抑制不期望的相干路径。
所述阈值可以定义一个最小的去相位梯度动量,其中这样进行扩散编码梯度的梯度动量的设置,使得由扩散编码梯度产生的梯度动量至少与最小的去相位梯度动量一样大。扩散编码梯度的梯度动量的设置例如可以这样进行,使得对于拍摄序列的其它相干路径(即,不期望的相干路径)的每一个,对于从扩散编码梯度产生的各个相干路径的梯度动量至少与最小的去相位梯度动量一样大。对于不同的相干路径可以产生(wirken)不同的扩散编码梯度,使得相干路径的分开考察和扩散编码梯度的相应设置是具有优势的。
这样确定所述阈值,使得MR信号的、相应于拍摄序列的其它相干路径的k空间的中央位置,位于利用该拍摄序列采样的k空间外部。通过基于所述阈值设置扩散编码梯度,由此可以将不期望的MR信号的中央的k空间位置移动到k空间的区域外部的如下位置,对于该位置利用拍摄序列拍摄MR信号。
所述阈值可以是最小的去相位梯度动量Mspoil并且按照等式Mspoil=N·Mreadout来确定,其中N≥1。在此,读出梯度动量Mreadout可以按照Mreadout=RES/(γFOV)来确定,其中,RES表示分辨率并且FOV表示读出模块的视野,并且γ表示回转磁比(gyromagnetische)。
扩散编码梯度的梯度动量例如可以借助梯度脉冲长度来设置。扩散模块包括至少四个具有至少四个梯度脉冲长度作为参数的扩散编码梯度。这些参数中的三个可以通过三个由拍摄序列确定的条件来确定,并且至少一个剩下的参数可以基于所述阈值来确定。例如所述至少一个剩下的参数可以从对于其它相干路径产生的扩散编码梯度的梯度动量与阈值的比较来确定。对于梯度脉冲长度的设置例如可以考虑以下条件:由拍摄序列预先给出的扩散编码梯度的接通的总持续时间;由扩散编码梯度要达到的磁化的去相位,其相应于预先确定的相干路径;和梯度脉冲持续时间的这样的匹配:使得通过拍摄序列要产生的自旋回波的去相位在读出模块的预先确定的拍摄时刻进行。最后的条件例如可以相应于在双自旋回波(DSE)序列的情况下用于回波时间的条件,其中DSE序列的去相位脉冲之间的时间相应于在激励脉冲和第一去相位脉冲和第二去相位脉冲和拍摄时刻之间的时间的和。
在一个实施方式中,扩散模块包括至少四个扩散梯度,其中分别在相对的方向上接通两个梯度。两个在相对的方向上接通的梯度还可以被称为双极的梯度。利用这样的梯度可以有效抑制涡流,并能够最佳利用MR设备的梯度系统。
读出模块可以包括平面回波成像(EPI)读出序列或分割的平面回波成像(EPI)读出序列。
在另一种实施方式中,扩散模块还包括多个去相位梯度,其产生相应于其它相干路径的MR信号的减小。由此附加的去相位梯度可以支持不期望的信号相干路径的抑制。去相位梯度还可以成对地具有相同的梯度,其中一对的一个梯度可以被设置在扩散模块的重聚焦脉冲之前并且另一个被设置在之后。
去相位梯度的梯度动量和/或去相位梯度的极化可以根据阈值来设置。如果阈值取决于读出模块,则由此可以确保,即使在小的b值的情况下也能够有效抑制不期望的信号路径。
去相位梯度可以与扩散梯度重叠。由此例如在使用附加的去相位梯度的情况下也能避免TE时间的延长,由此可以实现改善的信噪比。
优选地,当扩散编码梯度的振幅(GD)低于预先确定的振幅阈值时,接通去相位梯度。对于具有小的b值的扩散成像例如使用具有小的振幅的扩散编码梯度,从而可能仅不充分地进行了不期望的相干路径的抑制。此时可以这样确定振幅阈值,使得在通过扩散编码梯度不充分抑制的情况下接通附加的去相位梯度,使得即使对于小的b值也能实现有效抑制。
在此,去相位梯度的梯度动量可以根据扩散编码梯度的振幅来设置。
阈值例如可以定义一个最小的去相位梯度动量,其中可以这样进行扩散编码梯度和去相位梯度的梯度动量的设置,使得对于拍摄序列的其它相干路径的每个,由扩散编码梯度和去相位梯度对于相应相干路径产生的梯度动量至少与最小的去相位梯度动量(Mspoil)一样大。由此又可以对于各个相干路径设置如下的条件:基于这些条件既可以设置扩散编码梯度的参数也可以设置去相位梯度的参数。
此外,在至少两个扩散编码梯度之间的扩散模块可以包括具有可变持续时间的暂停,用于提供对于扩散编码梯度的梯度动量的选择的其它自由度。
拍摄序列可以包括具有两个重聚焦脉冲的双自旋回波序列。作为其它相干路径例如可以出现至少三个自由感应衰减,三个自旋回波和两个激励的回波。然后可以这样进行扩散编码梯度的梯度动量的设置,使得相应于这些相干路径的MR信号分量被减小。
前面描述的方法例如可以由控制单元或磁共振设备的计算单元自动地进行。
按照本发明的另一方面,提供一种用于利用如下的拍摄序列扩散加权地拍摄MR信号的磁共振设备,所述拍摄序列包括具有多个扩散编码梯度的扩散模块和具有用于拍摄MR信号的读出梯度的读出模块,其中拍摄序列被配置为用于拍摄相应于预先确定的信号相干路径的MR信号。该磁共振设备包括利用拍摄序列拍摄MR信号并且在拍摄序列期间以预先确定的梯度动量接通扩散编码梯度的拍摄单元,以及被构造为控制拍摄单元的控制单元。控制单元这样设置扩散编码梯度的梯度动量,使得相应于与预先确定的相干路径不同的相干路径的MR信号被减小,其中,在此基于阈值进行梯度动量的设置,以实现预先确定的减小。利用这样的磁共振设备能够实现类似于前面提到的优点。
按照一种实施方式,磁共振设备被构造为执行前面描述的方法之一。
按照本发明的另一方面,提供一种具有计算机程序的计算机程序产品,其在例如与磁共振设备在功能上相连的计算机系统中被执行时,执行前面提到的方法中的一种。这样的计算机程序例如可以在磁共振设备的计算单元或控制单元上被执行。此外提供了一种其上存储了电子可读的控制信息的电子可读的数据载体,其中这样构造所述控制信息,使得在例如与磁共振设备在功能上相连的计算机系统中使用该数据载体时,所述控制信息执行前面提到的方法的一种。
附图说明
以下借助附图详细解释本发明。其中,
图1按照本发明的磁共振设备的一种实施方式的示意图,
图2示出了按照本发明的方法的实施方式的流程图,
图3示出了按照本发明的方法的实施方式的时序图,
图4示出了按照本发明的方法的另一个实施方式的时序图,
图5示出了根据涡流时间常数在k空间中心中剩下的涡流的振幅,
图6示出了利用通常的方法拍摄的磁共振数据与利用按照本发明的方法的一种实施方式拍摄的磁共振数据的比较。
具体实施方式
图1示出了磁共振设备,为了拍摄磁共振数据其被配置为用于执行扩散加权的成像。这样的磁共振设备具有用于产生极化场B0的磁铁10。检查对象、此处是受检人11可以置于卧榻13上被驶入磁铁10中,如示意性通过箭头表示的。MR设备还具有用于产生为了成像和位置编码而被使用的磁场梯度的梯度系统14。为了扩散加权的成像还可以利用梯度系统14产生扩散编码梯度。
为了激励在主磁场中产生的极化,设置高频线圈装置15,其将高频(HF)场入射到被受检人11中,以便将磁化从均衡位置偏转。借助HF线圈装置15例如可以入射诸如90°-Sinc-脉冲的激励脉冲,或去相位脉冲、例如180°脉冲。为了控制磁场梯度设置梯度单元17,并且为了控制入射的HF脉冲设置HF单元16。梯度系统14和高频线圈装置15以及HF单元16和梯度单元17可以一起被称为拍摄单元21。
磁共振设备由控制单元18中央地控制,该控制单元可以包括计算单元,例如用于重建图像数据的计算单元。控制单元18控制HF脉冲的入射、梯度的施加和产生的MR信号的拍摄。操作人员通过输入单元19可以选择序列草案(Sequenzprotokoll)并且输入和改变在显示器20上显示的成像参数。例如可以选择要用来执行扩散加权的成像序列的b值。
在图1中示意性示出的磁共振设备自然可以包括磁共振设备通常具有的其它组件。MR设备的一般工作原理对于专业人员是公知的,从而不必对一般的组件进行详细描述。
此时这样构造控制单元18,使得其借助拍摄单元21执行例如象图3或图4中示出的那样的拍摄序列。该拍摄序列可以包括不同的读出模块,例如平面回波成像读出序列。在接通扩散编码梯度时,控制单元18可以改变振幅和/或梯度脉冲长度,以便设置特定的梯度动量。在一种实施方式中,控制单元18这样设置扩散编码梯度的梯度动量,使得由不期望的相干路径形成的MR信号被抑制。以下将详细参考图2-6描述可以如何进行这样的设置。例如可以考虑双自旋回波拍摄序列,其中仅应该拍摄相应于双自旋回波相干路径的MR信号。其它的相干路径,例如180°脉冲的自由感应衰减(FID)的信号以及各个自旋回波(SE)或激励的自旋回波(STE)应该被抑制。
控制单元18还可以接通特别是能够重叠于扩散编码梯度的附加的去相位梯度或扰相梯度。由此即使在小的b值的情况下也能够有效抑制不期望的相干路径。
图2示出了例如可以在图1中示出的磁共振设备上被执行的按照本发明的方法的一种实施方式。在第一步骤30中对于待使用的拍摄序列确定读出模块。用于产生扩散对比度的不同的拍摄序列对于专业人员是公知的,并且可以利用按照本实施方式的方法使用这些拍摄序列。例如,此处可以使用具有双极的扩散编码梯度的双自旋回波序列。这样的拍摄序列例如在图3中示出。图3在时间轴上示出了拍摄序列的时序,然而应该清楚的是,例如可以在不同的空间方向上接通仅示意性示出的梯度。拍摄序列50包括激励脉冲53以及两个去相位脉冲或重聚焦脉冲54和55,用于产生具有回波时间TE的双自旋回波,其中回波在k空间中心61中出现。拍摄序列50的扩散模块51包括具有梯度脉冲长度ti(i=0-3)和梯度动量Mi=GD·ti的多个扩散编码梯度56-59,其中GD表示梯度56-59的振幅。借助这些梯度在时间段T内进行扩散编码。如可以看出的,使用双极的梯度脉冲,即,分别在相对的方向上接通两个梯度。
拍摄序列50还包括用于采样具有合适的轨迹的k空间的读出模块52。这又通过接通不同的梯度、例如示意性通过梯度60表示的相位和频率编码梯度,来进行。例如可以进行蜿蜒形的平面回波的读出。为此例如可以使用平面回波成像(EPI)或分割的平面回波成像的读出序列。
在扩散编码之前或之后需要其它的时间间隔tprep和tadc,例如对于导航回波、层重聚焦和对k空间的第一部分的MR数据的拍摄。示出的拍摄序列50的扩散编码强度可以通过与编码梯度的振幅的平方成比例的b值来描述。在假定可忽略RF脉冲持续时间的条件下可以根据如下估计b值
b=2/3γ2GD 2(t0+t1)3 (1)
这样的双极的拍摄方案具有如下优点,梯度负载在正极和负极之间分布,由此可以最佳地利用梯度系统的硬件,并且一般地降低涡流。
仍然参考图2在下一步骤31中基于读出模块的读出梯度的梯度动量确定最小的去相位梯度动量(Mspoil)。最小的去相位梯度动量或扰相动量Mspoil表示一个阈值,其中利用超过该阈值的梯度可以进行不期望的信号相干路径的有效抑制。该阈值通过使用的读出模块来确定。对于EPI模块例如需要读出动量Mreadout的至少N倍的扰相动量。特别地使用阈值用于确定扩散编码梯度的梯度动量。这样设置梯度,使得例如不期望的回波在采样的k空间区域外部出现。在此N=1表示k空间的边界,并且在N值更大的情况下不期望的回波进一步通过采样的区域的边界向外移动。对于EPI读出模块产生阈值Mspoil=N·RES/(γFOV),其中RES表示像素的数量并且FOV表示视野。
在步骤32中对于拍摄序列确定不期望的相干路径。这又以图3的拍摄序列50为例来说明。利用在图3中示出的拍摄序列50产生新的相干路径。这是HF脉冲53-55的三个自由感应衰减(FID)、三个自旋回波(SE)、一个激励的回波(STE)、一个反激励的回波(ASTE)和一个双自旋回波(DSE)。预先确定的信号相干路径(按照该信号相干路径要以拍摄序列50拍摄MR信号),是双自旋回波(DSE)。所有其它相干路径会产生在重建的图像数据中的不期望的伪影。按照本实施方式,应该通过合适地设置扩散编码梯度的动量抑制这些其它相干路径。为此至少与预先确定的阈值(即,扰相动量Mspoil)一样大的梯度动量应该作用于该信号路径上。
为此,在步骤33中确定如下条件,这些条件确保,对于每个不期望的相干路径,对于用于扩散编码梯度的相干路径产生的梯度动量至少与最小的去相位梯度动量Mspoil一样大。对于在图3中所举例子可以构造如下等式:
1.FID |M0-M1+M2-M3| ≥Mspoil (5a)
2.FID |-M1+M2-M3| ≥Mspoil (5b)
3.FID |-M3| ≥Mspoil (5c)
1.SE |M0+M1-M2+M3| ≥Mspoil (5d)
2.SE |M0-M1+M2+M3| ≥Mspoil (5e)
3.SE |-M1+M2+M3| ≥Mspoil (5f)
STE |M0+M3| ≥Mspoil (5g)
ASTE |M0-M3| ≥Mspoil (5h)
DSE M0+M1-M2-M3 =0 (5i)
在等式中M0-M3表示扩散编码梯度56-59的梯度动量。以下假定,所有的动量M0-M3都是正的,其中即使在符号反向的情况下也能进行后续计算。在使用等式(5i)的条件下可以确定M0:M0=-M1+M2+M3,并且等式可以简化为
1.FID 2|M1-M2| ≥Mspoil (6a)
2.FID |-M1+M2-M3| ≥Mspoil (6b)
3.FID |M3| ≥Mspoil (6c)
1.SE 2|M3| ≥Mspoil (6d)
2.SE 2|M1-M2-M3| ≥Mspoil (6e)
3.SE |M1-M2-M3| ≥Mspoil (6f)
STE |M1-M2-2M3| ≥Mspoil (6g)
ASTE |M1-M2| ≥Mspoil (6h)
因为2|x|≥|x|总是成立,所以考虑以下不等式就足够了:
|M1-M2| ≥Mspoil (7a)
|M1-M2+M3| ≥Mspoil (7b)
|M3| ≥Mspoil (7c)
|M1-M2-M3| ≥Mspoil (7d)
|M1-M2-2M3| ≥Mspoil (7e)
这一般地可以通过如下来满足:使得设置
M3≥Mspoil+|M1-M2| (8a)
|M1-M2|≥Mspoil (8b)
通过使用特殊解
M3≥2|M1-M2| (9a)
|M1-M2|≥Mspoil (9b)
可以简化后续计算。然而应该注意,该特殊解提出了更高的要求和所需的梯度动量M3。
此时可以区分两种情况,以下用“+”“-”表示。第一种情况是M1-M2≥0(“+”),并且第二种情况是M1-M2<0(“-”)。因为通过GD·ti给出梯度动量Mi,所以可以将条件(9)简单地转换为时间不等式:
“+” “-”
t3≥2(t1-t2) t3≥2(t2-t1) (10a)
t1-t2≥Mspoil/GD t2-t1≥Mspoil/GD (10b)
在此,第二个条件取决于各个扩散梯度振幅GD。该振幅越高,则该不等式越容易满足。利用预先确定的时间段可能难以对于每个GD值找到有效的结果。特别是在使用小的b值的情况下要考虑这点。以下参考图4就此详细讨论。通过使用对于t1-t2最大允许的值可以松动(gelockert)第二条件:
“+” “-”
t3=2(t1-t2) t3=2(t2-t1) (11a)
t1-t2≥Mspoil/GD t2-t1≥Mspoil/GD (11b)
还可以在选择时间t3的情况下为了降低剩余的涡流而使用剩余的自由度。
基于此时确定的条件并基于通过拍摄序列预先给出的边界条件,在图2的步骤34中确定对于扩散编码梯度的梯度动量。在使用恒定的梯度振幅的情况下其相应于确定梯度脉冲持续时间ti。
对于在图3中解释的双自旋回波序列可以构造以下边界条件:
t0+t1=t2+t3 (2)
tprep+t0+t3+tadc=t1+t2 (3)
等式(2)在此确保进行梯度动量的完整的去相位。等式(3)是通过自旋回波的出现的时刻定义的条件。在脉冲53和54和脉冲55和自旋回波61的出现之间的时间段的和应该正好与重聚焦脉冲54和55之间的时间段一样大。此外,根据如下确定对于扩散编码梯度的入射的总持续时间T
T=t0+t1+t2+t3 (4)
这些等式从梯度的相同的梯度振幅GD以及从可忽略的梯度边沿和HF脉冲持续时间出发。回波时间TE和时间段T通过成像草案(Bildgebungsprotokoll)被预先给出,特别是通过期望的扩散灵敏度b。因此利用三个等式和四个参数剩下一个自由度,在按照本发明的实施方式中这样选择该自由度,使得不期望的相干路径被抑制。通过合适地选择该自由度可以隐含地通过扩散编码梯度而无需其它扰相梯度进行所述抑制。该自由的参数通过预先确定的用于抑制不期望的相干路径的条件来确定。
为此首先如下地修改等式(2)-(4):
(4):t0=T-t1-t2-t3 (12a)
(2):t0+t1=t2+t3
(3):TE-T+t0+t3=t1+t2
由此得到
t0=T-TE/2+T/2-t3-T/2+t3-t3=T-TE/2-t3 (13a)
t1=TE/2-T/2+t3 (13b)
t2=T/2-t3 (13c)
这又可以进一步被简化为:
t1-t2=TE/2-T/2+t3-T/2+t3=TE/2-T+2t3
此时将这些通过拍摄序列预先给出的边界条件与通过用于抑制不期望的相干路径的阈值所确定的条件(等式11)组合。得出:
“+” “-”
t3=2(t1-t2)=1/2(T-TE/2+t1-t2) t3=2(t2-t1)=1/2(T-TE/2+t1-t2)
4(t2-t1)=T-TE/2+(t2-t1)
(t2-t1)=1/5(T-TE/2) (15a)
t1-t2≥Mspoil/GD t2-t1≥Mspoil/GD (15b)
此时可以如下地计算梯度脉冲持续时间(利用其可以满足提到的条件):
t0=1/6(2T-TE) t0=1/10(6T-3TE) (16a)
t1=1/6(T+TE) t1=1/10(-T+3TE) (16b)
t2=1/6(-T+2TE) t2=1/10(T+2TE) (16c)
t3=1/6(4T-2TE) t3=1/10(4T-2TE) (16d)
在此要注意以下附加的条件:
t1-t2≥Mspoil/GD t2-t1≥Mspoil/GD
1/5(T-TE/2)≥Mspoil/GD (17)
在此,所有时间都是正的并且成立:TE>T,从而当2T-TE>0时,对于两种情况存在解。因为可以根据T=TE-tprep-tadc来计算T,可以将该条件表达为T>2(tprep+tadc)或T>tprep+tadc。对于需要许多时间的读出模块,例如高分辨率的读出模块,并且在小的b值的情况下,该条件限制了最短可达到的回波时间。然而这是双极的拍摄方法的原则上的限制,其与不期望的相干路径的隐含抑制没有关系,如从图3和等式(2)-(4)可以看出的。
在借助梯度脉冲持续时间计算了梯度动量之后,在步骤(35)中利用拍摄序列拍摄MR信号。在此接通具有特定的梯度动量的扩散编码梯度。通过接通该扩散编码梯度,根据双极性进行涡流的有效抑制。通过以事先确定的梯度脉冲持续时间来接通,在双自旋回波的同时去相位的情况下有效抑制不期望的相干路径。如从图3可以看出的,在此回波时间TE没有延长。除了几何失真的最小化之外,由此还可以实现高的信噪比。然后从利用读出模块52拍摄的磁共振信号中可以进行图像数据的重建。
应该清楚的是,在图2中解释的方法还可以包括通常在扩散加权的成像中使用的其它步骤。例如可以在不接通扩散编码梯度的情况下利用类似的拍摄序列进行MR信号的其它拍摄,以便将相应的图像数据例如通过减影与扩散加权的图像数据进行比较。
如以下参考图4所描述的,对于MR信号的扩散加权的拍摄还可以利用小的b值修改在图3中示出的拍摄序列50。如在图3中示出的拍摄序列50,也可以利用在图1中示出的磁共振设备执行在图4中示出的拍摄序列70。如已经解释的,隐含抑制的效率取决于扩散编码梯度的各个梯度振幅。尽管回波时间TE的放大一般地松动了所述条件,但是这一方面一般是不期望的,并且此外不是对于所有的值GD都可能的(例如在b值为零并且由此GD=0的情况下)。因此在图4中示出的实施方式中使用附加的扰相和去相位梯度。按照该实施方式,在高的扩散梯度振幅的情况下,使用参考图3描述的隐含抑制,并且当振幅GD低于特定的振幅阈值时,使用附加的显式抑制。在此可以随着振幅GD的下降连续地提高显式抑制的强度。在图4的实施方式中直接地在两个重聚焦脉冲54和55之前和之后接通附加的去相位梯度71和72以及73和74。通过使用一对相同的梯度脉冲可以保持双自旋回波相干路径不变。这样的安排还导致附加的梯度对扩散加权的最小的贡献,即,对b值的尽可能小的效力。在此第一以及第二梯度对的梯度动量MA、MB根据扩散编码梯度的振幅GD改变。此外去相位梯度71-74重叠于扩散编码梯度56-59。由此可以实现远远更小的TE时间,这由于通过减少T2弛豫提高了信噪比而提高了图像质量。这是相对于常规方法的主要优点,在常规的方法中附加的扰相梯度总会导致回波时间的延长。
在使用这些附加的梯度的情况下确保,重叠的去相位和扩散梯度的振幅保持在梯度系统的限制内。通过修改等式(5)可以导出对于拍摄序列70的相应条件。在此,特别是通过M0+MA代替M0,通过M1-MA代替M1,通过M2+MB代替M2,通过M3-MB代替M3,这得出以下条件:
1.FID |M0+MA-M1+MA+M2+MB-M3+MB| ≥Mspoil (18a)
2.FID |-M1+MA+M2+MB-M3+MB| ≥Mspoil (18b)
3.FID |-M3+MB| ≥Mspoil (18c)
1.SE |M0+MA+M1-MA-M2-MB+M3-MB| ≥Mspoil (18d)
2.SE |M0+MA-M1+MA+M2+MB+M3-MB| ≥Mspoil (18e)
3.SE |-M1+MA+M2+MB+M3-MB| ≥Mspoil (18f)
STE |M0+MA+M3-MB| ≥Mspoil (18g)
ASTE |M0+MA-M3+MB| ≥Mspoil (18h)
DSE M0+MA+M1-MA-M2-MB-M3+MB=0 (18i)
MA和MB的符号可以是正的也可以是负的。指定等式M0=-M1+M2+M3(18i)得到
1.FID 2|-(M1-M2)-MA+MB| ≥Mspoil (19a)
2.FID |-(M1-M2)-M3+MA+2MB| ≥Mspoil (19b)
3.FID |M3-MB| ≥Mspoil (19c)
1.SE 2|M3-MB| ≥Mspoil (19d)
2.SE 2|-(M1-M2)+M3+MA| ≥Mspoil (19e)
3.SE |-(M1-M2)+M3+MA| ≥Mspoil (19f)
STE |-(M1-M2)+2M3+MA-MB| ≥Mspoil (19g)
ASTE |-(M1-M2)+MA+MB| ≥Mspoil (19h)
又成立2|x|≥|x|,从而可以简化为
|-(M1-M2)-M3+MA+2MB| ≥Mspoil (20a)
|M3-MB| ≥Mspoil (20b)
|-(M1-M2)+M3+MA| ≥Mspoil (20c)
|-(M1-M2)+2M3+MA-MB| ≥Mspoil (20d)
|-(M1-M2)+MA+MB| ≥Mspoil (20e)
除了MA和MB的出现之外,这些不等式与不等式(7)是相同的。这意味着,对于M3≥|M1-M2|和|M1-M2|≥Mspoil,利用MA=MB=0满足这些不等式。
以下又区分“+”和“-”的情况。因为不应该取决于b值改变拍摄序列的时序,以避免不期望的信号改变,所以作出如下选择(比较等式(11)):
“+” “-”
M3=2(M1-M2) M3=2(M2-M1)(21a)
由此得到
|-3(M1-M2)+MA+2MB|≥Mspoil |-(M2-M1)+MA+2MB|≥Mspoil (22a)
|2(M1-M2)-MB|≥Mspoil |2(M2-M1)-MB|≥Mspoil (22b)
|(M1-M2)+MA|≥Mspoil |3(M2-M1)+MA|≥Mspoil (22c)
|3(M1-M2)+MA-MB|≥Mspoil |5(M2-M1)+MA-MB|≥Mspoil (22d)
|-(M1-M2)+MA+MB|≥Mspoil |(M2-M1)+MA+MB|≥Mspoil (22e)
如果
M1-M2<Mspoil M2-M1<Mspoil (23)
则在本实施方式中可以采用显式抑制(参见等式(11))。
在考虑这些等式的条件下并且在考虑具有去相位梯度的连续提高的梯度动量的解的条件下(当出现条件(23)时其为零),可以选择对于MA和MB的以下解:
MA=Mspoil-(M1-M2) MA=3Mspoil-3(M2-M1) (24a)
MB=-2MA MB=-2MA (24b)
此外在以下构造一个为此的条件,使得去相位梯度与扩散编码梯度的重叠低于通过硬件给出的限制。对于给出的去相位动量所需的去相位梯度的振幅取决于可用的时间。参考图4,可用于接通去相位动量MA的时间tmin A,等于时间段t0和t1的最小值。可用于接通去相位动量MB的时间tmin A相应于时间段t2和t3的最小值。如下计算去相位梯度的所需振幅Gspoil的数值
Gspoil A=(Mspoil-GD(t1-t2))/tmin A Gspoil A=3(Mspoil-GD(t2-t1))/tmin A (25a)
Gspoil B=2(Mspoil-GD(t1-t2))/tmin B Gspoil B=6(Mspoil-GD(t2-t1))/tmin B (25b)
其中,振幅落入如下范围
GD∈[0,Mspoil/(t1-t2)] GD∈[0,Mspoil/(t2-t1)] (26)
并且振幅落入如下范围
中。去相位梯度和扩散梯度的和应该低于通过使用的梯度系统限制的最大振幅Gmax:
Gspoil A+GD≤Gmax Gspoil A+GD≤Gmax (28a)
GD(1-3(t2-t1)/tmin A)+3Mspoil/tmin A≤Gmax
Gspoil B+GD≤Gmax Gspoil B+GD≤Gmax (28b)
GD(1-6(t2-t1)tmin B)+6Mspoil/tmin B≤Gmax
因为该不等式的左边随着GD线性增加,所以,分析极限值,以便如下地检验不等式(28)的满足就够了:
Min:Mspoil/tmin A≤Gmax 3Mspoil/tmin A≤Gmax (29a)
Max:Mspoil/(t1-t2)≤Gmax Mspoil/(t2-t1)≤Gmax
Min:Mspoil/tmin B≤Gmax 6Mspoil/tmin B≤Gmax (29b)
Max:Mspoil/(t1-t2)≤Gmax Mspoil/(t2-t1)≤Gmax
如果满足该条件,则所需的去相位梯度可以叠加于扩散编码梯度,而不会超过通过梯度系统给出的限制。在此可以确定,“+”的情况比“-”的情况限制小。
通过使用附加的去相位梯度71-74,由此在图4中示出的拍摄序列70的情况下对于小的b值,即,对于扩散编码梯度的小的振幅也能有效抑制不期望的相干路径。特别地没有延长回波时间TE地实现了这点。
此外,还可以通过如下插入附加的自由度,使得记录在扩散模块51中的暂停。尽管这会减小扩散编码的效率,但是可能可以扩大可用的解的范围。
按照第一变形,在梯度57和58之间即,在t0和t1之间插入附加的暂停tp。利用该附加的暂停tp可以匹配等式(2)-(4)。因为附加的暂停不影响磁化的演变,所以利用相应的匹配可以使用事先描述的计算。对于“+”和“-”的情况在此又得到对于扩散模块的梯度的下列时间段:
“+” “-”
t3=2(t1-t2)=2/3T-1/3TE t3=2(t2-t1)=2/5T-1/5TE (37)
t2=T/2-tp/2-t3=-1/6T+1/3TE-tp/2 t2=T/2-tp/2-t3=1/10T+1/5TE-tp/2
t1=TE/2-T/2-tp/2+t3=1/6T+1/6TE-tp/2 t1=TE/2-T/2-tp/2+t3=-1/10T+3/10TE-tp/2
t0=T-TE/2-t3=1/3T-1/6TE t0=T-TE/2-t3=3/5T-3/10TE
由此可以看出,t0和t3保持不变,而t1和t2被减小tp/2。条件2T-TE>0没有改变。
在另一种替换方案中,在t0之后在t3之前插入暂停tp/2。又可以相应地匹配等式(2)-(4)。对于该情况也可以类似地应用上面示出的用于计算梯度脉冲持续时间ti的途径。得到时间段
“+” “-”
t0=1/6(2T-TE-2tp) t0=1/10(6T-4TE-6tp) (47)
t1=1/6(T+TE-tp) t1=1/10(-T+3TE+tp)
t2=1/6(-T+2TE+tp) t2=1/10(T+2TE-tp)
t3=1/6(4T-2TE-4tp) t3=1/10(4T-2TE-4tp)
如从这些等式可以看出的,对于t0的条件变得尖锐。条件2T-TE>0没有被松动。在两种变形中应该满足该条件,以便确保全面抑制不期望的信号相干路径。如已经提到的,这是双极的扩散编码序列的一般的限制。
按照另一个实施方式,建议使用附加的重聚焦脉冲,以便松动该条件。然而在此缺陷是,产生附加的相干路径和用于入射HF脉冲的附加所需的时间。
此外,还可以最大化时间段T,以便进一步改善对于按照前面描述的实施方式的时序的使用的条件。这例如可以通过使用短的准备时间tprep和短的读出时间tadc来实现。这可以通过短的回波串(例如小的分辨率、小的回波间隔、局部k空间采样、iPAT等)来实现。
对于其它改进可以将在其上满足自旋回波条件(3)的点移动到前面的时间(即,自旋回波与k空间中心的拍摄不一致)。由此有效减小了时间tadc。这会导致在k空间中心中的小的分量T2 *弛豫,而不是纯的T2弛豫。然而取决于使用的序列,回波时间t3可以被减小到直到原来的2Δtadc倍。
以下说明利用按照本发明的方法和按照本发明的磁共振设备的一种实施方式达到的结果。表1示出了举例对于不同的b值和对于三个不同的拍摄序列(单极的Stejskal-Tanner序列,具有附加的单独的扰相梯度的常规的双极的序列,和按照本发明的拍摄序列“双极+”的一种实施方式)的最小可实现的回波时间TE。使用具有不同的分辨率的两个成像草案。对于标准分辨率使用以下值:
视野FOV=230×230mm,分辨率RES=180×128,层厚SL=5mm,局部傅里叶拍摄PF=6/8,带宽BW=1502Hz/Pixel,回波间隔ES=0.73ms。由此,对于三个拍摄序列达到下列最小回波时间:
对于具有高分辨率的成像草案,使用下列值:FOV=230×230mm,RES=192×192,iPAT2,SL=5mm,PF=6/8,BW=1184Hz/Pixel,ES=0.94ms。由此得到下列最小回波时间:
可以看出,对于b值的整个考察范围,按照本发明的一种实施方式的双极+序列实现比常规的双极方法更短的回波时间。对于1000s/mm2的b值,双极+方法需要比单极的方法长很少的回波时间(≤5ms区别),并且甚至对于超过5000s/mm2的b值达到比单极方法更小的回波时间。
然而,在此利用按照本发明的方法可以比单极方法远远有效地抑制涡流。这在图5中示出。在采用矩形的扩散梯度形状(其利用单个时间常数产生指数衰减的涡流)的条件下,估计剩下的动态场失真。图5示出了在k空间中心中的剩下的涡流振幅与涡流的时间常数的依赖关系。曲线81示出了对于单极方法(在该方法中在时间段T期间接通恒定的梯度)的振幅,曲线80示出了在使用双极+方法的情况下振幅的估计。对于两种情况采用具有最大梯度振幅的扩散编码和预修正,该预修长将涡流场降低到最初振幅的0.01%。1μT/m的剩下的涡流场相应于在230mm的视野内部的半个像素的最大错误配准。尽管双极+方法的时序不是对于涡流的最大补偿被优化,但是与单极的方法相比,双极的梯度脉冲的使用还是将剩下的涡流场以5-10的系数降低。这表明,双极+方法也能有效抑制涡流并且避免相应的图像干扰,而同时达到短的TE时间。
为了相应检查MR信号、即不期望的相干路径的有效抑制,以多个扩散方向和权重拍摄多个扩散加权的图像。例如对于6、10、12、20、32和46个方向拍摄图像,以及以0.50s/mm2、500s/mm2和1000s/mm2的b值拍摄图像。此外,对于256个方向以0和1000s/mm2的b值以及以正常的以及相反的扩散矢量取向拍摄图像。利用双极+方法拍摄的图像数据没有不期望的相干路径的信号分量的迹象。特别是不可见任何干涉(“条纹伪影”)。
对于20个扩散方向利用常规的双极方法(90)、利用新的双极+方法(91)和利用标准单极方法(92)拍摄图像数据。图6示出了对于三种方法的每一种的标准偏差图。它们解释了通过取决于取向的涡流(亮的轮廓)的空间的失真。利用双极+方法(91),涡流感应的失真仅比利用常规的双极方法(90)的稍高。然而其远远小于利用单极方法(92)的。注意到,常规的方法在此对于该涡流的减小是最佳的。对相同的时间,对于此处使用的拍摄序列,与常规的双极方法相比,双极+方法使得可以减少回波时间大约10%,在该例子中是8ms。
可以组合前面描述的实施方式的特征。前面描述的实施方式使用了具有四个扩散编码梯度和一个EPI或分割的EPI读出模块的双自旋回波拍摄序列。在按照本发明的方法的和按照本发明的磁共振设备的其它实施方式中还可以使用其它拍摄序列,例如具有附加的或具有很少扩散编码梯度的附加的HF脉冲和其它读出模块。
通过根据读出模块优化扩散编码梯度的梯度动量,可以有效抑制不期望的相干路径。该方法组合常规的双极方法的优点(小的涡流感应的干扰、梯度负载均匀分布到两个GPA极),和常规的单极方法的优点(短的回波时间)。此外,对于以小的b值的成像,前面描述的接通不会延长回波时间的附加的扰相梯度是具有优势的。对模体的测量已经表明了该方法的可应用性和有效性。本发明例如可以被用于改进对于活体扩散成像的图像质量。
Claims (22)
1.一种用于利用拍摄序列(50)扩散加权地拍摄磁共振信号的方法,该拍摄序列包括具有多个扩散编码梯度(56-59)的扩散模块(51)和具有用于拍摄磁共振信号的读出梯度(60)的读出模块(52),其中,所述拍摄序列被配置为用于拍摄相应于预先确定的信号相干路径的磁共振信号,
其中,所述方法包括:
-利用所述拍摄序列(50)拍摄磁共振信号,其中,在所述拍摄序列(50)期间以预先确定的梯度动量(M0-M3)接通扩散编码梯度(56-59),
其特征在于,这样设置所述扩散编码梯度(56-59)的梯度动量(M0-M3),使得相应于与预先确定的相干路径不同的相干路径的磁共振信号被减小,其中,基于阈值进行所述梯度动量(M0-M3)的设置,以便实现预先确定的减小。
2.根据权利要求1所述的方法,其特征在于,根据所使用的读出模块来(52)确定所述阈值。
3.根据权利要求2所述的方法,其特征在于,基于所述读出模块的读出梯度的产生的梯度动量(Mreadout)来确定所述阈值。
4.根据权利要求1所述的方法,其特征在于,所述阈值定义一个最小的去相位梯度动量(Mspoil),其中,这样进行扩散编码梯度的梯度动量的设置,使得对于拍摄序列的其它相干路径的每一个,由扩散编码梯度对于相应相干路径产生的梯度动量至少与最小的去相位梯度动量(Mspoil)一样大。
5.根据权利要求1所述的方法,其特征在于,这样确定所述阈值,使得磁共振信号的、相应于拍摄序列的其它相干路径的k空间的中央位置,位于利用该拍摄序列采样的k空间外部。
6.根据权利要求1所述的方法,其特征在于,所述阈值是最小的去相位梯度动量Mspoil并且按照等式
Mspoil=N·Mreadout
来确定,其中N≥1,在此读出梯度动量Mreadout能够按照Mreadout=RES/(γ·FOV)来确定,其中RES表示分辨率并且FOV表示读出模块(52)的视野,并且γ表示回转磁比。
7.根据权利要求1所述的方法,其特征在于,借助梯度脉冲长度来设置所述扩散编码梯度的梯度动量,其中,所述扩散模块(51)包括至少四个扩散编码梯度(56-59),所述至少四个扩散编码梯度的梯度脉冲长度作为参数分别为t0-t3,其中,能够通过三个由所述拍摄序列(50)确定的条件来确定这些参数中的三个,并且其中,基于所述阈值来确定至少一个剩下的参数。
8.根据权利要求7所述的方法,其特征在于,从对于其它相干路径产生的扩散编码梯度的梯度动量与阈值的比较来确定所述至少一个剩下的参数。
9.根据权利要求1所述的方法,其特征在于,借助梯度脉冲长度(t0-t3)设置扩散编码梯度的梯度动量,其中该设置考虑以下条件:
-由拍摄序列预先给出的扩散编码梯度的接通的总持续时间(T);
-由扩散编码梯度要达到的磁化的去相位,其相应于预先确定的相干路径;
-梯度脉冲持续时间的如下的匹配:使得通过拍摄序列要产生的自旋回波的去相位在读出模块的预先确定的拍摄时刻(61)进行。
10.根据权利要求1所述的方法,其特征在于,所述扩散模块包括至少四个扩散编码梯度,其中分别在相对的方向上接通两个梯度。
11.根据权利要求1所述的方法,其特征在于,所述读出模块包括平面回波成像(EPI)读出序列或分割的平面回波成像(EPI)读出序列。
12.根据权利要求1所述的方法,其特征在于,所述扩散模块(51)还包括多个去相位梯度(71-74),这些去相位梯度产生相应于其它相干路径的磁共振信号的减小。
13.根据权利要求12所述的方法,其特征在于,能够根据所述阈值来设置所述去相位梯度(71-74)的梯度动量和/或所述去相位梯度(71-74)的极化。
14.根据权利要求12或13所述的方法,其特征在于,所述去相位梯度(71-74)与所述扩散梯度(56-59)重叠。
15.根据权利要求12或13所述的方法,其特征在于,当扩散编码梯度(56-59)的振幅(GD)低于预先确定的振幅阈值时,接通所述去相位梯度。
16.根据权利要求12或13所述的方法,其特征在于,根据扩散编码梯度的振幅(GD)来设置所述去相位梯度的梯度动量(MA,MB)。
17.根据权利要求12或13所述的方法,其特征在于,所述阈值定义一个最小的去相位梯度动量(Mspoil),其中,这样进行扩散编码梯度和去相位梯度的梯度动量的设置,使得对于拍摄序列的其它相干路径的每个,由扩散编码梯度和去相位梯度对于相应相干路径产生的梯度动量至少与最小的去相位梯度动量(Mspoil)一样大。
18.根据权利要求1所述的方法,其特征在于,在至少两个扩散编码梯度之间的扩散模块(51)包括具有可变持续时间的暂停,用于获得对于扩散编码梯度的梯度动量的选择的其它自由度。
19.根据权利要求1所述的方法,其特征在于,所述拍摄序列包括具有两个重聚焦脉冲的双自旋回波序列。
20.根据权利要求19所述的方法,其特征在于,作为其它相干路径出现至少三个自由感应衰减、三个自旋回波和两个激励的回波,其中,这样进行扩散编码梯度的梯度动量的设置,使得相应于这些相干路径的磁共振信号分量被减小。
21.一种用于利用如下的拍摄序列(50)扩散加权地拍摄磁共振信号的磁共振设备,所述拍摄序列包括具有多个扩散编码梯度(56-59)的扩散模块(51)和具有用于拍摄磁共振信号的读出梯度的读出模块(52),其中所述拍摄序列被配置为用于拍摄相应于预先确定的信号相干路径的磁共振信号,该磁共振设备包括:
-利用拍摄序列(50)拍摄磁共振信号并且在所述拍摄序列期间以预先确定的梯度动量接通扩散编码梯度(56-59)的拍摄单元(21),
-被构造为用于控制所述拍摄单元(21)的控制单元(18),其中,所述控制单元(18)这样设置扩散编码梯度的梯度动量,使得相应于与预先确定的相干路径不同的相干路径的磁共振信号被减小,其中,基于阈值进行梯度动量的设置,以实现预先确定的减小。
22.根据权利要求21所述的磁共振设备,其中,所述磁共振设备被构造为执行权利要求1-20中任一项所述方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102009019895.4 | 2009-05-04 | ||
DE102009019895A DE102009019895B4 (de) | 2009-05-04 | 2009-05-04 | Verfahren und Magnetresonanzanlage zur diffusionsgewichteten Aufnahme von MR-Signalen |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101879063A CN101879063A (zh) | 2010-11-10 |
CN101879063B true CN101879063B (zh) | 2014-07-16 |
Family
ID=42978975
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201010169936.9A Expired - Fee Related CN101879063B (zh) | 2009-05-04 | 2010-05-04 | 用于扩散加权地拍摄磁共振信号的方法和磁共振设备 |
Country Status (3)
Country | Link |
---|---|
US (1) | US8436613B2 (zh) |
CN (1) | CN101879063B (zh) |
DE (1) | DE102009019895B4 (zh) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011034004A1 (ja) * | 2009-09-17 | 2011-03-24 | 株式会社 日立メディコ | 磁気共鳴イメージング装置及び傾斜磁場印加方法 |
CN102141603B (zh) * | 2010-01-28 | 2013-07-31 | 西门子(深圳)磁共振有限公司 | 一种平面回波成像方法和系统 |
US8970217B1 (en) | 2010-04-14 | 2015-03-03 | Hypres, Inc. | System and method for noise reduction in magnetic resonance imaging |
DE102010029932B4 (de) * | 2010-06-10 | 2012-10-04 | Siemens Aktiengesellschaft | Korrektur von Bildverzerrungen in diffusionsgewichteten Magnetresonanz-Bildern |
DE102011005084B3 (de) * | 2011-03-04 | 2012-06-28 | Siemens Aktiengesellschaft | Reduktion von Artefakten in der diffusionsgewichteten Bildgebung |
DE102011081411B4 (de) * | 2011-08-23 | 2013-04-11 | Friedrich-Alexander-Universität Erlangen-Nürnberg | Abtastmuster für iterative MR-Rekonstruktionsverfahren |
DE102011084289B4 (de) * | 2011-10-11 | 2013-06-27 | Siemens Aktiengesellschaft | Magnetresonanzbildgebung mit Präparation der Magnetisierung durch örtlich veränderliche Spoiler-Gradientenfelder |
DE102012217992A1 (de) * | 2012-10-02 | 2014-04-03 | Siemens Aktiengesellschaft | Verfahren und Magnetresonanzanlage zur Erfassung von MR-Daten mit Diffusionsinformation |
KR101967244B1 (ko) * | 2012-12-05 | 2019-04-09 | 삼성전자주식회사 | 자기 공명 영상 방법 및 장치 |
DE102013201616B3 (de) * | 2013-01-31 | 2014-07-17 | Siemens Aktiengesellschaft | TSE-basierte, gegen lokale B0-Feldvariationen unempfindliche MR-Mulitschicht-Anregung |
WO2014163164A1 (ja) * | 2013-04-04 | 2014-10-09 | 株式会社東芝 | 磁気共鳴イメージング装置 |
JP6510173B2 (ja) * | 2014-01-31 | 2019-05-08 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置 |
DE102014206929B4 (de) | 2014-04-10 | 2015-10-29 | Siemens Aktiengesellschaft | Unterdrückung unerwünschter Kohärenzpfade in der MR-Bildgebung |
DE102014220776B4 (de) * | 2014-10-14 | 2016-06-30 | Siemens Healthcare Gmbh | Verfahren zur Magnetresonanz-Bildgebung |
CN107615080A (zh) | 2015-05-29 | 2018-01-19 | 梅约医学教育与研究基金会 | 用于在具有非对称梯度的磁共振成像中进行伴随场校正的系统和方法 |
CN108020796B (zh) * | 2016-10-31 | 2021-10-29 | 西门子(深圳)磁共振有限公司 | 一种磁共振扩散加权成像方法和装置 |
CN109923426B (zh) | 2016-11-09 | 2021-09-03 | Cr发展公司 | 对样本执行扩散加权磁共振测量的方法 |
CN109143134B (zh) * | 2017-06-13 | 2021-01-29 | 西门子(深圳)磁共振有限公司 | 一种分段读出扩散加权成像方法、装置及可存储介质 |
AU2018329432A1 (en) * | 2017-09-07 | 2020-03-12 | Random Walk Imaging Ab | Method for performing diffusion weighted magnetic resonance measurements |
DE102017219382A1 (de) * | 2017-10-27 | 2019-05-02 | Siemens Healthcare Gmbh | Verfahren zur Steuerung eines Magnetresonanztomographie-Systems zur Diffusions-Bildgebung |
SE543292C2 (en) * | 2019-04-26 | 2020-11-17 | Cr Dev Ab | A method of performing diffusion weighted magnetic resonance measurements |
DE102020212281B4 (de) * | 2020-09-29 | 2022-05-12 | Siemens Healthcare Gmbh | Verfahren zur zeitsparenden Erzeugung einer B0-Karte basierend auf einer Doppelecho-Sequenz mit stimulierten Echos und Magnetresonanzvorrichtung |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10109511C2 (de) * | 2001-02-28 | 2003-03-27 | Max Planck Gesellschaft | Verfahren und Gerät zum Gewinnen von Daten für diffusionsgewichtete Magnetresonanz-Bildgebung |
US7202662B2 (en) * | 2002-12-11 | 2007-04-10 | The Board Of Trustees Of The Leland Stanford Junior University | Correction of the effect of gradient field non-linearities in phase contrast MRI |
DE102005053269B3 (de) * | 2005-11-08 | 2007-04-12 | Siemens Ag | Verfahren und Vorrichtung zur Ermittlung eines diffusionsgewichteten Bildes |
JP5269342B2 (ja) * | 2006-04-25 | 2013-08-21 | 株式会社東芝 | 磁気共鳴イメージング装置および磁気共鳴イメージング装置における撮影条件設定方法 |
US7622919B2 (en) * | 2006-07-31 | 2009-11-24 | Schlumberger Technology Corporation | Nuclear magnetic resonance measurement techniques in non-uniform fields |
US8125223B2 (en) * | 2006-10-03 | 2012-02-28 | Singapore Agency For Science, Technology And Research Act | Segmenting infarct in diffusion-weighted imaging volume |
GB0708655D0 (en) * | 2007-05-04 | 2007-06-13 | Imp Innovations Ltd | Imaging of moving objects |
US7804299B2 (en) * | 2007-05-04 | 2010-09-28 | Wisconsin Alumni Research Foundation | Diffusion weighted preparatory sequence for magnetic resonance imaging pulse sequence |
US8134363B2 (en) * | 2007-07-20 | 2012-03-13 | Medical College Of Georgia Research Institute, Inc. | Test object for use with diffusion MRI and system and method of synthesizing complex diffusive geometries using novel gradient directions |
US7994784B2 (en) * | 2007-07-20 | 2011-08-09 | Medical College Of Georgia Research Institute, Inc. | Systems and methods for rescaling image intensities with diffusion MRI to synthesize complex diffusive geometries |
-
2009
- 2009-05-04 DE DE102009019895A patent/DE102009019895B4/de not_active Expired - Fee Related
-
2010
- 2010-05-04 CN CN201010169936.9A patent/CN101879063B/zh not_active Expired - Fee Related
- 2010-05-04 US US12/773,121 patent/US8436613B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN101879063A (zh) | 2010-11-10 |
DE102009019895B4 (de) | 2011-05-12 |
US8436613B2 (en) | 2013-05-07 |
US20100277169A1 (en) | 2010-11-04 |
DE102009019895A1 (de) | 2010-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101879063B (zh) | 用于扩散加权地拍摄磁共振信号的方法和磁共振设备 | |
US8228063B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
CN102959388B (zh) | 利用压缩感测重建的动态对比度增强mr成像 | |
Yang et al. | Optimizing abdominal MR imaging: approaches to common problems | |
CN101632584B (zh) | 磁共振成像装置 | |
US7800367B2 (en) | Method and apparatus for generating T2* weighted magnetic resonance images | |
CN103238082B (zh) | 使用多点Dixon技术和低分辨率校准的MR成像 | |
KR100671092B1 (ko) | 자기 공명 촬영 장치 및 방법 | |
RU2605524C2 (ru) | Магнитно-резонансная визуализация с подавлением артефактов потока | |
US9488711B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
US7952356B2 (en) | System and method for reconstructing multi-spectral 3D MR images | |
JP5127841B2 (ja) | 磁気共鳴イメージング装置及び磁化率強調画像撮影方法 | |
CN104204838B (zh) | 具有Dixon类型水/脂肪分离和关于主磁场的不均匀性的先验知识的MRI | |
US20120074940A1 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
KR101453297B1 (ko) | 복수 타입의 자기 공명 영상들을 동시에 생성하는 장치 및 방법 | |
US10302729B2 (en) | Method and magnetic resonance apparatus for speed-compensated diffusion-based diffusion imaging | |
CN102713658B (zh) | 磁化率梯度绘图 | |
US10203387B2 (en) | MR imaging with enhanced susceptibility contrast | |
CN103282790B (zh) | 快速双对比度mr成像 | |
US11119175B2 (en) | Systems and methods for suppressing Nyquist ghost for diffusion weighted magnetic resonance imaging | |
CN107607896A (zh) | 在磁共振技术中减少伪影 | |
JP2002000581A (ja) | 磁気共鳴画像化を使用して得られた複素画像の合成 | |
JPH0241967B2 (zh) | ||
KR102190023B1 (ko) | 자기 공명 영상 장치 및 그 제어 방법 | |
US20220057467A1 (en) | Epi mr imaging with distortion correction |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220118 Address after: Erlangen Patentee after: Siemens Healthineers AG Address before: Munich, Germany Patentee before: SIEMENS AG |
|
TR01 | Transfer of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140716 |
|
CF01 | Termination of patent right due to non-payment of annual fee |