CN101263398A - 用于欠采样mri的约束逆投影重建方法 - Google Patents
用于欠采样mri的约束逆投影重建方法 Download PDFInfo
- Publication number
- CN101263398A CN101263398A CNA2006800246594A CN200680024659A CN101263398A CN 101263398 A CN101263398 A CN 101263398A CN A2006800246594 A CNA2006800246594 A CN A2006800246594A CN 200680024659 A CN200680024659 A CN 200680024659A CN 101263398 A CN101263398 A CN 101263398A
- Authority
- CN
- China
- Prior art keywords
- composograph
- projection
- space
- value
- picture frame
- 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
- 238000000034 method Methods 0.000 title claims abstract description 147
- 230000002093 peripheral effect Effects 0.000 claims description 16
- 230000000007 visual effect Effects 0.000 claims description 15
- 238000010606 normalization Methods 0.000 claims description 14
- 238000002360 preparation method Methods 0.000 claims description 11
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 9
- 230000004044 response Effects 0.000 claims description 7
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims 2
- 241001269238 Data Species 0.000 claims 2
- 230000035790 physiological processes and functions Effects 0.000 claims 1
- 239000002131 composite material Substances 0.000 abstract 2
- 230000008569 process Effects 0.000 description 74
- 238000005481 NMR spectroscopy Methods 0.000 description 31
- 238000003384 imaging method Methods 0.000 description 25
- 238000010586 diagram Methods 0.000 description 19
- 238000005070 sampling Methods 0.000 description 18
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 17
- 229910052704 radon Inorganic materials 0.000 description 12
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 12
- 230000008859 change Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 238000011160 research Methods 0.000 description 7
- 238000001914 filtration Methods 0.000 description 6
- 230000002792 vascular Effects 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 238000002583 angiography Methods 0.000 description 4
- 238000013500 data storage Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000012549 training Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 230000033001 locomotion Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 230000010287 polarization Effects 0.000 description 3
- 238000009877 rendering Methods 0.000 description 3
- 230000029058 respiratory gaseous exchange Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- NCGICGYLBXGBGN-UHFFFAOYSA-N 3-morpholin-4-yl-1-oxa-3-azonia-2-azanidacyclopent-3-en-5-imine;hydrochloride Chemical compound Cl.[N-]1OC(=N)C=[N+]1N1CCOCC1 NCGICGYLBXGBGN-UHFFFAOYSA-N 0.000 description 2
- 210000001367 artery Anatomy 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 229940079593 drug Drugs 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 238000007429 general method Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 210000005166 vasculature Anatomy 0.000 description 2
- 101100129500 Caenorhabditis elegans max-2 gene Proteins 0.000 description 1
- DGAQECJNVWCQMB-PUAWFVPOSA-M Ilexoside XXIX Chemical compound C[C@@H]1CC[C@@]2(CC[C@@]3(C(=CC[C@H]4[C@]3(CC[C@@H]5[C@@]4(CC[C@@H](C5(C)C)OS(=O)(=O)[O-])C)C)[C@@H]2[C@]1(C)O)C)C(=O)O[C@H]6[C@@H]([C@H]([C@@H]([C@H](O6)CO)O)O)O.[Na+] DGAQECJNVWCQMB-PUAWFVPOSA-M 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 210000003734 kidney Anatomy 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 229910052708 sodium Inorganic materials 0.000 description 1
- 239000011734 sodium Substances 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
在动态研究目标期间,采集两维或三维时间分辨MR帧图像。合成MR图像被生成并用于通过针对该图像帧采集的每个投影视图的逆投影(约束逆投影)进行加权来重建每个图像帧。该合成图像可根据单独采集的视图来重建,或者可通过组合在动态研究的过程期间采集的视图来生成。描述了该方法的大量不同临床应用。
Description
联邦资助研究的相关声明
本发明是在国家健康研究所(National Institute of Health)所赠予的基金No.1R01HL72260-01下由政府支持而做出的。美国政府对本发明拥有特定权利。
相关申请的交叉引用
本申请要求2005年7月8日提交并且题为“Backprojection ReconstructionMethod For Undersampled Time-Resolved MR Imaging(用于欠采样时间分辨MR成像的逆投影重建方法)”的美国临时专利申请No.60/697,607、以及2005年9月22日提交并且题为“Highly Constrained Image Reconstruction Method(高约束图像重建方法)”的美国临时专利申请No.60/719,445的优先权。
发明背景
本发明的领域是磁共振成像(“MRI”),尤其涉及MR图像的重建。
磁共振成像使用核磁共振(NMR)现象来生成图像。当诸如人体组织的物质处于均匀磁场(极化场B0)中时,组织中单独自旋(spin)磁矩试图与此极化场对齐,但是在其特征拉莫尔频率下围绕该极化场以随机次序进动。如果物质或组织处于在x-y平面中并接近拉莫尔频率的磁场(激发场B1)中,则净对齐力矩MZ可被旋转或“倾斜”到x-y平面以生成净横向磁矩Mt。由受激自旋发射一信号,并且在激发信号B1结束之后,此信号可被接收并处理以形成图像。
当使用这些信号来生成图像时,可使用磁场梯度(Gx、Gy和Gz)。通常,通过其中这些梯度根据所使用的具体定位方法变化的一系列测量周期来扫描要被成像的区域。在本领域中,每次测量称为“视图(veiw)”,并且视图的数量决定图像的分辨率。所得到的一组接收到的NMR信号、或视图、或k空间样本被数字化并处理以使用许多众所周知的重建技术之一来重建图像。总的扫描时间部分地由图像所需的测量周期、或视图的数量来确定,因此,扫描时间可在以图像分辨率或图像信噪比(“SNR”)为代价的前提下通过减少所需视图的数量来减少。
用于重建采集图像所用的NMR数据集的最普遍的方法称为“傅立叶变换”成像技术或“自旋卷绕(spin-warp)技术”。在由W.A.Edelstein等人在Physicsin Medicine and Biology(医学和生物学中的物理学)的第25卷第751-756页(1980)中的题为“Spin-Warp NMR Imaging and Application to HumanWhole-Body Imaging(自旋卷绕成像以及其对人体全身成像的应用)”的文章中讨论了这种技术。它在采集NMR信号之前使用可变幅相编码磁场梯度脉冲来在此梯度方向上对空间信息进行相位编码。例如在两维实现(2DFT)中,通过沿一个方向施加相位编码梯度(Gy)来在该方向上对空间信息进行编码,然后在与相位编码方向正交的方向上存在读出磁场梯度(Gx)的情况下采集信号。在自旋回波(spin-echo)的采集过程中出现的读出梯度在正交方向上对空间信息进行编码。在典型2DFT脉冲序列中,在扫描期间所采集的视图序列中相位编码梯度脉冲Gy的幅度递增(Gy)。在三维实现(3DFT)中,在每个信号被读出之前施加第三梯度(GZ),以沿第三轴进行相位编码。在扫描期间,此第二相位编码梯度脉冲Gz的幅度也按值来步进。这些2DFT和3DFT方法以诸如图2A中所示的直线图案对k空间进行采样,并且它们需要相当多的扫描时间以便充分地对k空间进行采样。
近来的大量研究使用了多个接收器线圈阵列来缩短成像扫描时间。在由Griswold等人在1999年6月的Magnetci Resonance In Medicine(医学中的磁共 振)的41(6):1235-45的“Simultaneous Acquisition Of Spatial Harmonics(SMASH)(空间谐波的同步采集(SMASH))”描述的SMASH技术中,多个线圈被谨慎地放置在傅立叶相位编码方向之一上。使用线圈灵敏度知识,可综合非采集相位编码,从而增加可获得给定分辨率的图像的速率,或者增加在同一速率下所获得的图像的分辨率。由Pruessmann等人在MRM 42:952-962(1999)的“CoilSensitivity Encoding For Fast MRI(用于快速MRI的线圈灵敏度解码)”中描述的SENSE技术是另一种减少扫描时间的多个接收信道方法。SMASH和SENSE方法由表示对于给定分辨率超出常规方法2到3个量级的速度增加的因子“R”来表征。它们还可由因子“g”来表征,该因子表示对于给定成像时间超过所期望的1-1.2量级的噪声增加。
近来的研究还使用了如美国专利No.6,487,435中公开的用于采集MRI数据的投影重建方法。投影重建方法从磁共振成像开始就已众所周知。与如傅立叶成像中所进行并在图2A中所示的以直线扫描图案对k空间采样不同,投影重建方法使用一系列视图对k空间进行采样,这些视图对如图2B中所示的从k空间的中心向外延伸的径向线进行采样。采样k空间所需的视图的数量决定了扫描的长度,并且如果未获得足够的视图数量,则在经重建的图像中产生条纹伪像。在专利No.6,487,435中公开的技术通过使用交错视图采集连续欠采样图像以及在连续图像之间共享***k空间数据来减少这种条纹。共享采集到的***k空间数据的方法在本领域中公知为缩写“TRICKS”。
存在例如在美国专利No.6,710,686中描述的、用于从采集到的k空间投影视图集重建图像的两种方法。最普遍的方法是将k空间样本从它们的在径向采样轨迹上的位置重新栅格化(regrid)到笛卡尔栅格。然后,通过以常规方式对经重新栅格化的k空间样本进行2D或3D傅立叶变换来重建图像。
用于重建图像的第二种方法是通过傅立叶变换每个投影视图来将径向k空间投影视图变换成Radon空间。通过滤波并将其逆投影到视场(FOV)来从这些信号投影重建图像。如本领域众所周知的,如果采集到的信号投影在数量上不足以满足奈奎斯特(Nyquist)采样定理,则将在经重建的图像中产生条纹伪像。
图3中示出了标准的逆投影方法。通过沿如箭头17所示的投影路径将轮廓11中的每个信号样本15通过FOV 13投影,来将每个Radon空间信号投影轮廓11逆投影到视场13上。在将每个信号样本15投影到FOV 13时,不具备学科的先验知识,并且作出以下假设:FOV 13中的NMR信号是均匀的并且信号样本15应当平均地分布在投影路径所穿过的每个像素中。例如,图3中示出了一个信号投影轮廓11中单个信号样本15的投影路径8,它穿过FOV 13中的N个像素。此信号样本15的信号值(P)以常规逆投影方式在此N个像素之间平均分割:
Sn=(P×1)/N (1)
其中Sn是分配到通过FOV 13的具有N个像素的投影路径上的第N个像素的NMR信号值。
显然,FOV 13中的NMR信号均匀的假设并不正确。然而,如本领域中众所周知的,如果对每个信号轮廓11作出特定校正,并且在相应数量的投影角下采集数量足够的轮廓,则将使由这种不完善的假设导致的误差将被最小化,并且图像伪像被抑制。在图像重建的典型滤波逆投影方法中,256×256个像素的2D图像需要400个投影,而256×256×256个像素的3D图像需要203,000个投影。如果使用上述美国专利No.6,487,435中描述的方法,则针对同样的图像所需的投影视图的数量可减小到100(2D)和2000(3D)。
由Tsao J.、Besinger P.、Pruessman KP在2003年11月的Magn.Reson.Med.的50(5):1031-43的“kt-Blast and k-t Sense:Dynamic MRI with High Frame RateExploiting Spatiotemporal Correlations(kt-Blast和k-tSense:采用时空相关的具有高帧速率的动态MRI)”、以及Hansen MS.,Tsao J.、Kozerke S.和Eggers H.在Miami Florida的2005年的ISMRM的文摘684页的“k-t Blast ReconstructionFrom Arbitrary k-t Sampling:Application to Dynamic Radial Imaging(根据任意k-t采样的k-t Blast重建:应用到动态径向成像)”公开的kt-blast技术认识到:在所采集的时间序列中,在与所采集的时间帧集相关联的k空间数据中存在大量相关。在已被应用到径向采集的kt-blast中,需要较低的空间频率训练数据集来消除当在空域和时域中执行欠采样时发生的混叠。使用迭代图像重建,可显著减少所需的数据。
Huang Y.、Gurr D.和Wright G.在Miami Florida的2005年的ISMRM的文摘1707页的“Time-Resolved 3D MR Angiography By Interleaved BiplaneProjections(根据交错双面投影的时间分辨的3D MR血管造影术)”描述了结合了使用训练数据集来引导使用正交2D投影图像的图像重建的思想的血管造影技术。在此方法中,迭代图像重建使用根据训练数据的相关数据分析来引导,该训练数据集由所有采集到的正交2D投影图像构成。
本发明的概述
本发明是一种用于重建磁共振图像的新方法,尤其是一种经改进的逆投影方法。合成图像被采集作为MRI扫描的一部分,并且对其进行重建来提供正被成像的目标的先验知识。在高度欠采样图像帧的重建期间,此合成图像用于对每个图像帧中的逆投影视图的分布进行加权。结果,需要采集更少的投影视图,从而导致更短的扫描时间。从10到100的速度增加因子是有可能的,这取决于图像的细节。
本发明的发现在于,如果FOV中的信号轮廓的先验知识被用在重建过程中,则可使用少得多的投影信号轮廓来生成高质量图像。参看图4,例如,可已知FOV 13中的信号轮廓包括诸如血管19和21的结构。情况是,当逆投影路径8穿过这些结构时,每个像素中的信号样本15的更精确的分布通过对作为该像素位置处的已知信号轮廓的函数的分布进行加权来实现。使用这种加权,在图4的示例中,大多数信号样本15将分布在与已知结构19和21相交的像素处。对于具有N个像素的逆投影路径8,这种高约束逆投影可表示如下:
其中:Sn=在正被重建的图像帧中的第n个像素处的逆投影信号的大小。
P=正被逆投影的投影剖面中的信号样本值;以及
Cn=在沿逆投影路径的第n个像素处的先验合成图像的信号值。合成图像根据在扫描期间所采集的数据来重建,并且可包括用于重建图像帧的数据以及描绘视图区域中的结构的其它采集到的数据。等式(2)中的分子使用合成图像中的相应信号值来对每个像素进行加权,而分母对该值进行归一化,使得所有逆投影信号样本反映该图形帧的投影总和,并且不与合成图像的总和相乘。应当注意:虽然可在执行了逆投影之后,分别对每个像素进行归一化,但是在许多临床应用中,在逆投影之前对投影P进行归一化要容易得多。在此情况中,投影P通过除以整个合成图像中处于同一视角的投影中的相应值Pc来归一化。接着,归一化投影P/Pc被逆投影,然后,所得图像与合成图像相乘。
图5中图示了由视角θ和φ表征的单个3D投影视图的高约束逆投影的3D实施例。此投影视图沿轴17逆投影并沿逆投影轴17在距离r处扩展成Radon平面25。与其中对投影信号值进行滤波并均匀地分布到连续的Radon平面的滤波逆投影不同,投影信号值使用合成图像中的信息来沿轴17分布在Radon平面25中。图5示例中的合成图像包含脉管19和21。基于合成图像中的相应位置x、y、z的强度,在Radon平面25的图像位置x、y、z处沉积加权信号轮廓值。这是逆投影信号轮廓值P与相应合成图像体素(voxel)值的简单乘法。然后,此乘积通过将该乘积除以来自由合成图像形成的对应图像空间投影轮廓的投影轮廓值来归一化。3D高约束重建的公式是:
I(x,y,z)=∑(P(r,θ,φ)*C(x,y,z)(r,θ,φ)/Pc(r,θ,φ) (2a)
其中总和(∑)覆盖正被重建的图像帧中的所有投影,并且特定Radon平面中的x、y、z值通过使用在该平面的适当r,θ,φ值处的投影轮廓值P(r,θ,φ)来计算。Pc(r,θ,φ)是来自合成图像的对应投影轮廓值,而C(x,y,z)(t,θ,φ)是在(r,θ,φ)的合成图像值。
本发明的另一个发现在于,存在大量临床MR应用,其中先验信息可用并且合成图像可被重建并用于加强欠采样图像的重建。当在动态研究中采集一系列时间分辨图像时,每个图像帧可使用极有限的采集视图集来重建。然而,每个这种视图集与对其它图像帧所采集的视图交错,在采集了大量图形帧之后,足够数量的不同视图可用于重建高质量合成图像以根据本发明来使用。
本发明的另一个方面是将高约束图像重建方法应用到对比增强磁共振血管造影术(“CEMRA”)。使用CEMRA,可在施加造影剂之前和之后采集目标脉管***的图像。前者图像用作从对比增强图像所扣除的掩模(mask),以从最后血管造影照片移除所有固定组织。使用本发明的图像重建方法,可以附加方式扣除掩模以提高其有效性。首先,在将合成图像用于重建最后图像之前,可从合成图像扣除掩模图像。或者,所采集的数据集中的每个k空间投影可使在掩模数据集中对应k空间投影在其高约束逆投影之前被扣除。或者,在同一重建期间,可使用上述两种掩模减法。
本发明的另一个方面是将高约束图像重建方法应用到其中采集一系列图像帧的动态研究中。在在这种动态研究期间,目标发生变化并且在整个研究中,单个合成图像可能未精确地示出该目标。例如,当造影剂流入感兴趣的脉管***中时,可采集CEMRA图像帧。为了更好地查看在研究期间发生的变化,多个合成图像用于重建图像帧。更具体地,用于重建图像帧的合成图像由用于重建图像帧本身的投影视图加上在周围时间窗中采集的交错投影视图构成。时间窗越窄,合成图像将在采集图像帧时的变化目标反映得越精确。
本发明的又一个方面是将高约束图像重建方法应用到其中一系列图像帧被采集并使用TRICKS视图共享方法重建的动态研究。在此应用中,可对k空间的中心部分和k空间的***部分生成独立的合成图像,以备在图像重建期间使用,或者可针对所有区域重建单个合成图像。
根据以下描述,本发明的前述以及其它目的和优点将变得显而易见。在描述中,对形成其一部分并作为说明而示出本发明的优选实施例的示例的附图进行参照。然而,此实施例不一定表示本发明的整个范围,并且由此对权利要求和本文作出参考以解释本发明的范围。
附图的简要描述
图1是其中使用本发明的MRI***的框图;
图2A是其中在使用图1的MRI***的典型傅立叶、或自旋卷绕图像采集期间对k空间进行采样的方式的图示;
图2B是其中在典型投影重建图像采集期间对k空间进行采样的方式的图示;
图3是图像重建过程中的常规逆投影步骤的图示表示;
图4是与根据本发明实现的相同的步骤的图示表示;
图5是根据本发明的逆投影步骤的3D实施例的图示表示;
图6A是用于引导图1的MRI***来采集3D投影重建图像的优选脉冲序列的图示;
图6B是用于实践本发明的另一优选脉冲序列的图示;
图6C是使用图6B的脉冲序列执行的k空间采样图案的图示表示;
图7是用于在MR成像应用中使用本发明的第一优选方法的流程图;
图8A和8B是用于在使用TRICKS视图共享方法的采集中使用本发明的另一优选方法的流程图;
图9是用在图8A和8B的TRICKS实施例中的k空间采样序列的图示表示;
图10是使用图8A和8B的TRICKS实施例生成的k空间数据集的图示表示;
图11是用在图8A和8B的方法中的合成图像的重新投影的图示表示;
图12是用于在CEMRA成像应用中使用本发明的优选方法的流程图;
图13是用在图12的CEMRA方法中的合成图像更新程序的图示表示;
图14是用于在CEMRA成像应用中使用本发明的又一方法的流程图;以及
图15是可用于从正被成像的FOV中消除非希望物体的流程图变体。
优选实施例的详细描述
具体参看图1,本发明的优选实施例被用在MRI***中。MRI***包括具有显示器112和键盘114的工作站110。工作站110包括处理器116,该处理器是运行商用操作***的商用可编程机器。工作站110提供使扫描命令能够被输入到MRI***中的操作员接口。
工作站110被耦合到四个服务器:脉冲序列服务器118;数据采集服务器120;数据处理器服务器122;以及数据存储服务器23。在优选实施例中,数据存储服务器123由工作站处理器116执行并与盘驱动器接口电路相关联。剩余三个服务器118、120和122由安装在单个外壳中独立的处理器来执行并使用64位底板总线互连。脉冲序列服务器118使用商用微处理器和商用四通道通信控制器。数据采集服务器120和数据处理器服务器122都使用同一商用微处理器,并且该数据处理服务器122还包括基于商用并行向量处理器的一个或多个阵列处理器。
工作站110以及服务器18、20和22的各个处理器被连接到串行通信网络。此串行通信网络传送从工作站110下载到服务器118、120和122的数据,并且它传送在服务器之间以及工作站与服务器之间传递的标签(tag)数据。另外,在数据处理器服务器122与工作站110之间提供高速数据链路,以便向数据存储服务器123传送图像数据。
脉冲序列服务器118响应于从工作站110下载的程序元件而发挥功能,以操作梯度***124和RF***126。需要执行规定扫描的梯度波形被生成并施加到梯度***,该***激励组件128中的线圈以生成用于位置编码NMR信号的磁场梯度Gx、Gy和Gz。梯度线圈组件128形成包括极化磁体132和整体RF线圈134的磁体组件130的一部分。
RF激励波形由RF***126施加到RF线圈134以执行规定磁共振脉冲序列。由RF线圈134检测到的响应NMR信号由RF***126在由脉冲序列服务器118生成的命令的指导下进行接收、放大、解调、滤波以及数字化。RF***126包括用于生成用在MR脉冲序列中的多种RF脉冲的RF发射器。该RF发射器响应于来自脉冲序列服务器118扫描命令和指导,以生成所需频率、相位和脉冲幅度波形的RF脉冲。将所生成的RF脉冲施加到整体RF线圈134,或者一个或多个局域线圈或线圈阵列。
RF***126还包括一个或多个RF接收器信道。每个RF接收器信道包括对与其连接的线圈接收到的NMR信号进行放大的RF放大器,以及检测和数字化接收到的NMR信号的I和Q正交分量的正交检测器。因而,根据I和Q分量的平方和的平方根来确定任何采样点处的接收到的NMR信号的大小:
而且接收到的NMR信号的相位也可被确定:
φ=tan-1Q/I
脉冲序列服务器118还任选地从生理采集控制器136接收病人数据。控制器136从连接到病人的许多不同传感器接收信号,诸如来自电极的ECG信号或来自风箱的呼吸信号。这些信号通常被脉冲序列服务器118用来同步或“选通(gate)”扫描执行与目标的呼吸或心搏。
脉冲序列服务器118还连接到从与病人的状况以及磁体***相关联的各种传感器接收信号的扫描室接口电路138。病人定位***140还通过扫描室接口电路138接收命令以便在扫描期间将病人移动到期望位置。
显而易见的是,脉冲序列服务器118在扫描期间对MRI***元件执行实时控制。结果,有必要使用以及时方式由运行时程序执行的程序指令操作其硬件元件。以对象的形式从工作站110下载扫描命令的描述组件。脉冲序列服务器118包含接收这些对象并将它们转换成可由运行时程序使用的对象的程序。
由RF***126生成的数字化NMR信号样本通过数据采集服务器120来接收。数据采集服务器120响应于从工作站110下载的描述组件进行操作,以接收实时NMR数据并提供缓冲存储,使得不会由于数据溢出而丢失数据。在一些扫描中,数据采集服务器120只将所采集的NMR数据传递到数据处理器服务器122。然而,在需要从所采集的NMR数据导出信息以控制扫描的进一步执行的扫描中,数据采集服务器120被编程为生成这种信息并将其传送到脉冲序列服务器118。例如,在预扫描期间,NMR数据被采集并用于校准由脉冲序列服务器118所实现的脉冲序列。而且,在扫描期间,导航信号可被采集并用于调节RF或梯度***操作参数,或者用于控制对k空间进行采样的视图次序。而且,数据采集服务器120可用于处理用来在MRA扫描中检测造影剂的到达的NMR信号。在所有这些是示例中,数据采集服务器120采集NMR数据,并对其进行实时处理以生成用于控制扫描的信息。
数据处理服务器122从数据采集服务器120接收NMR数据,并根据从工作站110下载的描述组件来对其进行处理。这些处理可包括例如:对原始k空间NMR数据进行傅立叶变换以生成两维或三维图像;对经重建的图像使用滤波器;对所采集的NMR数据执行逆投影图像重建;以及计算功能MR图像;计算运动或流图像等。如以下将更详细描述的,响应于由数据处理器服务器122执行的程序,通过MRI***来实现本发明。
由数据处理服务器122重建的图像被传送回工作站110,并存储其中。实时图像被存储在数据基本高速缓冲存储器(未示出),并可从该存储器向位于磁体组件130附近的操作员显示器112或显示器142输出以供主治医师使用。批处理模式图像或所选实时图像被存储在盘存储144的主机数据库中。当这些图像已被重建或转移到存储时,数据处理服务器122通知工作站110上的数据存储服务器123。工作站110可***作员用来存档图像、生成胶片或经由网络向其它设备发送图像。
为了实践本发明的某些优选实施例,在3D k空间球坐标系中采集NMR数据,并且读出梯度方向由来自kz轴的角θ和来自ky轴的角φ来定义。采样方法由一系列等距投影构成,所有投影都通过k空间的中心。最大k空间半径值(kmax)决定所得图像的所有三个空间方向上的分辨率。径向样本间隔(Δkr)决定经重建的图像的整个视场(FOV)的直径(D)。如果满足奈奎斯特条件Δk、Δk<Δkr,则可在无伪像的情况下重建整个FOV图像。然而,如果不满足这个条件,则在小于整个FOV(D)的减小直径(d)的范围内仍存在无混叠重建。如果假设以均匀的间隔(Δk=Δk=Δkr)采集投影,则在kmax处与投影相关联的表面积A为:
其中Np是所采集视图或投影的数量。等式(3)确定了Δk,籍此可将由于投影视图的角间距而减小的FOV的直径(d)与整个FOV直径D关联如下:
其中NR是FOV上的矩阵大小(即,在信号读出期间的样本数量)。在图像域中,良好构建的经减小的FOV呈现为以每个目标为中心,即使不满足奈奎斯特条件。然而,径向条纹伪像可从外部进入局部FOV。k空间被完全采样或d=D的条件要求被采样的投影的数量为:
例如,如果在每个所采集的NMR信号的读出期间采集NR=256个样本,则满足奈奎斯特条件所需的投影Np的数量约为103,000。
图6中示出了用于采集数据作为3D投影的脉冲序列。在装备有高性能梯度子***(40mT/m最大幅度以及150T/m/sec最大压摆率)的上述MRI***上实现该序列。可在数据采集窗200期间执行全回波或部分回波读出。如果选择部分回波,则仅部分地采集k空间的下半区(kz<0)。因为在所有方向上的较大FOV,所以非选择性射频(RF)脉冲202可用于生成贯穿图像FOV的横向磁化。
梯度回聚NMR回波信号203由受激FOV中的自旋生成,并在存在三个读出梯度206、208和210出现时进行采集。由于不需要平坦选择(slab-select)梯度,所以读出梯度波形Gx、Gy和Gz具有类似波形。这种对称仅在需要扰乱序列时才被打断,这通过实施移相梯度波瓣204来实现。Gx和Gy读出梯度208和210通过对应的梯度脉冲212和214而反绕(rewind)以实现稳态。
在扫描期间,读出梯度波形Gx、Gy和Gz被调制以在不同角采样径向轨迹。角间距被选择成在被采样k空间球体的***边界(kmax)处均匀分布的k空间采样点。尽管计算分布的若干方法是公知的,但是采用在恒定路径速度和表面积覆盖范围的条件下,通过以螺旋轨迹对球面进行采样的方法。这种解决方案还具有生成连续的采样路径的优点,这减小了梯度切换和涡流。对于总共N个投影,作为投影数量n的函数的梯度幅度的等式为:
如果要执行完全采样图像采集,则将N设置为如以上在等式(4)中定义的Np,并且执行一系列N=Np个脉冲序列。此系列中的第n个脉冲序列的读出梯度幅度由等式(5)、(6)和(7)给出。虽然在扫描期间n可从1到N以单调次序变化,但是应当理解,其它次序是可能的。如以下将描述的,本发明使得能够以少得多的投影视图来对球形k空间进行采样,这导致了更短的扫描时间。
图6B中示出了用于实践本发明的多个实施例的另一脉冲序列。这是快速梯度回聚回波脉冲序列,其中在出现平坦选择梯度232的情况下生成选择性不对称截短sinc射频激励脉冲218。射频脉冲218的翻转角(flip angle)被设置成接近通常为30°到40°的T1缩短的血液Ernst角。
此脉冲序列可用于通过在单个k空间圆形切片(slice)中进行采样来采集单个2D切片,或者可用于对如图6C中的234、236和238处所示的对多个圆形k空间切片进行采样。当采集多个2D切片时,径向梯度生成跟随有相位编码梯度波瓣240以及相反极性的反绕梯度波瓣242的平坦选择梯度232。在扫描期间,此径向相位编码梯度240在多个值中步进以从2D k空间切片234、236和238的每一个进行采样。显而易见的是,通过施加对应数量的不同径向相位编码,可使用此脉冲序列来对任何数量的2D k空间切片进行采样。
两个平面内读出梯度224和226在NMR回波信号228的采集期间放出,以便沿径向轨迹在2D平面234、236和238中采样k空间。这些平面内梯度224和226垂直于轴向梯度,并且它们彼此垂直。在扫描期间,它们在一系列值中步进以旋转径向采样轨迹的视角。平面内读出梯度的每一个落后于预定相梯度波瓣220和222,并跟随有反绕梯度波瓣224和226。对于每个轴向相位编码240,采集径向投影视图的完全集以采样2D k空间切片。如以下将描述的,本发明使得能够以更少的径向投影视图采样2D k空间切片,这导致了更短的扫描时间。
对于本领域技术人员而言,可使用与从k空间***边界上的一点穿过k空间的中心到达k空间***边界上的一相对点的优选直线轨迹不同的采样轨迹。如上所述,一种变体是采集部分NMR回波信号228,该回波信号沿不延伸跨越所采样的k空间体积的整个范围的轨迹进行采样。等效于直线投影重建脉冲序列的另一种变体是沿弯曲路径而非直线进行采样。例如,在F.E.Boada等人在1997年的MRM的37:706-715的“Fast Three Dimensional Sodium Imaging(快速三维钠成像)”、以及K.V.Koladia等人在Proc.Intl.Soc.Magn.Reson.Med.13(2005)的“Rapid 3D PC-MRA Using Spiral Projection Imaging(使用螺旋投影成像的快速3D PC-MRA)”和J.G.Pipe和Koladia在Proc.Intl.Soc.Mag.Reson.Med.13(2005)的“Spiral Projection Imaging:a new fast 3D trajectory(螺旋投影成像:新的快速3D轨迹)”中描述了这些脉冲序列。还应当显而易见的是,本发明可随这些采样方法的2D以及3D版本一起使用,并且本文对在如下文中所用术语“像素”的引用旨在指2D或3D图像中的位置。
本发明是用于根据所采集的k空间投影重建图像的经改进的方法。该方法需要重建FOV的合成图像,并且采集此合成图像的方式以及使用其的方式将取决于具体临床应用。现在将描述本发明的不同较佳实施例。
具体参看图7,第一较佳实施例可应用于其中几乎实时地采集时间分辨图像系列的临床情况。例如,这可用在干涉MR程序中。该程序中的第一步骤是采集并重建规定视场(FOV)中的目标的合成图像,如进程框300所示。此图像可使用任何MR脉冲序列来采集,但是通常其使用用于采集时间分辨图像的同一投影重建脉冲序列,在此实施例中,该序列是以上所述以及图6A中所示的脉冲序列。然而,因为不考虑时间分辨率,所以将采集许多投影视图,并且较佳地,采集足够的数据来满足奈奎斯特准则并生成不具有显著条纹伪像的所需分辨率的图像。
如果检查的目标响应于目标的呼吸或心动周期而移动,则可对合成图像的采集进行选通,如302所示。取决于正被检查的目标,可通过呼吸选通信号或心动选通信号或者两者来触发选通。使用常规重建方法从所采集的k空间数据重建合成图像。重建合成图像指示贯穿2d或3D FOV中每个图像像素处的NMR信号的大小。合成图像还可通过将不超出最小阈值大小的所有像素或体素设置成零来进行滤波。这些阈值滤波器使合成图像中的背景变暗,并使如下所述的随后使用该合成图像重建的时间分辨图像中的背景变暗。合成图像提供了关于扫描的目标的先验信息。
在合成图像被采集并重建之后,***进入其中几乎实时地采集和重建一系列帧图像的循环。更具体地,如进程框304所示,使用图6A的脉冲序列来采集图像帧投影集。对2D采集,这种采集可少至2-20个投影视图,而对3D采集可少至250到500个视图。结果,这得以极快地完成。如306所示,如果在合成图像的采集期间使用选通,则同一选通被用于采集帧图像。为了使帧图像中的条纹伪像最小化,所采集的N个投影视图应当如上所述地按角度隔开以均匀地采样k空间。
在采集图像帧投影之后,对它们进行运动补偿,如进程框308所示。合成图像用作目标的基准位置,并且对图像帧投影进行相位校正,以有效地参照合成图像中所示的基准位置来记录所示目标。用于记录两个图像的多种方法在本领域中是公知,并且在较佳实施例中,使用了由威斯康星(Wisconsin)大学在2002年发表的Oliver Wieben的题为“Novel Acquisition Strategies For TimeResolved 3D,Magnetic Resonance Angiography(用于时间分辨3D的磁共振血管造影术的新采集策略)”的博士论文中所述的方法。
如进程框310所示,下一步骤是通过对其进行傅立叶变换来将帧图像k空间投影变换到radon空间。结果是如图4中所示的一组信号轮廓11。如图7中的进程框312所示,然后将这些信号轮廓的每一个逆投影到如图4中由路径8所示的VOA中。此逆投影通过合成图像来加权,如上参照等式(2A)所述。即,任何像素(x,y,z)处的逆投影值(P)如上所述地归一化(P/Pc),并且通过经先前重建的合成图像中的同一像素的大小(C (x,y,z))进行加权。
如进程框314所示,然后,将逆投影信号值添加到正在重建的帧图像。然后,***在判定框316返回以便如进程框318和312所示对下一信号轮廓11进行逆投影。将所有逆投影信号轮廓11的信号值添加到帧图像,然后如进程框320所示,显示完成的帧图像。
如进程框324所示,附加图像帧以同一方式被采集、重建并显示直至进程结束,如判定框322所示。然而,可对每个图像帧采集不同的投影视图集。因此,对连续帧图像所采集的投影视图在k空间中被交错。在较佳实施例中,等式5-7用于生成采集经交错的投影所需的梯度场。
采集交错投影视图使得能够在图像重建方法中进行变化,其中在连续图像帧之间共享***k空间数据。上述美国专利No.6,487,435在常规图像重建的上下文中描述了这种视图共享概念。这种成像方法可用于使用图6B的脉冲序列的本发明,其中从如图6C中所示的中心区域A以及两个***区域B和C采集k空间数据。在较佳实施例中,在每个区域A、B和C中使用10到30kz相位编码,从而可对每个区重建对应数量的切片。
具体参看图8A,执行扫描来在投影视角交错的情况下采集一系列欠采样图像帧,从而可通过组合来自多个所采集的图像帧数据来形成合成图像。进入循环,其中如进程框400所示,标记(flag)被设置成“1”以指示正采集新的图像帧。如进程框402所示,然后执行图6B的脉冲序列以在中心区域A的每个kz相位编码处采集一个投影视角θ。然后,取决于标记,***分流到三条路径之一。在标记设置成“1”的情况下,***在区域B中的每个kz相位编码处采集一个投影视角θ,如进程框404所示,并且在区域C中的每个kz相位编码处采集一个投影视角θ,如进程框406所示。然后在408,标记被设置成“2”,并在410递增投影视角θ以采集下一投影角。然后,***返回进程框402以在中心区域A中采集新投影视角θ。
因为现在标记被设置成“2”,所以在此视角θ下仅对区域B进行采样,如进程框412所示,并且在414,标记被设置成“3”。然后,在进程框410,再次递增视角θ,并且***再次返回以便在进程框402从区域A采集k空间数据。因为现在标记被设置成“3”,这次在新的视角θ下从区域C采集投影视图,如进程框416所示,并且在进程框418,标记被设置回“2”。
显而易见的是,随着采集连续的视角,标记在“2”与“3”之间切换,结果以图9中所示的模式对k空间数据进行采样。如判定框420所确定的,当最后的视角θn已被采集时,一个图像帧已被采集,并且如进程框422所示,区域A、B和C的k空间数据被保存。在较佳实施例中,在每个图像帧采集期间,采集15个等距投影视角θ。
仍参照图8A,扫描以同一方式继续以采集连续的图像帧。然而,如进程框414所示,不采集相同的投影视角,相反,投影视角与已采集视角交错。因而,在采集了两个图像帧之后,已在每个区域A、B和C中已对2×15=30个不同投影视角进行采样,在已采集了三个图像帧之后,已对45个不同的投影视角进行采样,如此类推,直至在扫描结束时采集15×n个不同的交错投影角,如判定框426所确定的。
显而易见的是,扫描因两个因素而显著缩短。首先,对每个图像帧,仅采集15个投影视图,而非常规用于避免图像伪像的通常400到800个视图。另外,虽然在每个视角θ对k空间的中心区域A进行采样,但是每隔一投影视角θ才对***区域B和C进行采样,如图9中所示。这大致将扫描时间较小了三分之一,但是它也意味着错过k空间数据。如图8A中的进程框428所示,这通过对错过的投影视角内插k空间数据来解决。图9中对一个视角θ2进行说明,其中内插数据C2通过在相邻投影视图C1和C3中的对应的k空间样本点进行线性内插来生成。结果,如图10中,生成完整但高度欠采样图像帧k空间数据集。
参看图8B和10,下一步骤是沿kz执行一维傅立叶逆变换,如进程框430中所示。此变换解决信号沿z轴的位置,以定义沿z轴设置的多个2D切片。
然后,对每个2D断面重建合成图像,如进程框432所示。较佳地,来自所有采集到的图像帧的所有投影视图用于重建每个2D合成图像,以便生成可能实现的最佳图像,尽管并不使用所有图像帧。例如,如果采集十个图像帧,则可在每个2D合成图像的重建中使用总共10×15=150个交错投影视图。执行每个2D切片的常规图像重建,并且在较佳实施例中,这通过将投影视图k空间样本重新栅格化到2D笛卡尔栅格并沿每个轴执行常规两维傅立叶变换来实现。
可存在许多替代的合成图像重建方法。与对每个2D切片生成2D合成图像不同,可对整个体积重建单个3D合成图形,或者对每个区域A、B和C重建单独的3D合成图像。
具体参看图8B,现在使用由本发明示教的方法来重建每个所采集的图像帧中的每个2D切片。根据以下描述将变得明显的是,对每个图像帧重建一个2D切片图像,然后对2D切片的每一个重复该过程直至图像帧被完全重建。结果是每个所采集的图像帧的3D图像。
在此重建进程中的下一步骤是对当前2D切片的合成图像进行重新投影,如进程框434所示以及图11中所示。这是如在Jiang Hsieh在SPIE Press 2003的第三章的“Computed Tomography Principles,Design,Artifacts and RecentAdvances(计算X线断层摄影术的原理、设计、制品及最新进展)”中所述的常规Radon变换,并且对扫描期间所采集的每个视角θ生成合成图像投影Pc。然后,进入循环,其中重建一个帧图像中的一个2D切片,如436概括所示的。
如进程框438所示的,2D切片重建中的第一步骤是对所采集的2D切片的投影视图P进行归一化。这通过首先执行k空间投影视图到Radon空间的傅立叶变换,然后将其中的值除以同一视角θ下的合成图像投影Pc中的对应值来完成。然后,将结果归一化投影P/Pc逆投影到FOV中,如进程框440所示。将结果逆投影值添加到由对当前帧图像切片中的所有逆投影归一化投影视图求和所形成的无约束投影视图,如442所示。按判定框444所确定的,对所采集的投影视图(在此实施例中为15个)的每一个重复此进程,然后将形成的所得2D无约束切片图像与对应的合成图像2D切片相乘,如进程框446所示。这是矩阵乘法,其中无约束图像阵列中的像素值与合成图像阵列中的对应像素值相乘。所得重建切片图像被存储为当前图像帧的一部分,如进程框448所示。
对下一图像帧中的同一断面重复上述进程436,如进程框450所示。当已按判定框452所确定的对所采集的图像帧的每一个重建当前切片时,对每个图像帧重建下一切片,如进程框454所示。当每个图像帧的最后的2D切片被重建时,重建进程完成,如判定框456所检测到的。
使用更完全采样的合成图像的每个高度欠采样图像帧的约束重建导致更少的图像伪像,同时对每个图像帧保持较短的扫描时间。
本发明的另一临床应用是对比增强磁共振血管造影术(CEMRA)。在动态CEMRA研究中,在高时间速率下采集图像帧,以便观察流入感兴趣的脉管***的造影剂。结果,在研究期间的不同时刻,脉管亮度将得到增强,并且用于根据本发明的示教重建连续帧图像的合成图像将不保持恒定。换言之,在动态研究期间,用于重建帧图像的合成图像应当进行更新,以便适当地反映脉管逐渐填充造影剂并且在所采集的图像中变得更亮的事实。
具体参看图12,较佳CEMRA程序中的第一步骤是采集预对比(pre-contrast)掩模图像和初始合成图像,如进程框330所示。使用图6A的成像脉冲序列,并且对掩模图像采集足够的视图以满足奈奎斯特准则。然后,如进程框332所示,造影剂被注入到目标中,并且进入循环,在其中尽可能快地采集帧图像。应当理解,如果在造影剂流入FOV之前进行采集,则可在造影剂施加之后采集掩模。
如进程框334所示,采集一个帧图像的k空间投影视图,并且如进程框336所示,重建该帧图像。如上述实施例中所示,选择一个帧图像的投影视图以尽可能均匀地对k空间进行采样,并且所采集视图的数量由规定的时间分辨率来确定。当然,本发明的优点是可减少视图的数量以实现较高的时间分辨率而不丧失图像分辨率并且不增加条纹伪像。图像帧重建336与以上参照图7中的进程框310、312和314所述的相同,其中合成图像中的先验信息用于对逆投影进行高度约束。
在重建图像帧之后,从其扣除预对比掩模图像,如进程框338所示。如此进行以从稍后被显示的所得帧图像中移除非脉管组织,如进程框340所示。
如果按判定框342确定的将要采集附加帧图像,则首先更新合成图像,如进程框344所示。如上所述,在动态研究期间,正被成像的目标脉管***持续变化,并且目的是尽可能实时地保持合成图像,从而对下一图像帧实现更精确的高约束逆投影。图13中示出了此合成更新步骤,其中框345-352表示在动态研究期间进行的连续图像帧采集。例如,如果要重建的下一图像帧由框349来表示,则已更新的合成图像通过对n=3个先前采集的图像帧与当前图像帧进行组合来形成。更具体地,针对n个先前图像帧加上当前图像帧的交错k空间投影视图被组合以形成单个k空间图像,如进程框352所示。从此组合数据集重建已更新的合成图像,如进程框354所示。如上所示,这是常规图像重建进程。在以上参照进程框312所述的随后图像帧的逆投影349中,使用已更新的合成图像。
仍参照图13,因而通过n个先前所采集的图像帧和当前图像帧的窗口形成已更新合成图像,该当前图像帧最精确地反映了正被检查的目标的当前状态。当目标中的变化相对较慢时,可增加n的值以包括更多先前所采集的图像帧。所得的更大量的投影视图改进了所得已更新合成图像的质量。另一方面,当快速变化时,可将n减小到少至n=1个图像帧,以便正确地反映正被检查的目标的真实状态。因此,在对n的选择导致的高SNR方面与更精确的动态变化描绘的另一方面之间存在平衡。
如果在动态扫描完成之后重建帧图像,则用于更新合成图像的所采集的图像帧窗口可扩展成包括在当前图像帧之后采集的图像帧。例如,正被重建的图像帧可处于窗口的中间,并且当前图像帧之前和之后所采集的其它图像帧的数量基本上相等。或者,可在窗口的起始处采集当前图像帧。在这种所采集的图像帧的后处理中,可重建许多不同的图像帧,其中可改变窗口大小和窗口相对于当前图像帧的定位以实现最佳结果。
再次参照图12,在合成图像被更新之后,***返回以采集下一帧图像,如进程框360所示。因为已更新的合成图像通过组合先前所采集的图像帧来形成,所以当来自n个连续图像帧的k空间数据被组合时,从一个帧到下一帧的投影视图应当如上所述地交错,从而基本上均匀地对k空间进行采样。
在图12中所示以及上述的CEMRA方法中,当执行动态研究时,几乎实时地重建并显示图像帧。图14中示出了替代CEMRA方法,其中在执行图像重建之前,在进程框380采集所有图像帧。在此实施例中,在进程框332处注入造影剂之前,在进程框330处预先采集掩模投影,并从对应的注入后k空间投影扣除该掩模投影,如进程框382所示。大量经扣除的投影可用于重建合成图像,如进程框384所示。使用标准图像重建方法,因为可从研究获得更完整的径向投影集。然而,为了使合成图像随图像帧与在进程386重建的图像帧保持更新,仅与所选时间帧同时采集的经扣除的投影通常用于形成合成图像,如图13所示。在进程框386通过将所选择的采集到的投影视图傅立叶变换到Radon空间并使用上述合成图像执行高约束逆投影来重建图像帧。
如进程框388和390所示,可通过更新从经扣除的投影产生的合成图像并使用正确选择的投影重复逆投影来重建附加图像帧。
本发明该实施例的优点是在执行高约束逆投影之前,从所采集的图像数据扣除掩模图像数据。这从FOV移除了没有临床价值的许多结构并导致“稀疏”数据集,该数据集使得逆投影过程更精确地聚焦在感兴趣的临床结构上。而且,因为在图像重建之前采集所有数据,所以可使用在采集当前图像帧之前和之后所采集的投影视图的窗口来更新合成图像。这使得更多投影视图能被组合。
尽管上述实施例的目的是使用尽可能精确地描述正被成像的目标的合成图像,但是存在其中有意改变合成图像是有益的临床情况。例如,一种此情况是当进行其中主动脉呈现在视场中的肾脏动脉的CEMRA研究时。来自大动脉的优势NMR信号可生成难以移除的条纹伪像。在此示例中,期望抑制该信号,因为其没有临床价值。这可通过改变用于对帧图像进行逆投影的合成图像来用本发明实现。
图15中示出了本发明中对合成图像作出这种改变以滤除视场中的目标物体的一实施例。此实施例是增强预先采集的欠采样及时间分辨图像帧的质量的后处理方法。更具体地,在目标的研究期间,扫描包括采集一系列图像帧投影,如进程框364所示。在扫描完成之后,此所采集的k空间投影数据被存储并处理。
如由进程框366所示的第一后处理步骤是通过从多个所采集的图像帧投影集组合交错k空间投影视图来重建一个或多个合成图像。这是常规图像重建,并且之后显示所得的合成图像,如进程框367所示。如进程框368所示,操作员编辑合成图像。在上述示例临床应用中,操作员可圈出描述大动脉的像素,并将它们的值设置成零。换言之,从合成图像移除大动脉。为了进一步抑制非期望物体,所圈出的待删除像素也被傅立叶变换回k空间,并从每个图像帧的对应k空间投影中扣除。
然后,使用这种经改变的合成图像和经改变的k空间投影来重建图像帧。如进程框370所示,一组图像帧k空间投影被傅立叶变换到Radon空间,并且使用本发明的经编辑的合成图像和约束逆投影来进行逆投影,以生成如进程框372所示进行显示的帧图像。此步骤与以上参照图7中的进程框310、312和314描述的重建相同,并且按判定框374所确定的对所有采集到的图像帧重复该步骤。显而易见的是,为了适应在研究期间目标中的突然变化,多个合成图像可被创建并用于重建如上所述系列中的帧图像。在此情况中,编辑每个已更新的合成图像以便在将其用于逆投影步骤之前移除所顾虑的物体。
本发明的又一实施例生成相位对比MRA图像。在此实施例中,所用的脉冲序列使用如题为“Phase Contrast Imaging Using Interleaved Projection Data(使用交错投影数据的相位对比成像)”的美国专利No.6,188,922所述的附加运动编码梯度,并且常常不使用造影剂。在使用同时采集的投影来形成合成图像的情况下,被重建的相位图像使用上述逆投影方法。
Claims (21)
1.一种用于对被定位在磁共振成像(MRI)***的视场(FOV)中的目标生成图像的方法,包括如下步骤:
a)使用所述MRI***对被定位在所述FOV中的所述目标采集投影视图集;
b)使用所述MRI***生成指示被定位在所述FOV中的所述目标的每个合成图像的像素处的值的合成图像;和
c)重建所述目标的图像,通过:
c)i)将所述集合中的每个投影视图逆投影到所述FOV并且通过所述合成图像中的对应像素的所述值对逆投影到每个图像像素的所述值加权;以及
c)ii)对每个图像像素的所述逆投影值进行求和。
2.如权利要求1所述的方法,其特征在于,如下计算步骤c)i)中的每个图像像素的逆投影值Sn:
其中:P=正被逆投影的所述投影视图值;
Cn=所述合成图像中的对应像素值;
Sn=沿逆投影路径的第n个像素的所述值;以及
N=沿所述逆投影路径的像素的总数。
3.如权利要求1所述的方法,其特征在于,步骤b)包括编辑所述合成图像以移除其中的物体以及由此基本上最小化该物体在所述重建图像中的出现。
4.如权利要求1所述的方法,其特征在于,步骤c)i)中的所述加权包括使用来自所述合成图像的对应投影视图对每个投影视图进行归一化,并且将所述逆投影值与所述合成图像中的所述对应像素的所述值相乘。
5.如权利要求4所述的方法,其特征在于,还包括:
d)在所述图像系列的重建期间周期性地更新所述合成图像以在其中描绘在检查期间所述目标中发生的变化。
6.如权利要求5所述的方法,其特征在于,所述合成图像的所述更新包括使用步骤a)中所采集的投影视图来重建所述合成图像。
7.如权利要求1所述的方法,其特征在于,步骤a)包括采集所述目标的k空间投影视图,而步骤c)包括对所述k空间投影视图进行傅立叶变换。
8.如权利要求1所述的方法,其特征在于,响应于表示所述目标中选定的生理事件的选通信号来执行步骤a)。
9.如权利要求1所述的方法,其特征在于,所述FOV是三维的,生成三维图像,并且步骤c)中重建的所述图像I(x,y,z)为:
I(x,y,z)=∑(P(r,θ,φ)*C(x,y,z)(r,θ,φ)/Pc(r,θ,φ)
其中求和(∑)覆盖所采集的集合中的所有投影视图;
I(x,y,z)是在FOV像素位置x,y,z处的所述图像值;P(r,θ,φ)是来自视角θ,φ处采集的所述视图的所述逆投影值;C(x,y,z)是在像素位置x,y,z处的所述合成图像值;以及Pc(r,θ,φ)是从所述视角θ,φ处的所述合成图像投影的轮廓值。
10.如权利要求1所述的方法,其特征在于,包括:
d)在施加造影剂之前,采集在其每个图像像素处描绘被定位在所述FOV中的所述目标的掩模图像;
e)在执行步骤a)和b)之前,将造影剂施加到所述目标;以及
f)在执行步骤c)之前,从所述合成图像扣除所述掩模图像。
11.如权利要求10所述的方法,其特征在于,所述掩模图像和所述合成图像被采集成投影视图集,并且步骤f)通过从所述合成图像集中的对应投影视图扣除所述掩模图像集中的投影视图来执行。
12.如权利要求10所述的方法,其特征在于,步骤f)通过从所述合成图像中的对应像素扣除所述掩模图像中的像素来执行。
13.如权利要求1所述的方法,其特征在于,步骤c)包括在执行步骤c)i)之前对每个投影视图进行傅立叶变换。
14.如权利要求1所述的方法,其特征在于,包括:
d)在用于执行步骤a)的视角下重新投影所述合成图像;以及
其中步骤c)i)中的所述加权包括通过将其值除以同一视角下的所述合成图像的所述投影视图中的对应值来归一化每个投影视图。
15.一种用于对被定位在磁共振成像(MRI)***的视场(FOV)中的目标生成一系列图像帧的方法,包括如下步骤:
a)采集一系列图像帧k空间数据集,每个图像帧k空间数据集包括来自k空间的中心区域的k空间样本以及来自多个***k空间区域的k空间样本,并且其中每个图像帧k空间数据集中的所述k空间样本与其它图像帧k空间数据集中的所述k空间样本交错;
b)根据来自多个所述图像帧k空间数据集的中心区域k空间样本来重建中心区域合成图像;
c)根据来自多个所述图像帧k空间数据集中的所述***k空间区域之一的k空间样本来重建第一***区域合成图像;
d)根据来自多个所述图像帧k空间数据集中的所述***k空间区域的另一个的k空间样本来重建第二***区域合成图像;以及
e)使用来自步骤b)、c)和d)的所述合成图像中的先验信息来根据其对应的k空间数据集重建所述图像帧系列的每一个,从而加权所述图像帧中的所述像素值。
16.如权利要求15所述的方法,其特征在于,采集所述k空间样本作为投影视图,并且步骤e)包括:
e)i)对图像帧k空间数据集中的所述投影视图进行傅立叶变换;以及
e)ii)将每个经傅立叶变换的投影视图逆投影到所述FOV中并通过所述合成图像之一中的所述对应像素的所述值对逆投影到每个图像帧像素的所述值加权。
17.如权利要求16所述的方法,其特征在于,步骤e)ii)中的所述加权包括使用来自所述一合成图像的对应投影视图来归一化每个经傅立叶变换的投影视图,并且将所述逆投影值与所述一合成图像中的所述对应像素的所述值相乘。
18.如权利要求15所述的方法,其特征在于,k空间的每个区域包括多个两维切片,并且每个区域合成图像包括对应的多个两维合成图像。
19.一种用于对被定位在磁共振成像(MRI)***的视场(FOV)中的目标生成一系列图像帧的方法,包括如下步骤:
a)采集一系列不完全图像帧k空间数据集,每个图像帧k空间数据集包括来自k空间的中心区域的k空间样本以及来自多个***k空间区域的k空间样本,并且其中每个图像帧k空间数据集中的所述k空间样本与其它图像帧k空间数据集中的所述k空间样本交错;
b)通过根据所采集的***区域k空间数据计算***区域k空间数据来生成一系列对应的完全图像帧k空间数据集;
c)根据来自多个所述完全图像帧k空间数据集的k空间样本来重建合成图像;
d)使用所述合成图像中的先验信息来根据其对应的完全图像帧k空间数据集重建所述图像帧系列的每一个,从而加权所述图像帧中的所述像素值。
20.如权利要求19所述的方法,其特征在于,采集所述k空间样本作为投影视图,并且步骤d)包括:
d)i)对完全图像帧k空间数据集中的所述投影视图进行傅立叶变换;以及
d)ii)将每个经傅立叶变换的投影视图逆投影到所述FOV中并通过所述合成图像中的所述对应像素的所述值对逆投影到每个图像帧像素的所述值加权。
21.如权利要求20所述的方法,其特征在于,步骤d)ii)中的所述加权包括使用来自所述合成图像的对应投影视图来归一化每个经傅立叶变换的投影视图,并且将所述逆投影值与所述合成图像中的所述对应像素的所述值相乘。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US69760705P | 2005-07-08 | 2005-07-08 | |
US60/697,607 | 2005-07-08 | ||
US71944505P | 2005-09-22 | 2005-09-22 | |
US60/719,445 | 2005-09-22 | ||
PCT/US2006/026201 WO2007008532A1 (en) | 2005-07-08 | 2006-07-06 | Constrained backprojection reconstruction method for undersampled mri |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101263398A true CN101263398A (zh) | 2008-09-10 |
CN101263398B CN101263398B (zh) | 2012-07-18 |
Family
ID=39947938
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2006800246594A Active CN101263398B (zh) | 2005-07-08 | 2006-07-06 | 用于欠采样mri的约束逆投影重建方法 |
CN200680024410.3A Expired - Fee Related CN101300600B (zh) | 2005-07-08 | 2006-07-06 | 用于ct成像的反投影重构方法 |
CN2006800155580A Active CN101248458B (zh) | 2005-07-08 | 2006-07-06 | 高约束图像重构方法 |
Family Applications After (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200680024410.3A Expired - Fee Related CN101300600B (zh) | 2005-07-08 | 2006-07-06 | 用于ct成像的反投影重构方法 |
CN2006800155580A Active CN101248458B (zh) | 2005-07-08 | 2006-07-06 | 高约束图像重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (3) | CN101263398B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102085097A (zh) * | 2010-12-16 | 2011-06-08 | 中国科学院深圳先进技术研究院 | 磁共振动态成像方法 |
CN102565738A (zh) * | 2010-11-25 | 2012-07-11 | 西门子公司 | 对拍摄的磁共振数据的补充 |
CN102906791A (zh) * | 2010-05-19 | 2013-01-30 | 威斯康星校友研究基金会 | 使用在先图像约束的图像重建的辐射计量减少的方法 |
CN104081437A (zh) * | 2011-09-20 | 2014-10-01 | 西门子公司 | 用于基于可调谐二极管激光吸收光谱的气体浓度重建的贝叶斯方法 |
CN104903740A (zh) * | 2012-12-18 | 2015-09-09 | 皇家飞利浦有限公司 | 对极化磁场的基于体模的mr场映射 |
CN105022009A (zh) * | 2014-02-20 | 2015-11-04 | 西门子公司 | 采集二维体积片段的磁共振数据的方法以及磁共振设备 |
CN105259525A (zh) * | 2015-10-28 | 2016-01-20 | 北京大学 | 基于邻域共享压缩感知的动态对比增强磁共振快速成像方法 |
CN105980875A (zh) * | 2013-10-08 | 2016-09-28 | 皇家飞利浦有限公司 | 经校正的多切片磁共振成像 |
CN109727297A (zh) * | 2018-12-29 | 2019-05-07 | 上海联影医疗科技有限公司 | 医学图像重建方法、***、可读存储介质和设备 |
CN110121293A (zh) * | 2017-11-24 | 2019-08-13 | 佳能医疗***株式会社 | 医用数据处理装置、磁共振成像装置及学习完毕模型生成方法 |
CN110333466A (zh) * | 2019-06-19 | 2019-10-15 | 东软医疗***股份有限公司 | 一种基于神经网络的磁共振成像方法和装置 |
CN111157935A (zh) * | 2019-12-31 | 2020-05-15 | 上海联影智能医疗科技有限公司 | 磁共振成像方法、装置、存储介质及计算机设备 |
CN111373445A (zh) * | 2017-10-23 | 2020-07-03 | 皇家飞利浦有限公司 | 对具有重叠的全身正电子发射断层摄影(pet)扫描重建图像,并且改变个体床位的曝光时间 |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8180130B2 (en) * | 2009-11-25 | 2012-05-15 | Imaging Sciences International Llc | Method for X-ray marker localization in 3D space in the presence of motion |
EP2632337A1 (en) * | 2010-10-27 | 2013-09-04 | Koninklijke Philips Electronics N.V. | Under-sampled, multi-energy computed tomography (ct) data acquisition data processing |
WO2013121379A2 (en) * | 2012-02-16 | 2013-08-22 | Koninklijke Philips N.V. | Spatially corrected nuclear image reconstruction |
CN103284749B (zh) * | 2012-02-24 | 2015-09-09 | 株式会社东芝 | 医用图像处理装置 |
CN104812306A (zh) * | 2012-11-29 | 2015-07-29 | 株式会社东芝 | 医用信息处理装置、医用图像诊断装置以及医用信息处理方法 |
CN103793890A (zh) * | 2014-03-05 | 2014-05-14 | 南方医科大学 | 一种能谱ct图像的恢复处理方法 |
CN105334479B (zh) * | 2014-06-10 | 2019-02-22 | 西门子(深圳)磁共振有限公司 | 一种磁共振成像方法和装置 |
EP2995254A1 (en) * | 2014-09-12 | 2016-03-16 | Koninklijke Philips N.V. | Imaging system |
US10456116B2 (en) | 2014-09-30 | 2019-10-29 | Siemens Medical Solutions Usa, Inc. | Shadow suppression in ultrasound imaging |
CN104751429B (zh) * | 2015-01-27 | 2018-02-02 | 南方医科大学 | 一种基于字典学习的低剂量能谱ct图像处理方法 |
DE102015208905A1 (de) | 2015-05-13 | 2016-11-17 | Siemens Healthcare Gmbh | Verfahren zum Erzeugen eines Bildes |
US10331936B2 (en) * | 2015-07-14 | 2019-06-25 | Idex Asa | Duplicate pattern reconstructions |
DE102016219709B3 (de) * | 2016-10-11 | 2018-03-01 | Siemens Healthcare Gmbh | Verfahren zur Ermittlung eines Perfusionsdatensatzes, sowie Röntgenvorrichtung, Computerprogramm und elektronisch lesbarer Datenträger |
US10624599B2 (en) * | 2017-07-31 | 2020-04-21 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for motion signal correction |
CN108491583B (zh) * | 2018-02-28 | 2021-09-10 | 哈尔滨工程大学 | 一种不规则形状放射源的自适应辐射剂量计算仿真方法 |
CN108957515B (zh) * | 2018-09-18 | 2020-09-08 | 上海联影医疗科技有限公司 | 确定探测器的能量响应函数的方法、装置和成像*** |
CN112233157B (zh) * | 2020-10-14 | 2022-03-08 | 俐玛精密测量技术(苏州)有限公司 | 微纳米ct***转台误差的频域校正方法 |
CN113466279B (zh) * | 2021-06-30 | 2024-03-29 | 中国科学技术大学 | 宽场超分辨自旋磁成像装置及方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5522390A (en) * | 1991-11-21 | 1996-06-04 | U.S. Philips Corporation | Magnetic resonance imaging method |
EP0627633A1 (en) * | 1993-05-18 | 1994-12-07 | Koninklijke Philips Electronics N.V. | Method and apparatus for magnetic resonance imaging |
DE4319539A1 (de) * | 1993-06-12 | 1994-12-15 | Philips Patentverwaltung | Verfahren zur Erzeugung einer MR-Bildfolge und Anordnung zur Durchführung des Verfahrens |
DE19647537A1 (de) * | 1996-11-16 | 1998-05-20 | Philips Patentverwaltung | MR-Verfahren zur Reduzierung von Bewegungsartefakten und Anordnung zur Durchführung des Verfahrens |
JP4008193B2 (ja) * | 2000-08-29 | 2007-11-14 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像処理方法、画像処理装置、mri装置および記録媒体 |
JP4130055B2 (ja) * | 2000-08-31 | 2008-08-06 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 加算断層画像作成方法およびx線ct装置 |
JP4175791B2 (ja) * | 2001-08-20 | 2008-11-05 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像生成方法およびx線ct装置 |
CN100403985C (zh) * | 2004-04-09 | 2008-07-23 | 清华大学 | 数字化近红外光漫射断层成像*** |
-
2006
- 2006-07-06 CN CN2006800246594A patent/CN101263398B/zh active Active
- 2006-07-06 CN CN200680024410.3A patent/CN101300600B/zh not_active Expired - Fee Related
- 2006-07-06 CN CN2006800155580A patent/CN101248458B/zh active Active
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102906791A (zh) * | 2010-05-19 | 2013-01-30 | 威斯康星校友研究基金会 | 使用在先图像约束的图像重建的辐射计量减少的方法 |
CN102906791B (zh) * | 2010-05-19 | 2015-02-11 | 威斯康星校友研究基金会 | 使用在先图像约束的图像重建的辐射计量减少的方法 |
CN102565738A (zh) * | 2010-11-25 | 2012-07-11 | 西门子公司 | 对拍摄的磁共振数据的补充 |
CN102565738B (zh) * | 2010-11-25 | 2015-10-28 | 西门子公司 | 对拍摄的磁共振数据的补充 |
CN102085097A (zh) * | 2010-12-16 | 2011-06-08 | 中国科学院深圳先进技术研究院 | 磁共振动态成像方法 |
CN102085097B (zh) * | 2010-12-16 | 2012-07-25 | 中国科学院深圳先进技术研究院 | 磁共振动态成像方法 |
CN104081437A (zh) * | 2011-09-20 | 2014-10-01 | 西门子公司 | 用于基于可调谐二极管激光吸收光谱的气体浓度重建的贝叶斯方法 |
CN104081437B (zh) * | 2011-09-20 | 2017-02-15 | 西门子公司 | 用于基于可调谐二极管激光吸收光谱的气体浓度重建的贝叶斯方法 |
CN104903740A (zh) * | 2012-12-18 | 2015-09-09 | 皇家飞利浦有限公司 | 对极化磁场的基于体模的mr场映射 |
CN105980875A (zh) * | 2013-10-08 | 2016-09-28 | 皇家飞利浦有限公司 | 经校正的多切片磁共振成像 |
CN105980875B (zh) * | 2013-10-08 | 2019-03-08 | 皇家飞利浦有限公司 | 经校正的多切片磁共振成像 |
CN105022009A (zh) * | 2014-02-20 | 2015-11-04 | 西门子公司 | 采集二维体积片段的磁共振数据的方法以及磁共振设备 |
CN105259525A (zh) * | 2015-10-28 | 2016-01-20 | 北京大学 | 基于邻域共享压缩感知的动态对比增强磁共振快速成像方法 |
CN111373445A (zh) * | 2017-10-23 | 2020-07-03 | 皇家飞利浦有限公司 | 对具有重叠的全身正电子发射断层摄影(pet)扫描重建图像,并且改变个体床位的曝光时间 |
CN110121293A (zh) * | 2017-11-24 | 2019-08-13 | 佳能医疗***株式会社 | 医用数据处理装置、磁共振成像装置及学习完毕模型生成方法 |
CN109727297A (zh) * | 2018-12-29 | 2019-05-07 | 上海联影医疗科技有限公司 | 医学图像重建方法、***、可读存储介质和设备 |
CN109727297B (zh) * | 2018-12-29 | 2023-08-08 | 上海联影医疗科技股份有限公司 | 医学图像重建方法、***、可读存储介质和设备 |
CN110333466A (zh) * | 2019-06-19 | 2019-10-15 | 东软医疗***股份有限公司 | 一种基于神经网络的磁共振成像方法和装置 |
CN111157935A (zh) * | 2019-12-31 | 2020-05-15 | 上海联影智能医疗科技有限公司 | 磁共振成像方法、装置、存储介质及计算机设备 |
CN111157935B (zh) * | 2019-12-31 | 2022-04-26 | 上海联影智能医疗科技有限公司 | 磁共振成像方法、装置、存储介质及计算机设备 |
Also Published As
Publication number | Publication date |
---|---|
CN101263398B (zh) | 2012-07-18 |
CN101248458B (zh) | 2012-08-15 |
CN101300600A (zh) | 2008-11-05 |
CN101300600B (zh) | 2016-01-20 |
CN101248458A (zh) | 2008-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101263398B (zh) | 用于欠采样mri的约束逆投影重建方法 | |
CN101573629B (zh) | 利用相衬合成图像进行高度约束反投影重建的对比增强磁共振血管造影 | |
JP6998218B2 (ja) | 動き検出を用いるmr撮像 | |
EP1902328B1 (en) | Constrained backprojection reconstruction method for undersampled mri | |
EP1444530B1 (en) | Three-dimensional phase contrast magnetic resonance imaging using interleaved projection-reconstruction data | |
US7245124B2 (en) | Under-sampled 3D MRI using a shells k-space sampling trajectory | |
US7865227B2 (en) | Image reconstruction method for cardiac gated magnetic resonance imaging | |
CN101297214B (zh) | 心脏选通mri中的高度约束的背投重构过程 | |
US7408347B2 (en) | Highly constrained magnetic resonance spectroscopy image reconstruction method | |
US8744551B2 (en) | Method for non-contrast enhanced magnetic resonance angiography | |
US8155419B2 (en) | MRI acquisition using sense and highly undersampled fourier space sampling | |
US7457655B2 (en) | Motion correction of magnetic resonance images using moments of spatial projections | |
JP2005503874A (ja) | 浮動式テーブル投影イメージングを使用した磁気共鳴血管造影法 | |
EP2699926A1 (en) | Spatially encoded phase-contrast mri | |
EP2350698B1 (en) | System and method for moving table mri | |
EP4057022A1 (en) | System and method for magnetization-prepared three-dimensional unbalanced steady-state free precession magnetic resonance imaging | |
US7167740B2 (en) | Measuring method in magnetic resonance imaging device and magnetic resonance imaging device | |
JP2009508654A (ja) | 高度に限定された磁気共鳴スペクトロスコピー画像再構成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 1124666 Country of ref document: HK |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |