CN101609132A - 确定控制高频发送线圈的脉冲序列的方法和磁共振设备 - Google Patents

确定控制高频发送线圈的脉冲序列的方法和磁共振设备 Download PDF

Info

Publication number
CN101609132A
CN101609132A CNA2009101468944A CN200910146894A CN101609132A CN 101609132 A CN101609132 A CN 101609132A CN A2009101468944 A CNA2009101468944 A CN A2009101468944A CN 200910146894 A CN200910146894 A CN 200910146894A CN 101609132 A CN101609132 A CN 101609132A
Authority
CN
China
Prior art keywords
equation
controlled variable
coil
equations
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.)
Granted
Application number
CNA2009101468944A
Other languages
English (en)
Other versions
CN101609132B (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 Healthineers AG
Original Assignee
Siemens AG
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Siemens AG filed Critical Siemens AG
Publication of CN101609132A publication Critical patent/CN101609132A/zh
Application granted granted Critical
Publication of CN101609132B publication Critical patent/CN101609132B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/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/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/288Provisions within MR facilities for enhancing safety during MR, e.g. reduction of the specific absorption rate [SAR], detection of ferromagnetic objects in the scanner room
    • 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
    • G01R33/5612Parallel RF transmission, i.e. RF pulse transmission using a plurality of independent transmission 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/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/5659Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the RF magnetic field, e.g. spatial inhomogeneities of the RF magnetic field
    • 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/483NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
    • G01R33/4833NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices

Abstract

本发明涉及一种确定用于控制在磁共振设备中具有一个或多个发送通道的高频发送线圈的脉冲序列的方法,其中执行如下步骤:规定用于磁共振设备的梯度线圈的对于脉冲序列的多个时间步骤确定的控制顺序,以便根据时间选择所选择的层的或者所选择的体积的区域;对于每个时间步骤在规定期望的目标磁化条件下和在考虑梯度线圈的控制顺序的条件下求解用于发送线圈的单个发送通道的控制参数的非线性方程组,其中,在空间和时间中离散的非线性方程组除了包括从布洛赫方程得到的按照控制参数的非线性方程之外至少包括另一个描述边界条件的附加方程。

Description

确定控制高频发送线圈的脉冲序列的方法和磁共振设备
技术领域
本发明涉及一种确定用于控制在磁共振设备中具有一个或多个发送通道的高频发送线圈的脉冲序列的方法。
背景技术
常规的磁共振设备对于控制用于从均衡磁化产生自旋的偏转的高频发送线圈使用这样构造的参数,产生的偏转场(通常也称为B1场)在不由人占据的测量区域内部尽可能均匀。事先指出,“高频发送线圈”的名称也完全包括发送和接收线圈。
然而,特别是在新的设备中、例如在具有大于或等于3特斯拉的基本磁场强度的设备中表明,该均匀性会由于患者的进入而被中断,因为在其中会产生涡流,该涡流本身又产生干扰场,该干扰场使得偏转场的和由此磁化的均匀性消失。涡流效应的这种上升是由使用更高的频率所引起的。
在上面提到的久已公知的情况下通常这样进行控制,在线圈的n个激励通道中分别错开2π/n地控制依次的通道,以便达到圆形极化场。该控制变形在空的空间中产生均匀的偏转场并且通常被称为“模式1”。
用于解决上述问题的第一种方案被称为“静态方案”。在此,对于激励的整个持续时间对每个通道保持相位和振幅为恒定,以便实现更均匀的磁化。然而,由此达到的改进还不够,并且特别是既不能优化患者的SAR负担也不能仅仅局部偏转自旋,而这通常是期望的。
为了解决最后提到的问题已经提出,在高频激励脉冲期间也使用在激励期间通常直到层选择梯度非激活的梯度线圈,以便时间分辨地选择待激励的层的特定区域。同时还提出,在激励的持续时间期间允许发送通道的相位和振幅的动态改变,从而产生脉冲序列。然后可以这样说,选择并且相应地激励特定的区域,由此在下一个时间步骤中处理另一个区域。为此规定用于梯度线圈的k空间轨迹,该k空间轨迹选择一种预定的方式的激励的层的体积的子空间(也称为体素)。
然而,在此出现的问题是,要解高复杂的方程组。该方程组直接来自于所谓的布洛赫方程
d M → dt = γ M → x B → 1 - - - ( 1 )
其中M是磁化,t表示时间并且B1表示偏转场。通常该方程还包含描述弛豫的项,但是可以忽略这些项,因为激励的持续时间比时间常数T1和T2短很多。如果下面提到布洛赫方程,则由此是指方程(1)。相应地作为在脉冲序列的持续时间(一个或多个毫秒)上的时间积分以及在时间离散的情况下作为在所有时间步骤上的总和而得到期望的均匀的或者在特定的区域中均匀的磁化。也就是在要确定,通过从其可以导出对于发送通道的控制参数的B1场和梯度场的哪个组合可以得出该磁化时,预先规定该期望的磁化。
然而,布洛赫方程由于其矢量结构和在空间和时间中要进行的离散化分裂为多个方程,它们形成方程组。如果考虑每个发送通道(例如8个)对于每个时间单元(通常在一个毫秒的整个激励脉冲的持续时间期间数百个)需要复数值的系数(也就是作为控制参数的振幅和相位),则在使用布洛赫方程的条件下得到具有数千个变量的方程组。在该方程组中讨论期望的磁化(例如均匀的)以及在所选择的层中的期望的磁化分布(例如在特定区域中均匀的)。
为此,公知的脉冲计算方法从可以对所产生的方程组线性化的小角度近似出发,这意味着,从小的偏转角出发。然后期望的磁化的实部和虚部被输入到该方程组中,其中磁化仅仅允许为均衡磁化的小的百分比(约<10%),由此该小角度近似是允许的。
该解决方案除了不允许大的偏转角度之外,还有另一个出现的问题不被注意。在HF入射期间叠加梯度场导致空间的选择性的激励。带来的总的损耗功率、也就是SAR(特殊吸收率)因此上升,这会导致超过预定的边界值并危害患者。此外,在单个激励的区域中出现所谓的热区,也就是在其中SAR局部特别高的区域。最后,不能或者仅仅很难观察到限制激励脉冲的其它参数,例如在高频放大器的被限制的输出功率情况下。
发明内容
因此,本发明要解决的技术问题是,提出一种相应改进的方法,利用该方法对于大的偏转角在考虑其它(线性和非线性)限制、特别是SAR限制的条件下,也可以确定用于达到期望的磁化的理想的激励脉冲。
为了解决该技术问题,在本文开头提到种类的方法中按照本发明设置如下步骤:
-规定用于磁共振设备的梯度线圈的对于脉冲序列的多个时间步骤确定的控制顺序,以便根据时间选择所选择的层的或者所选择的体积的区域,
-对于每个时间步骤在规定期望的目标磁化条件下和在考虑控制顺序的条件下求解用于发送线圈的单个发送通道的控制参数的非线性方程组,其中,在空间和时间中离散的非线性方程组除了包括从布洛赫方程得到的按照控制参数的非线性方程之外至少包括另一个描述边界条件的附加方程。
在按照本发明的方法中,因此首先规定用于梯度线圈的控制顺序。除了层选择梯度因此还给其余的梯度线圈按照该控制顺序通电流,以便在不同的时间步骤中选择在所选择的层或者所选择的体积内部的不同区域。在此,还可以考虑不是相关的层区域或者子体积。也就是可以在任意的区域中规定磁化。控制顺序同样已经是时间上和空间上离散的。空间上的离散化表达为所选择的层分布在不同的体素中。
用于为梯度线圈确定这样的控制顺序的方法在现有技术中已经公知。通常选择特别是三维的k空间轨迹,该k空间轨迹基本上遵循螺旋变化,并且通过位置频率的间隔,以便最后在所选择的区域中建立磁化。随着控制顺序的规定相应地已知,在哪个时间步骤中哪个区域受到高频发送线圈的脉冲的影响。
借助该规定和对期望的目标磁化的附加规定,此时按照本发明对于每个时间步骤应该求解用于发送线圈的单个发送通道的控制参数的非线性方程组。控制参数是已经提到的复数值的系数,其再现了用来进行至对应于发送通道的馈入点的馈入的振幅和相位。非线性方程组在空间和时间中如已经提到的那样离散化,并且除了包括由布洛赫方程(1)得出的、按照控制参数的非线性方程之外还包括至少另一个描述边界条件的方程。边界条件在此是可以在解中设置的任何的附加条件(也称为“标准”)。相应地不是只提出,求解作为以其非线性形式的方程组的布洛赫方程,而是利用产生的可能性,以便以附加的方程的形式将其它边界条件带入方程组。
因此,不进行(如在现有技术中进行的)方程组的线性化,而是优选地采用用于非线性方程组的求解方法,特别是非线性最小平方方法,以便近似地确定控制参数。相对现有技术仅仅这点就带来许多优点。这样,此时可以对任意的磁化或者说偏转角度计算高频脉冲。在通常的方法中由于线性化仅仅能有意义地计算小角度激励,而现在还可以处理大角度激励。
此外,非线性的解决方案还可以只规定期望的目标磁化的量值(Betrag),这对于磁共振图像拍摄过程是足够的。在使用线性化的方程组时这是不可能的,因为量值形成是非线性的过程。由于现在量值的说明就够了,由此提供用于实现其它目标、例如最小化患者负担的自由度。
另一个优点是,许多非线性求解方法允许对变量的限制,这目前相应于高频放大器的最大输出功率的限制。由此脉冲序列的确定可以与可用的高频硬件匹配。然而非线性解决方案的最大优点是,可以引入其它、特别也是非线性附加方程,这些非线性附加方程在解中可以直接地、也就是与期望的目标磁化同时地被考虑。
在按照本发明方法的优选实施方式中作为附加方程使用至少一个描述SAR负担的、按照控制参数的非线性方程。例如通过直接地在脉冲计算中讨论用于SAR的特定的边界值,非线性解决方案特别是可以将描述SAR负担的非线性方程直接地列入方程组的解中。由此可以避免患者暴露于太高的SAR中。在规定恒定的、也就是最大均匀目标磁化(每个体素和时间步骤一个方程)和同时非常小的损耗功率目标(例如0,提供在每个时间步骤中用于总损耗功率的另一个方程,或者如果在所有的时间步骤上求和,提供一个方程)的情况下非线性方程组的近似解提供脉冲序列,该脉冲序列导致磁化的较好直到非常好的均匀性和同时小的总损耗功率、也就是小的患者负担。
为了带来描述SAR相关的损耗功率的方程,基本上有两种可能性。一种是作为附加方程使用根据直接取决于控制参数的矢量值的散射变量和散射参数矩阵给出患者中的损耗功率的方程。损耗功率与散射参数矩阵的该关系实际上是普遍公知的并且通过下式表达
PL ≈ Σ N a + ( I - S + S ) a - - - ( 2 )
其中,PL表示总的损耗功率,N表示时间步骤,a表示对应于对于一个时间步骤的所有通道的复数值的系数(也就是控制参数)的散射参数,S表示散射参数矩阵以及I表示单元矩阵并且其中“+”表示“复数共轭的和转置的”。如通常的那样通过除以对应的权重得到SAR。
一方面可以借助三维的电磁模拟来确定散射参数矩阵,其中例如可以基于人模型,但是优选地还可以测量散射参数矩阵。在此,可以利用现有技术中的所有已知方法。使用包括散射参数矩阵的附加方程具有如下优点,即其可以被测量,然而缺陷是,只能观察总损耗功率以及每个时间步骤的损耗功率,使得由此不能定位热区。
因此,描述SAR负担的附加方程的另一种变形也是可能的。按照本发明的另一种实施方式,作为附加方程使用根据取决于控制参数的电场分布和传导率分布(
Figure G2009101468944D00051
)给出患者中损耗功率的方程。从作为基础的物理关系中产生的方程也可以表示为如下公式:
PL ≈ Σ v = 1 V σ V Σ n = 1 N | Σ ch R E → ( v , t n , ch ) | 2 - - - ( 3 )
其中,PL仍表示总的损耗功率,σ表示传导率分布,V和所属的连续指标v表示体素或者体积单元,N利用连续变量n给出时间步骤的数量,t表示时间,E表示电场分布,R利用连续变量ch表示发送通道。该公式具有如下优点,即通过观察单个v加数而不是整个方程,用该公式还可以考虑局部SAR。也就是优选地可以位置分辨地观察附加方程(3)。因此通过在方程组中的附加的非线性方程可以将最高局部SAR热区的可预定的数量一起考虑到脉冲计算中,以便将局部SAR值保持在限制之下或者说保持在与全局SAR值的特定比例之下。也就是,位置分辨的附加方程可以描述局部SAR值。如果限制在最高的局部SAR值(热区)出现的定义的区域上观察位置分辨的附加方程,则在此可以大大降低计算花费。在此5-10个这样的区域就已经足够了。按照本发明的方法使得可以首次有针对地致力于防止出现这样的SAR热区。
按照本发明,例如可以从事先的患者图像拍摄中、特别是磁共振和/或计算机断层成像拍摄确定组织分布,然后因为不同组织类型的取决于频率的传导率是已知的,从中可以确定传导率分布;而用于电场分布的测量方法是未知的。然而,电场分布可以特别具有优点地在给出在患者中的传导率分布的情况下通过三维电磁模拟来确定。在此,为了确定组织分布优选地可以使用人模型,该人模型特别是借助关于患者的输入数据确定组织分布。因此,例如可以使用电子的患者档案获得的关于患者的数据,如其性别、体重、身高、年龄等,以便在人模型中产生人的组织分布。然后它们形成用于三维电磁模拟的基础,其结果是每个单个发送通道的控制中的电场分布。以这种方式可以在按照本发明的方法中还使用方程(3)作为附加方程。
除了从来自事先的患者图像拍摄中的组织分布间接确定传导率分布之外,同样还可以借助电的阻抗断层成像直接确定传导率分布。作为另一种替换,还要提到的是借助在表格中存储的传导率图集确定传导率分布。
一般地,非线性方程组的解在使用用于描述边界条件的附加方程的条件下提供具有优势地实现所述方法的一系列其它可能性。如已经提到的,可以规定用于控制参数的、相应于高频硬件(特别是高频放大器)的限制的最小和/或最大值。
在按照本发明的方法的另一种优选实施方式中,为了调整解还使用对方程组的方程加权的加权系数。由此可以将方程组的方程对应一个值,该值给出该方程必须满足的重要性。为了将确定的解、也就是确定的控制参数,在部分地矛盾的目标(例如磁化的高均匀性和小的患者负担)之间朝向一个或另一个目标调整,这是特别具有优势的。例如,可以为从布洛赫方程得出的方程提供一个加权系数,为一个或多个附加方程提供另一个加权系数。
为了检验SAR边界值的保持,还可以在确定控制参数之后从这些控制参数确定SAR负担。这样的附加的检验对于患者安全性是具有优势的。如果SAR负担例如还是太高,则可以改变加权系数或者例如在使用其它边界值的条件下形成附加方程。利用按照本发明的方法原则上还可以考虑,执行迭代的优化过程,其中按照本发明的方法进行多个过程,以便最后获得在尽可能小的患者负担和最大的均匀性之间的理想的折中。
在此要再次指出的是,利用按照本发明的方法以改进的方式还可以选择性地利用目标磁化仅激励所选择的层的或者所选择的体积的特定区域。
如本文开头所示,特别是还可以观察不一定是一层,而是任意的体积。在按照本发明的方法的特殊实施方式中,由此可以根本不使用层选择梯度,这意味着,通过控制顺序仅仅控制不是用于层选择的梯度线圈。例如仅仅控制在x和y方向上的梯度线圈,从而进行在z方向上延伸的体积激励。
最后使用所确定的激励参数,以便(与梯度线圈的控制顺序同步地)控制高频发送线圈的不同发送通道并且由此获得相应的磁化。
按照本发明的方法使用的按照本发明的磁共振设备,例如可以包括输入单元以及控制单元。操作者可以经过输入单元输入期望的磁化,紧接着,用于实施按照本发明的方法的控制单元,在规定的用于梯度线圈的控制顺序条件下,自动地根据按照本发明的方法确定用于高频发送线圈的控制参数并且相应地控制高频发送线圈。
按照本发明的方法在此还可以被用于具有低的通道数的高频发送线圈。特别是提出如下实施方式,在这些实施方式中由高频放大器提供的信号可以被分配到两个或多个调制单元上,其调制参数相应于用于本发明的控制参数。同样对于这些情况可以优选地采用本方法。
附图说明
本发明的其它优点和细节在以下描述的实施例中以及结合附图给出。在此附图中:
图1示出了按照本发明的磁共振设备的示意图,
图2示出了按照本发明的方法的流程图,
图3示出了用于梯度线圈的控制顺序的k空间轨迹,
图4示出了借助按照本发明的方法实现的磁化分布,以及
图5示出了根据现有技术借助模式1激励实现的磁化分布。
具体实施方式
图1以示意图示出了可用以实施按照本发明的方法的磁共振设备1。磁共振设备1包括基本磁铁2,其产生具有强度大于或等于3特斯拉的静态磁场。此外设置三个梯度线圈3、4、5,其中梯度线圈3对应于z方向,在该方向上激励脉冲序列形成层选择梯度。梯度线圈4和5产生在x以及y方向上的梯度。此外设置了高频发送线圈的阵列6,在这种情况下是也可以用作接收线圈的鸟笼线圈。其具有八个发送通道,经过该八个发送通道发出通过控制参数确定的脉冲序列,以便使得待检查的患者的核自旋从均衡磁化偏转。为此所必要的高频场通过多个高频放大器7产生。在此一个放大器7馈入阵列6的一个发送天线。作为替换一个放大器7可以馈入多个发送天线,其中每个发送天线配备一个调制器和移相器,以允许与振幅和相位独立地提供电流。即使只有一个发送线圈,自然也可以应用按照本发明的方法。
对于磁共振设备1还设置了中央控制设备8,通过该控制设备可以控制磁共振设备1的不同组件、特别是线圈3-6。控制设备8被构造为用于实施按照本发明的方法,即在规定用于梯度线圈的控制顺序和期望的目标磁化的条件下计算用于脉冲序列的控制参数并且然后相应地控制梯度线圈3、4、5以及特别是阵列6的高频发送线圈的发送通道。此外,控制设备8还配备有输入单元9,经过该输入单元用户可以输入用于磁共振设备1的相应的参数,特别是用于按照本发明的方法的输入数据。
图2示出了按照本发明的方法的流程图。
在此,首先在步骤10中确定或者说找出用于梯度线圈的预定的控制顺序。在此作为待观察的参数输入对应的硬件和梯度线圈3、4、5本身的限制,框11。由梯度线圈3产生的层选择梯度在脉冲序列期间可以是连续地激活的;而在该脉冲序列期间在单个时间步骤中也接通梯度线圈4和5,以便选择所选择的层的特定区域。如在图3在示例性示出的k空间轨迹已经证明作为控制顺序是特别合适的。该k空间轨迹螺旋形地延伸并且由此采集位置频率的特定的间隔。每个单个的、在图3中示出的点在此相应于对于一个时间步骤的梯度线圈参数。
如果已知用于梯度线圈3、4、5的控制顺序,则在步骤12中借助非线性最小平方求解方法求解非线性方程组。该方程组由在时间和空间中离散的布洛赫方程(1)导出的非线性方程(其中控制参数经过相应的B1值进入这些方程中)以及由附加方程组成。这些附加方程(也可能只有一个)描述边界条件、本例中是SAR负担,并且在该实施例中本身是非线性的。在此,可以考虑描述SAR负担的两个方程,即方程(2)和(3),在以下还要详细讨论。在采用最大损耗功率(并且由此最大SAR)的条件下形成一个或多个附加方程,从而SAR边界值作为边界条件通过附加方程巧妙地进入方程组中,使得在该解中可以直接地考虑它而不必再进行费事的计算过程。
除了控制顺序,作为输入参数首先还需要用于所有体素或者体积单元的期望的目标磁化的量值,框13。期望的目标磁化是在布洛赫方程(1)上的积分。作为目标磁化可以设置在整个所选择的层中的均匀的目标磁化,然而还也可以设置任意的磁化分布,例如仅仅在特定区域中的磁化。
此外,对于步骤12规定用于控制参数的最小和最大值作为输入参数,框14。在使用最小平方求解方法时可以限制控制参数的大小。该最小和/或最大值通过硬件特性特别是高频放大器7来确定。
为了选择附加方程基本上有两种可能性,即方程(2)和(3)。两种可能性具有其优点和缺点并且相应地根据需要被选择。
用方程(2)、即散射参数矩阵的入口,关于散射参数的控制参数直接包含在方程中。此外,优选地可以测量散射参数矩阵。当然还可以,通过相应的三维电磁模拟、例如用人模型来确定散射参数矩阵。然而,方程(2)仅提供整个损耗功率以及用于一个时间步骤的整个损耗功率的入口。由此不能有针对地排除局部的热区。
在使用方程(3)的条件下情况不同。可以在对于每个体素或者说每个子体积的单个方程中求解方程(3),从而还可以局部预先给出SAR限制。然而缺陷是,迄今为止还没有公知在测量技术上如何测量包含控制参数的电场分布,从而必须通过模拟来确定该电场分布。为此必须已知患者的组织分布。然而这可以通过使用可调整的人模型来实现,例如可以基于患者参数如身高、体重、年龄、性别等调整该人模型。然后可以通过三维电磁模拟来确定电场分布。同样在方程(3)中讨论的传导率分布可在按照本发明的方法中间接地从来自事先的图像拍摄、特别是磁共振或者计算机断层成像拍摄的组织分布,或者直接从阻抗断层成像或者结合在表格中存储的传导率图集来确定。
如果在步骤12中解得方程组并且得出用于控制参数的近似值,则可以在可选步骤15中从中计算SAR负担、特别是局部SAR负担。以这种方式可以再次检验,是否所有限制值都保持,也就是呈现最佳的脉冲序列。如果不是这样,则可以如通过箭头16所示,例如相应地改变方程组的方程或者其加权系数,以便进行重新的计算。因此可以考虑优化过程,在该优化过程中多次重复步骤12和14,直到达到优化标准。
此外,非线性方程组的方程在按照图2的方法中具有给出相应评估的加权系数。通过设置该加权系数、例如用于由布洛赫方程导出的方程的加权系数和另一个用于附加方程的加权系数,可以方便地设置,在矛盾的要求情况下解遵循哪个规定,例如当均匀性的要求和小的SAR负担的要求矛盾时。
最后,在步骤17中控制设备8根据控制顺序和确定的控制参数控制梯度线圈3、4、5和高频发送线圈6的通道,从而得到最佳的脉冲序列,利用该脉冲序列既可以实现均匀性的改进也可以实现SAR负担的改进。
这点将针对在图4和5中的例子详细描述。在该例子中,考虑具有八个发送通道的身体线圈,其中在小角度激励(偏转角0.57°)的情况下要实现均衡磁化的1%的均匀磁化。在此,图4示出用按照本发明的方法的磁化分布的结果,其中各个阴影范围表示偏差,其中阴影部分的强度给出偏差的强度。图5示出在相同患者情况下在模式1激励时会达到的磁化分布。可以看出,存在明显较强的偏差。实际上均匀性利用按照本发明的方法相对于模式1激励以因数5得到改进(考虑磁化相对于在不含手臂的胸部的横截面中的磁化的平均值的标准偏差)。在此,总损耗功率仅为在相同脉冲持续时间情况下导致相同平均磁化的模式1激励的1.5倍。
按照本发明的方法还可以被用于大角度激励、例如48.6°的偏转角度。对于具有八个发送通道的身体线圈的相应的实验得出,均匀性利用按照本发明的方法可以按照几乎因数3被改进,其中总损耗功率利用按照本发明的方法仅为模式1激励的75%。
由此,利用按照本发明的方法在考虑SAR边界值的条件下可以达到极度的改进。

Claims (18)

1.一种确定用于控制在磁共振设备中具有一个或多个发送通道的高频发送线圈的脉冲序列的方法,其中执行如下步骤:
-规定用于磁共振设备的梯度线圈的对于脉冲序列的多个时间步骤确定的控制顺序,以便根据时间选择所选择的层的或者所选择的体积的区域,
-对于每个时间步骤在规定期望的目标磁化条件下和在考虑梯度线圈的控制顺序的条件下求解用于发送线圈的单个发送通道的控制参数的非线性方程组,
其中,在空间和时间中离散的非线性方程组除了包括从布洛赫方程得到的按照控制参数的非线性方程之外至少包括另一个描述边界条件的附加方程。
2.根据权利要求1所述的方法,其中,作为附加方程使用至少一个描述SAR负担的、按照控制参数的非线性方程。
3.根据权利要求2所述的方法,其中,作为附加方程使用根据直接取决于控制参数的矢量值的散射变量和散射参数矩阵给出患者中的损耗功率的方程。
4.根据权利要求3所述的方法,其中,借助三维的电磁模拟来确定所述散射参数矩阵。
5.根据权利要求3所述的方法,其中,测量所述散射参数矩阵。
6.根据权利要求2所述的方法,其中,作为附加方程使用根据取决于控制参数的电场分布和传导率分布给出患者中损耗功率的方程。
7.根据权利要求6所述的方法,其中,从由事先的患者图像拍摄中、特别是磁共振和/或计算机断层成像拍摄中得到的组织分布,间接确定所述传导率分布。
8.根据权利要求6或7所述的方法,其中,在给出在患者中的传导率分布的情况下通过三维电磁模拟来确定所述电场分布。
9.根据权利要求8所述的方法,其中,为了确定所述组织分布使用人模型,该人模型特别是借助关于患者的输入数据确定所述组织分布。
10.根据权利要求6所述的方法,其中,借助电的阻抗断层成像确定所述传导率分布。
11.根据权利要求6至10中任一项所述的方法,其中,位置分辨地考虑所述附加方程。
12.根据权利要求11所述的方法,其中,所述位置分辨的附加方程描述了局部SAR值。
13.根据权利要求1至12中任一项所述的方法,其中,为了解方程组使用非线性最小平方方法。
14.根据权利要求1至13中任一项所述的方法,其中,规定用于控制参数的最小和/或最大值。
15.根据权利要求1至14中任一项所述的方法,其中,为了调整解使用对方程组的方程加权的加权系数。
16.根据权利要求1至15中任一项所述的方法,其中,为了检验SAR边界值的保持在确定控制参数之后从这些控制参数中确定SAR负担。
17.根据权利要求1至16中任一项所述的方法,其中,通过所述控制顺序只控制未用于层选择的梯度线圈。
18.一种磁共振设备(1),被构造为用于实施按照权利要求1至17中任一项所述的方法。
CN200910146894.4A 2008-06-19 2009-06-19 确定控制高频发送线圈的脉冲序列的方法和磁共振设备 Active CN101609132B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102008029175A DE102008029175B4 (de) 2008-06-19 2008-06-19 Verfahren zur Ermittlung einer Pulssequenz zur Ansteuerung einer Hochfrequenz-Sendespule
DE102008029175.7 2008-06-19

Publications (2)

Publication Number Publication Date
CN101609132A true CN101609132A (zh) 2009-12-23
CN101609132B CN101609132B (zh) 2013-08-21

Family

ID=41334896

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200910146894.4A Active CN101609132B (zh) 2008-06-19 2009-06-19 确定控制高频发送线圈的脉冲序列的方法和磁共振设备

Country Status (4)

Country Link
US (1) US8120359B2 (zh)
JP (1) JP5641713B2 (zh)
CN (1) CN101609132B (zh)
DE (1) DE102008029175B4 (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102210587A (zh) * 2010-04-01 2011-10-12 西门子公司 用于确定磁共振系统控制序列的方法和装置
CN102288929A (zh) * 2011-05-11 2011-12-21 浙江大学 用于磁共振系统的激发脉冲序列生成器
CN102955143A (zh) * 2011-08-24 2013-03-06 西门子公司 磁共振信号的流量敏感的、取决于频率的激励
CN103364749A (zh) * 2012-03-30 2013-10-23 西门子公司 磁共振系统的控制
CN104573308A (zh) * 2013-10-16 2015-04-29 布鲁克碧奥斯平股份公司 包括通过维度缩减的参数空间内的优化确定目标场的匀场过程
CN106569158A (zh) * 2016-10-24 2017-04-19 浙江大学 一种降低特殊吸收率的压缩型螺旋梯度射频脉冲设计方法
CN107874757A (zh) * 2016-09-30 2018-04-06 通用电气公司 磁共振射频子系统及用于其的线圈去耦合装置及方法
CN111308405A (zh) * 2018-12-12 2020-06-19 西门子医疗有限公司 用于监控局部线圈的方法和设备

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8120358B2 (en) * 2006-04-13 2012-02-21 The Regents Of The University Of California Magnetic resonance imaging with high spatial and temporal resolution
DE102008061455B4 (de) * 2008-12-10 2011-03-17 Siemens Aktiengesellschaft Verfahren und Vorrichtung zur Bestimmung einer vorbestimmten Signalamplitude bei MR-Messungen
DE102008063630B4 (de) * 2008-12-18 2011-04-14 Siemens Aktiengesellschaft Verfahren zur Kontrolle einer Hochfrequenzsendeeinrichtung, Hochfrequenzsendeeinrichtung, Hochfrequenz-Kontrolleinrichtung und Magnetresonanztomographiesystem
DE102009024077B4 (de) * 2009-06-05 2012-09-13 Siemens Aktiengesellschaft Verfahren und Vorrichtung zur SAR-Überwachung bei Transmit-Array-Sendesystemen
DE102009025476B4 (de) * 2009-06-18 2012-03-08 Siemens Aktiengesellschaft Verfahren zur Ermittlung von Parametern zur Ansteuerung der Gradientenspulen und Hochfrequenzspulen einer Magnetresonanzeinrichtung
DE102010015066B4 (de) * 2010-04-15 2014-12-11 Siemens Aktiengesellschaft Verfahren und Einrichtung zur Ermittlung einer Magnetresonanzsystem-Ansteuersequenz und Verfahren zum Betrieb eines Magnetresonanzsystems
DE102010015044A1 (de) * 2010-04-15 2011-10-20 Siemens Aktiengesellschaft Verfahren und Einrichtung zur Ermittlung einer Magnetresonanzsystem-Ansteuersequenz und Verfahren zum Betrieb eines Magnetresonanzsystems
DE102011004432A1 (de) * 2011-02-20 2012-08-23 Forschungszentrum Jülich GmbH Magnetresonanzverfahren und Kernspintomograph zur Durchführung des Magnetresonanzverfahrens
DE102011005174A1 (de) * 2011-03-07 2012-09-13 Friedrich-Alexander-Universität Erlangen-Nürnberg Verfahren und Einrichtung zur Ermittlung einer Magnetresonanzsystem-Ansteuersequenz
DE102011006151B4 (de) 2011-03-25 2014-04-30 Siemens Aktiengesellschaft Ermittlung einer hinsichtlich eines HF-Energieparameters optimierten HF-Pulslänge eines Magnetresonanzsystems
DE102011006149B4 (de) * 2011-03-25 2012-10-31 Siemens Aktiengesellschaft Ermittlung einer Magnetresonanzsystem-Ansteuersequenz unter Optimierung der Sendebandbreite
DE102011006573B4 (de) * 2011-03-31 2014-09-04 Siemens Aktiengesellschaft Verfahren zur Ermittlung einer optimalen Struktur eines Hochfrequenzschirms einer Hochfrequenzantenne einer Magnetresonanzanordnung
DE102011007823B4 (de) * 2011-04-20 2012-11-15 Bruker Biospin Mri Gmbh Verfahren zur Erzeugung eines gewünschten zeitlichen Verlaufs des Magnetisierungszustandes in einem Untersuchungsobjekt im Rahmen eines Experimentes der magnetischen Resonanz
DE102011087915B4 (de) * 2011-12-07 2015-04-09 Siemens Aktiengesellschaft Verfahren zur Überwachung der Patientenbelastung in einer Magnetresonanzeinrichtung
DE102012212402B3 (de) * 2012-03-20 2013-10-17 Siemens Aktiengesellschaft Ermittlung einer MR-Sequenz unter Berücksichtigung der Energieverteilungsfunktion im k-Raum
DE102012205297B4 (de) * 2012-03-30 2015-02-12 Siemens Aktiengesellschaft Ermittlung einer Magnetresonanzsystem-Ansteuersequenz und Verfahren zum Betrieb eines Magnetresonanzsystems
JP6101352B2 (ja) * 2013-08-27 2017-03-22 株式会社日立製作所 磁気共鳴撮像装置および撮像パラメータ決定方法
JP6410633B2 (ja) * 2015-02-20 2018-10-24 株式会社日立製作所 磁気共鳴イメージング装置およびsar算出方法
JP6965073B2 (ja) * 2017-09-19 2021-11-10 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置
EP3644085A1 (de) * 2018-10-25 2020-04-29 Siemens Healthcare GmbH Auswahl von messspulen bei der magnetresonanz-bildgebung

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1433737A (zh) * 2002-01-23 2003-08-06 Ge医疗系统环球技术有限公司 磁共振成像装置
US6989673B2 (en) * 2003-11-26 2006-01-24 General Electric Company Method and apparatus to reduce RF power deposition during MR data acquisition
CN1851494A (zh) * 2005-04-22 2006-10-25 西门子公司 改进发送侧加速的体积选择磁共振成像的方法和设备
WO2007092740A3 (en) * 2006-02-07 2008-01-17 Univ California Method and system for the amplification of nuclear magnetic resonance imaging
US7336145B1 (en) * 2006-11-15 2008-02-26 Siemens Aktiengesellschaft Method for designing RF excitation pulses in magnetic resonance tomography
DE102006058162B3 (de) * 2006-12-09 2008-06-12 Bruker Biospin Mri Gmbh Verfahren zur Gewinnung von Amplituden- und Phasenverläufen von HF-Pulsen für die räumlich-selektive Anregung

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4709212A (en) * 1986-01-03 1987-11-24 General Electric Company Method of enhancing image signal-to-noise ratio by combining NMR images of differing pulse sequence timing
WO1999009428A1 (en) * 1997-08-13 1999-02-25 Btg International Ltd. Apparatus for and method of determining values of relaxation parameters
WO2007098011A2 (en) * 2006-02-17 2007-08-30 Regents Of The University Of Minnesota High field magnetic resonance
US7786729B2 (en) * 2006-10-31 2010-08-31 Regents Of The University Of Minnesota Method for magnetic resonance imaging

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1433737A (zh) * 2002-01-23 2003-08-06 Ge医疗系统环球技术有限公司 磁共振成像装置
US6989673B2 (en) * 2003-11-26 2006-01-24 General Electric Company Method and apparatus to reduce RF power deposition during MR data acquisition
CN1851494A (zh) * 2005-04-22 2006-10-25 西门子公司 改进发送侧加速的体积选择磁共振成像的方法和设备
WO2007092740A3 (en) * 2006-02-07 2008-01-17 Univ California Method and system for the amplification of nuclear magnetic resonance imaging
US7336145B1 (en) * 2006-11-15 2008-02-26 Siemens Aktiengesellschaft Method for designing RF excitation pulses in magnetic resonance tomography
DE102006058162B3 (de) * 2006-12-09 2008-06-12 Bruker Biospin Mri Gmbh Verfahren zur Gewinnung von Amplituden- und Phasenverläufen von HF-Pulsen für die räumlich-selektive Anregung

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102210587A (zh) * 2010-04-01 2011-10-12 西门子公司 用于确定磁共振系统控制序列的方法和装置
CN102210587B (zh) * 2010-04-01 2014-12-17 西门子公司 用于确定磁共振系统控制序列的方法和装置
CN102288929A (zh) * 2011-05-11 2011-12-21 浙江大学 用于磁共振系统的激发脉冲序列生成器
CN102288929B (zh) * 2011-05-11 2013-04-17 浙江大学 用于磁共振系统的激发脉冲序列生成器
CN102955143A (zh) * 2011-08-24 2013-03-06 西门子公司 磁共振信号的流量敏感的、取决于频率的激励
US9651643B2 (en) 2011-08-24 2017-05-16 Siemens Aktiengesellschaft Flow-sensitive, frequency-independent excitation of magnetic resonance signals
US9594137B2 (en) 2012-03-30 2017-03-14 Siemens Aktiengesellschaft Controlling magnetic resonance systems
CN103364749B (zh) * 2012-03-30 2016-12-28 西门子公司 磁共振系统的控制
CN103364749A (zh) * 2012-03-30 2013-10-23 西门子公司 磁共振系统的控制
CN104573308A (zh) * 2013-10-16 2015-04-29 布鲁克碧奥斯平股份公司 包括通过维度缩减的参数空间内的优化确定目标场的匀场过程
CN107874757A (zh) * 2016-09-30 2018-04-06 通用电气公司 磁共振射频子系统及用于其的线圈去耦合装置及方法
CN107874757B (zh) * 2016-09-30 2023-06-09 通用电气公司 磁共振射频子系统及用于其的线圈去耦合装置及方法
CN106569158A (zh) * 2016-10-24 2017-04-19 浙江大学 一种降低特殊吸收率的压缩型螺旋梯度射频脉冲设计方法
CN106569158B (zh) * 2016-10-24 2019-02-05 浙江大学 一种降低特殊吸收率的压缩型螺旋梯度射频脉冲设计方法
CN111308405A (zh) * 2018-12-12 2020-06-19 西门子医疗有限公司 用于监控局部线圈的方法和设备
US11226384B2 (en) 2018-12-12 2022-01-18 Siemens Healthcare Gmbh Method and device for monitoring local coils

Also Published As

Publication number Publication date
JP5641713B2 (ja) 2014-12-17
US20090322329A1 (en) 2009-12-31
DE102008029175A1 (de) 2009-12-24
CN101609132B (zh) 2013-08-21
DE102008029175B4 (de) 2010-09-30
JP2010000358A (ja) 2010-01-07
US8120359B2 (en) 2012-02-21

Similar Documents

Publication Publication Date Title
CN101609132B (zh) 确定控制高频发送线圈的脉冲序列的方法和磁共振设备
CN101067649B (zh) 磁共振成像射频发射线圈的有效解耦
CN103250049B (zh) 用于测量固体和浆体含水量的低场nmr设备
EP1651974B1 (en) Multicoil nmr data acquisition and processing methods
CN103154761B (zh) 在并行发射mri中的虚拟线圈仿真
CN103649767B (zh) 具有多通道阻抗匹配网络的磁共振成像系统
CN102283649B (zh) 用于磁共振成像的三维空间选择性激励的序列设计方法
CN102159965B (zh) 用于mri的b1映射和b1l匀场
CN102640010B (zh) 用于产生均匀磁场的方法和装置
CN101055308B (zh) 用于控制磁共振系统的方法
CN101236237B (zh) 在磁共振断层摄影中设计射频激励脉冲的方法
CN101529268B (zh) 采用多个发射线圈的mri rf编码
CN100520435C (zh) 在磁共振设备的检查空间内产生均匀高频磁场的方法
Sbrizzi et al. Fast design of local N‐gram‐specific absorption rate–optimized radiofrequency pulses for parallel transmit systems
US8115482B2 (en) Magnetic resonance anatomical image generating method and system
CN103364747A (zh) 用于控制磁共振系统的方法和控制设备
CN104755951B (zh) 基于功率需要的并行rf发射中的rf放大器控制
CN102288929B (zh) 用于磁共振系统的激发脉冲序列生成器
CN102369452B (zh) 耦合天线阵列中的噪声匹配
Dar et al. Transfer learning for reconstruction of accelerated MRI acquisitions via neural networks
Zelinski et al. Designing RF pulses with optimal specific absorption rate (SAR) characteristics and exploring excitation fidelity, SAR and pulse duration tradeoffs
Blümler et al. Hardware developments: Single-sided magnets
Song et al. Computing transient electromagnetic responses of a metallic object using a spheroidal excitation approach
JP2013520263A (ja) 核磁気共鳴イメージング装置用高周波アンテナを制御するためのモジュールの作成方法
Zhang et al. Validation of transmit SENSE with reciprocity

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
TR01 Transfer of patent right

Effective date of registration: 20220111

Address after: Erlangen

Patentee after: Siemens Healthineers AG

Address before: Munich, Germany

Patentee before: SIEMENS AG