CN106030329A - 使用多回波分段k空间采集的MR成像 - Google Patents

使用多回波分段k空间采集的MR成像 Download PDF

Info

Publication number
CN106030329A
CN106030329A CN201480075316.5A CN201480075316A CN106030329A CN 106030329 A CN106030329 A CN 106030329A CN 201480075316 A CN201480075316 A CN 201480075316A CN 106030329 A CN106030329 A CN 106030329A
Authority
CN
China
Prior art keywords
space
echo
peripheral
sampling
signal
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
CN201480075316.5A
Other languages
English (en)
Other versions
CN106030329B (zh
Inventor
G·M·贝克
G·H·范艾佩伦
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
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 Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN106030329A publication Critical patent/CN106030329A/zh
Application granted granted Critical
Publication of CN106030329B publication Critical patent/CN106030329B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
    • 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/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/482—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/4824—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5616—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using gradient refocusing, e.g. EPI
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5617—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using RF refocusing, e.g. RARE

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设备(1)的检查体积中的至少对象(10)进行MR成像的方法。本发明的目的是使用相对于运动鲁棒的多回波成像序列来实现快速MR成像。本发明的方法包括如下步骤：‑使所述对象(10)经受多回波成像序列的多次发射，通过每次发射生成回波信号的串；采集所述回波信号，其中，每个回波信号表示k空间分布，其中，k空间(20)被分成中心k空间部分(21)和一个或多个外围k空间部分(22、23)，其中，所述中心k空间部分(21)是通过对所述多回波成像序列的单次发射而被采样的，并且其中，所述外围k空间部分(22、23)是通过对所述多回波序列的一次或多次另外的发射(25、26、27、28)而被采样的，并且根据所述k空间分布来重建MR图像。此外，本发明涉及用于执行这种方法的MR设备以及涉及用于在MR设备上运行的计算机程序。

Description

使用多回波分段k空间采集的MR成像
技术领域
本发明涉及磁共振(MR)成像的领域。本发明涉及对对象进行MR成像的方法。本发明还涉及MR设备以及在MR设备上运行的计算机程序。
背景技术
目前，特别是在医学诊断的领域中，广泛使用利用磁场与核自旋之间的相互作用以形成二维或三维图像的图像形成MR方法，因为对于软组织的成像而言，它们在许多方面优于其他成像方法，不需要电离辐射并且通常是无创的。
根据一般的MR方法，要被检查的对象(例如，患者的身体)被布置在强的均匀磁场中，其方向同时定义了测量所基于的坐标系的轴(通常为z轴)。磁场根据磁场强度针对个体核自旋产生不同的能量水平，能够通过应用定义频率(所谓的拉莫尔频率或MR频率)的电磁交变场(RF场)来激励个体核自旋(自旋共振)。从宏观的角度讲，个体核自旋的分布产生总体磁化，能够通过应用合适频率的电磁脉冲(RF脉冲)使其偏离平衡状态，从而所述磁化执行关于z轴的进动运动。进动运动描绘锥形的表面，所述锥形的孔径角被称为翻转角。翻转角的幅度取决于所应用的电磁脉冲的强度和持续时间。在所谓的90°脉冲的情况下，所述自旋从z轴偏转到横向平面(翻转角为90°)。
在RF脉冲终止之后，所述磁化弛豫返回到平衡的原始状态，其中，在z方向的磁化以第一时间常量T1(自旋晶格或纵向驰豫时间)再次建立，并且在垂直于z方向的方向上的磁化以第二时间常量T2(自旋-自旋或横向弛豫时间)弛豫。借助接收RF线圈能够探测磁化的变化，所述接收RF线圈以这样的方式在MR设备的检查体积之内进行布置和取向，使得在垂直于z轴的方向中测量磁化的变化。在应用例如90°的脉冲之后，横向磁化的衰减伴随有核自旋(由局部磁场不均匀感生的)从具有相同相位的有序状态到所有相位角均匀分布(失相)的状态的转变。所述失相能够借助重聚焦脉冲(例如，150°的脉冲)来补偿。这在接收线圈中产生回波信号(自旋回波)。
为了在身体中实现空间分辨，沿三个主轴延伸的恒定磁场梯度被叠加在均匀磁场上，得到自旋共振频率的线性空间相关性。在接收线圈中拾取的信号然后包含不同频率的分量，其与身体中的不同位置相关联。经由所述接收线圈获得的信号数据对应于空间频率域并且被称为k空间数据。所述k空间数据通常包括利用不同的相位编码采集的多条线。每条线通过收集若干样本来数字化。k空间数据的集合借助图像重建算法被转换为MR图像。
所谓的“快速自旋回波”(TSE)序列是周知的多回波成像序列。TSE序列的一次“发射(shot)”包括初始RF脉冲，其用于激励磁共振，跟随有多个快速施加的(通常为180°)重聚焦RF脉冲，其生成不同相位编码的自旋回波信号的串。采集所述回波信号，其中，每个回波信号表示k空间分布(profile)，即，k空间的一维样本，在k空间中的所述k空间分布的位置是通过所应用的频率编码的和相位编码的切换的磁场梯度的序列来确定的。所谓的“快速因子”(TF)是在每次激励之后采集的回波的数量。通常，应用所述多回波序列的多次发射以对k空间进行完全采样，以便能够根据所采集的k空间分布来重建MR图像。k空间通常被分段，使得每个段归属于由所述序列的每次发射所生成的回波的串中的一个回波索引数。换言之，每次发射通常从k空间的所有分段采集k空间分布，其中k空间中的分段的数量等于TF。
目前，在几乎所有应用中应用TSE序列。尽管其广泛使用，已知其由于分段的k空间分布的采集而对运动非常敏感。被成像对象的小的移动，例如皮下脂肪的小的移动，能够引入伪影，所述伪影与感兴趣部分交叠，干扰扫描的总体诊断值。
上文提到的k空间分段原理和问题适用于自旋回波(TSE)和梯度回波序列(TFE)两者以及‘回波平面成像’-EPI-序列。
美国专利US2006/0273791公开了一种磁共振成像方法，其中，在每个子集中所采集的k空间分布的次序通常不同于其他分布子集的次序。
发明内容
根据上述内容，容易认识到，存在对改进的MR成像技术的需求。因此，本发明的目的是使用相对于运动更为鲁棒的多回波成像技术来实现快速MR成像。
根据本发明，公开了一种对被放置在MR设备的检查体积中的对象进行MR成像的方法。所述方法包括如下步骤：
-使所述对象经受多回波成像序列的多次发射，通过每次发射生成回波信号的串，
-采集所述回波信号，其中，每个回波信号表示k空间分布，其中，k空间被分成中心k空间部分以及一个或多个外围k空间部分，其中，所述中心k空间部分是通过对所述多回波成像序列的单次发射而被采样的，并且其中，所述外围k空间部分是通过对所述多回波序列的一次或多次另外的发射而被采样的，并且
-根据所述k空间分布来重建MR图像。
本发明提出了在单次发射中采集中心k空间。因此，利用在基本单个时间点处对中心k空间部分进行采样，有效降低了对运动的敏感度。
必须指出，表述“在单次发射中进行采集”包括这样的变量，即：中心k空间是通过仅单次发射的回波信号的串的子集来采样的，其中，相应发射的剩余回波信号被用于对外围k空间位置进行采样。
根据本发明的优选实施例，在k空间中对所述中心k空间部分的所述采样开始的位置处开始对所述外围k空间部分的采样，并且/或者在k空间中对所述中心k空间部分的所述采样结束的位置处结束对所述外围k空间部分的所述采样。外围采集分别在k空间中在对中心采集的开始点和结束点处开始和结束，以避免信号幅度和/或相位中的不连续，所述不连续将呈现为图像伪影。
更为一般地说，优选根据本发明执行k空间采样，使得从k空间中的相邻位置采集的k空间分布具有相邻(即，连贯)或相同的回波索引数，其中，所述连贯回波索引数指示在回波信号串内的每个回波信号的时间位置。
在三维(3D)采集中，k空间分段可以是二维的，例如，在ky和kz中。在本发明的优选实施例中，不仅仅对外围和中心k空间采样的开始点和结束点要被对齐，而且可能中心和外围k空间部分的所有k空间分段应当被对齐，以便避免在信号幅度和/或相位中的不连续。
本发明的方法提出能够针对各自的应用根据需要来调整针对中心和外围k空间部分的采样策略。
根据本发明的优选实施例，本发明的多回波成像序列是自旋回波序列，所述自旋回波序列包括：
-至少一个RF脉冲，其用于激励磁共振，
-多个重聚焦RF脉冲，
-用于对每个回波信号进行频率编码的多个切换的磁场梯度，以及
-用于对每个回波信号进行相位编码的多个切换的磁场梯度。换言之，本发明的方法可以与自身已知的‘快速’自旋回波序列(例如TSE或自旋回波型EPI)组合应用。
备选地，所述多回波成像序列是梯度回波序列，所述梯度回波序列包括：
-至少一个RF脉冲，其用于激励磁共振，
-用于对磁共振进行重聚焦并且用于对每个回波信号进行频率编码的多个切换的磁场梯度，以及
-用于对每个回波信号进行相位编码的多个切换的磁场梯度。换言之，本发明的方法可以与自身已知的‘快速’自旋回波序列(例如TFE或EPI)组合应用。
本发明的方法尤其非常适于与笛卡尔k空间采样方案进行组合。然而，所述方法也可以与非笛卡尔采样(诸如PROPELLER(螺旋桨)或螺旋扫描)组合应用。针对k空间采样的所有类型，术语“中心k空间部分”在本发明的意义内是指接近k空间原点的k空间分布。
关于螺旋桨扫描，已知单发射螺旋桨叶采集与多发射螺旋桨方法相比相对于运动通常是非常鲁棒的。本发明与螺旋桨的组合意味着，每个螺旋桨叶的中心k空间分布是在单发射中采集的，而每个螺旋桨叶的外围k空间分布是使用所应用的多回波序列的另外的发射分段的。这造成抗运动的多发射螺旋桨采集。
根据本发明的另一优选实施例，k空间的至少部分是不完全采样的。这意味着，例如，本发明的方法能够与诸如本身已知的‘部分矩阵’技术或者‘压缩感测’(CS)技术的技术进行组合。
‘部分矩阵’是扫描时间减少的方法，其利用k空间的复共轭属性。由于来自负k空间位置的k空间分布与来自对应的正k空间位置的复共轭k空间分布大致相同，因而根据本发明可以仅对外围正或负部分进行采样，以便复制对整个k空间的采样。这造成扫描的时间显著减少，同时保持完整的分辨率。
根据本发明的又一优选实施例，k空间是利用可变密度采样的。换言之，这意味着，从外围k空间部分采集的k空间分布可以在相位编码的方向上不规则地分布。这可以被采用，以便通过运动优化的可变密度压缩感测来进一步减少扫描时间。CS的理论已知有很大的潜力用于显著的信号数据减少。在CS理论中，在变换域中具有稀疏表示的信号数据集能够通过应用合适的正则化算法来从欠采样的测量结果恢复。欠采样的可能性导致显著减少的采集时间。作为用于信号采样和重建的数学架构，CS规定了这样的状况，在所述状况下，信号数据集能够即使在k空间采样密度远低于奈奎斯特准则的情况下也能够准确地或者至少以高的图像质量来进行重建，并且其也提供用于这样的重建的方法。在大多数现有基于CS的MR采集和重建方案中，使用基本的CS构想，其仅利用这样的先决条件，即，MR信号数据在变换域中是稀疏的。例如，M.Lustig等人已经提出了针对快速MR成像应用CS(M.Lustig等人：“Sparse MRI:The Application of CompressedSensing for Rapid MR Imaging”，Magnetic Resonance in Medicine，58，1182-1195，2007)。也已知的是，由于关于未知MR图像的额外先验信息在特定应用中是可用的，有利的是将该先验信息并入到CS重建中。
根据本发明的另一优选实施例，k空间的至少中心部分是用于信号平均的目的利用若干采集来采样的，以便补偿FID和运动伪影。这意味着，本发明的方法能够与部分平均化进行组合，其中，分别在中心和外围k空间中应用不同的信号平均化。
根据又一优选实施例，在对回波信号的采集期间执行对所述对象的运动的检测，其中，在所述对象的运动期间从外围k空间部分采集的k空间分布被舍弃。
运动例如可以通过自身已知的导航技术进行检测。出于该目的，导航回波的采集可以被并入到本发明的方法中。备选地，可以通过使用k空间一致性标准或者通过内部和/或外部运动传感器(例如，呼吸传感器)来检测运动。
如果运动是在k空间分布的采集期间从外围部分在特定点处采集的，则该k空间分布(或者通过各自的发射而采集的所有k空间分布)可以被舍弃，并且已知的迭代重建算法或POCS类方法可以被用于重建无运动的图像。
‘凸集投影(POCS)’是有力的数学工具，其用于重建不完整和/或不一致的MR数据。其已知非常适于从部分k空间数据进行MR图像重建。
检测发射中的运动的另一种可能是识别在不同发射之间的k空间不一致性。例如，已知的COCOA方法(‘Data Convolution and CombinationOperation for Motion Ghost Artefacts Reduction’)可以与本发明的方法组合使用。COCOA目前被用于执行对k空间数据的后处理，以减少运动伪影。来自k空间核区域的样本被用于检测和校正运动污染的k空间样本。基于所采集的k空间样本，能够对相邻k空间样本进行估计(参见Feng Huang等人，Magnetic Resonance in Medicine，64，157-166，2010)。根据本发明，在COCOA中使用的k空间核可以根据从中心k空间部分采集的k空间样本中导出。
相比于在常规多回波技术中的k空间分布次序，根据本发明从外围k空间部分移除(运动污染的)k空间分布不影响中心k空间，因此所提到的重建技术将更为成功。
至此所描述的本发明的方法能够借助MR设备来执行，所述MR设备包括：至少一个主磁体线圈，其用于在检查体积之内生成均匀的、静态的磁场；若干梯度线圈，其用于在检查体积之内的不同空间方向上生成切换的磁场梯度；至少一个RF线圈，其用于在检查体积之内生成RF脉冲和/或用于从定位在检查体积之内的对象接收MR信号；控制单元，其用于控制RF脉冲和切换的磁场梯度的时间演替；以及重建单元。本发明的方法能够例如通过对MR设备的重建单元和/或控制单元的对应编程来实施。
本发明的方法能够有利地在当前临床使用的大多数MR设备中执行。最后，仅仅需要利用计算机程序，通过所述计算机程序，MR设备被控制为使得其执行本发明的以上解释说明的方法的步骤。所述计算机程序能够存在于数据载体上或者能够存在于数据网络上，从而能够被下载以供在MR设备的控制单元中进行安装。
附图说明
随附的附图公开了本发明的优选实施例。然而，应当理解，附图仅仅是出于图示说明的目的而设计的，并且不是对本发明的限制范围的限定。在附图中：
图1示出了用于执行本发明的方法的MR设备；
图2示出了示意性图示本发明的采样方案的第一实施例的k空间的图表；
图3示出了示意性图示本发明的采样方案的第二实施例的k空间的图表；
图4示出了作为从本发明的k空间采样方案得到的k空间位置的函数信号幅度图表；
图5示出了示意性图示本发明的采样方案的第三实施例的k空间的图表；
图6示出了根据本发明的第四实施例的归因于k空间位置的发射数量和回波索引数；
图7示出了示意性图示利用螺旋采样的本发明的第四实施例的k空间的图表。
具体实施方式
参考图1，示出了MR设备1。所述设备包括超导或电阻主磁体线圈2，使得沿通过检查体积的z轴创建充分均匀的、时间上恒定的主磁场。
磁共振生成和操控系统应用一系列RF脉冲和切变的磁场梯度，从而：反转或激励核磁自旋、诱发磁共振、重聚焦磁共振、操纵磁共振、空间地或者以其他方式对磁共振进行编码、使自旋饱和等等，以执行MR成像。
更为具体地，梯度脉冲放大器3向沿检查体积的x、y和z轴的全身梯度线圈4、5和6中的选定的那些应用电流脉冲。数字RF频率发射器7经由发送/接收开关8向全身体积RF线圈9发射RF脉冲或脉冲包，以将RF脉冲发射到检查体积中。典型的MR成像序列包括彼此一起进行的短持续时间的RF脉冲段的包，并且任何施加的磁场梯度实现对核磁共振的选定操控。RF脉冲被用于使共振饱和、激励共振、反转共振、重聚焦共振或者操纵共振，并且选择被定位在检查体积中的身体10的部分。MR信号也由全身体积RF线圈9拾取。
为了生成身体10的限定区域的MR图像，邻近选定用于成像的区域放置一组局部阵列RF线圈11、12、13。阵列线圈11、12、13能够被用于接收由身体线圈RF发射所感生的MR信号。
所得到的MR信号由全身体积RF线圈9和/或由阵列RF线圈11、12、13拾取，并且由优选包括前置放大器(未示出)的接收器14进行解调。接收器14经由发送/接收开关8被连接到RF线圈9、11、12和13。
主计算机15控制梯度脉冲放大器3和发射器7，以生成多个MR成像序列中的任一个，例如回波平面成像(EPI)、回波体积成像、梯度和自旋回波成像、快速自旋回波(TSE)成像等。针对选定的序列，接收器14在每个RF激励脉冲之后快速相继的接收单条或多条MR数据线。数据采集系统16执行对所接收的信号的模数转换，并且将每条MR数据线转换为适于进一步处理的数字格式。在现代MR设备中，数据采集系统16是单独的计算机，所述计算机被专用于原始图像数据的采集。
最后，数字原始图像数据被重建处理器17重建为图像表示，重建处理器17应用傅里叶变换或其他适当的重建算法。MR图像可以表示通过患者的平面切片、平行平面切片的阵列、三维体积等。所述图像然后被存储在图像存储器中，其中，所述图像存储器可以被访问，用于将切片、投影或图像表示的其他部分转换为合适的格式以供例如经由视频监视器18进行可视化，所述视频监视器18提供所得到的MR图像的人类可读的显示。
继续参考图1并且进一步参考图2-7，解释了本发明的成像方法的实施例。
如在图2中所示，k空间20被分成中心k空间部分21、正外围k空间部分22以及负外围k空间部分23。中心k空间部分21通过由箭头24指示的TSE序列的单次发射进行完全采样。箭头24指示k空间采样的方向，即，由所施加的多回波序列的单次发射所生成的回波的次序。外围k空间部分22、23中的每个是通过多回波序列中的至少两次另外的发射进行采样的。这是由箭头25、26和27、28来指示的。发射25、26和27、27以交错的方式对各自的外围k空间部分22、23进行采样。箭头25、26和27、28再次指示由多回波序列的各自发射所生成的回波的次序。如能够从箭头的方向看到的，在对中心k空间部分的采样24开始的位置处开始对外围k空间部分23的采样27、28，而在k空间中在对中心k空间部分21的采样24结束的位置处结束对外围k空间部分22的采样25、26。对于在例如发射26中检测到运动的情况，该发射的k空间分布能够被舍弃并且合适的迭代重建、COCOA或POCS类方法能够被用于重建无运动的MR图像。移除发射26的分布不影响中心k空间部分21，因此对k空间分布的舍弃能够被期望仅诱发不显著量的伪影。
在图3中所示的实施例中，外围k空间部分22、23通过多回波发射25、26、27、28进行采样，再次在开始对中心k空间部分21的采样的位置处开始采集k空间分布，并且在中心k空间扫描的结束位置处解释采集，其中，由连接弧29指示相位变化从最小ky值到最大ky值的跳跃。图4示出了作为k空间位置(在相位编码方向ky)的函数的回波信号幅度的对应表现。在图表中的每个点指示单个回波信号。如能够在图表中看到的，任何的信号不连续被有效地避免。图4的图表也示出了通过虚线对外围k空间部分22、23的分段以及根据本发明对中心k空间部分21的“未分段的”采集。
根据又一实施例，单次发射24的仅部分被用于对中心k空间部分21的采样，如在图5中所示。这对于从开始到结束的具有非常大(T2-感生)的信号变化的非常长的发射可能是有用的。在图5中所示的实施例中，k空间20被分成中心k空间部分21和外围k空间部分22、23、30、31，其中，k空间部分30、31能够被称为中间k空间部分。中间k空间部分30、31通过施加发射24和25以多发射方式进行采样。中间k空间部分30、31可以被不同地分在正和负k空间中。外围k空间部分22、23通过发射26、27、28进行采样。在该实施例中，在对“扩展的”中心k空间部分(由部分21、30、31构成)的采样24、25开始的位置处开始对外围k空间部分23的采样26、27、28，而在k空间中对部分21、30、31的采样24、25结束的位置处结束对外围k空间部分22的采样26、27、28。
在图2-5中所示的实施例的所有采样24、25、26、27、28应用相对于相位编码的k空间分布的线性次序。
图6图示了图5的采样方案的特定实施方案。所述图表示出了发射数S和回波索引数I，后者指示每个回波信号在由TSE序列的各自发射所生成的连贯回波信号串内的序列号。如能够看到的，从k空间中的相邻位置采集的k空间分布具有相继或相同的回波索引数。在所描绘的范例中，应用使用快速因子为六的线性分布次序的TSE序列的五次发射。发射1(S＝1)在中间部分处开始采集并且跳跃到中心k空间部分21，其在单次发射中被采集。发射2-6(S＝2-6)采集外围k空间部分。
在图5和图6中所示的采样方案允许在k空间上调谐T2相关的对比度改变。对于典型的指数T2衰减，在正外围k空间部分22、30中的较强的衰减能够通过调整分段的尺寸来进行调谐。快速因子和发射的总数量能够保持相同。T2信号分布仅仅通过将较大的T2差异进一步偏移到外部k空间来重新排序。
本发明的方法可以与翻转角扫掠进行组合，因为其通常被用在TSE方法中，以便在每次发射中尽可能长地维持磁化。在TSE读出期间的翻转角变化可以针对k空间分段调整以便优化成像结果。
本发明的方法可以与现有方法进行组合，以补偿T2模糊，诸如，例如反转k空间分布次序。同样地，与部分平均、k空间一致性导航(‘ONAV’)和/或压缩感测的组合在本发明的上下文中是相关的。具体而言，本发明实现如针对压缩感测所需要的运动鲁棒的可变k空间密度方法。
将部分平均、压缩感测和本发明的k空间分布排序原理进行组合允许容易地采集k空间的中心部分21两次或者更多次，用于足够的信噪比(SNR)和FID伪影抑制，而可变密度随机方法可以被应用以利用例如泊松盘分布来实现接近中心部分21的更密的k空间采样。能够利用来自压缩感测重建的数据舍弃来自可变密度采集的被运动污染的数据。
同样地，呼吸触发、呼吸保持方法以及还在中心k空间部分中的更为关键的k空间一致性拒绝机制是在本发明的上下文中能够想像的，以便考虑被成像对象的运动。
如在图2-6中所展示的，本发明的方法尤其非常适于与笛卡尔k空间采样方案进行组合。然而，所述方法也可以与非笛卡尔采样(诸如螺旋扫描)组合应用，如在图7中所示。中心k空间部分21通过利用从k空间原点朝外的螺旋k空间轨迹进展的多回波序列的单次发射来扫描。沿着从k空间向内的外部边界进展的相反旋转的螺旋k空间轨线来扫描外围k空间部分22，其中，在k空间中对中心k空间部分21的采样结束的位置处结束对外围k空间部分22的螺旋采样。
尽管上文在二维成像的背景下图示了本发明的方法，其也能够被应用于三维(3D)成像，而不脱离本发明的原理。在3D成像中，如上文所提到的，k空间分段能够是二维的，在ky和kz中。因此，在本发明的优选实施例中，不仅仅对外围和中心k空间采集的开始点和结束点被对齐，而且中心和外围k空间部分的k空间分段和各自的回波信号索引数也被对齐。在示范性实施例中，中心k空间能够被定义为椭圆中心ky-kz部分，其在单次发射(或者单次发射的部分)中被采集。椭圆中心k空间采集发射的开始点和结束点与对外围k空间采集的发射的开始点和结束点对齐。

Claims (15)

1.一种对被放置在MR设备(1)的检查体积中的对象(10)进行MR成像的方法，所述方法包括如下步骤：
-使所述对象(10)经受多回波成像序列的多次发射，通过每次发射生成回波信号的串；
-采集所述回波信号，其中，每个回波信号表示k空间分布，其中，k空间(20)被分成中心k空间部分(21)和一个或多个外围k空间部分(22、23)，其中，所述中心k空间部分(21)通过对所述多回波成像序列的单次发射(24)而被采样，并且其中，所述外围k空间部分(22、23)通过对所述多回波序列的一次或多次另外的发射(25、26、27、28)而被采样，其中，对所述外围k空间部分(22、23)的采样开始于在k空间(20)中对所述中心k空间部分(21)的所述采样开始的位置处，并且/或者对所述外围k空间部分(22、23)的所述采样结束于在k空间(20)中对所述中心k空间部分(21)的所述采样结束的位置处；并且
-根据所述k空间分布来重建MR图像。
2.根据权利要求1所述的方法，其中，所述多回波成像序列是自旋回波序列，所述自旋回波序列包括：
-至少一个RF脉冲，其用于激励磁共振，
-多个重聚焦RF脉冲，
-用于对每个回波信号进行频率编码的多个切换的磁场梯度，以及
-用于对每个回波信号进行相位编码的多个切换的磁场梯度。
3.根据权利要求1或2所述的方法，其中，所述多回波成像序列是梯度回波序列，所述梯度回波序列包括：
-至少一个RF脉冲，其用于激励磁共振，
-用于对磁共振进行重聚焦并且用于对每个回波信号进行频率编码的多个切换的磁场梯度，以及
-用于对每个回波信号进行相位编码的多个切换的磁场梯度。
4.根据权利要求2或3所述的方法，其中，所述多回波成像序列是TFE序列或TSE序列。
5.根据权利要求1至4中的任一项所述的方法，其中，所述中心k空间部分(21)和所述外围k空间部分(22、23)在笛卡尔采样网格上被采样。
6.根据权利要求1至4中的任一项所述的方法，其中，所述中心k空间部分(21)和所述外围k空间部分(22、23)是在螺旋桨采集中的个体k空间桨叶的部分。
7.根据权利要求1至5中的任一项所述的方法，其中，所述中心k空间部分(21)和所述外围k空间部分(22、23)是沿着在k空间中的螺旋轨迹被采样的。
8.根据权利要求1-7中的任一项所述的方法，其中，从k空间(20)中的相邻位置采集的所述k空间分布具有相邻或者相同的回波索引数(I)，每个回波索引数(I)指示每个回波信号在通过所述多回波成像序列的每次发射所生成的所述回波信号的串内的位置。
9.根据权利要求1-8中的任一项所述的方法，其中，所述外围k空间部分(22、23)中的至少部分是不完整采样的。
10.根据权利要求1-9中的任一项所述的方法，其中，k空间(20)是利用可变密度采样的。
11.根据权利要求1-10中的任一项所述的方法，其中，所述中心k空间部分(21)被采集两次或更多次以进行信号平均。
12.根据权利要求1-11中的任一项所述的方法，其中，执行对所述对象(10)的运动的检测，其中，来自所述外围k空间部分(22、23)的在所述对象(10)的运动期间所采集的k空间分布被舍弃。
13.根据权利要求1-12中的任一项所述的方法，其中，所述MR图像是使用压缩感测来重建的。
14.一种MR设备，包括：至少一个主磁体线圈(2)，其用于在检查体积之内生成均匀的、静态的磁场；若干梯度线圈(4、5、6)，其用于在所述检查体积内的不同空间方向上生成切换的磁场梯度；至少一个RF线圈(9)，其用于在所述检查体积之内生成RF脉冲和/或用于从被定位在所述检查体积中的对象(10)接收MR信号；控制单元(15)，其用于控制RF脉冲和切换的磁场梯度的时间演替；以及重建单元(17)，其中，所述MR设备(1)被布置为执行如下步骤：
-使所述对象(10)经受多回波成像序列的多次发射，通过每次发射生成回波信号的串；
-采集所述回波信号，其中，每个回波信号表示k空间分布，其中，k空间(20)被分成中心k空间部分(21)和一个或多个外围k空间部分(22、23)，其中，所述中心k空间部分(21)是通过对所述多回波成像序列的单次发射(24)而被采样的，并且其中，所述外围k空间部分(22、23)是通过对所述多回波序列的一次或多次另外的发射(25、26、27、28)而被采样的，
-在k空间(20)中对所述中心k空间部分(21)的所述采样开始的位置处开始对所述外围k空间部分(22、23)的采样，并且/或者在k空间(20)中对所述中心k空间部分(21)的所述采样结束的位置处结束对所述外围k空间部分(22、23)的所述采样，并且
-根据所述k空间分布来重建MR图像。
15.一种用于在MR设备上运行的计算机程序，所述计算机程序包括指令用于：
-产生多回波成像序列的多次发射，通过每次发射生成回波信号的串，
-采集所述回波信号，其中，每个回波信号表示k空间分布，其中，k空间(20)被分成中心k空间部分(21)和一个或多个外围k空间部分(22、23)，其中，所述中心k空间部分(21)是通过对所述多回波成像序列的单次发射(24)而被采样的，并且其中，所述外围k空间部分(22、23)是通过对所述多回波序列的一次或多次另外的发射(25、26、27、28)而被采样的，
-在k空间(20)中对所述中心k空间部分(21)的所述采样开始的位置处开始对所述外围k空间部分(22、23)的采样，并且/或者在k空间(20)中对所述中心k空间部分(21)的所述采样结束的位置处结束对所述外围k空间部分(22、23)的所述采样，并且
-根据所述k空间分布来重建MR图像。
CN201480075316.5A 2013-12-12 2014-12-04 使用多回波分段k空间采集的MR成像方法 Active CN106030329B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP13196861.2 2013-12-12
EP13196861 2013-12-12
PCT/EP2014/076478 WO2015086415A1 (en) 2013-12-12 2014-12-04 Mr imaging using multi-echo segmented k-space acquisition

Publications (2)

Publication Number Publication Date
CN106030329A true CN106030329A (zh) 2016-10-12
CN106030329B CN106030329B (zh) 2019-09-27

Family

ID=49759133

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201480075316.5A Active CN106030329B (zh) 2013-12-12 2014-12-04 使用多回波分段k空间采集的MR成像方法

Country Status (5)

Country Link
US (1) US10191133B2 (zh)
EP (1) EP3080633A1 (zh)
JP (1) JP2016539722A (zh)
CN (1) CN106030329B (zh)
WO (1) WO2015086415A1 (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109814058A (zh) * 2017-11-22 2019-05-28 西门子保健有限责任公司 Caipirinha波space
CN110389310A (zh) * 2018-04-16 2019-10-29 西门子医疗有限公司 用于产生磁共振记录的方法和控制装置
CN110412493A (zh) * 2018-04-26 2019-11-05 西门子医疗有限公司 记录磁共振数据组的方法、数据载体以及磁共振设备
CN110687488A (zh) * 2019-09-18 2020-01-14 上海联影医疗科技有限公司 磁共振扫描方法和磁共振成像装置
CN110730912A (zh) * 2017-06-09 2020-01-24 皇家飞利浦有限公司 并行多切片mr成像
CN111257809A (zh) * 2020-01-23 2020-06-09 上海东软医疗科技有限公司 磁共振成像方法、装置、存储介质及医疗设备
CN114041778A (zh) * 2022-01-13 2022-02-15 康达洲际医疗器械有限公司 一种应用于骨关节磁共振成像的融合成像方法与系统
US11543482B2 (en) * 2018-10-16 2023-01-03 Koninklijke Philips N.V. Magnetic resonance imaging using motion-compensated image reconstruction

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6998218B2 (ja) * 2015-07-15 2022-01-18 コーニンクレッカ フィリップス エヌ ヴェ 動き検出を用いるｍｒ撮像
EP3324844A4 (en) 2015-07-23 2018-08-22 Shanghai United Imaging Healthcare Co., Ltd. System and method for magnetic resonance imaging
KR101802336B1 (ko) 2016-02-19 2017-11-28 삼성전자주식회사 다중 여기 rf 펄스를 이용한 자기공명영상 획득 방법 및 이를 위한 자기공명영상 장치
WO2018015298A1 (en) 2016-07-21 2018-01-25 Koninklijke Philips N.V. Motion-corrected compressed sensing magnetic resonance imaging
JP7023954B2 (ja) * 2016-11-24 2022-02-22 コーニンクレッカ フィリップス エヌ ヴェ プロペラｍｒイメージング
US10557907B2 (en) * 2017-09-15 2020-02-11 Uih America, Inc. System and method for reducing nyquist ghost artifact
EP3540453A1 (en) 2018-03-13 2019-09-18 Koninklijke Philips N.V. Mr imaging with spiral acquisition
GB2573557B (en) 2018-05-10 2020-06-10 Siemens Healthcare Gmbh Method of reconstructing magnetic resonance image data
EP3599476A1 (de) 2018-07-24 2020-01-29 Siemens Healthcare GmbH Erfassen von magnetresonanz-daten mittels echozügen
US10928473B2 (en) * 2018-09-17 2021-02-23 General Electric Company Methods and systems for reduced shading and blurring in magnetic resonance imaging
DE102019102799A1 (de) * 2019-02-05 2020-08-06 Albert-Ludwigs-Universität Freiburg Verfahren der Magnetresonanztomographie
JP7350568B2 (ja) 2019-08-28 2023-09-26 キヤノンメディカルシステムズ株式会社 磁気共鳴データ収集方法及び磁気共鳴イメージング装置
DE102020202576B4 (de) 2020-02-28 2022-05-25 Bruker Biospin Mri Gmbh Verfahren zum Erzeugen eines Magnetresonanzbildes
JP7374822B2 (ja) * 2020-03-09 2023-11-07 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング方法及び装置
EP3913387A1 (en) 2020-05-19 2021-11-24 Koninklijke Philips N.V. Motion estimation and correction in magnetic resonance imaging
CN113835058A (zh) 2020-06-24 2021-12-24 通用电气精准医疗有限责任公司 采集和处理mr数据的方法、mri系统和方法、存储介质
US20230169704A1 (en) * 2021-11-29 2023-06-01 Canon Medical Systems Corporation Advanced signal combination for improved overlapping image correction

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6255820B1 (en) * 1999-03-08 2001-07-03 Picker International, Inc. Variable bandwidth MRI data collection
CN101144853A (zh) * 2006-09-13 2008-03-19 西门子公司 用于三维磁共振成像的测量序列以及磁共振设备
US20120112745A1 (en) * 2009-07-16 2012-05-10 Hitachi Medical Corporation Magnetic resonance imaging apparatus and magnetic resonance imaging method
CN102772210A (zh) * 2011-05-10 2012-11-14 西门子公司 弥散加权磁共振成像
US20130307537A1 (en) * 2011-01-28 2013-11-21 Hitachi Medical Corporation Magnetic resonance imaging device and angiography using same

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5402067A (en) * 1993-08-04 1995-03-28 Board Of Trustees Of The Leland Stanford Junior University Apparatus and method for rare echo imaging using k-space spiral coverage
WO1999014616A1 (en) 1997-09-15 1999-03-25 Ge Medical Systems Israel, Ltd. Keyhole mri
JP3847512B2 (ja) 2000-02-07 2006-11-22 株式会社日立メディコ 磁気共鳴イメージング装置
JP3853585B2 (ja) 2000-10-11 2006-12-06 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｍｒｉ装置
US6617850B2 (en) 2001-12-03 2003-09-09 Mayo Foundation For Medical Education And Research Motion correction of magnetic resonance images using phase difference of two orthogonal acquisitions
US7362099B2 (en) 2003-09-08 2008-04-22 Koninklijke Philips Electronics N.V. Randomized ordered k-space sub-sets for shared pre-pulses in MRI
US10073160B2 (en) * 2011-04-21 2018-09-11 Koninklijke Philips N.V. Magnetic resonance imaging of object in motion
US8768034B2 (en) * 2012-03-20 2014-07-01 Siemens Medical Solutions Usa, Inc. Motion compensated MR imaging system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6255820B1 (en) * 1999-03-08 2001-07-03 Picker International, Inc. Variable bandwidth MRI data collection
CN101144853A (zh) * 2006-09-13 2008-03-19 西门子公司 用于三维磁共振成像的测量序列以及磁共振设备
US20120112745A1 (en) * 2009-07-16 2012-05-10 Hitachi Medical Corporation Magnetic resonance imaging apparatus and magnetic resonance imaging method
US20130307537A1 (en) * 2011-01-28 2013-11-21 Hitachi Medical Corporation Magnetic resonance imaging device and angiography using same
CN102772210A (zh) * 2011-05-10 2012-11-14 西门子公司 弥散加权磁共振成像

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110730912A (zh) * 2017-06-09 2020-01-24 皇家飞利浦有限公司 并行多切片mr成像
CN110730912B (zh) * 2017-06-09 2022-06-14 皇家飞利浦有限公司 并行多切片mr成像
CN109814058A (zh) * 2017-11-22 2019-05-28 西门子保健有限责任公司 Caipirinha波space
CN109814058B (zh) * 2017-11-22 2021-09-17 西门子保健有限责任公司 磁共振系统和在其中生成检查对象的磁共振图像的方法
CN110389310A (zh) * 2018-04-16 2019-10-29 西门子医疗有限公司 用于产生磁共振记录的方法和控制装置
CN110412493A (zh) * 2018-04-26 2019-11-05 西门子医疗有限公司 记录磁共振数据组的方法、数据载体以及磁共振设备
US11016155B2 (en) 2018-04-26 2021-05-25 Siemens Healthcare Gmbh Method, apparatus and storage medium for recording a magnetic resonance dataset
US11543482B2 (en) * 2018-10-16 2023-01-03 Koninklijke Philips N.V. Magnetic resonance imaging using motion-compensated image reconstruction
CN110687488A (zh) * 2019-09-18 2020-01-14 上海联影医疗科技有限公司 磁共振扫描方法和磁共振成像装置
CN110687488B (zh) * 2019-09-18 2022-04-26 上海联影医疗科技股份有限公司 磁共振扫描方法和磁共振成像装置
CN111257809A (zh) * 2020-01-23 2020-06-09 上海东软医疗科技有限公司 磁共振成像方法、装置、存储介质及医疗设备
CN114041778A (zh) * 2022-01-13 2022-02-15 康达洲际医疗器械有限公司 一种应用于骨关节磁共振成像的融合成像方法与系统

Also Published As

Publication number Publication date
WO2015086415A1 (en) 2015-06-18
CN106030329B (zh) 2019-09-27
JP2016539722A (ja) 2016-12-22
US20160313431A1 (en) 2016-10-27
US10191133B2 (en) 2019-01-29
EP3080633A1 (en) 2016-10-19

Similar Documents

Publication Publication Date Title
CN106030329A (zh) 使用多回波分段k空间采集的MR成像
US10466322B2 (en) MR imaging using multi-echo k-space acquisition
CN103477238B (zh) 采用来自先验采集的约束的压缩感测mr图像重建
JP6998218B2 (ja) 動き検出を用いるｍｒ撮像
CN107750338A (zh) 使用星形叠层采集的mr成像
CN105074491B (zh) 具有使用压缩感知的图像重建的动态mri
JP6243522B2 (ja) 正則化された検出再構成を使用するマルチエコーディクソン水−脂肪分離及びｂ０歪補正による並列ｍｒｉ
JP5547800B2 (ja) 並列信号収集を利用したｍｒ撮像
US10401456B2 (en) Parallel MR imaging with Nyquist ghost correction for EPI
CN103502831B (zh) 对运动中目标的磁共振成像
US20150077112A1 (en) System, method and computer-accessible medium for highly-accelerated dynamic magnetic resonance imaging using golden-angle radial sampling and compressed sensing
CN105143906B (zh) 金属抗mr成像
US11937908B2 (en) Method and apparatus for accelerated magnetic resonance imaging
CN105531597B (zh) 耐金属的mr成像
CN105814449B (zh) 具有水/脂肪分离的零回波时间mr成像
CN109814058A (zh) Caipirinha波space
CN106133545B (zh) 利用对k空间中心的采样的零回波时间MR成像
CN107076818A (zh) 零回波时间mr成像
JP7446328B2 (ja) ソフト動きゲーティングによる３ｄラジアル又はスパイラル収集を用いたｍｒ画像
JP2021512671A (ja) 固有運動補正を有するスタック・オブ・スター収集を用いたｍｒ撮像
KR20160011012A (ko) 자기 공명 영상 장치 및 그 제어 방법
EP3640662A1 (en) Magnetic resonance imaging using motion-compensated image reconstruction
US20240069136A1 (en) Mr mammography
US20220229134A1 (en) Spiral mr imaging with off-resonance artefact correction
SPARSIFIED Fuleah A. Razzaq, Shady Mohamed, Asim Bhatti, Saeid Nahavandi

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