CN102007423A - 用于患者安全性和经改进的扫描性能的实时局部和全身sar估计 - Google Patents
用于患者安全性和经改进的扫描性能的实时局部和全身sar估计 Download PDFInfo
- Publication number
- CN102007423A CN102007423A CN2009801130917A CN200980113091A CN102007423A CN 102007423 A CN102007423 A CN 102007423A CN 2009801130917 A CN2009801130917 A CN 2009801130917A CN 200980113091 A CN200980113091 A CN 200980113091A CN 102007423 A CN102007423 A CN 102007423A
- Authority
- CN
- China
- Prior art keywords
- sar
- pulse
- absorption rate
- specific absorption
- matrix
- 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
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/58—Calibration of imaging systems, e.g. using test probes, Phantoms; Calibration objects or fiducial markers such as active or passive RF coils surrounding an MR active material
- G01R33/583—Calibration of signal excitation or detection systems, e.g. for optimal RF excitation power or frequency
-
- 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
- G01R33/5612—Parallel RF transmission, i.e. RF pulse transmission using a plurality of independent transmission channels
-
- 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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/288—Provisions 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
-
- 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/5659—Correction 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
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
在使能在高场强MR扫描中的经提高的RF占空比的一种方法或装置中,比吸收率(SAR)计算处理器(36)计算局部或全身SAR或者甚至空间SAR绘图。通过使用预平均的数据(基于E-场)的有效实施使得快速和精确的SAR估计成为可能。通过整合例如患者位置的附加信息,以及通过使用更多的患者特异性预计算信息(例如基于不同的生物网格),所谓的Q-矩阵,能够提高SAR计算精度。任选地,序列控制器(24)创建全身SAR最优RF脉冲。在应用所述最优RF脉冲之后,确定所述SAR以及其空间分布。还确定SAR热点。对绕所述热点的适当的半径内的Q-矩阵进行平均,并以加权的形式将其添加到全身Q-矩阵中。在更新所述全身Q-矩阵之后,创建了新的最优RF脉冲。能够重复所述步骤中的一个或多个,从而将所述SAR收敛为最小值。
Description
技术领域
本申请涉及诊断成像技术。其具体应用于患者安全性以及相关的经改进的扫描性能(在射频(RF)占空比方面)的领域,并且将具体参考其进行描述。此外,与高场强磁共振成像(MRI)有关的局部比吸收率(SAR)热点的估计和抑制。然而,应当认识到,其还可应用于对其他信息的优化和处理,并且不必受限于上述应用。
背景技术
对于许多较高场强的MR应用而言,局部SAR是限制性的因素。SAR的沉积随着较高的场强而增加并限制RF功率、占空比和可用的翻转角,从而导致扫描采集时间的延长以满足指定的SAR限值。对于单发送器系统而言,SAR相对容易计算,因为所有的天线元件以相同的幅度和其间固定的相移发送。此外,实验所需的RF脉冲的形状是已知的并将其连同它们的SAR存储于形状库中。随着其中每个线圈元件都具有潜力来独立发送其特有的幅度和相位的多发送系统的出现,基于每个通道计算SAR必须同时考虑并行RF发送脉冲,该并行RF发送脉冲仅仅能够基于例如B1绘图的附加信息进行计算,并且因此是实验/患者特异性的。
RF安全性是在活体内并行发送MRI扫描的先决条件,换言之,必须保证使用多通道RF发送线圈在SAR限值内进行扫描。除非扫描是“SAR安全的”,否则它们不能够开始。在具有多个发送通道的MR系统中,能够通过将电场信息整合到RF脉冲设计中来计算SAR降低的RF脉冲。在过去,已经使用了对于每个个体(例如,眼)通用的考虑到已知SAR热点来构建RF脉冲的方法。这通常不能满足于全身成像,因为SAR热点能够在患者之间、以及在RF脉冲之间在位置和幅度两方面有变化。因此,在一个患者体内将SAR限制为可接受水平的RF脉冲序列可能不能对另一患者也那样的限制。此外,适合已知静态热点的RF序列可能无意间加重在其他位置的未知的患者特异性的热点。
一种可能的解决方案是开发对所有患者安全的SAR的最差情况方案评估。然而,该解决方案将相当大地限制所允许的RF占空比,使得变得严重危及MRI系统在结合在活体内的并行发送扫描中的使用。针对患者定制SAR计算的能力将比使用一揽子方案或已知针对所有患者的条件更为有益。
当前RF序列不是基于逐个患者构建的一个具体的原因是对于临床相关空间RF脉冲(例如局部激发和放大成像)而言,需要能够加速这些类型的RF脉冲(TxSENSE)的并行发送系统。对于RF序列而言,基本的先决条件是对SAR进行有效估计的可获得性。此外,精确的SAR估计(包括全身和局部SAR值以及任选SAR绘图)高度需要患者相关的E-场和患者位置的可获得性。由仿真获得的场数据某种程度上不同于扫描仪中的实际场。用于E-场仿真而非实际患者的生物网格模型的使用导致难以特征化的系统误差。对于标准的单通道鸟笼线圈RF发送组件而言,RF波形对于每个Tx线圈元件是相同的,并且仅存在相位的增加(例如8个元件45°)。对于多个Tx线圈元件而言,该计算更为复杂,因为每个通道可能具有不同但固定的幅度和相位。在更为复杂的扫描中,诸如对于2D/3D空间选择性脉冲,每个通道可能具有动态改变的幅度和相位。
为了计算在标准中指定的所有SAR类型(局部和全身的)以及任选地计算针对多通道RF发送系统(例如,八个发送通道)的患者的SAR绘图,该系统执行大量的计算(例如TeraFLOPs):高达1010或更高的计算量,其取决于用于计算的模型和单元的分辨率。该过程将持续若干分钟并且实际不能够随着患者在扫描仪内等待实际诊断扫描开始而实时执行。
本申请提供了一种新的经改进的磁共振系统,其克服了以上涉及的问题和其他问题。
发明内容
根据一个方面,提供了一种磁共振系统。主磁体在检查区域生成大体上均匀的主磁场。射频组件在检查区域内的受检者的选定偶极子中诱发磁共振,并接收该磁共振。比吸收率计算处理器计算比吸收率并确定局部比吸收率热点。序列控制器设计RF激励脉冲,其负责局部比吸收率热点并保持传输给热点的能量处在可接受的水平之下。
根据另一方面,提供了一种磁共振系统。主磁体在检查区域生成大体上均匀的主磁场。射频组件在检查区域内的受检者的选定偶极子中诱发磁共振,并接收该磁共振。比吸收率计算处理器计算比吸收率并确定局部比吸收率热点。图形卡并行处理非图形信息。
根据另一方面,提供了一种磁共振的方法。在检查区域生成大体上均匀的主磁场。在检查区域内的受检者的选定偶极子中诱发磁共振,并接收该磁共振。确定在检查区域中的受检者的位置。计算比吸收率。计算负责所计算的比吸收率的全身安全RF脉冲波形。如果例如第一RF脉冲估计不满足SAR限值,或如TR(重复时间)的其他系统参数被延长,这可以迭代地完成。
一个优点在于有效验证并行发送扫描不违反现有FDA或国际电工委员会(International Electrotechnical Commission)的限制的能力。
另一优点在于对SAR值、热点和空间SAR分布的经提高的计算速度。
另一优点在于针对个体患者自定义SAR计算的能力。
另一优点在于基于患者的SAR图创建最优RF脉冲序列的能力。
另一优点在于基于患者与患者确定诸如E-场和患者位置的特定信息的能力。
另一优点在于检测用于SAR模型适应性调整的异物或外科移植物的能力。
本领域的普通技术人员通过阅读和理解下面的详细描述将认识到本发明的更进一步优点。
附图说明
本发明可以具体化为不同的部件或部件布置,以及具体化为不同的步骤和步骤安排。附图仅用于图示说明不同的方面,而不应解释为是对本发明的限制。
图1是根据本申请的磁共振成像装置的图解说明;
图2包括在RF激励期间用于PUC采样的示例性波形;
图3是热点抑制对经验性确定的权重因子的依赖性的曲线图;
图4是并未考虑所有可获得的信息的备选方法与考虑所有信息的实施例的比较;
图5是计算SAR的最差情况方案方法与考虑所有信息的实施例的比较。
具体实施方式
参考图1,磁共振扫描仪10被图示为包括螺线管主磁体组件12的闭合孔径系统,尽管也涵盖开放的和其他磁体构造。主磁体组件12产生沿成像区域的水平轴取向的大体上恒定的主磁场B0。应当理解,也涵盖诸如垂直的和其他构造的其他磁体布置。孔径型系统内的主磁体12通常可以具有大约0.5T到7.0T或更大的场强。
梯度线圈组件14在成像区域内产生磁场梯度,用于对主磁场进行空间编码。优选地,磁场梯度线圈组件14包括配置为在通常为纵向或z、横向或x和垂直方向或y的三个正交方向上产生磁场梯度的线圈段。
包括n个线圈元件161、162、…16n的射频线圈组件16生成用于在受检者的偶极子中激发共振的射频脉冲。射频线圈组件16发送的信号通常被称为B1场。射频线圈组件16还用于检测从成像区域中发射出的共振信号。所图示的射频线圈组件16是对整个成像区域成像的发送/接收线圈,然而,也涵盖了局部发送/接收线圈、局部专用接收线圈或专用发送线圈。在一个实施例中,射频线圈组件16包括8通道发送/接收天线。
梯度脉冲放大器18将受控的电流传输到磁场梯度组件14,以产生选定磁场梯度。包括n个发送器201、202、…20n的优选是数字式的射频发送器阵列20向射频线圈组件16施加射频脉冲或脉冲分组,以激发选定共振。在所图示的实施例中,线圈元件的数目和发送器的数目是相同的。然而,能够将多于一个线圈元件与每个发送通道相关联。在图示的实施例中,包括n个接收器221、222、…22n的射频接收器阵列22与线圈组件16或分离的接收线圈阵列耦合以接收并调制诱发的共振信号。
为了采集受检者的共振成像数据,受检者被放置在成像区域内。序列控制器24与梯度放大器18和射频发送器201、202、…20n通信以激发并操纵感兴趣区域内的磁共振。例如,序列控制器24产生选定的重复回波稳态或其他响应序列,对此类共振空间编码,选择性地操纵或扰乱共振,或以其他方式生成表征受检者的选定磁共振信号。所生成的共振信号通过RF线圈组件16或局部线圈组件(未示出)检测,传送给射频接收器22、被调制,并存储于k-空间存储器26中。成像数据通过重建处理器28重建,以产生一幅或多幅存储在图像存储器30中的图像表示。在一个合适的实施例中,重建处理器28执行逆傅里叶变换重建。
所生成的图像表示通过视频处理器32处理,并显示在配备有人可读显示器的用户接口34上。接口34优选是个人电脑或工作站。除了产生视频图像,图像表示可以由打印机驱动器处理和打印、通过计算机网络或因特网传送等等。优选地,用户接口34还允许技术人员或其他操作人员与序列控制器24通信,从而选择磁共振成像序列、修改成像序列、执行成像序列等等。在用户接口34,用户能够选择SAR模型,并且可以利用用户交互和反馈来确定所有或部分其余参数。
比吸收率(SAR)处理器36计算受检者在线圈组件16内的部分的SAR。在一个实施例中,SAR计算处理器36创建包括增加的SAR或热点的区域的全身的SAR绘图。对于其中标准RF脉冲的幅度和相位恒定变化的标准扫描而言,能够非常快地计算SAR,因为仅仅与相位/幅度关系有关。因此,单个RF样本的计算是足够的,因为幅度/相位关系不改变脉冲。可以通过使用恒定相位和幅度的常规MRI系统和并行发送系统实施标准的扫描。
假设受检者内的RF场线性响应于驱动该场的电流,SAR可以在脉冲样本中表述为二次形式,其中,表示共轭转置,b是RF波形样本,并且Q是由麦克斯韦(Maxwell)等式的解导出的并对应于特定的受检者体积的埃尔米特(Hermitian)正定矩阵。通过考虑包括轨线、B1场绘图、靶激发模式和全身Q矩阵的若干输入创建SAR绘图。现有的SAR优化算法通常仅约束特定的已知静态局部区域,诸如眼睛。正如先前所提到的,这对于全身成像是不充分的,因为可能存在随受检者不同而变化的其他热点。从统计学角度讲,约束热点出现的空间区域可能导致在其他位置出现新的热点。存储器35能够存储一个或多个患者位置的预计算的数据,从而防止当不必要时对Q-矩阵的重新计算。此外,存储器35能够存储SAR值,这样就不需要对同一脉冲重复计算SAR。唯一的ID可以用于识别脉冲。
对于局部SAR计算,患者模型的每个体积要素的SAR被平均直到达到了预期量。体积要素的SAR值是对沿体积要素的边缘的SAR的指示,并且该数据被内插,从而采集每个体素的中心处的SAR值。
某些信息可以预计算并存储在查找表(LUT)37中。诸如电场和B1场绘图的扫描仪特异性信息存储在LUT 37中。能够从LUT 37的身体模型存储器中通过获知患者的身高、体重、性别和在MR扫描仪中的位置来选择适当的开始生物网格。而身高、性别和体重能够由操作人员在扫描前输入,患者的位置由受检者位置处理器39确定。获得患者位置以及细化患者模型的一种方式是使用移动床的方法。在将患者向扫描仪的孔径内移动时采集图像,这导致低分辨率的3D体积数据集。备选地,一旦患者处于孔径内的最终位置时可以执行短预扫描。可以通过例如阈值分割或一些其他处理方法对该数据进行分割。随后,可以利用现有模型通过相关性方法例如从患者的横向切片或对界标的检测中获得患者的位置。同时,可以估计患者的体积和大小。还能够检测到诸如移植物或缺失器官的异常。以这种方式,针对当前患者自定义了开始身体模型。
一旦已经确定了患者的位置,SAR计算处理器36根据输入参数(体重、性别、位置,等等)使用其适当的组合查阅LUT 37以得到对应的E-场数据。如果与存储在LUT 37中的所有模型都偏离过大,那么可以使用非常保守的SAR估计。在具体对移植物扫描的情况下,能够从LUT 37中检索该设备的合适的SAR限值。在已经确定了初始位置之后,任何桌台运动可以通过受检者位置处理器39进行监测并能够用于精确地确定患者的新的位置。
在备选实施例中,粗略分割的身体模型可以从移动床成像数据或预扫描中获得,其能够用于对相似患者的现有E-场的快速适应性调整或对例如同类模型的快速估计。同类模型的使用在传导率和电容率方面引入相对小的误差。使用来自同类模型的数据和实际数据之间的差别是可忍受的,从而模型的使用是可行的备选。
在另一备选实施例中,患者位置可以通过使用拾取线圈(PUC)确定。多通道发送线圈的每个发送元件配备有用于监测在每个元件中的电流的PUC,以确保患者安全并有利于系统的调整。一般而言,患者的存在影响线圈的性能。因此,RF线圈元件的负载在患者通过磁体孔径的移动期间改变。可以将该移动检测为相位改变。可以将这种改变转换为患者在MR系统中的近似位置。这是可能的,因为在RF脉冲期间对线圈元件中的电流进行了采样,如图2所示。提供了示例性的RF和梯度波形。点线表示RF激发波形,虚线表示MR信号采样波形,而实线表示PUC采样波形。此外,PUC能够用于感测RF线圈通道中的异常电流,并且如果超过了安全参数就启动扫描终止。
一旦创建了SAR绘图,序列控制器24就设计针对当前受检者的SAR绘图定制的RF脉冲序列。这还可以通过计算RF脉冲的主重建器或分离的图形卡来完成。序列控制器24引入了指定不同热点区域和全身SAR之间的折衷的权重因子。例如,根据相对全身SAR进行优化的RF脉冲的空间SAR分布,经由Q1=Q全身+∑qiQ临界_区域(i)的热点降低是可能的,其中,Q1是经调整的Q-矩阵,Q全身是原始全身Q-矩阵,qi是权重因子,并且Q临 界_区域(i)是直接绕热点的体积(例如3x3x3的体素体积)的Q-矩阵。序列控制器24迭代地处理SAR绘图,以找到最佳的权重因子qi,从而满足现有的SAR限值并降低最大限度的SAR值。
更具体而言,序列控制器24引导梯度组件14和RF组件16以应用新近设计的RF脉冲序列。SAR计算处理器36随后重新计算SAR绘图。再次确定局部热点的位置。随后,序列控制器24体积平均热点的Q-矩阵(Q临 界_区域(i))并对其进行加权。权重因子已经基于热点与磁体的等中心之间的距离(z)经验性地确定。现在参考图3,曲线的波谷代表对该距离最优的权重因子。曲线40代表z=20cm的热点抑制。曲线42代表z=40cm的热点抑制。曲线44代表z=60cm的热点抑制。曲线46代表z=80cm的热点抑制。最后,曲线48代表z=100cm的热点抑制。
经加权、体积平均的Q-矩阵被添加到全身Q-矩阵。当重新计算经更新的Q-矩阵Q1时,绕每个热点空间平均的半径、热点位置、局部Q-矩阵和选定的权重因子都纳入考虑。一旦已经对Q-矩阵进行更新,序列控制器24基于经更新的Q-矩阵Q1设计新的SAR经优化的RF脉冲序列。此前,已经在临界区域获得了较低的SAR值。序列控制器24和SAR计算处理器36能够迭代地应用一个或多个以上步骤,直到SAR在热点收敛为最小值,或者备选地,直到达到了预期的安全SAR水平。在某些情况下,如果SAR在其收敛之前就达到了安全水平,可能没有必要应用迭代直到SAR收敛为最小值。备选地,还可以延长重复时间TR,或者能够降低翻转角,或者这两者的组合。同时,如果患者移动,可以重新优化RF脉冲。
该迭代过程是计算密集型的,每次更新全身Q-矩阵时需要大量数据处理能力。利用现有系统,每次迭代将持续若干分钟的时间,其对于在扫描仪中等待的患者而言是不切实际的。每个Q-矩阵计算负责每个相关通道的正确的幅度信息和正确的相位信息。在一个实施例中,分离地计算身体的每个体素,给予可能的最高分辨率。生物网格中体素大小为5mm的体素的平均数量处于750,000的量级。当针对作用于每个体素的每个RF通道处理相位和幅度信息时,需要大量的计算(例如TeraFLOPs)以计算全身和局部SAR以及产生SAR绘图。正如先前所提到的,一个实施例包括具有八个通道的RF组件16,但应当理解,具有更多通道的组件也是可能的,具有操作在任何给定时间的通道的任何组合。相应地针对这些情况计算SAR。
在图1的实施例中,SAR计算处理器36为诸如高性能图形卡的子处理器38分派任务。子处理器38可以位于SAR计算处理器36自身中、位于主计算机中、或位于分光计中。因为各个体素的SAR计算并不取决于其他体素,不必要对它们一个接一个的处理,换言之,能够对它们进行并行处理。诸如图形卡的子处理器38提供许多并行处理通道(例如128、256等)以加速SAR的计算。例如,通过使用具有128个处理通道的图形卡,SAR的计算比单独使用3GHz的处理器计算SAR加速了100倍。结果,针对单个生物网格的RF脉冲的SAR的计算在几秒内而不是几分钟内执行。这使得能够在实际可接受的时间内执行上述将SAR热点收敛为最小的迭代处理。在一个实施例中,如果子处理器38是不可获得的(例如,如果图形卡已损坏),那么SAR计算处理器36能够完成计算,这样扫描仍然是可能的。
在备选实施例中,能够按照体素的接近程度对它们进行分组并求平均,将体积要素的数目从大约750,000降低到例如100,000。这进一步降低了SAR的计算时间,但牺牲了所计算的SAR绘图中的分辨率和精度。结果,添加额外的安全余量以获得用于扫描的所估计的SAR值。
在另一备选实施例中,考虑幅度和通道信息,但不考虑相位信息。这也能够加速计算,但将计算较低精度的SAR绘图,这倾向于谨慎。在该实施例中过估计SAR值。
在另一备选实施例中,在每个对应的通道中将幅度设置为最大值。该方法再次削减计算量,因为仅仅考虑每个通道的最大幅度,但牺牲了所生成的计算的质量,这再次倾向于谨慎。
在另一备选实施例中,最差情况方案实施例,仅仅考虑最大幅度,而不管通道。这导致仅仅对实际SAR绘图的粗略估计。图4和图5图示了与考虑了所有可获得信息的实施例相比的一些并未考虑所有信息的备选实施例。在图4中,曲线50代表考虑正确幅度而不考虑相位的实施例。曲线52代表考虑在正确通道内的最大幅度的实施例。曲线54代表最差情况方案的实施例,其中,甚至不考虑通道信息。曲线50、52、54将针对实际SAR所计算的比率描绘为降低因子的函数。正如所见到的,随着更多的信息被考虑,SAR的估计值更接近实际的SAR。然而,如果这些误差比率是可接受的,能够通过使用备选方法中的一种来节省计算时间。
图5图示了使用与实际计算60对比的最差情况方案方法58执行的计算。图5是仿效标准的一通道体线圈、具有45°的相位和在所有通道上幅度为一的8-通道体线圈的局部躯体SAR的位置依赖性。显然,最差情况方案方法显著地过估计了SAR,尤其是在患者的上腹部,这导致较低精度的SAR计算。还显然的是SAR对位置的依赖性。
在另一备选实施例中,SAR计算处理器36、子处理器38或任何其他的部件可以位于远程服务器中。多个客户端可以同时由该服务器服务。当同时出现多个对SAR值的请求时,该服务器能够基于到达的顺序、或基于其他优先级别区分它们的优先顺序。
本文已经参考多个实施例对创新进行了描述。他人通过阅读和理解前面的详细描述,可以做出修改和变型。只要这些修改和变型在权利要求或其等同的范围内,意欲将本发明解释为包括所有这些修改和变型。
Claims (15)
1.一种磁共振系统,包括:
主磁体(12),其在检查区域中生成大体上均匀的主磁场;
包括多个线圈元件(161、162、…16n)的多通道射频线圈组件(16)和包括多个发送器(201、202、…20n)的发送器阵列(20),所述发送器阵列在所述检查区域中的受检者的选定偶极子中诱发磁共振;
接收阵列(22),其接收磁共振信号;
比吸收率计算处理器(36),其针对所述多通道射频线圈组件(16)计算比吸收率;以及
序列控制器(24),其重构针对所述局部比吸收率热点的RF激发脉冲。
2.如权利要求1所述的磁共振系统,还包括图形处理单元(38),其包括辅助处理所述比吸收率和能量吸收率热点的所述计算的多个并行处理通道。
3.如权利要求1所述的磁共振系统,还包括:
存储器(37),其包括针对电场、Q-矩阵和模型患者数据中的至少一个预确定的值。
4.如权利要求1所述的磁共振系统,还包括:
受检者位置计算处理器(39),其确定在所述主磁体(12)内的所述受检者的位置,并为所述比吸收率计算处理器(36)提供所述位置。
5.一种磁共振方法,包括:
在检查区域中生成大体上均匀的主磁场;
在所述检查区域中的受检者的选定偶极子中诱发磁共振,并接收所述磁共振;
确定所述受检者在所述检查区域中的位置;
计算比吸收率;以及
计算针对所计算的比吸收率的全身安全RF脉冲波形。
6.如权利要求5所述的方法,还包括:
确定至少一个局部比吸收率热点在所述受检者体内的位置。
7.如权利要求6所述的方法,还包括:
施加所述全身安全RF脉冲;以及
计算由所述全身安全RF脉冲在所述受检者体内诱发的空间比吸收率分布。
8.如权利要求6所述的方法,还包括:
通过添加所述至少一个热点的经加权体积平均的Q-矩阵来更新全身Q-矩阵。
9.如权利要求8所述的方法,其中,通过考虑绕所述至少一个热点空间平均的半径、所述至少一个热点的位置、至少一个局部Q-矩阵以及权重因子来计算所述经加权体积平均的Q-矩阵。
10.如权利要求8所述的方法,还包括:
根据经更新的全身Q-矩阵计算经优化的RF脉冲波形;
施加所述经优化的RF脉冲波形。
11.如权利要求10所述的方法,还包括:
至少重复更新所述全身Q-矩阵以及计算经优化的RF脉冲波形的所述步骤,直到所述比吸收率收敛为预期值为止。
12.如权利要求5所述的方法,其中,计算所述全身安全RF脉冲波形的所述步骤包括考虑所述RF脉冲的轨线、B1场绘图,靶激发模式和全身Q-矩阵。
13.如权利要求5所述的方法,还包括:
基于所述受检者的所述位置更新所述全身安全RF脉冲。
14.一种计算机可读介质,其存储有控制磁共振系统执行如权利要求5所述的方法的程序。
15.如权利要求5所述的方法,其中,计算比吸收率的所述步骤包括基于标准扫描和空间选择性扫描中的至少一种来计算所述比吸收率,所述标准扫描包括具有恒定幅度和相位的RF脉冲,所述空间选择性扫描使用用于患者安全性和热点抑制中至少一个的3D RF脉冲。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP08103561 | 2008-04-16 | ||
EP08103561.0 | 2008-04-16 | ||
PCT/IB2009/051531 WO2009128013A1 (en) | 2008-04-16 | 2009-04-13 | Real-time local and global sar estimation for patient safety and improved scanning performance |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102007423A true CN102007423A (zh) | 2011-04-06 |
CN102007423B CN102007423B (zh) | 2015-11-25 |
Family
ID=40765597
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200980113091.7A Active CN102007423B (zh) | 2008-04-16 | 2009-04-13 | 磁共振系统和方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8941380B2 (zh) |
EP (1) | EP2269085A1 (zh) |
JP (1) | JP5731373B2 (zh) |
CN (1) | CN102007423B (zh) |
WO (1) | WO2009128013A1 (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102736045A (zh) * | 2011-04-08 | 2012-10-17 | 西门子公司 | 具有局部sar 约束的并行发射射频脉冲设计 |
CN103703384A (zh) * | 2011-07-20 | 2014-04-02 | 皇家飞利浦有限公司 | 用于mri的感应耦合的局部并行发射线圈,其元件的每个都包括可变阻抗 |
CN103705237A (zh) * | 2012-09-28 | 2014-04-09 | 上海联影医疗科技有限公司 | 磁共振成像系统的保护方法、装置及自关断测试方法 |
CN105190342A (zh) * | 2013-03-13 | 2015-12-23 | 皇家飞利浦有限公司 | 用于磁共振成像的多元件rf发射线圈 |
CN105740534A (zh) * | 2016-01-28 | 2016-07-06 | 珠海和佳医疗设备股份有限公司 | 一种轴对称物理场的计算方法 |
CN105792749A (zh) * | 2013-12-18 | 2016-07-20 | 株式会社日立制作所 | 磁共振成像装置及其控制方法 |
CN110088642A (zh) * | 2016-09-30 | 2019-08-02 | 维宁尔美国公司 | 用于检测目标对象的前进方向和速度的系统和方法 |
CN107495967B (zh) * | 2017-08-24 | 2020-06-19 | 上海联影医疗科技有限公司 | 射频能量沉积预测及控制方法、装置、系统及存储介质 |
Families Citing this family (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2260318A1 (en) * | 2008-03-26 | 2010-12-15 | Koninklijke Philips Electronics N.V. | Determination of local sar in vivo and electrical conductivity mapping |
DE102009030721B4 (de) * | 2009-06-26 | 2013-04-04 | Siemens Aktiengesellschaft | SAR-Berechnung für mehrkanalige MR-Sendesysteme |
EP2296000B1 (en) * | 2009-09-10 | 2014-03-12 | Commissariat à l'Énergie Atomique et aux Énergies Alternatives | Parallel excitation of nuclear spins with local SAR control |
JP5542591B2 (ja) * | 2009-11-12 | 2014-07-09 | 株式会社東芝 | 磁気共鳴イメージング装置、および、磁気共鳴イメージング方法 |
DE102010011588B4 (de) * | 2010-03-12 | 2012-04-05 | Friedrich-Alexander-Universität Erlangen-Nürnberg | Verfahren zur Bestimmung von Sensitivitätsmatrizen für kritische Hotspots |
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 |
US9297870B2 (en) | 2010-05-27 | 2016-03-29 | Koninklijke Philips N.V. | Magnetic resonance examination with instrument detection |
US9689945B2 (en) | 2010-08-20 | 2017-06-27 | Koninklijke Philips Electronics N.V. | Virtual coil emulation in parallel transmission MRI |
EP2745133B1 (en) | 2011-08-17 | 2022-03-09 | Koninklijke Philips N.V. | Reducing the radio-frequency transmit field in a predetermined volume during magnetic resonance imaging |
WO2014034722A1 (ja) * | 2012-08-29 | 2014-03-06 | 株式会社東芝 | 磁気共鳴イメージング装置及びその比吸収率の演算方法 |
JP6395388B2 (ja) * | 2013-06-26 | 2018-09-26 | 株式会社日立製作所 | 磁気共鳴イメージング装置およびその制御方法 |
US10661091B2 (en) * | 2014-03-26 | 2020-05-26 | The General Hospital Corporation | System and method for hyperthermia treatment using radiofrequency phased arrays |
DE102014209783B4 (de) | 2014-05-22 | 2018-06-14 | Siemens Healthcare Gmbh | Erstellung eines MR-Bildes mit Hilfe einer Kombination aus einer Ganzkörperspule und einer lokalen Sendespule |
RU2700468C2 (ru) * | 2014-06-04 | 2019-09-17 | Конинклейке Филипс Н.В. | Удельная скорость поглощения, модулируемая пространственной близостью к пациенту |
DE102014211574A1 (de) * | 2014-06-17 | 2015-12-17 | Siemens Aktiengesellschaft | Anpassung der Tischposition bei der MR-Bildgebung |
US20160091583A1 (en) * | 2014-09-30 | 2016-03-31 | Siemens Aktiengesellschaft | Patient-Specific Estimation of Specific Absorption Rate |
US9939501B2 (en) * | 2014-11-12 | 2018-04-10 | The Board Of Trustees Of The Leland Stanford Junior University | Iterative minimization procedure with uncompressed local SAR estimate |
WO2016170177A1 (en) | 2015-04-24 | 2016-10-27 | Koninklijke Philips N.V. | A multi-channel transmit/receive radio frequency (rf) system |
DE102016200611A1 (de) * | 2016-01-19 | 2017-07-20 | Siemens Healthcare Gmbh | Lagerungsbewertung eines Untersuchungsobjekts |
JP2017164211A (ja) * | 2016-03-15 | 2017-09-21 | 株式会社日立製作所 | 磁気共鳴イメージング装置、q値算出方法及び比吸収率管理方法 |
US10444305B2 (en) * | 2017-06-21 | 2019-10-15 | Canon Medical Systems Corporation | Magnetic resonance imaging apparatus and specific absorption rate calculation method |
DE102018212091A1 (de) * | 2018-07-19 | 2020-01-23 | Siemens Healthcare Gmbh | Ermittlung einer belastungsoptimierten MR-Sequenz |
DE102021208408A1 (de) | 2021-08-03 | 2023-02-09 | Siemens Healthcare Gmbh | Verfahren zur Überwachung einer Absorption, Magnetresonanzvorrichtung und Computerprogrammprodukt |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005001502A1 (en) * | 2003-06-30 | 2005-01-06 | Koninklijke Philips Electronics N.V. | Control of specific absorption rate (asr) in mri |
US20070108976A1 (en) * | 2003-10-31 | 2007-05-17 | Koninkijke Hpilips Electronics Nv | B1 field control in magnetic resonance imaging |
US20070241753A1 (en) * | 2006-02-21 | 2007-10-18 | Beth Israel Deaconess Medical Center, Inc. | Magnetic resonance imaging and radio frequency impedance mapping methods and apparatus |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6249121B1 (en) * | 1999-05-17 | 2001-06-19 | General Electric Company | RF body coil |
DE10150138B4 (de) * | 2001-10-11 | 2009-10-08 | Siemens Ag | Verfahren zur Magnetresonanz-Bildgebung |
DE10150137B4 (de) * | 2001-10-11 | 2006-08-17 | Siemens Ag | Verfahren und Vorrichtung zur Magnetresonanz-Bildgebung |
US6989673B2 (en) * | 2003-11-26 | 2006-01-24 | General Electric Company | Method and apparatus to reduce RF power deposition during MR data acquisition |
DE102005019859A1 (de) * | 2004-05-11 | 2005-12-08 | Siemens Ag | Verfahren zur Durchführung einer Magnetresonanzuntersuchung eines Patienten |
CN101171761B (zh) * | 2005-04-29 | 2012-06-20 | 皇家飞利浦电子股份有限公司 | 在多通道发送/接收天线设备中用于功率监视的方法和电路装置 |
US8238678B2 (en) * | 2006-08-30 | 2012-08-07 | Siemens Medical Solutions Usa, Inc. | Providing representative image information |
DE102006052711A1 (de) * | 2006-11-08 | 2008-05-15 | Siemens Ag | Verfahren und Vorrichtung zur Oberflächenerfassung einer Untersuchungsperson |
DE102007013060A1 (de) * | 2007-03-19 | 2008-09-25 | Siemens Ag | Verfahren zur Einstellung eines B1-Felds |
-
2009
- 2009-04-13 EP EP09732227A patent/EP2269085A1/en not_active Ceased
- 2009-04-13 JP JP2011504587A patent/JP5731373B2/ja not_active Expired - Fee Related
- 2009-04-13 US US12/936,540 patent/US8941380B2/en active Active
- 2009-04-13 CN CN200980113091.7A patent/CN102007423B/zh active Active
- 2009-04-13 WO PCT/IB2009/051531 patent/WO2009128013A1/en active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005001502A1 (en) * | 2003-06-30 | 2005-01-06 | Koninklijke Philips Electronics N.V. | Control of specific absorption rate (asr) in mri |
US20070108976A1 (en) * | 2003-10-31 | 2007-05-17 | Koninkijke Hpilips Electronics Nv | B1 field control in magnetic resonance imaging |
US20070241753A1 (en) * | 2006-02-21 | 2007-10-18 | Beth Israel Deaconess Medical Center, Inc. | Magnetic resonance imaging and radio frequency impedance mapping methods and apparatus |
Non-Patent Citations (6)
Title |
---|
A.C.ZELINSKI ET AL: "Designing RF Pulses With Optimal Specific Absorption Rate (SAR) Characteristics And Exploring Excitation Fidelity, SAR And Pulse Duration Tradeoffs", 《PROC.INTL.SOC.MAG.RESON.MED.》 * |
A.C.ZELINSKI ET AL: "Pulse Design Methods For Reduction Of Specific Absorption Rate In Parallel RF Excitation", 《PROC.INTL.SOC.MAG.RESON.MED.》 * |
ANDREAS KEIL ET AL: "Patient Position Detection for SAR Optimization in Magnetic Resonance Imaging", 《MICCAI 2006》 * |
CORNELIS A.T.VAN DEN BERG ET AL: "Simultaneous B1+ Homogenization and Specific Absorption Rate Hotspot Suppression Using a Magnetic Resonance Phased Array Transmit Coil", 《MAGNETIC RESONANCE IN MEDICINE》 * |
D.O.BRUNNER ET AL: "Enforcing strict constraints in multiple-channel RF pulse optimization", <PROC.INTL.SOC.MAG.RESON.MED.> * |
YUDONG ZHU: "Parallel Excitation With an Array of Transmit Coils", 《MAGNETIC RESONANCE IN MEDICINE》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102736045A (zh) * | 2011-04-08 | 2012-10-17 | 西门子公司 | 具有局部sar 约束的并行发射射频脉冲设计 |
CN102736045B (zh) * | 2011-04-08 | 2014-09-03 | 西门子公司 | 具有局部sar 约束的并行发射射频脉冲设计 |
CN103703384A (zh) * | 2011-07-20 | 2014-04-02 | 皇家飞利浦有限公司 | 用于mri的感应耦合的局部并行发射线圈,其元件的每个都包括可变阻抗 |
CN103703384B (zh) * | 2011-07-20 | 2016-12-14 | 皇家飞利浦有限公司 | 用于mri的感应耦合的局部并行发射线圈,其元件的每个都包括可变阻抗 |
CN103705237A (zh) * | 2012-09-28 | 2014-04-09 | 上海联影医疗科技有限公司 | 磁共振成像系统的保护方法、装置及自关断测试方法 |
CN105190342A (zh) * | 2013-03-13 | 2015-12-23 | 皇家飞利浦有限公司 | 用于磁共振成像的多元件rf发射线圈 |
CN105190342B (zh) * | 2013-03-13 | 2018-01-30 | 皇家飞利浦有限公司 | 用于磁共振成像的多元件rf发射线圈 |
CN105792749A (zh) * | 2013-12-18 | 2016-07-20 | 株式会社日立制作所 | 磁共振成像装置及其控制方法 |
CN105740534A (zh) * | 2016-01-28 | 2016-07-06 | 珠海和佳医疗设备股份有限公司 | 一种轴对称物理场的计算方法 |
CN110088642A (zh) * | 2016-09-30 | 2019-08-02 | 维宁尔美国公司 | 用于检测目标对象的前进方向和速度的系统和方法 |
CN110088642B (zh) * | 2016-09-30 | 2022-11-29 | 安致尔软件有限责任公司 | 用于检测目标对象的前进方向和速度的系统和方法 |
CN107495967B (zh) * | 2017-08-24 | 2020-06-19 | 上海联影医疗科技有限公司 | 射频能量沉积预测及控制方法、装置、系统及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
US20110043205A1 (en) | 2011-02-24 |
EP2269085A1 (en) | 2011-01-05 |
US8941380B2 (en) | 2015-01-27 |
WO2009128013A1 (en) | 2009-10-22 |
JP2011517983A (ja) | 2011-06-23 |
CN102007423B (zh) | 2015-11-25 |
JP5731373B2 (ja) | 2015-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102007423A (zh) | 用于患者安全性和经改进的扫描性能的实时局部和全身sar估计 | |
CN102736045B (zh) | 具有局部sar 约束的并行发射射频脉冲设计 | |
CN103430037B (zh) | 将针对mri的成像区域限制在不均匀的磁场中的医疗装置及方法 | |
US10145929B2 (en) | Method and magnetic resonance apparatus for generating a parameter map to compensate local inhomogeneities in the basic magnetic field | |
EP3160339B1 (en) | Positioning of a magnetic resonance imaging antenna within the homogeneous field zone | |
CN102743171A (zh) | 用于在mr成像中使用的患者支撑台控制系统 | |
CN101858965B (zh) | 确定k空间位置的方法以及执行该方法的磁共振设备 | |
CN106918794A (zh) | 磁共振系统及成像的方法 | |
US20160245882A1 (en) | Method and magnetic resonance apparatus for determining absolute receive sensitivity maps for reception coils | |
KR101687177B1 (ko) | 자기 공명 시스템 제어 | |
US20200367765A1 (en) | Method and system for measuring blood flow | |
CN105308473A (zh) | 平行成像加速参数的自动优化 | |
CN103782184A (zh) | Rf阵列线圈/天线阻抗的动态修改 | |
KR20200006003A (ko) | 자기 공명 이미징을 위한 펄스 시퀀스들의 데이터베이스를 구축하는 컴퓨터 구현 방법 및 이러한 데이터베이스를 이용하여 자기 공명 이미징을 수행하는 방법 | |
US10162033B2 (en) | Magnetic resonance imaging method and apparatus | |
CN102483450A (zh) | 沿弯曲辐条k空间轨线的RF匀场的MRI切片激励 | |
US11428767B2 (en) | Computer-implemented method for providing an actuation sequence and method for acquiring measurement data by a magnetic resonance device | |
US7956611B2 (en) | Magnetic resonance imaging apparatus and method | |
US11733325B2 (en) | Method for controlling an MR apparatus | |
US11754651B2 (en) | Patient-model-based determination of a stimulation of a magnetic resonance imaging | |
EP2741097A1 (en) | Method and apparatus for acquiring B1 magnetic field phase information | |
CN111090067B (zh) | 用于在磁共振成像中对k空间数据进行采样的方法和系统 | |
US20180011160A1 (en) | Reduction of eddy currents during flow encoded magnetic resonance imaging | |
EP3832336B1 (en) | A method of designing a pulse sequence for parallel-transmission mri, and a method of performing parallel-transmission mri using such a pulse sequence | |
Hardy et al. | Bench to bore ramifications of inter-subject head differences on RF shimming and specific absorption rates at 7T |
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 |