CN111615637B - 用于执行扩散加权磁共振测量的方法 - Google Patents
用于执行扩散加权磁共振测量的方法 Download PDFInfo
- Publication number
- CN111615637B CN111615637B CN201880055673.3A CN201880055673A CN111615637B CN 111615637 B CN111615637 B CN 111615637B CN 201880055673 A CN201880055673 A CN 201880055673A CN 111615637 B CN111615637 B CN 111615637B
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- gradient
- time
- field gradient
- refocusing pulse
- 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
Links
Images
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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/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
- G01R33/482—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 using a Cartesian trajectory
- G01R33/4822—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 using a Cartesian trajectory in three dimensions
-
- 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/50—NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences
-
- 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/5611—Parallel 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
- G01R33/5617—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using RF refocusing, e.g. RARE
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56572—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of a gradient magnetic field, e.g. non-linearity of a gradient magnetic field
- G01R33/56581—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of a gradient magnetic field, e.g. non-linearity of a gradient magnetic field due to Maxwell fields, i.e. concomitant fields
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)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Vascular Medicine (AREA)
- Nonlinear Science (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
Description
技术领域
本发明概念涉及用于执行扩散加权磁共振测量的方法。
背景
扩散加权磁共振成像(dMRI)可用于探测组织微结构(Nilsson等人,2013年,每个参考文献的详细列表在参考文献章节中找到),且具有科学和临床应用。在所谓的b-张量编码中的最新发展已经证明,当使用线性和球形b-张量编码的组合时,可以从成像数据收回微结构特征的更大特异性(Lasic等人,2014年;Westin等人,2016年;Szczepankiewicz等人,2016年)。通过使用超越传统脉冲波形范式的定制梯度波形来实现必要的编码(等人,2015年)。然而,并非所有波形都被创建为相同的。一些波形不太有效,需要较长的编码时间,因此降低了信噪比,其它波形由于所谓的伴随磁场而引起图像伪影(Bernstein等人,1998年)。先前的方法描述了对于给定的b张量形状和编码时间的给定分布如何最大化波形的功效(等人,2015),但没有考虑伴随梯度场的影响。
发明概念的概述
本发明概念的目的是提供用于执行扩散加权磁共振测量的方法,其最小化或至少降低了伴随磁场效应。
根据本发明概念的方面,提供了一种用于扩散加权磁共振成像的方法,该方法包括:
由磁共振成像扫描仪的梯度线圈生成时间相关的磁场梯度G(t)=[Gx(t)Gy(t)Gz(t)]T,其中梯度G相对于重聚焦脉冲在时间上是不对称的,并且其中梯度G使得为零,其中TE是回波时间以及h(t)是时间的函数,其在重聚焦脉冲之前的间隔期间为正而在重聚焦脉冲之后的时间间隔期间为负。
根据本发明概念的第二方面,提供了一种用于扩散加权磁共振成像的方法,该方法包括:
由磁共振成像扫描仪的梯度线圈生成时间相关的磁场梯度G=[Gx(t)Gy(t)Gz(t)]T,其中梯度G相对于重聚焦脉冲在时间上是不对称的,并且其中梯度G使得归因于T2*弛豫的衰减因子为1,其中根据t是回波位移时间,
其中kp=k·np,其中k是残余梯度矩(residual gradient moment),以及np定义了相位编码方向,
其中,Δk是由Δk=n/FOV给出的两个所获取的k空间线之间的距离,其中n是平行成像因子且FOV是视场,以及
其中,Δk是两个连续k空间线的中心回波的获取之间的时间。
根据本发明概念的第三方面,提供了一种用于扩散加权磁共振成像的方法,该方法包括:
由磁共振成像扫描仪的梯度线圈生成时间相关的磁场梯度G(t)=[Gx(t)Gy(t)Gz(t)]T,其中梯度G相对于重聚焦脉冲在时间上是不对称的,并且其中梯度G使得m=(Tr[AA])1/2小于阈值,
其中A=∫P1G(t)G(t)Tdt-∫P2G(t)G(t)Tdt,其中P1和P2表示在重聚焦脉冲之前和之后的时间间隔。
根据本发明概念的第四方面,提供了一种用于为扩散加权磁共振成像设计非对称时间相关磁场梯度G(t)=[Gx(t) Gy(t) Gz(t)]T的方法,该梯度G相对于扩散编码序列的重聚焦脉冲在时间上是不对称的,该方法包括:
执行确定磁场梯度G的波形分量[Gx(t) Gy(t) Gz(t)]T的波形优化过程,其中优化过程被约束以:
其中TE是回波时间以及h(t)是时间的函数,其在重聚焦脉冲的时刻之前的间隔期间为正以及在重聚焦脉冲的时刻之后的时间间隔期间为负,或最小化m=(Tr[AA])1/2,
其中A=∫P1G(t)G(t)Tdt-∫P2G(t)G(t)Tdt,其中P1和P2表示在重聚焦脉冲的时刻之前和之后的时间间隔。
本发明概念基于下面的见解:关于重聚焦脉冲不对称的梯度波形可以被设计成使得它们将伴随场的影响降低到可忽略的水平。因此,这样设计的磁场梯度可以被称为所谓的“麦克斯韦补偿”。此外,如从下文可以进一步理解的,这可以针对任何形状的b张量实现,例如平面的、扁长的或球形的b张量编码,同时保持任意定时分布、高功效(对于短回波时间的大b值)和旋转不变性。
尽管存在用于降低由于伴随梯度场而引起的图像伪影的现有技术方法,但它们涉及在扩散编码块之后施加包括位置相关校正梯度的校正块,并且通常仅允许来自成像体积的中心区域的信号的校正。校正块的详细形式也取决于成像体积的定向,因此需要针对每个定向进行计算。
相反,本发明概念涉及以导致伴随场梯度的最小化或至少被抑制的影响的方式来谨慎地设计磁场梯度。因此,由于伴随场梯度引起的信号衰减——其原本会产生测量伪影——可以在不使用位置相关的后施加校正梯度的情况下被减轻或避免。
进一步的优点是,在可以提供伴随场伪影的最小化或减少而不管扩散编码b张量的形状或定向的意义上,麦克斯韦补偿磁场梯度可以是旋转不变的。
又一个优点是,伴随场伪影不仅在来自样本的中心部分而且在来自外围区域的信号贡献中都被最小化或减少。这又使得该方法变得适于供激发多于一个薄切割(例如,同时多切割成像(SMS)技术)的成像读出策略使用。
在上文中,参数x、y和z表示轴,梯度G沿着该轴被施加。
此外,G(t)G(t)T应被理解为G(t)与自身的外积。因此,G(t)G(t)T可以形成3x3大小的矩阵。
相对于(即关于)重聚焦脉冲(即,180°RF脉冲)在时间上“不对称”的时间相关的磁场梯度G在重聚焦脉冲之前和之后可以不呈现同一性(在重聚焦脉冲之前的G的一部分可以不是在重聚焦脉冲之后的G的一部分的副本),梯度G也可以在重聚焦脉冲之前和之后不呈现时间反转同一性(关于重聚焦脉冲的对称性,即,在重聚焦脉冲之前的G的一部分可以不是在重聚焦脉冲之后的G的时间反转副本)。
根据一些实施例,h(t)可以是时间的函数,其对于0<t<TE/2为正而对于TE/2<t<TE为负。
根据一些实施例,磁场梯度G可以形成自旋回波序列的一部分。因此,该方法还可以包括在生成磁场梯度之后获取回波。获取可以包括回波平面成像(EPI)读出。
根据一些实施例,编码梯度G可以被配置成使得扩散编码张量表示b具有至少两个非零特征值,其中:
根据一些实施例,编码梯度G可以适用于产生各向同性扩散编码。
根据一些实施例,梯度G进一步使得衰减因子AFp为1,以及衰减因子AFs=sinc(ks·ST)为1,其中ks=k·ns,其中ns是切割平面的法线方向,并且其中ST是切割厚度。
第三方面的参数m表示由磁梯度G提供的指示伴随场效应的补偿度的指数。该参数m可被称为“麦克斯韦指数”。麦克斯韦指数m的优点在于它是旋转不变的标量。这可以便于设计梯度G以实现关于伴随场效应的期望效果。根据第三方面的一些实施例,梯度G使得m=(Tr[AA])1/2为零。
附图简述
参考附图通过对本发明概念的优选实施例的以下说明性且非限制性详细描述,将更好地理解本发明概念的以上以及附加的目的、特征和优点。在附图中,相似的参考数字将用于相似的元件,除非另有说明。
图1-图4示出了在没有施加和施加附加积分限制条件的情况下优化的示例梯度波形。
图5是MRI扫描仪的示意图。
图6是编码和检测序列的示意图。
详细描述
为了便于理解本发明概念,现在将提供一些理论概念的讨论。
理论
麦克斯韦方程表明,由MRI扫描仪的梯度线圈产生的线性场梯度伴随着在空间上相关的附加场(伴随场)(Norris和Hutchison,1990年;Bernstein等人,1998年;周等人,1998年;Meier等人,2008年)。它们对场的贡献可以用在笛卡尔坐标x、y和z(其中,原点在等角点处)中磁场的展开来近似:
其中B0是主磁场,Gx x、Gy y和Gz z是由梯度线圈产生的时间相关的场梯度,以及其他项是所谓的麦克斯韦项。
与“预期”梯度的偏差可能引起不期望的信号衰减,这对于扩散MRI来说尤其是一个问题(Baron等人,2012年)。不期望的衰减由两种效应引起:穿透平面相位色散和回波形成的延迟。这两种效应都可以根据残余梯度矩k进行预测,对于自旋回波序列,残余梯度矩k由下式给出:
k=kpre-kpost 方程2
以及
其中,t0和t1是关于kpre的从激励到180脉冲的时间以及关于kpost的从180脉冲到回波时间的时间。在自旋回波序列中,如果g(t)在180脉冲之前和之后的时间周期内是相同的,则该矩可以为零,因为那时kpre=kpost。在其他情况下,需要附加的分析。
我们可以根据下式从方程1的空间导数的伴随梯度计算有效梯度g(t):
g(t)=G(t)+Mr 方程4
其中,G(t)=[Gx(t)Gy(t)Gz(t)]T是由梯度线圈产生的期望的时变梯度,r=[x yz]T是在空间中的位置,以及
穿透平面相位色散由沿着切割编码方向的非零矩引起,由下式给出:
ks=k·ns 方程6
其中,ns是切割平面的法向量。得到的衰减由切割选择剖面定义。在这里,我们假设它是矩形函数,因此衰减因子(AF)由sinc函数给出(Baron等人,2012年):
AFs=sinc(ks·ST) 方程7
其中,ST是切割厚度。
如发明人所认识到的,存在引起信号衰减的之前尚未描述的第二潜在效应。在这里,回波时间中的延迟从沿着相位编码方向的残余梯度矩给出(假设使用了EPI读出),由下式给出:
kp=k·np 方程8
这个梯度矩根据下式被转换成回波位移时间t(其可替代地可以被标记为t’):
其中,Δk是在两个获取的k空间线之间的距离,由Δk=n/FOV给出,其中n是平行成像因子以及FOV是视场,以及Δt是在获取两个连续的k空间线的中心回波之间的时间。回波位移可以由于T2*弛豫根据下式得出附加的衰减:
如何最小化伴随梯度的影响
如发明人所认识到的,ks和kp都必须被最小化,以便除去不期望的信号衰减。为此实现的充分(但非必要)条件是下列条件的满足:
其中,TE是回波时间,h(t)是函数,其在0<t<TE/2时为正而在TE/2<t<TE时为负,以考虑到在自旋回波实验中的180脉冲的影响,以及G(f)G(f)T是G(f)与自身的外积,即,大小3x3的矩阵。更一般地,h(t)可以是函数,其在重聚焦脉冲之前的间隔期间为正而在重聚焦脉冲之后的时间间隔期间为负。在考虑到180脉冲之后,A的对角元素将包含Gx 2、Gy 2和Gz 2项的积分,而非对角元素将包含Gi Gj项。如果所有这些项都为零,那么从激励到回波形成的M的平均值将为零,因此k将为零,使得不应出现不想要的信号衰减。A=0的条件可以被实现为在等人(2015年)中描述的波形优化程序的附加条件。
在方程11中的问题表述的一个重要结果是满足条件的任何梯度波形都可以自由地旋转和重新按比例缩放,因为当A=0时,波形旋转和重新按比例缩放了R得出A’=R ART=0。因此,保证了伴随场伪影的最小化,而不考虑张量形状。
此后,接下来是利用上述推导的波形优化器的示例实现。
方法:在波形优化器中的伴随梯度限制条件的实现
可以使用NOW工具箱来实现波形优化器,该工具箱可以从https://github.com/isjol/NQW(其中,提交ID 6c69462)下载。可以进行以下修改以产生期望的梯度波形:
1.优化限制条件
该限制条件用于最小化梯度波形与自身的外积的分量。此外,通过改变第二部分的积分的符号来考虑重聚焦脉冲的影响。两个间隔、即在重聚焦之前和之后的梯度波形由前面和后面的索引表示,以及轴——沿着该轴施加梯度——由x、y和z表示。为了最小化每个分量的大小,在优化中考虑平方和,降到被表示为tolMaxwell的公差水平。因此,优化器的新部分根据下式呈现限制条件(const)的形式:
const=(sum(g(pre,x).*g(pre,x))-sum(g(post,x).*g(post,x))).^2+(sum(g(pre,x).*g(pre,y))-sum(g(post,x).*g(post,y))).^2+(sum(g(pre,x).*g(pre,z))-sum(g(post,x).*g(post,z))).^2+(sum(g(pre,y).*g(pre,y))-sum(g(post,y).*g(post,y))).^2+(sum(g(pre,y).*g(pre,z))-sum(g(post,y).*g(post,z))).^2+(sum(g(pre,z).*g(pre,z))-sum(g(post,z).*g(post,z))).^2-tolMaxwell
其中,g是沿着三个正交轴的时间相关的梯度波形。
下面示出了NOW框架的修改(用Matlab代码)createConstraintGradientFunction.m的实现。
createConstraintGradientFunction.m
2.公差包含
为了保证限制条件,需要向最小化函数提供公差。因此,根据下面的代码用problem.tolMaxwell来扩展对fmincon的调用:
[x,fval,exitflag,output,lambda,grad]=fmincon(@(x)objFun(x),x0,A,b,Aeq,beq,[],[],@(x)feval(nonlconFileName,x,problem.tolIsotropy,problem.gMaxConstraint,Problem.integralConstraint,Problem.targetTensor,problem.tolMaxwell),options);
结果
在波形优化器中实现了附加的积分限制条件。生成了四个波形,包括对梯度振幅的欧几里德范数或最大范数的使用或者在有和没有伴随梯度限制条件的情况。在图1-图4中示出了这些波形。图1-图2示出了在没有保证新限制条件(“无补偿(non-comp)”)的情况下优化的波形。作为结果,GTG的积分很大,大约103到104的量级(在这里,G表示Nx3矩阵,其中N表示时间点的数量)。图3-图4示出了限制条件被保证的波形,且作为结果GTG项明显更小并且总是低于100。
对于图1:
对于图2:
对于图3:
对于图4:
结果证明有可能获得用于各向同性扩散编码的有效梯度波形,而不会由于伴随场引起不期望的附加梯度累积。
可替代波形优化程序
还可以使用“麦克斯韦指数”m=(Tr[AA])1/2应小于预定阈值(即在公差内)作为波形优化过程的附加条件来设计不对称磁场梯度波形G(t)=[Gx(t)Gy(t)Gz(t)]T,使伴随场的影响降低到可忽略的水平。这种限制条件可以通过提供公差在NOW框架中得到保证,类似于上面针对A=0所概述的内容。
实施例的描述
根据本发明概念,提供了一种用于扩散磁共振成像(dMRI)的方法。可以使用在图5中示意性示出的MRI扫描仪1来执行该方法。参考图6,dMRI测量可以包括编码块,该编码块包括一个或更多个时间相关的磁场梯度G(t)=[Gx(t)Gy(t)Gz(t)]T。磁梯度G可以由MRI扫描仪的梯度线圈12生成。梯度线圈12可以包括用于生成梯度G的每个相应分量的线圈部分。可以通过磁梯度分量的相对定向和由扫描仪1的主磁体10生成的静态主磁场B0来控制梯度的定向。
编码块还可以包括被配置成提供弛豫编码(即,由于纵向和/或横向弛豫引起的衰减)的射频(RF)序列。RF序列可以由扫描仪1的RF发射系统14生成。RF序列和磁梯度G可以被生成并施加于样本或测量对象“S”,其例如可以是包括水的生物样本,诸如脑组织或任何器官细胞(的悬浮液)的活组织检查样本。更一般地,样本包括核自旋系统,其性质可以通过核磁共振技术来进行测量。如在图6中所指示的,编码块后面可以是检测块,在该检测块期间,由编码块产生的回波衰减信号可以被获取。检测块可以包括EPI读出。信号可以由扫描仪1的RF接收系统16获取。所获取的信号可以被采样和数字化,并作为测量数据存储在扫描仪1的存储器18中。测量数据可以由扫描仪1的图像处理器20处理以生成样本的数字图像,该数字图像例如可以显示在连接到扫描仪1的监视器22上。
MRI扫描仪可以包括控制器24,其用于控制扫描仪1、特别是磁体10、梯度线圈12、RF发射系统14和接收系统16以及信号获取等的操作。控制器24可以在扫描仪1的一个或更多个处理器上实现,其中用于生成磁梯度和RF序列的控制数据可以使用软件指令来实现,该软件指令可以存储在计算机可读介质(例如,非暂时性计算机可读存储介质)上并且由设备的一个或更多个处理器执行。软件指令可以例如存储在设备的存储器的程序/控制部分中,设备的一个或更多个处理器可以访问该程序/控制部分。然而,控制器方法也可以代替使用软件指令而以设备/计算机的专用电路的形式实现,诸如仅举几个示例,在一个或更多个集成电路中、在一个或更多个专用集成电路(ASIC)或现场可编程门阵列(FPGA)中。
根据该方法,由梯度线圈12产生的时间相关的磁场梯度G(t)可以被配置为呈现关于RF序列的180°重聚焦脉冲的不对称性。还可以以A=为零的方式配置/设计梯度G,其中TE是回波时间,并且h(t)是时间的函数,其在重聚焦脉冲之前的间隔期间为正而在重聚焦脉冲之后的时间间隔期间为负。
梯度G也可以被定义为使得归因于T2*弛豫的衰减因子为1,其中t是根据kp=k·np的回波位移时间,其中k是残余梯度矩,以及np定义了相位编码方向。Δk是由Δk=n/FOV给出的在两条所获取的k空间线之间的距离,其中n是平行成像因子以及FOV是视场。Δt是获取两个连续的k空间线的中心回波之间的时间。
梯度G还可以被配置成使得衰减因子AFs=sinc(ks·ST)为1,其中ks=k·ns,其中ns是切割平面的法线方向,并且其中ST是切割厚度。
梯度G也可以被定义成使得m=(Tr[AA])1/2为零或至少小于阈值。阈值可以是预定阈值集,使得伴随场的影响可以降低到可忽略的水平。
可以例如使用波形优化器(诸如,以上所述的所修改的波形优化器)来确定满足这些条件中的一个或更多个的梯度波形G(t)=[Gx(t)Gy(t)Gz(t)]T。更一般地,可以通过执行被配置为确定磁场梯度G的波形分量[Gx(t)Gy(t)Gz(t)]T的波形优化过程来设计不对称梯度波形G(t)。可以在下面的限制条件之一下执行优化程序:
其中TE是回波时间,并且h(t)是时间的函数,其在重聚焦脉冲的时刻之前的间隔期间为正而在重聚焦脉冲的时刻之后的时间间隔期间为负,或者
最小化m=(Tr[AA])1/2
其中,A=∫P1G(t)G(t)Tdt-∫P2G(t)G(t)Tdt,其中P1和P2表示在重聚焦脉冲的时刻之前和之后的时间间隔。波形优化器可以配置有最小化容差。换句话说,波形优化器可以被配置成将A或m最小化为在预定公差内(即,在从零开始的公差内),该预定公差可以根据设备灵敏度和测量要求来确定。波形优化器可以被实现为包括具有指令的非暂时性计算机可读存储介质的计算机程序产品,该指令适于在由(计算机的)处理器执行时执行用于确定梯度G的波形分量的方法。如本文使用的术语“非暂时性”可以涵盖任何有形存储介质,例如非易失性存储器(诸如ROM(旨在提供数据的持久存储))以及易失性存储器(诸如RAM(不旨在永久地存储数据))。因此,“非暂时性”应该被考虑为对存储介质本身的限制而不是对数据存储持久性的限制。
一旦梯度波形G(t)已经被确定,定义编码序列和梯度波形的数据就可以被提供到扫描器1,其中dMRI测量可以被执行。
本发明的设计梯度波形G的方式与平面的、扁长的以及球形的b张量编码兼容。例如,编码梯度G可以具有使得扩散编码张量表示b具有至少两个非零特征值的形式。如在本领域中已知的,张量b由给出,其中是失相向量。编码梯度G可以被配置成使得张量b确切地具有三个非零特征值,从而提供3D扩散编码。各向异性和各向同性3D扩散编码都是可能的,其中各向异性3D编码暗示张量b的特征值不是都相等的,以及各向同性3D编码暗示张量b的特征值是都相等的。
在上文中,主要参考有限数量的示例来描述本发明概念。然而,如本领域中的技术人员容易认识到的,除了以上公开的示例之外的其他示例在如由所附权利要求限定的本发明概念的范围内同样是可能的。
参考文献
接下来是在上文中提到的参考文献的列表:
Baron,C.A.,Lebel,R.M.,Wilman,A.H.&Beaulieu,C.2012。The effect ofconcomitant gradient fields on diffusion tensor imaging。Magnetic Resonance inMedicine,68,1 190-201。
Bernstein,M.A.,Zhou,X.J.,Polzin,J.A.,King,K.F.,Ganin,A.,Pelc,N.J.,&Glover,G.H.(1998)。Concomitant gradient terms in phase contrast MR:analysisand correction。Magn Reson Med,39(2),300-308。
Lasic,S.,Szczepankiewicz,F.,Eriksson,S.,Nilsson,M.&Topgaard,D.2014。Microanisotropy imaging:quantification of microscopic diffusion anisotropyand orientational order parameter by diffusion MRI with magic-angle spinningof the q-vector。Frontiers in Physics,2,11。
Meier,C,Zwanger,M.,Feiweier,T.,&Porter,D.(2008)。Concomitant fieldterms for asymmetric gradient coils:consequences for diffusion,flow,and echo-planar imaging.Magn Reson Med,60(1),128-134。http://doi.org/10.1002/mrm.21615
Nilsson,M.,van Westen,D.,Stahlberg,F.,Sundgren,P.C,&Latt,J.(2013)。Therole of tissue microstructure and water exchange in biophysical modelling ofdiffusion in white matter。Magnetic Resonance Materials in Physics,Biology andMedicine,26(4),345-370。
Norris,D.G.&Hutchison,J.M.S.1990年。Concomitant magnetic fieldgradients and their effects on imaging at low magnetic field strengths。Magnetic Resonance Imaging,8,33-37。
J.,Szczepankiewicz,F.,Nilsson,M.,Topgaard,D.,Westin,C.F.&Knutsson,H.2015。Constrained optimization of gradient waveforms forgeneralized diffusion encoding。Journal of Magnetic Resonance,261,157-168。
Szczepankiewicz,F.,Van Westen,D.,Englund,E.,Westin,C.F.,Stahlberg,F.,Latt,J.,Sundgren,P.C.&Nilsson,M.2016年。The link between diffusion MRI andtumor heterogeneity:Mapping cell eccentricity and density by diffusionalvariance decomposition(DIVIDE)。Neuroimage,142,522-532。
Westin,C.-F.,Knutsson,H.,Pasternak,O.,Szczepankiewicz,F.,Ozarslan,E.,van Westen,D.等人(2016)。Q-space trajectory imaging for multidimensionaldiffusion MRI of the human brain。Neuroimage,135(C),345-362。
Zhou,X.J.,Du,Y.P.,Bernstein,M.A.,Reynolds,H.G.,Maier,J.K.&Polzin,J.A.1998。Concomitant magnetic-field-induced artifacts in axial echo planarimaging。Magnetic Resonance in Medicine,39,596-605。
Claims (20)
2.根据权利要求1所述的方法,其中h(t)是时间的函数,该函数对于0<t<TE/2为正而对于TE/2<t<TE为负。
3.根据权利要求1所述的方法,其中所述磁场梯度G形成自旋回波序列的一部分。
4.根据权利要求2所述的方法,其中所述磁场梯度G形成自旋回波序列的一部分。
6.根据权利要求1-4中任一项所述的方法,其中所述磁场梯度G适用于产生各向同性扩散编码。
7.根据权利要求5所述的方法,其中所述磁场梯度G适用于产生各向同性扩散编码。
9.根据权利要求8所述的方法,其中所述磁场梯度G进一步使得衰减因子AFs=sinc(ks·ST)为1,其中ks=k·ns,其中ns是切割平面的法线方向,并且其中ST是切割厚度。
10.根据权利要求8所述的方法,其中所述磁场梯度G形成自旋回波序列的一部分。
11.根据权利要求9所述的方法,其中所述磁场梯度G形成自旋回波序列的一部分。
13.根据权利要求8-11中任一项所述的方法,其中所述磁场梯度G适用于产生各向同性扩散编码。
14.根据权利要求12所述的方法,其中所述磁场梯度G适用于产生各向同性扩散编码。
15.一种用于扩散加权磁共振成像的方法,包括:
由磁共振成像扫描仪的梯度线圈生成时间相关的磁场梯度G(t)=[Gx(t)Gy(t)Gz(t)]T,其中所述磁场梯度G相对于重聚焦脉冲在时间上是不对称的,并且其中所述磁场梯度G使得m=(Tr[AA])1/2小于阈值,其中m是麦克斯韦指数,表示由所述磁场梯度G提供的指示伴随场效应的补偿度,
其中A=∫P1G(t)G(t)Tdt-∫P2G(t)G(t)Tdt,其中P1和P2表示在所述重聚焦脉冲之前和之后的时间间隔。
16.根据权利要求15所述的方法,其中所述磁场梯度G使得m为零。
19.根据权利要求15-18中任一项所述的方法,其中所述磁场梯度G适用于产生各向同性扩散编码。
20.一种用于为扩散加权磁共振成像设计非对称时间相关的磁场梯度G(t)=[Gx(t)Gy(t)Gz(t)]T的方法,所述方法包括:
执行确定所述磁场梯度G的波形分量[Gx(t)Gy(t)Gz(t)]T的波形优化过程,其中所述优化过程被约束以:
其中TE是回波时间以及h(t)是时间的函数,该函数在重聚焦脉冲的时刻之前的间隔期间为正而在所述重聚焦脉冲的时刻之后的时间间隔期间为负,或
最小化m=(Tr[AA])1/2,
其中m是麦克斯韦指数,表示由所述磁场梯度G提供的指示伴随场效应的补偿度,其中A=∫P1G(t)G(t)Tdt-∫P2G(t)G(t)Tdt,其中P1和P2表示在所述重聚焦脉冲的时刻之前和之后的时间间隔。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201762555165P | 2017-09-07 | 2017-09-07 | |
US62/555,165 | 2017-09-07 | ||
PCT/SE2018/050901 WO2019050462A1 (en) | 2017-09-07 | 2018-09-07 | METHOD FOR MAKING DIFFUSION WEIGHTED MAGNETIC RESONANCE MEASUREMENTS |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111615637A CN111615637A (zh) | 2020-09-01 |
CN111615637B true CN111615637B (zh) | 2022-06-14 |
Family
ID=65634235
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201880055673.3A Active CN111615637B (zh) | 2017-09-07 | 2018-09-07 | 用于执行扩散加权磁共振测量的方法 |
Country Status (9)
Country | Link |
---|---|
US (1) | US11525880B2 (zh) |
EP (1) | EP3679388A4 (zh) |
JP (1) | JP7512190B2 (zh) |
KR (1) | KR102701873B1 (zh) |
CN (1) | CN111615637B (zh) |
AU (1) | AU2018329432A1 (zh) |
BR (1) | BR112020004488A8 (zh) |
CA (1) | CA3074649A1 (zh) |
WO (1) | WO2019050462A1 (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SE543292C2 (en) * | 2019-04-26 | 2020-11-17 | Cr Dev Ab | A method of performing diffusion weighted magnetic resonance measurements |
KR102701517B1 (ko) | 2020-07-15 | 2024-08-30 | 주식회사 엘지에너지솔루션 | 배터리 관리 장치, 배터리 팩, 에너지 저장 시스템 및 배터리 관리 방법 |
US11768264B2 (en) * | 2021-04-30 | 2023-09-26 | Case Western Reserve University | System and method for magnetic resonance fingerprinting with relaxation and diffusion data acquisition |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1214622A (zh) * | 1996-03-28 | 1999-04-21 | 株式会社日立医药 | 磁共振成像装置和方法 |
CN101879063A (zh) * | 2009-05-04 | 2010-11-10 | 西门子公司 | 用于扩散加权地拍摄磁共振信号的方法和磁共振设备 |
CN102264291A (zh) * | 2008-12-26 | 2011-11-30 | 株式会社日立医疗器械 | 磁共振成像装置及读出梯度磁场误差修正方法 |
CN102428383A (zh) * | 2009-05-22 | 2012-04-25 | Cr发展公司 | 磁共振成像的方法和系统及其用途 |
CN104363829A (zh) * | 2012-05-23 | 2015-02-18 | 皇家飞利浦有限公司 | 多回波presto |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3369688B2 (ja) * | 1993-12-27 | 2003-01-20 | 株式会社日立製作所 | 核磁気共鳴を用いた検査装置 |
JPH09285457A (ja) * | 1996-04-25 | 1997-11-04 | Toshiba Corp | 磁気共鳴イメージング装置 |
US6150815A (en) * | 1997-04-10 | 2000-11-21 | The University Of British Columbia | Method of correcting for magnetic field inhomogeneity in magnetic resonance imaging |
WO2013025487A1 (en) * | 2011-08-12 | 2013-02-21 | The United States Of America, As Represented By The Secretary, Dpt Of Health And Human Services | Spin echo sequences for diffusion weighted imaging of moving media |
DE102011082009B3 (de) | 2011-09-01 | 2013-02-28 | Siemens Aktiengesellschaft | Schichtselektive Detektion und Korrektur von fehlerhaften Magnetresonanz-Bilddaten in Schicht-Multiplexing-Messsequenzen |
JP6176825B2 (ja) | 2013-02-28 | 2017-08-09 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 磁気共鳴装置 |
WO2014164963A1 (en) * | 2013-03-12 | 2014-10-09 | Levin Pavel | Method and apparatus for magnetic resonance imaging |
US9897678B2 (en) | 2013-04-19 | 2018-02-20 | General Electric Company | Magnetic resonance imaging data correction methods and systems |
US9903930B2 (en) * | 2013-11-06 | 2018-02-27 | King's College London | RF pulses for magnetic resonance |
CN105980876B (zh) | 2014-02-10 | 2019-05-17 | Cr发展公司 | 用于定量样品中的各向同性弥散和/或各向异性弥散的方法 |
EP2955536B1 (en) | 2014-06-12 | 2020-08-12 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | MRI method to quantify iron amount in tissues using diffusion magnetic resonance imaging |
CN105433944B (zh) | 2014-07-31 | 2018-07-03 | 西门子公司 | 用于获取对象的磁共振数据的方法及装置 |
DE102015205693A1 (de) | 2015-03-30 | 2016-10-06 | Siemens Healthcare Gmbh | Geschwindigkeitskompensierte diffusionssensibilisierte Diffusionsbildgebung |
WO2016180983A1 (en) | 2015-05-13 | 2016-11-17 | Koninklijke Philips N.V. | Multi-echo mri with correction of concomitant gradient-induced phase errors |
US10613174B2 (en) | 2015-10-29 | 2020-04-07 | Siemens Healthcare Gmbh | Method and magnetic resonance apparatus for maxwell compensation in simultaneous multislice data acquisitions |
SE538834C2 (sv) | 2015-12-29 | 2016-12-20 | Cr Dev Ab | Method of extracting information about a sample by nuclear magnetic resonance measurements |
EP3408679A4 (en) | 2016-01-25 | 2019-09-11 | The Regents of The University of California | OPTIMIZED WEIGHTED DIFFUSION IMAGING SYSTEM AND METHOD |
DE102016200975A1 (de) | 2016-01-25 | 2017-07-27 | Siemens Healthcare Gmbh | Diffusionsgewichtete Magnetresonanz-Bildgebung |
EP3441781A1 (de) * | 2017-08-11 | 2019-02-13 | Siemens Healthcare GmbH | Beschleunigte magnetresonanz-messung |
-
2018
- 2018-09-07 WO PCT/SE2018/050901 patent/WO2019050462A1/en unknown
- 2018-09-07 BR BR112020004488A patent/BR112020004488A8/pt not_active Application Discontinuation
- 2018-09-07 CN CN201880055673.3A patent/CN111615637B/zh active Active
- 2018-09-07 JP JP2020512717A patent/JP7512190B2/ja active Active
- 2018-09-07 EP EP18854163.5A patent/EP3679388A4/en active Pending
- 2018-09-07 KR KR1020207009466A patent/KR102701873B1/ko active IP Right Grant
- 2018-09-07 US US16/645,569 patent/US11525880B2/en active Active
- 2018-09-07 AU AU2018329432A patent/AU2018329432A1/en not_active Abandoned
- 2018-09-07 CA CA3074649A patent/CA3074649A1/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1214622A (zh) * | 1996-03-28 | 1999-04-21 | 株式会社日立医药 | 磁共振成像装置和方法 |
CN102264291A (zh) * | 2008-12-26 | 2011-11-30 | 株式会社日立医疗器械 | 磁共振成像装置及读出梯度磁场误差修正方法 |
CN101879063A (zh) * | 2009-05-04 | 2010-11-10 | 西门子公司 | 用于扩散加权地拍摄磁共振信号的方法和磁共振设备 |
CN102428383A (zh) * | 2009-05-22 | 2012-04-25 | Cr发展公司 | 磁共振成像的方法和系统及其用途 |
CN104363829A (zh) * | 2012-05-23 | 2015-02-18 | 皇家飞利浦有限公司 | 多回波presto |
Non-Patent Citations (3)
Title |
---|
"An Optimal Dimensionality Sampling Scheme on the Sphere with Accurate and Efficient Spherical Harmonic Transform for Diffusion MRI";IEEE Signal Processing Letters;《IEEE Signal Processing Letters》;20161231;全文 * |
"MRI不同序列在膝关节骨损伤中的应用";龚美琳;《放射学实践》;20071231;第22卷(第10期);全文 * |
"磁共振扩散加权成像量化诊断肝脏局灶性病变的研究进展";张海兵;《中华临床医师杂志》;20111231;第5卷(第23期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
US11525880B2 (en) | 2022-12-13 |
EP3679388A1 (en) | 2020-07-15 |
JP2020533051A (ja) | 2020-11-19 |
WO2019050462A1 (en) | 2019-03-14 |
JP7512190B2 (ja) | 2024-07-08 |
KR102701873B1 (ko) | 2024-08-30 |
AU2018329432A1 (en) | 2020-03-12 |
US20200284865A1 (en) | 2020-09-10 |
BR112020004488A2 (pt) | 2020-09-08 |
BR112020004488A8 (pt) | 2023-01-24 |
EP3679388A4 (en) | 2021-09-29 |
KR20200087749A (ko) | 2020-07-21 |
CA3074649A1 (en) | 2019-03-14 |
CN111615637A (zh) | 2020-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Setsompop et al. | Pushing the limits of in vivo diffusion MRI for the Human Connectome Project | |
Tan et al. | High slew‐rate head‐only gradient for improving distortion in echo planar imaging: preliminary experience | |
Wilm et al. | Single‐shot spiral imaging enabled by an expanded encoding model: D emonstration in diffusion MRI | |
Gallichan et al. | Retrospective correction of involuntary microscopic head movement using highly accelerated fat image navigators (3D FatNavs) at 7T | |
Koch et al. | Optimization of static magnetic field homogeneity in the human and animal brain in vivo | |
Wilm et al. | Higher order reconstruction for MRI in the presence of spatiotemporal field perturbations | |
US10488479B2 (en) | MRI apparatus and image processing apparatus | |
CN111615637B (zh) | 用于执行扩散加权磁共振测量的方法 | |
CN105676155B (zh) | 磁共振成像系统的磁场不均匀性值获取方法和装置和失真校正方法和装置 | |
Stirnberg et al. | Rapid fat suppression for three‐dimensional echo planar imaging with minimized specific absorption rate | |
Eggenschwiler et al. | Improving T2‐weighted imaging at high field through the use of kT‐points | |
US20170146624A1 (en) | Simultaneous multi-slice mri measurement | |
Liu et al. | Imaging neural architecture of the brain based on its multipole magnetic response | |
US20160116560A1 (en) | Mr imaging with enhanced susceptibility contrast | |
CA2872348C (en) | Pulse sequence method for mri | |
Mueller et al. | Eddy current compensated double diffusion encoded (DDE) MRI | |
US9753108B2 (en) | Magnetic resonance imaging method for at least two separate radio-frequency transmit coils with time-delayed slice-selective excitation pulses | |
US10578696B2 (en) | System, method and computer-accessible medium for spectroscopic localization using simultaneous acquisition of double spin and stimulated echoes | |
Sandgaard et al. | Incorporating the effect of white matter microstructure in the estimation of magnetic susceptibility in ex vivo mouse brain | |
Blasche et al. | Magnet homogeneity and shimming | |
Wu et al. | Velocity selective spin labeling using parallel transmission | |
US11959986B2 (en) | MR imaging with spiral acquisition | |
Kopanoglu et al. | Specific absorption rate reduction using nonlinear gradient fields | |
Lowen et al. | Design of calibration‐free RF pulses for T 2 _2‐weighted single‐slab 3D turbo‐spin‐echo sequences at 7T utilizing parallel transmission | |
CN113625211A (zh) | 借助磁共振装置进行扩散成像的方法、磁共振装置、计算机程序和电子可读的数据载体 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20221124 Address after: Okap, Sweden Patentee after: Landmwork Imaging Co.,Ltd. Address before: Lund, Sweden Patentee before: Douglick Effective date of registration: 20221124 Address after: Lund, Sweden Patentee after: Douglick Address before: Lund, Sweden Patentee before: CR DEVELOPMENT AB |