CN105143904A - 用于多信道发射机的多频带射频/磁共振成像脉冲设计 - Google Patents
用于多信道发射机的多频带射频/磁共振成像脉冲设计 Download PDFInfo
- Publication number
- CN105143904A CN105143904A CN201480017259.5A CN201480017259A CN105143904A CN 105143904 A CN105143904 A CN 105143904A CN 201480017259 A CN201480017259 A CN 201480017259A CN 105143904 A CN105143904 A CN 105143904A
- Authority
- CN
- China
- Prior art keywords
- pulse
- channel
- band
- modulation
- phase
- 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
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/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
-
- 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/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
- G01R33/4833—NMR 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
- G01R33/4835—NMR 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 of multiple slices
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
- A61B5/0042—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
-
- 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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/32—Excitation or detection systems, e.g. using radio frequency signals
- G01R33/36—Electrical details, e.g. matching or coupling of the coil to the receiver
-
- 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
-
- 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/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/246—Spatial mapping of the RF magnetic field B1
-
- 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/32—Excitation or detection systems, e.g. using radio frequency signals
- G01R33/36—Electrical details, e.g. matching or coupling of the coil to the receiver
- G01R33/3607—RF waveform generators, e.g. frequency generators, amplitude-, frequency- or phase modulators or shifters, pulse programmers, digital to analog converters for the RF signal, means for filtering or attenuating of the RF signal
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Optics & Photonics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Neurology (AREA)
- Signal Processing (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
提供了一种用于磁共振成像(“MRI”)***的、用于设计多信道、多频带射频(“RF”)脉冲的方法。该方法包括:通过对包含复数值向量的原函数进行最小化,来为多信道RF线圈中的每一个信道确定多个RF幅度调制和多个RF相位调制。原函数还包含***矩阵,所述***矩阵说明多信道RF线圈中每一个信道的空间敏感度轮廓以及多频带RF脉冲中每一个激发频带的磁场图。
Description
相关申请的交叉引用
本申请要求以下两件美国临时专利申请的权益:2013年1月25日提交的第61/756,775号、题为“SYSTEMANDMETHODFORMULTIBANDRADIOFREQUENCYMAGNETICRESONANCEIMAGINGPULSEDESIGNFORUSEWITHAMULTICHANNELTRANSMITTER(用于多信道发射机的多频带射频磁共振成像脉冲设计的***和方法)”以及2013年12月11日提交的第61/914,775号、题为“SYSTEMANDMETHODFORMULTIBANDRADIOFREQUENCYMAGNETICRESONANCEIMAGINGPULSEDESIGNFORUSEWITHAMULTICHANNELTRANSMITTER(用于多信道发射机的多频带射频磁共振成像脉冲设计的***和方法)”。
关于联邦赞助研究的声明
本发明通过政府支持、依据国家健康协会获奖的RR008079、EB015894和MH091657而作出。政府对本发明有特定的权利。
发明背景
发明领域是用于磁共振成像(“MRI”)的***和方法。具体而言,本发明涉及用于设计、优化和填补用于MRI中的多频带射频(“RF”)脉冲的***和方法。
MRI中的高及超高磁场(诸如3特斯拉(3T)和7特斯拉(7T)的场)的最近可用性已经能以不断增加的空间分辨率来进行人脑的解剖学及功能性MRI(“fMRI”)研究,且在fMRI应用的情况下,能对神经元活性升高的站点启用保真度。然而,即使在采用单脉冲切片采集法、诸如回波成像(“EPI”)时,在整个人脑上获得这种高分辨率fMRI数据会遇到长空间重复时间(“TRs”)的不期望结果。已经使用如下方法来解决上述问题:即,通过使用多切片、二维EPI策略的并行成像原理,进行多切片的同时多频带(“MB”)RF激发,其中对同时采集的切片进行后续去混叠。该方法允许空间TR(volumeTR)直接降低一个因子(称为MB因子),该因子等于同时激发的切片的数目。
MB方法的使用已经在基于任务的fMRI和休息状态MRI中取得了显著成功,导致休息状态网络的改进检测以及揭露这种网络的临时动态的新分析策略。采用该方法还显著降低在扩散成像技术中通常会很长的成像采集时间,所述扩散成像技术诸如高角扩散加权成像(“HARDI”)和扩散频谱成像(“DSI”)。如此,MBRF脉冲的使用已经使这些成像技术的使用对于像人类连接体项目(HumanConnectomeProject)这样的工作来说变得切实可行和不可或缺。在这些应用中,去混叠过程所容忍的可实现MB因子可以通过使用EPI中的梯度标志将同时采集的切片沿相位编码维度部分地移位而显著改进。注意到,MB方法可与除EPI以外的空间编码策略一同使用,诸如一次(例如,闪光FLASH)一个k空间行地采集的梯度回叫回波、快速自旋回波等等。此外,MB方法可应用于除fMRI以外的成像应用,一般包括其中基本上同时采集多个切片的任何成像应用。
然而,尽管有这些收获,但是高(3T和4T)磁场和超高(7T及更高)磁场下的切片加速多频带方法的优选使用被发射B1(“B1+”)不同质性和功率沉积约束条件所妨碍。信噪比(“SNR”)和图像对比度变得空间不均匀,且在一些位置处由于不均匀B1+而变得次优,所述不均匀B1+由于在这些场强下存在行波行为而产生的破坏性干扰的结果从而产生。这些B1+不同质性已在4T处特别是7T处很好地记录,但在3T处足够强以引起大脑中央和***部分之间的SNR差异,特别是在基于自旋回波(“SE”)的序列中。类似地,特别是在超高磁场下和/或在采用SE序列时,最大可实现的切片加速度因子可受功率沉积所限制。当切片数量和空间TR被保持相同时,多频带方法与常规的单切片激发相比不沉积更多功率,即使峰值功率会随MB因子呈二次地递增;然而,通过MB因子加速导致功率沉积中的MB倍增加,这对可实现的加速度产生限制。
这些限制在临床诊断中人体躯干及四肢的成像中也是最重要的。特别是,使用诸如涡轮自旋回波(也称为快速自旋回波及相关衍生物)等序列的成像在用多频带方法实现时会经受功率沉积和B1+不同质性,其中所述序列是大脑、躯干及四肢成像中大量临床扫描的基础。
因此,会期望提供一种用于为多信道传输应用、设计和提供多频带RF脉冲的方法,其中所设计的多频带RF脉冲与当前可用的多信道、多频带RF技术相比已经减少了B1+不同质性且降低了功率沉积(包括降低的全局SAR、降低的峰值局部SAR和/或降低的峰值RF功率)。
发明概述
本发明通过提供一种指示磁共振成像(“MRI”)***来生成多个多频带射频(“RF”)脉冲的方法,克服了上述缺陷。该方法包括预计算阶段,在此阶段选择对要被每一个多频带RF脉冲所激发的切片数量进行定义的多频带因子,确定要设计的多频带RF脉冲的数量,以及确定要基于所选择的多频带因子和所确定的多频带RF脉冲数目而激发的多个板。该方法包括通过为多频带RF线圈中的每一个信道确定RF幅度调制数量和RF相位调制数量来设计多频带RF脉冲。这些值通过使包括复数值向量的原函数最小化来确定,所述原函数包含每一个脉冲的RF脉冲幅度调制和RF脉冲相位调制。所述原函数还包含***矩阵,所述***矩阵说明多信道RF线圈中每一个信道的空间敏感度轮廓、多频带RF脉冲中每一个激发频带的磁场图、以及每一个多频带RF脉冲的磁场梯度数量。MRI***针对使用通过这一最小化确定的RF幅度调制和RF相位调制来产生多频带RF脉冲。
本发明另一方面是提供一种指示磁共振成像(“MRI”)***来产生多频带射频(“RF”)脉冲的方法。该方法包括通过为多频带RF线圈中的每一个信道确定RF幅度调制和RF相位调制来设计多频带RF脉冲。这些值通过使包括复数值向量的原函数最小化来确定,所述原函数包含RF脉冲幅度调制和RF相位调制。所述原函数还包含***矩阵,所述***矩阵说明多信道RF线圈中每一个信道的空间敏感度轮廓以及多频带RF脉冲中每一个激发频带的磁场图。MRI***针对使用通过这一最小化确定的RF幅度调制和RF相位调制来产生多频带RF脉冲。
本发明的以上及其他方面和优点将从以下说明书中显而易见。在说明书中,参考形成说明书一部分的附图,其中通过图示示出本发明的优选实施例。然而,这种实施例不必要表示本发明的完全范围,因此参考权利要求书且在此解释本发明的范围。
附图简述
图1A是用于激发两个板(Nslab=2)的一个多频带(“MB”)射频(“RF”)脉冲(K=1)的板指示的示例;
图1B示出以下板指示的示例:用于激发三个板(Nslab=3)的一个MBRF脉冲(K=1)、用于激发四个板(Nslab=4)的两个MBRF脉冲(K=2)、用于激发六个板(Nslab=6)的两个MBRF脉冲(K=2)、用于激发六个板(Nslab=6)的三个MBRF脉冲(K=3)、以及用于激发九个板(Nslab=9)的三个MBRF脉冲(K=3);
图2A和2B示出用于激发四个板(Nslab=4)的两个MBRF脉冲(K=2)的不同板指示,其中板位置未空间交织;
图3A和3B描述了板内的切片位置的示例;
图4是磁共振成像(“MRI”)***的示例的框图;以及
图5是射频(“RF”)***的示例的框图,所述射频***形成图4的MRI***的一部分并且包括多信道发射机和接收机。
发明的详细描述
提供了一种用于多信道发射(“pTx”)多频带(“MB”)脉冲设计的方法,所述方法解决了发射B1(“B1+”)不同质性和RF功率沉积的上述问题。所述方法能够产生pTxMB脉冲设计,该设计能同时在多个离散切片中提供显著改进的B1+同质化而无需使用相对于单信道圆极化(“CP”)模式应用较高的RF功率。脉冲设计也能够实现功率沉积的节约。因此,本发明能够改进B1同质性而得到改进的SNR和对比噪声比(“CNR”),而同时减少目标功率沉积,这可以包括全局SAR、峰值局部SAR和/或峰值RF功率。
通常,在多信道MB激发中可以采用两种策略来查找基础RF脉冲的射频(“RF”)大小和相位调制(即,RF填充值)以便减轻B1+不同质性。一种策略,在此称为频带联合设计,是为全部M个单频带(“SB”)脉冲获得一组公共的信道专用RF填充值,所述全部M个单频带脉冲被组合以同时激发M个切片。该解决方案可以使用携带MB脉冲形状的单个信道发射机来应用,所述所述MB脉冲形状然后被分成多个不同信道并且经受信道专用的相位和幅度变化。在该方法中,MB脉冲的M个SB脉冲中的每一个SB脉冲经历相同的信道专用相位和幅度变化。另一种策略,在此称为频带专用设计,是为多信道发射机的每一个信道计算用于M个SB脉冲的每一个SB脉冲的一组不同的填充值。以下将依次讨论这些方法的每一种。
在针对均匀激发(即,均匀|B1+|)的频带联合方法中,M-频带Q-信道激发的RF填充值可以通过求解以下幅值最小平方问题来获得:
其中
以及
在公式(1)-(3)中,A是组合的复数值矩阵,其中Am(m=1,2,...,M)是包括第m个频带的个别信道的B1+空间敏感度轮廓和ΔB0图的***矩阵;w是复数值向量,其中wq(q=1,2,...,Q)是第q个信道的RF填充值;d是表示期望的横向磁化的标量;以及λ是正则化参数。
在频带专用的方法中,可以求解:
其中
以及
在公式(4)-(6)中,Afull是块对角矩阵,其中Am是对角元素矩阵,而wfull是级联向量,其中wm(m=1,2,...,M)是第m个频带的个别信道的RF填充值的向量。
值得注意,构成最终加总MB脉冲的M个基础RF脉冲随着引入的时间而具有不同的相位演化(即,不同斜率的线性相位坡度)以便定向同时激发的M个切片的每一个切片的频移。因此,根据Parseval的理论,给定时间段内的功率沉积由所采集的切片数量来确定,无论切片是使用MB方法还是使用常规的单切片激发过程来采集。结果,项的使用是对于频带专用pTxMB脉冲设计中的总RF功率的最优约束条件。类似地,以下项可用作全局SAR的最优约束条件:
其中S0是全局SAR矩阵。峰值局部SAR管理也可以在频带专用pTxMB脉冲设计中使用将最终加总脉冲形状结合到约束条件中的约束条件、并同时优化一个以上参数(诸如全局SAR和峰值局部SAR)来实现。
还值得注意,在频带专用pTxMB脉冲设计中,最终RF脉冲的应用要求对个别信道的独立控制RF脉冲。这是因为每信道的最终RF脉冲是M个不同基础脉冲分别与不同权重相乘后的加总,因此不能仅仅用信道专用权重w′q乘以一个公共脉冲来表示。
例如,MBRF脉冲可以被设计为单辐条(即,在横向平面中没有梯度编码),或者被设计为更复杂的构造(例如,多辐条)以便在所关注空间(volume)中同时激发期望的、具有|B1+|同质化的谐振频带。这一脉冲设计的统一目标可以被定义为期望切片中的同质|B1+|。该目标可以通过手动创建空间掩码来产生,该空间掩码仅覆盖频带中所关注空间的部分。RF幅值和相位调制可以用可变交换算法来计算,该可变交换算法具有在脉冲计算时未考虑的所关注空间之外的体元。目标也可以是统一(即,均匀激发)之外的其他目标,并且可以被裁剪以实现B1+在所关注区域上的具体的、预定义的空间分布。例如,目标可以可以是聚焦或“缩放”以便在一个子区域或多个子区域中实现高B1+以及在其他子区域中实现低B1+。
设计有单辐条的RF脉冲也可以通过将RF幅值按比例增加以达180度翻转角而同时维持相同级别的翻转角同质化而用于重新聚焦。然而,当多辐条被期望用于进一步改进翻转角同质性时,应当考虑大尖角脉冲设计方法(诸如最优控制方法)来维持激发保真度。
以此方式设计的MBRF脉冲可以被构造成一次用多个切片来覆盖所关注空间,诸如人的头部或躯干的分段,其中每组切片自身的优化取决于在所关注空间中操纵的切片的具***置。MB脉冲的上述设计也可以结合静态磁场B0不同质性以得到更好的脉冲性能。
在多频带成像方法中,多个切片被同时(或基本同时)激发,随后一起采集,且最终被去混叠以产生对应于个别切片的图像。该过程用一组不同的切片进行重复以覆盖所关注空间(“VOI”)。例如,如果期望分辨率为1mm各向同性,则切片“厚度”为1mm。因此,如果VOI在切片方向具有120mm的尺寸,则会需要120个切片、以1mm各向同性分辨率来覆盖VOI。
在常规的多切片方法中,每一个切片分别在二维采集中被激发和成像,例如使用读出梯度和相位编码梯度来编码两个正交的面方向。因此,在常规的多切片方法中,对于上述示例会需要120个分开的激发和数据采集。对于多频带方法,使用例如MB因子4,则可以仅使用30个激发和后续数据采集来覆盖相同的VOI。在该情况下,此处描述的脉冲优化会对于30组四离散切片而执行,而不考虑沿切片厚度的B1和B0的潜在变化,因为它很小。
替代的“混合”方法会是激发非1mm切片,而是较厚的板,诸如10mm厚的板。该实例中大于1的MB因子会对应于基本上同时激发多个板。例如,对于MB因子4,会使用三组四板激发来覆盖上例中120mm厚的VOI。在该情况下,为了获得1mm各向同性分辨率,会使用三维采集。在这一采集中,沿板方向添加相位编码维度。在上述示例中,会需要十个这样的相位编码来获得沿10mm厚板的名义1mm分辨率。
在混合方法中,使用三维傅立叶变换来将板沿板方向以及沿其他两个正交方向解析为1mm的“切片”。然后,同时或基本同时激发的“板”使用并行成像原理彼此去混叠。在该采集方案中,用于“板”选择性脉冲的脉冲优化可以如上所述实现。任选地,沿板维度的B1或B0变化可以被忽略,它在板的位置处平均优化,如同采集个别薄切片的情况一样。然而,通过使用混合成像方法,优化也可以考虑B1的不同质性,或者B1和B0一起的不同质性,所述不同质性不仅在面方向上(如同薄切片的情况),也可以沿板的厚度。
使用pTx脉冲设计进行MB激发可以实现以下好处。首先,B1+不同质性的影响可以最小化。其次,可以实现其他约束条件,诸如最小化局部或全局功率沉积,所述局部或全局功率沉积可以被测量为比吸收率(“SAR”)的减少。在一示例中,这些优点可以通过利用具有仅平面中梯度标志(gradientblips)的所谓的3D辐条切片选择性RF脉冲来提供。该方法允许平面中的B1同质化,同时允许锐利切片轮廓的生成。而且,全局和局部SAR的管理以及平均和峰值RF功率也可以被结合至脉冲设计中,以便符合技术和安全约束条件两者。
当获得大的所关注空间(“VOI”)的MB数据时,诸如整个脑部,即使在沿被成像的对象的最短形态维度进行采集时,也需要激发和采集许多切片。例如,为了覆盖整个人脑,包括大脑、小脑和脑干,会要求总共104个切片来填充130mm的典型维度,沿切片维度具有1.25mm的空间分辨率。对于这么多切片,对于医疗应用来说,为每一个个别成像切片的个别发射信道获得B1+轮廓以及执行后续大规模优化问题以使用上述按切片过程来设计多个pTxMBRF脉冲所需要的时间会不切实际地长。
例如,具有用于104个切片的典型脉冲设计规范(例如,约5mm的平面中分辨率和用于脉冲优化的单辐条)的16-信道脉冲设计会在脉冲设计问题中得到166400×1664个矩阵。求解这一大优化问题是有挑战的,特别是在功率沉积包括附加约束条件时。查找这一情况下令人满意的解决方案所需要的计算时间以当前可用的计算资源来说可能过长。
而且,根据MB脉冲序列的现有技术实现方式的当前状态,由于有限的最大允许存储器分配,因此作为切片专用的脉冲设计的结果而维持大量不同的pTxSB和pTxMB脉冲可能是困难的。例如,应用切片专用的pTxMB脉冲来以MB因子4(即,M=4)采集104个切片可能要求该序列创建和准备总共130个pTx脉冲(104个pTxSB脉冲加上26个pTxMB脉冲)用于RF激发。如果手头的脉冲序列是自旋回波序列,则附加的用于重新聚焦的130个脉冲会是必要的。这一数目的pTx脉冲在现有技术MRI***的当前状态下会难以维持,现有技术MRI***可能受限于它们对于最大数量RF脉冲(诸如最大数量128个RF脉冲)被允许的存储器分配。
即使对RF脉冲的最大数量没有限制,应用大量切片专用的pTx脉冲可能使得时间效率较差,因为脉冲序列准备时间与要应用的pTx脉冲的数量成比例。注意到,较高分辨率的使用会导致大量的切片,大量切片进而可以产生进一步的困难。因此,以逐切片的方式对于整个脑部覆盖范围执行pTxMB脉冲设计可能造成时间成本和计算要求成本非常高,并且实际上在用当前仪器的医疗时间约束条件内可能不可行。
为了解决这些问题,上述脉冲设计方法可适于提供一种用于覆盖大VOI(诸如整个脑部覆盖范围)的pTxMBRF脉冲设计的有效、高效且实际的策略。以下讨论为说明目的使用整个脑部的示例,但可用于任何大VOI覆盖范围,如可应用于躯干成像或整个人体成像。
在上述方法的此适配中,对于要为数据采集的切片不实施任何逐切片的pTx脉冲设计。相反,对于为覆盖整个VOI而规定的少量板设计出构成pTxMB脉冲的pTxSB脉冲。板的数量由所期望的pTxMB脉冲的数量来确定,其中MB因子有较小的限制。此外,每一个板的脉冲设计可以通过仅考虑来自所述板内的少量有间隙切片的B1+信息来确定。
使用上述方法的此适配将显著减少脉冲设计问题的规模。例如,考虑具有四个切片的四个板,每一个都用于设计一个16-信道、单辐条pTxMB4脉冲(即,MB因子M=4)将会使传输矩阵的尺寸减小为25600×64个,这与使用相同平面中分辨率的逐切片设计时166400×1664个的尺寸相比有169倍的减少。
设计问题的尺寸的显著减少将显著降低找到令人满意的脉冲解决方案所需的计算时间。此外,在B1+映射中仅使用16个切片以及在序列中仅使用五个pTx脉冲(4个pTxSB脉冲加上1个pTxMB脉冲)会大大减少扫描时间。
上述方法描述了在同时的多切片多频带采集中使用pTxRF脉冲以便解决在非常高的场下的增加的B1+不同质性和RF功率沉积的问题。该pTxMB脉冲设计可以在人脑中同时在多个切片中提供显著改进的B1+均匀性(即,同质性),而不需要使用相对于单信道圆极化(“CP”)模式应用更高的RF功率,以及/或者实现功率沉积的节省。
然而,如上所述,该方法对于整个脑部覆盖范围需要长采集时间和计算时间。尽管这种计算限制在将来随着不断增长的计算功率和网络速度会毫无疑问地减轻,且用于映射发射B1场的采集速度可以增加,但是当前的实际实现方式可通过适配上述方法来实现。
在上述pTx脉冲设计方法的适配中,可以设计单个或多个pTxMB脉冲来在整个VOI上在每一个个别采集切片内创建均匀的激发。脉冲设计不是基于从要在最终MB采集中采集的每一个切片获得的B1+信息;相反,它是基于仅从VOI的独立于采集切片的少量稀疏分布的切片获得的B1+信息。该设计方法的实现方式包括两个一般阶段:预计算阶段和计算阶段,以下分别讨论。
预计算阶段的主要目的是规定少量相邻的成像板,所述成像板包含与要设计的个别pTxMB脉冲相对应的VOI。然后,预计算阶段定义每一个板内的切片位置,从所述切片位置将采集B1+图用于后续的脉冲设计。为此,提供了两步骤的过程。
第一步骤是确定要设计的pTxMB脉冲的数量K,所述数量K与MB因子M一起得出需要规定的板的数量Nslab,,
Nslab=K×M(8)。
板的数量Nslab应当被视为K个组,每一个组具有与不同的pTxMB脉冲相对应的M个板。在数据采集期间,每一个pTxMB脉冲将如同在常规MBRF激发中那样被应用,以便通过一次激发M个切片来连续地激发驻留在其M个板中的全部成像切片,且K个pTxMB脉冲将被依次展开以激发覆盖整个VOI的全部Nslab个板。
图1A和1B图示在从整个脑部采集同时的多切片多频带数据时用于pTxMB2和MB3(即,MB因子M=2和3)RF脉冲设计的示例的板指示,其中轴向切片在头到脚趾的方向上跨过所述脑部。如图所示,对于给定的MB因子,可以考虑不同数量的板,这取决于要设计多少pTxMB脉冲。较大数量的pTxMB脉冲将得到较大数量的板。注意到,板厚度仅需要对于同一MB脉冲的M个板相同,而不需要在不同的pTxMB脉冲之间相同。
还应当注意,图1A和1B中图示的全部板指示以交织方式布置,其中最大可能的切片间距离跨全部MB脉冲是恒定的,以便实现后续MB重构中的全局最优切片去混叠。然而,板位置通常不需要被交织并且可以为任何任意次序。例如,在指示用于设计两个pTxMB2脉冲的四个板时,一个pTxMB脉冲的两个板可以被另一个MB脉冲的两个板夹在中间,或者可以堆叠在另一个MB脉冲的两个板顶部,如图2A和2B分别示出。注意到在前一情况下,切片间距离最好在两个MB脉冲之间不同,无论使用什么板厚度。另一方面,在堆叠板的后一情况下,如果同一板厚度被指派给全部四个板,则切片间距离可以对于两个MB脉冲相同。
预计算阶段中的第二步骤是设计个别板内的B1+切片位置。该过程包括确定每一个板的切片数量以及切片间的距离。图3A和3B描述了一个板内的切片位置的两个示例。在图3A所示的示例中,总共14个切片位于已被指示用于设计两个pTxMB3脉冲的六个板内部,其中切片间距离恒定且边界切片被两个相邻板共享。在图3B所示的示例中,不同数量的切片以随机方式被置于六个板内。应当注意,最优切片位置可以通过基于电磁模拟或试验采集的B1+图沿着切片选择性方向分析B1+变化来获得。
重要的是要注意:随着要设计的pTxMB脉冲的数量增加,要被指示的板的数量将会增加且平均板厚度会增加。当pTxMB脉冲的数量达到采集切片总数除以MB因子时,按板的脉冲设计会返回至上述的按切片的设计。例如,按板方式设计26个pTxMB4脉冲来采集104个切片实际上是按切片的设计,因为需要被指示的104个板将减少至与104个采集切片相一致的104个切片。这也意味着上述的按切片的脉冲设计可以被视为更一般化的按板设计的特殊情况。
注意到,板可以从所关注对象的模型得出,所关注对象可以是人脑。例如,板可以从诸如单个人类头部模型和/或被导出为许多不同受试者的复合或平均的模型等模型中得出。后一模型可以表示人类头部的共同平均RF属性。通过结合这样的模型数据,可能获得一些预设计的pTxRF脉冲,这些脉冲说明受试者间的变异并且可用于不同的扫描,而对于B1+映射或脉冲计算没有进一步的个别扫描。
一旦完成了预计算阶段,就可以实施计算阶段来获得可减少B1+不同质性的RF填充值。这里,目的是设计彼此不同的K个pTxMB脉冲。首先,提供对RF功率和SAR具有明确约束条件的脉冲设计问题的广义公式。然后,提供用于设计个别pTxMB脉冲的两个策略并且概括了小尖角区中的相应设计问题。
广义化的脉冲设计问题可以如下写出:
minwf0(w)(9);
符合
(1)f1(w)≤gSAR(10);
(2)f2(w)≤lSAR(11);
(3)f3(w)≤Pave(12);
(4)f4(w)≤Ppeak(13)。
这里,原函数f0一般表示激发误差,未知数w是将个别pTxMB辐条脉冲的RF填充值级联的复数值向量。关于w的四个约束条件函数f1-4表示所关心的四个最相关的RF功率相关量,并且以不等式形式使用以便对以下施加明确约束条件:(1)全局SAR、(2)局部SAR、(3)总RF功率、以及(4)瞬时峰值RF功率,其中,gSAR、lSAR、Pave和Ppeak是相应的SAR和功率限制。求解这一受约束的最小化问题导致尝试找到会提供最低可实现激发误差(即,最佳可实现激发保真度)的一组RF填充值,而同时确保所得到的RF功率和SAR数量在相应的规定限制内。
对于小尖角脉冲设计,原函数f0通常被列公式为平方误差或幅值平方误差以便对所得激发模式距规定目标模式的偏差进行量化。这里A是包括相对于相应信道的输入RF功率归一化的多信道B1场分布的复数值***矩阵。***矩阵A也具有在规定的全部切片内的主要场不同质性,以及按梯度标志的平面中相位调制(对应于辐条的k空间位置)。参数d是表示期望的激发模式的复数值向量(例如,它可以是给定翻转角的均匀激发的单一性值)。而且,***矩阵A具有块对角的特征结构并且可以被构造为,
其中K是要设计的pTxMB脉冲的数量,对角元素矩阵Ak是对于第k个pTxMB脉冲设计的***矩阵。相应地,具有未知数的向量w和目标的向量d分别由以下给出:
以及
原函数中以或的形式使用A、w和d的这种定义,与脉冲设计的目标相一致,即为要设计的K个pTxMB脉冲的每一个找到一组不同的RF填充值。然而,个别***矩阵Akk=1,2,...,K的结构以及个别RF填充向量wk高度依赖于如何为每一个pTxMB脉冲计算个别的RF填充组wk。如以下将呈现的,对于wk的计算存在利用不同级别的自由度的两种可能的设计策略,每一个策略对于脉冲设计问题的公式化都具有一个不同的处方。
对于每一个个别pTxMB脉冲的设计可以考虑两种设计策略。一种策略在此定义为频带专用设计,其用于为组成pTxMB脉冲的M个pTxSB脉冲的每一个找到不同的RF填充组。另一种策略在此称为频带联合设计,其用于为构成pTxMB脉冲的全部M个pTxSB脉冲找到公共RF填充组。以下,首先提供用于频带专用设计的设计问题,此后提供用于频带联合设计的设计问题。
考虑受约束的幅值最小平方最小化,从用于频带专用设计的广义化脉冲设计公式化中衍生的脉冲设计问题可以写成:
符合,
这里,第k个pTxMB脉冲的***矩阵Ak和RF填充向量wk被定义为,
和
注意到矩阵Ak也被构建为块对角,该结构连同wk的相应定义反映出在为对应于给定MB脉冲的M个pTxSB脉冲的每一个找到不同RF填充组的频带专用设计中的目标。矩阵Ak的个别对角元素矩阵以及向量wk的个别元素向量进一步由下式给出,
Ak,m=[Ak,m,1,Ak,m,2,...,Ak,m,C]以及(23);
其中
Ak,m,c=[Ak,m,c,1,Ak,m,c,2,...,Ak,m,c,S](25);
且
c=1,2,...,C;s=1,2,...,S,(28);
其中m、c和s分别对频带、发射信道和辐条计数;矩阵Ak,m,c,s,l(l=1,2,...,L(k,m))包括在辐条s处从所应用的平面中梯度标志得到的空间依存的相位累加、以及在已经为设计第k个pTxMB脉冲的第m个频带而预定义的全部L(k,m)个切片的第l个切片内的信道c的B1+场图,而w(k,m,c,s)是表示MB脉冲k、频带m、信道c和辐条s的填充值的复数值。注意到确定在搜索最佳解决方案时可用的自由度的未知数的总数因此是(K×M×C×S),由总填充向量w的维度(即,向量元素的数量)给出。
相应地,目标向量dk写成,
其中
其中dk,m,l是表示切片l、频带m和pTxMB脉冲k的目标激发轮廓的向量。应当注意,上述的矩阵Ak,m以及向量wk,m和dk,m在常规的pTx辐条脉冲设计中使用以便基于单个切片(L=1)或者同时基于多个切片(L>1)来计算单个多信道多辐条RF填充组。
在四个RF功率和SAR相关约束条件的公式化中,w(k,m,s)是施加于MB脉冲k、辐条s和频带m的全部信道上的填充值的垂直级联(即,w(k,m,s)=[w(k,m,1,s),w(k,m,2,s),...,w(k,m,C,s)]T),而w(k,c,s)是仅施加于辐条s的信道c上、但用于全部频带的填充值的垂直级联(即,w(k,c,s)=[w(k,1,c,s),w(k,2,c,s),...,w(k,M,c,s)]T);Q是全局SAR矩阵,而{Sv}v=1,2,...,V表示为估计局部SAR(例如,10克SAR)得出的一组V个虚拟观察点(VOP);应当注意,这些SAR矩阵已结合了SAR计算中独立于RF脉冲的全部因子,包括组织特性和空间平均。这里,Nk是MB脉冲k被应用于MB数据采集中的次数(例如,在MB单次EPI中,将给出所采集的切片总数),T是用于RF功率平均的时间窗(例如,它可以是采集中使用的重复时间),Δt是RF脉冲的停留时间,而b(t)是用N个时间点定义的基础脉冲波形(例如,常规的sinc脉冲)。矩阵M是由M个常规的未填充SBRF脉冲组成的矩阵,am(t),m=1,2,...,M,并且可以被构造为
其中am(t)=b(t)·exp[-iγG·(m-1)z0·t],其中γ是旋磁比,G是切片选择性梯度的最大幅度,而z0是SB脉冲间的切片间距离。
重要的是注意到矩阵向量乘积Mw(k,c,s)得到在施加于信道c上的辐条s处的pTxMB脉冲k的最终加总的MB脉冲形状,其中其最大绝对值被平方(即,)得到相应的瞬时峰值RF功率,且因此,满足在全部K个pTxMB脉冲、全部C个信道和全部S个辐条上的第四不等式从而确保了全部所设计的pTxMB脉冲的瞬时峰值RF功率将在RF功率放大器所允许的最大峰值功率内。此外,脉冲设计形式可以有效地求解,例如经由内点方法。
还重要的是注意到频带专用pTxMB脉冲的应用依赖于具有完全独立的多个发射信道的pTx***(例如,独立的RF合成器、独立的RF放大器),因为信道专用的RF波形彼此完全不同。
类似地,频带联合设计的脉冲设计问题可以被列公式为:
符合
由于为对应于给定MB脉冲的全部M个pTxSB脉冲计算一组公共RF填充,因此每一个个别***矩阵Ak现在被构造为多个较低级别矩阵Ak,m的堆叠并且由下式给出,
其中个别的矩阵Ak,m被定义有与频带专用设计中相同的层级结构。相应地,每一个个别向量wk现在分解至m个维度并且由下式给出,
其中
这里w(k,c,s),s=1,2,...,S是表示在信道c上和辐条s处施加的MB脉冲k的全部M个SB脉冲的填充值的复数。目标向量dk具有与频带专用设计中相同的定义。注意到在该情况下,通过总填充向量w的维度确定的自由度仅为(K×C×S),与其中自由度由(K×M×C×S)给出的频带专用设计相比有M倍的减少。由于自由度的这种退化,因此在实现激发模式并同时满足所有约束条件时,联合设计可能不如分开的设计有效。
总体上,频带联合设计中的四个RF功率和SAR约束条件的公式与频带专用设计中的那些相当类似,大多数参数都相同。然而这里引入了两个新的参数。一个是向量w(k,s),它是施加于全部M个频带的辐条s的、MB脉冲k的所有信道的填充值的垂直级联(即,w(k,s)=[w(k,1,s),w(k,2,s),...,w(k,C,s)]T),另一个是标量A0,它是M个SB基础脉冲的加总的最大幅度(即,)。
尽管相似,相较于频带专用设计可以标识两个差异。第一个差异存在于全局SAR、局部SAR和总RF功率的前三个约束条件中。在频带联合脉冲设计中,SB脉冲的相应的量(由m索引)的加总减少至与数量M的简单乘法,因为全部M个SB脉冲对于给定的MB脉冲、给定的信道和给定的辐条现在共享一个相同的填充值。这也导致第二差异,该第二差异可以在施加于瞬时峰值RF功率上的最后一个约束条件中所见。对于频带联合脉冲设计方法,在该约束条件中,给定MB脉冲、给定信道和给定辐条的峰值功率现在可以简单地通过两个标量的乘积来提供、而不是如频带专用设计中包括M个基础脉冲形状的矩阵向量乘法来提供。
不像那些频带专用的pTxMB脉冲,从频带联合设计获得的pTxMB脉冲不必要依赖于完全pTx能力并且可以仅通过在低功率阶段把扫描仪的单个RF信道分成多个信道的较简单且较便宜的硬件来应用,并且向所得到的分开信道施加信道专用的相位和幅度变化。
对于频带专用和频带联合脉冲设计方法两者,可以进一步优化经由平面中梯度标志的辐条的k空间位置,以便实现更好的RF性能;然而,这种优化最好按MB脉冲基础来执行,这与在设计用于顺序多切片激发的pTx脉冲时的按切片基础相反。此外,尽管以上在仅有四个约束条件的受约束最小化的上下文中描述,上述用于pTxMBRF脉冲设计的方法也可以被公式化为任何其他可应用的优化,诸如正则化的最小化。脉冲设计方法也可以包括任何其他期望的约束条件,诸如约束每信道MB脉冲包络的变化率的约束条件。
上述用于MBRF脉冲设计的方法也可用于设计时间移位的MBRF脉冲,诸如那些在共同待批专利申请PCT/US2012/58541中描述的RF脉冲,该专利申请通过引用完全结合于此;或者用于设计具有较低功率沉积的MBRF脉冲,诸如那些在共同待批专利申请PCT/US2011/066081中描述的RF脉冲,该专利申请通过引用完全结合于此;或者用于设计使用例如优化控制方法的大尖角MBRF脉冲;或者用于设计可用于控制混叠条件的RF脉冲,诸如那些可在具有改进的成像性能的CAIPIRINHA(“较高加速度的并行成像结果中的受控混叠”)技术中采用的RF脉冲。
通常,对于有M倍加速的MB试验,可以使用M个单频带RF波形来生成最终的MBRF波形。在最简单的情况下,为每一个频带使用相同的基础RF形状(例如,sinc波形),但是通过调制波形的相位而对每一个频带应用不同的频率偏移。该频率偏移定义了频带间隔。计算M个单频带波形的复数和,以生成一个多频带波形。可以使用单频带波形来采集基准扫描数据,且多频带波形可用于加速的数据采集。在一些实施例中,对于多切片和多频带采集,这些M+1个RF波形以N个不同的发射机频率偏移依次应用(主发射频率由扫描仪动态调节)以采集总共N×M个切片。对于pTxMB脉冲的实现方式,有几种可能的方式来扩展该基本方案以利用多个发射信道的优势,如下所述。
在一示例中,使用激发k空间中的单个辐条轨迹来设计单个pTxMB脉冲。在该示例中,对于每一个频带使用相同的基础脉冲形状(例如,sinc波形),以如上所述相同的方式对每一个频带应用不同的频率偏移。然后在M个单频带RF波形的每一个的顶部应用发射信道专用的全局幅值和相位。计算这些波形的复数和以产生发射信道专用的多频带波形。现在对于每一个发射信道有M+1个波形,即总共T×(M+1)个波形,T是发射信道的数量。注意到可以并行地计算每组M+1个波形,这在实践中节约时间。同样,对于多切片和多频带采集,这些波形以N个不同的发射机频率偏移依次应用以采集总共N×M个切片。
在另一示例中,可以使用激发k空间中的单个辐条轨迹来设计多个不同的pTxMB脉冲。在这种情况下,取代仅计算T×(M+1)个RF波形并且以不同的发射机频率将它们重复使用N次,可以确定和使用多达N×T组发射信道专用幅值和相位以计算N×T×(M+1)个波形来激发N×M个切片。实践中,由于每一个发射信道的N组单频带波形仅在全局幅值和相位方面不同,因此仅需要保留一个组,而可以使扫描仪动态地修改主发射频率、RF幅值和RF相位来重复其余组。因此,仅有必要为每一个发射信道加载M+N个唯一波形,即总共T×(M+N)个波形。这在考虑可用于在扫描仪上存***形的有限存储器时是重要的。
上述方案可以在使用激发k空间中的多个辐条轨迹来设计pTxMB脉冲时容易地扩展,该情况中可能将整个多辐条RF形状细分成几个分段以形成个别的“辐条”脉冲,对于每一个辐条具有发射信道专用的幅值和相位。此外,也可能对于不同的激发频带或切片位置不使用不同的脉冲波形形状。
可能影响脉冲设计策略的选择的实际考虑因素包括计算速度以及可用于存储脉冲波形的***存储器。尽管可能为每一个信道、为每一个频带、为每一个切片位置以及为每一个k空间辐条等等计算唯一的波形,但是可以在标准医疗扫描仪上的序列中预计算和使用的唯一波形的数量可能受限制(例如,每发射信道128个)。在一些情况下,计算和上载波形也可能是相对慢的过程。对于实际的实现方式,可能优选的是最小化唯一波形的数量,这是以上方案被设计成要做的事。
现在特别参考图4,图示磁共振成像(“MRI”)***400的示例。MRI***400包括:操作者工作站402,其一般会包括显示器404;一个或多个输入设备406,诸如键盘和鼠标;以及处理器408。处理器408可以包括运行商业可用操作***的商业可用可编程机器。操作者工作站402提供操作者界面,该操作者界面使扫描指示被输入至MRI***400中。通常,操作者工作站402可以耦合至四个服务器:脉冲序列服务器410;数据采集服务器412;数据处理服务器414;以及数据存储服务器416。操作者工作站402和每一个服务器410、412、414和416相连以便彼此通信。例如,服务器410、412、414和416可以经由通信***440相连,通信***440可以包括任何适当的网络连接,无论有线、无线还是有线和无线的组合。例如,通信***440可以包括专属或专用网络以及开放式网络,诸如因特网。
脉冲序列服务器410响应于从操作者工作站402下载的指令而工作,以操作梯度***418和射频(“RF”)***420。执行所指示扫描必须的梯度波形被产生并且被应用于梯度***418,梯度***418激发组件422中的梯度线圈以产生用于对磁共振信号进行位置编码的磁场梯度Gx、Gy和Gz。梯度线圈组件422形成磁性组件424的一部分,磁性组件424包括极化磁铁426和全身RF线圈428。全身RF线圈428可以包括能实现本发明方法的单信道发射线圈或多信道发射线圈阵列。例如,全身RF线圈428可以包括双信道发射线圈阵列。
RF波形由RF***420应用于RF线圈428或者分开的局部线圈(图4未示出),以便执行所指示的磁共振脉冲序列。RF线圈428或者分开的局部线圈(图4未示出)所检测到的响应性磁共振信号被RF***420接收,在RF***420中,在脉冲序列服务器410所产生的命令的指示下对上述信号进行放大、解调、滤波和数字化。RF***420包括用于产生MRI脉冲序列中所使用的各种RF脉冲的一个或多个RF发射机。RF发射机响应于来自脉冲序列服务器410的扫描指示和指引,以产生具有期望频率、相位和脉冲幅度波形的RF脉冲。所产生的RF脉冲可以被应用于全身RF线圈428或者被应用于一个或多个局部线圈或线圈阵列(图4未示出)。
RF***420也包括一个或多个RF接收机信道。每一个RF接收机信道包括RF预放大器和检测器,所述RF预放大器对其所连接的线圈428接收到的磁共振信号进行放大,所述检测器检测并数字化接收到的磁共振信号的I和Q正交分量。因此,接收到的磁共振信号的幅值可以在任一采样点通过I和Q分量的平方之和的平方根来确定:
且接收到的磁共振信号的相位也可以根据以下关系来确定:
脉冲序列服务器410还任选地接收来自生理采集控制器430的病人数据。例如,生理采集控制器430可以接收来自连接至病人的多个不同传感器的信号,诸如来自电极的心电图(“ECG”)信号、或者来自呼吸风箱或其他呼吸监控设备的呼吸信号。这种信号一般由脉冲序列服务器410用来将扫描的性能与受试者的心跳或呼吸同步,或“门控”。
脉冲序列服务器410还连接至扫描室接口电路432,扫描室接口电路432从与病人和磁性***的状态相关联的各种传感器接收信号。病人定位***434还通过扫描室接口电路432接收命令以在扫描期间将病人移至期望位置。
RF***420所产生的数字化磁共振信号样本被数据采集服务器412接收。数据采集服务器412响应于从操作者工作站402下载的指令而工作,以接收实时磁共振数据并且提供缓冲存储器,使得不会因数据溢出而丢失任何数据。在一些扫描中,数据采集服务器412基本上仅将所采集的磁共振数据传递至数据处理器服务器414。然而,在需要自所采集的磁共振数据得到的信息来控制扫描的进一步性能的扫描中,数据采集服务器412被编程以生成这种信息并将其传送至脉冲序列服务器410。例如,在预扫描期间,采集和使用磁共振数据来校准脉冲序列服务器410所执行的脉冲序列。举另一个例子,可以采集和使用导航信号来调节RF***420或梯度***418的运行参数,或者控制k空间被采样的视图次序。在还有另一示例中,数据采集服务器412也可用于处理用来检测磁共振血管造影(“MRA”)扫描中造影剂的出现的磁共振信号。例如,数据采集服务器412采集磁共振数据并且实时处理它以产生用于控制扫描的信息。
数据处理服务器414接收来自数据采集服务器412的磁共振数据并且按照从操作者工作站402下载的指令来处理所述磁共振数据。这种处理可以例如包括以下的一或多者:通过执行原始k空间数据的傅立叶变换来重构二维或三维图像;执行其他图像重构算法,诸如迭代或反投影重构算法;向原始k空间数据或向经重构的图像应用过滤器;生成功能性磁共振图像;计算运动或流动图像;等等。
数据处理服务器414所重构的图像被传送回操作者工作站402并保存。实时图像被存储于数据库存储器高速缓存(图4未示出)中,并且它们可以从中被输出至操作者显示器412或位于磁性组件424附近的显示器436以便供参与的医师使用。批量模式图像或所选择的实时图像被存储于磁盘存储器438上的主机数据库中。当这种图像已经被重构和传输至存储器时,数据处理服务器414就通知操作者工作站402上的数据存储服务器416。操作者工作站402可由操作者用来存档图像、产生影像、或者经由网络将图像发送至其他设施。
MRI***400还可以包括一个或多个联网工作站442。例如,联网工作站442可以包括:显示器444;一个或多个输入设备446,诸如键盘和鼠标;以及处理器448。联网工作站442可以与操作者工作站402位于同一设施中,或者位于不同设施中,诸如不同的医疗机构或诊所。
联网工作站442无论与操作者工作站402位于同一设施还是不同设施中都可以经由通信***440获得对数据处理服务器414或数据存储服务器416的远程访问。相应地,多个联网工作站442可具有对数据处理服务器414和数据存储服务器416的访问。以此方式,可以在数据处理服务器414或数据存储服务器416与联网工作站442之间交换磁共振数据、经重构图像或其他数据,以使所述数据或图像可被联网工作站442远程地处理。该数据可以以任何适当格式被交换,诸如按照传输控制协议(“TCP”)、因特网协议(“IP”)或者其他已知的或适当的协议。
如图4所示,RF***420可以连接至全身RF线圈428,或者如图5所示,RF***420的发射部分可以连接至RF线圈阵列504的一个或多个发射信道502,而RF***420的接收机部分可以连接至RF线圈阵列504的一个或多个接收机信道506。发射信道502和接收机信道506通过一个或多个发射/接收(“T/R”)转换器508连接至RF线圈阵列504。在RF***428的替代配置中,其中接收线圈是与发射线圈分开的线圈集合,此配置中不需要且不使用T/R转换器508。取而代之,在这一配置中,接收阵列在发射期间被“失谐”以使它不耦合至发射机。同样,在接收期间,发射机被失谐。以此方式,发射路径和接收路径不发生混合。
特别参照图5,RF***420包括产生规定的RF激发场的一个或多个发射信道502。该RF激发场的基础频率、或载波频率在频率合成器510的控制下产生,频率合成器510接收来自脉冲序列服务器410的一组数字信号。这些数字信号指示在输出端512处产生的RF载波信号的频率和相位。RF载波被应用于调制器和上变频器514,在那里,RF载波的幅度响应于也从脉冲序列服务器410接收到的信号R(t)而被调制。信号R(t)定义了要产生的RF激发脉冲的包络,并且通过依次读出一系列所存储的数字值而产生。这些所存储的数字值可以变化以使任何期望的RF脉冲包络得以产生。
在输出端516处产生的RF激发脉冲的幅值可被励磁机衰减器电路518所衰减,励磁机衰减器电路518接收来自脉冲序列服务器410的数字命令。然后,经衰减的RF激发脉冲被应用于功率放大器520,功率放大器520驱动RF线圈阵列504。
受试者所产生的MR信号由RF线圈阵列502拾取并被应用于该组接收机信道506的输入端。每一个接收机信道506中的预放大器522放大信号,信号然后被接收机衰减器524衰减一个量,该量由接收自脉冲序列服务器410的数字衰减信号所确定。接收的信号处于或约为Larmor频率,这一高频信号由下变频器526在两步骤过程中下变频。下变频器526首先将MR信号与线路512上的载波信号混合,然后将所得差异信号与线路528上的基准信号混合,所述基准信号由基准频率发生器530所产生。经下变频的MR信号被应用于模数(“A/D”)转换器532的输入端,模数转换器532对模拟信号进行采样和数字化。作为高频信号的下变频的一种替代方式,接收到的模拟信号也可以用适当快速的模数(“A/D”)转换器和/或用适当欠采样直接检测。然后,经采样和数字化的信号被应用于数字检测器和信号处理器534,数字检测器和信号处理器534产生与接收信号相对应的16位同相(I)值和16位正交(Q)值。接收信号的经数字化的I和Q值的所得流被输出至数据采集服务器412。除了在线路528上生成基准信号以外,基准频率发生器530还在线路536上生成被应用于A/D转换器532的采样信号。
本发明已经就一个或多个优选实施例进行描述,应当理解,除了那些明确阐述的实施例以外,许多等价物、替代、变化和修改都是可能的并且落入本发明的范围内。
Claims (17)
1.一种指示磁共振成像(MRI)***来生成多频带射频(RF)脉冲的方法,所述方法的步骤包括:
a)通过使原函数最小化来确定多信道RF线圈中的每一个信道的RF幅度调制和RF相位调制,所述原函数包括:
包含所述RF幅度调制和RF相位调制的复数值向量;以及
***矩阵,所述***矩阵说明所述多信道RF线圈中每一个信道的空间敏感度轮廓以及所述多频带RF脉冲中每一个激发频带的磁场图;以及
b)指示所述MRI***使用在步骤a)中确定的所述RF幅度调制和RF相位调制来生成至少一个多频带RF脉冲。
2.如权利要求1所述的方法,其特征在于,步骤a)包括:为所述多信道RF线圈中的每一个信道确定不同的RF幅度调制和RF相位调制。
3.如权利要求1所述的方法,其特征在于,步骤b)包括:在将所述多频带RF脉冲波形分成所述多信道RF线圈的多个不同信道之前、使用单信道发射机来发射多频带RF脉冲波形,此时每个多频带RF脉冲波形都经受由步骤a)中确定的所述RF幅度调制和RF相位调制所确定的信道专用的幅度和相位变化。
4.如权利要求1所述的方法,其特征在于,步骤b)包括:独立地控制所述多信道RF线圈中的每一个信道以便使用在步骤a)中为该信道确定的RF幅度调制和RF相位调制来生成多频带RF脉冲。
5.如权利要求1所述的方法,其特征在于,步骤a)包括:
a)i)选择对要被每一个多频带RF脉冲激发的切片数量进行定义的多频带因子;
a)ii)确定要设计的多频带RF脉冲的数量;
a)iii)基于所选择的多频带因子和所确定的多频带RF脉冲数量来定义要激发的多个板;以及
其中步骤b)包括:生成多个多频带RF脉冲、指示所述MRI***使用步骤a)中确定的所述RF幅度调制和RF相位调制来生成多个多频带RF脉冲,其中所述多个多频带RF脉冲激发步骤a)iii)中定义的所述多个板。
6.如权利要求1所述的方法,其特征在于,步骤a)包括:使用k空间中的单个辐条轨迹来确定每一个RF幅度调制和每一个RF相位调制。
7.如权利要求1所述的方法,其特征在于,步骤a)包括:为k空间中的多个不同的辐条轨迹的每一个确定RF幅度调制和RF相位调制。
8.一种指示磁共振成像(MRI)***生成多个多频带射频(RF)脉冲的方法,所述方法的步骤包括:
a)选择对要被每一个多频带RF脉冲激发的切片数量进行定义的多频带因子;
b)确定要设计的多频带RF脉冲的数量;
c)基于所选择的多频带因子和所确定的多频带RF脉冲数量来定义要激发的多个板;
d)通过使原函数最小化来确定多信道RF线圈中的每一个信道的RF幅度调制数量和RF相位调制数量,所述原函数包括:
包含每一个多频带RF脉冲的所述RF幅度调制和RF相位调制的复数值向量;以及
***矩阵,所述***矩阵说明所述多信道RF线圈中每一个信道的空间敏感度轮廓、每一个多频带RF脉冲中每一个激发频带的磁场图、以及每一个多频带RF脉冲的磁场梯度数量;以及
e)指示所述MRI***使用在步骤d)中确定的所述RF幅度调制和RF相位调制来生成多个多频带RF脉冲,其中所述多个多频带RF脉冲激发步骤c)中定义的多个板。
9.如权利要求8所述的方法,其特征在于,步骤c)包括:定义每一个板的空间位置以及每一个板要被激发的次序。
10.如权利要求9所述的方法,其特征在于,步骤c)包括:将所述次序定义为交织的次序。
11.如权利要求9所述的方法,其特征在于,步骤c)还包括:对位于每一个板内的多个切片的每一个定义空间位置。
12.如权利要求10所述的方法,其特征在于,所述多个切片的空间位置相等地分布在每一个板中。
13.一种使用磁共振成像(MRI)***同时激发受试者中的多个不同切片位置的方法,所述方法的步骤包括:
a)计算要应用于多信道RF线圈中每一个信道的多频带射频(RF)脉冲波形;
b)通过使原函数最小化,同时对多信道发射机中的每一个信道计算发射信道专用的RF幅度调制以及对多信道发射机中的每一个信道计算发射信道专用的RF相位调制,所述原函数包括:
包含每一个发射信道的所述发射信道专用的RF幅度调制和所述发射信道专用的RF相位调制的复数值向量;以及
***矩阵,所述***矩阵说明所述多信道RF线圈中每一个信道的空间敏感度轮廓、每一个多频带RF脉冲中每一个激发频带的磁场图、以及每一个多频带RF脉冲的磁场梯度;
c)指示所述MRI***通过为每一个发射信道将给定发射信道的发射信道幅度和发射信道相位应用于所述给定发射信道的所计算的多频带RF脉冲波形来生成RF激发场,从而激发所述受试者中多个不同切片位置中的自旋。
14.如权利要求13所述的方法,其特征在于,每一个发射信道专用的RF幅度调制是相同的全局RF幅度调制,且每一个发射信道专用的相位调制是相同的全局RF相位调制。
15.如权利要求13所述的方法,其特征在于,每一个多频带RF脉冲的磁场梯度定义了相同的k空间轨迹。
16.如权利要求15所述的方法,其特征在于,所述相同的k空间轨迹是辐条轨迹。
17.如权利要求13所述的方法,其特征在于,每一个多频带RF脉冲的磁场梯度定义了不同的k空间轨迹。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201361756775P | 2013-01-25 | 2013-01-25 | |
US61/756,775 | 2013-01-25 | ||
US201361914775P | 2013-12-11 | 2013-12-11 | |
US61/914,775 | 2013-12-11 | ||
PCT/US2014/013004 WO2014116986A1 (en) | 2013-01-25 | 2014-01-24 | Multiband rf/mri pulse design for multichannel transmitter |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105143904A true CN105143904A (zh) | 2015-12-09 |
CN105143904B CN105143904B (zh) | 2018-11-09 |
Family
ID=51228074
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201480017259.5A Active CN105143904B (zh) | 2013-01-25 | 2014-01-24 | 用于多信道发射机的多频带射频/磁共振成像脉冲设计 |
Country Status (4)
Country | Link |
---|---|
US (1) | US10684337B2 (zh) |
KR (1) | KR20150110685A (zh) |
CN (1) | CN105143904B (zh) |
WO (1) | WO2014116986A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109923427A (zh) * | 2016-11-09 | 2019-06-21 | Cr发展公司 | 对样本执行扩散加权磁共振测量的方法 |
Families Citing this family (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20150110685A (ko) * | 2013-01-25 | 2015-10-02 | 리전츠 오브 더 유니버스티 오브 미네소타 | 다중 채널 송신기를 위한 다중 대역 rf/mri 펄스 디자인 |
US9869737B2 (en) * | 2013-03-04 | 2018-01-16 | Siemens Healthcare Gmbh | Methods and systems for accelerated MR imaging |
DE102013205785B4 (de) * | 2013-04-02 | 2014-10-09 | Siemens Aktiengesellschaft | Ermittlung einer Magnetresonanzsystem-Ansteuersequenz auf Basis einer reduzierten Anzahl an Feldverteilungskarten |
CN105334479B (zh) * | 2014-06-10 | 2019-02-22 | 西门子(深圳)磁共振有限公司 | 一种磁共振成像方法和装置 |
DE102014219784B4 (de) | 2014-09-30 | 2016-06-02 | Siemens Aktiengesellschaft | Betrieb eines bildgebenden medizinischen Untersuchungsgeräts mit einer Mehrzahl an Teilsystemen |
DE102014219782A1 (de) * | 2014-09-30 | 2016-03-31 | Siemens Aktiengesellschaft | Betrieb eines bildgebenden medizinischen Untersuchungsgeräts mit einer Mehrzahl an Teilsystemen |
DE102014219786A1 (de) | 2014-09-30 | 2016-03-31 | Siemens Aktiengesellschaft | Betrieb eines bildgebenden medizinischen Untersuchungsgeräts |
DE102014219779B4 (de) * | 2014-09-30 | 2016-06-02 | Siemens Aktiengesellschaft | Betrieb eines bildgebenden medizinischen Untersuchungsgeräts mit einer Mehrzahl an Teilsystemen |
DE102014219778B4 (de) | 2014-09-30 | 2016-09-01 | Siemens Healthcare Gmbh | Betrieb eines bildgebenden medizinischen Untersuchungsgeräts mit einer Mehrzahl an Teilsystemen |
US11241187B2 (en) * | 2014-11-03 | 2022-02-08 | Qian Xia | Electromagnetic wave sensing and modulating of neuronal activities |
US11058361B2 (en) * | 2014-11-06 | 2021-07-13 | Canon Medical Systems Corporation | Signal processing apparatus, imaging apparatus, and signal processing method |
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 |
US10495709B2 (en) * | 2016-07-25 | 2019-12-03 | Siemens Healthcare Gmbh | Multi-contrast simultaneous multislice magnetic resonance imaging with binomial radio-frequency pulses |
DE102016214088A1 (de) * | 2016-07-29 | 2018-02-01 | Siemens Healthcare Gmbh | Verbessertes MR-Schichtmultiplexing-Verfahren |
DE102016218955B4 (de) | 2016-09-30 | 2019-02-14 | Siemens Healthcare Gmbh | Optimiertes Magnetresonanz-Schichtmultiplexing-Verfahren |
US10234524B2 (en) * | 2016-10-18 | 2019-03-19 | Siemens Healthcare Gmbh | Shifted pulses for simultaneous multi-slice imaging |
US10641855B2 (en) * | 2016-11-02 | 2020-05-05 | The Board Of Trustees Of The Leland Stanford Junior University | Methods for mitigating local SAR hotspots and flip angle uniformity in ultra-high field simultaneous multislice imaging using parallel transmission |
GB201705577D0 (en) * | 2017-04-06 | 2017-05-24 | King's College London | Controlled excitiation and saturation of magnetization transfer systems |
EP3692381B1 (en) | 2017-10-06 | 2023-12-06 | University of Cincinnati | Systems and methods for estimating complex radiofrequency fields in a magnetic resonance imaging |
CN109714143B (zh) | 2017-10-26 | 2021-12-31 | 通用电气公司 | 基于无线或有限线缆互联的同步***和方法及mr*** |
EP3893403B1 (en) | 2018-12-28 | 2023-07-19 | Huawei Technologies Co., Ltd. | Signal processing apparatus and method and access network device |
US11874353B2 (en) | 2019-03-27 | 2024-01-16 | The General Hospital Corporation | Multi-shot echo planar imaging using reordered segments and recursive radio frequency pulse design giving matched slice profiles across segments |
US11733332B2 (en) * | 2019-12-09 | 2023-08-22 | Nous Imaging, Inc. | Systems and method of precision functional mapping-guided interventional planning |
US11980457B2 (en) * | 2021-04-30 | 2024-05-14 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for simultaneous multi-slice multitasking imaging |
US11874354B2 (en) * | 2021-05-13 | 2024-01-16 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for magnetic resonance imaging |
WO2023049524A1 (en) * | 2021-09-27 | 2023-03-30 | Regents Of The University Of Minnesota | Parallel transmit radio frequency pulse design with deep learning |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1870481A (zh) * | 2006-04-20 | 2006-11-29 | 华东师范大学 | F_tdm多通道接收方法 |
US7420371B2 (en) * | 2006-01-04 | 2008-09-02 | Enh Research Institute | Slab-selective RF coil for MR system |
CN101675354A (zh) * | 2007-05-04 | 2010-03-17 | 皇家飞利浦电子股份有限公司 | 用于mri的具有数字反馈的rf发射器 |
US7800368B2 (en) * | 2006-02-17 | 2010-09-21 | Regents Of The University Of Minnesota | High field magnetic resonance |
CN101860391A (zh) * | 1998-08-18 | 2010-10-13 | 直视集团公司 | 多层载波离散多音通信技术 |
WO2012088065A1 (en) * | 2010-12-21 | 2012-06-28 | Regents Of The University Of Minnesota | Method for reducing power deposition in magnetic resonance imaging using multiband pulses and multichannel transmission |
WO2012088060A1 (en) * | 2010-12-20 | 2012-06-28 | Regents Of The University Of Minnesota | Method for rapid whole brain magnetic resonance imaging with contrast preparation |
US20120268130A1 (en) * | 2011-03-25 | 2012-10-25 | Hans-Peter Fautz | Method and device for determining a magnetic resonance system activation sequence |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
US8154289B2 (en) * | 2008-04-11 | 2012-04-10 | The General Hospital Corporation | Method for joint sparsity-enforced k-space trajectory and radiofrequency pulse design |
US8148985B2 (en) * | 2008-10-15 | 2012-04-03 | Massachusetts Institute Of Technology | Method for reducing maximum local specific absorption rate in magnetic resonance imaging |
US8941381B2 (en) * | 2010-05-28 | 2015-01-27 | David Feinberg | Multiplicative increase in MRI data acquisition with multi-band RF excitation pulses in a simultaneous image refocusing pulse sequence |
KR20150110685A (ko) * | 2013-01-25 | 2015-10-02 | 리전츠 오브 더 유니버스티 오브 미네소타 | 다중 채널 송신기를 위한 다중 대역 rf/mri 펄스 디자인 |
US10180475B2 (en) * | 2014-04-24 | 2019-01-15 | Regents Of The University Of Minnesota | Systems and methods for multiband sweep imaging with fourier transformation |
-
2014
- 2014-01-24 KR KR1020157022565A patent/KR20150110685A/ko not_active Application Discontinuation
- 2014-01-24 CN CN201480017259.5A patent/CN105143904B/zh active Active
- 2014-01-24 WO PCT/US2014/013004 patent/WO2014116986A1/en active Application Filing
- 2014-01-24 US US14/762,093 patent/US10684337B2/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101860391A (zh) * | 1998-08-18 | 2010-10-13 | 直视集团公司 | 多层载波离散多音通信技术 |
US7420371B2 (en) * | 2006-01-04 | 2008-09-02 | Enh Research Institute | Slab-selective RF coil for MR system |
US7800368B2 (en) * | 2006-02-17 | 2010-09-21 | Regents Of The University Of Minnesota | High field magnetic resonance |
CN1870481A (zh) * | 2006-04-20 | 2006-11-29 | 华东师范大学 | F_tdm多通道接收方法 |
CN101675354A (zh) * | 2007-05-04 | 2010-03-17 | 皇家飞利浦电子股份有限公司 | 用于mri的具有数字反馈的rf发射器 |
WO2012088060A1 (en) * | 2010-12-20 | 2012-06-28 | Regents Of The University Of Minnesota | Method for rapid whole brain magnetic resonance imaging with contrast preparation |
WO2012088065A1 (en) * | 2010-12-21 | 2012-06-28 | Regents Of The University Of Minnesota | Method for reducing power deposition in magnetic resonance imaging using multiband pulses and multichannel transmission |
US20120268130A1 (en) * | 2011-03-25 | 2012-10-25 | Hans-Peter Fautz | Method and device for determining a magnetic resonance system activation sequence |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109923427A (zh) * | 2016-11-09 | 2019-06-21 | Cr发展公司 | 对样本执行扩散加权磁共振测量的方法 |
USRE49978E1 (en) | 2016-11-09 | 2024-05-21 | The Brigham And Women's Hospital, Inc. | Method of performing diffusion weighted magnetic resonance measurements on a sample |
Also Published As
Publication number | Publication date |
---|---|
US20150362574A1 (en) | 2015-12-17 |
WO2014116986A1 (en) | 2014-07-31 |
KR20150110685A (ko) | 2015-10-02 |
US10684337B2 (en) | 2020-06-16 |
CN105143904B (zh) | 2018-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105143904A (zh) | 用于多信道发射机的多频带射频/磁共振成像脉冲设计 | |
CN102870000B (zh) | 使用单个和多个通道接收器线圈的用于同时多切片磁共振成象的方法 | |
US9396562B2 (en) | MRI reconstruction with incoherent sampling and redundant haar wavelets | |
CN101669038B (zh) | 用于propeller mri的磁共振设备和方法 | |
US8154289B2 (en) | Method for joint sparsity-enforced k-space trajectory and radiofrequency pulse design | |
US9417298B2 (en) | Local SAR reduction in multi-slice pTx via SAR-hopping between excitations | |
US10215822B2 (en) | Image processing apparatus and magnetic-resonance imaging apparatus | |
US10345409B2 (en) | System and method for simultaneous multislice excitation using combined multiband and periodic slice excitation | |
WO2014054540A1 (ja) | 磁気共鳴イメージング装置及び画像処理装置 | |
CN102736045A (zh) | 具有局部sar 约束的并行发射射频脉冲设计 | |
WO2007013423A1 (ja) | 磁気共鳴イメージング装置 | |
RU2702859C2 (ru) | Параллельная мультисрезовая мр-визуализация с подавлением артефактов боковой полосы частот | |
EP3100067B1 (en) | Simultaneous multislice mri with random gradient encoding | |
US8473028B2 (en) | K-space sample density compensation for magnetic resonance image reconstruction | |
US10222439B2 (en) | System and method for spiral multislab magnetic resonance imaging | |
US11320506B2 (en) | Multiband spiral cardiac MRI with non-cartesian reconstruction methods | |
JP2015128585A (ja) | 磁気共鳴生データの再構成方法および装置 | |
US10120048B2 (en) | Systems and methods for efficient radial magnetic resonance imaging with azimuthal equidistant projections | |
US10466329B2 (en) | Partial Fourier acquisition and reconstruction for k-space shells based magnetic resonance imaging | |
US11266324B2 (en) | System and methods for fast multi-contrast magnetic resonance imaging | |
US10324149B2 (en) | Systems and methods for generalized slice dithered enhanced resolution magnetic resonance imaging | |
US11231476B2 (en) | Accelerated magnetic resonance imaging acquisition using two-dimensional pulse segments as virtual receivers | |
WO2018187739A1 (en) | System and method for simultaneous multi-slice magnetic resonance fingerprinting imaging using a parallel transmit radio frequency coil array |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |