CN104181483A - 一种基于定-变角优化激发的超极化磁共振成像方法 - Google Patents
一种基于定-变角优化激发的超极化磁共振成像方法 Download PDFInfo
- Publication number
- CN104181483A CN104181483A CN201410438269.8A CN201410438269A CN104181483A CN 104181483 A CN104181483 A CN 104181483A CN 201410438269 A CN201410438269 A CN 201410438269A CN 104181483 A CN104181483 A CN 104181483A
- Authority
- CN
- China
- Prior art keywords
- angle
- excites
- hyperpolarization
- excite
- excitation
- 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
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种基于定-变角优化激发的超极化磁共振成像方法,涉及磁共振成像领域。本发明综合定角激发和变角激发的优点,将激发分为定角激发和变角激发两个阶段,即前n*次激发采用定角激发方法,而随后N-n*次激发采用变角激发方法(其中N为总的激发次数),同时选取合适的定变角激发参数,使得重建的MR图像具有更高的SNR和更好的图像质量,算法鲁棒性强,且噪声、伪影等与现有技术相比明显减小。
Description
技术领域
本发明涉及磁共振成像技术领域,尤其是涉及一种基于定-变角优化激发的超极化磁共振成像方法。适用于超极化129Xe或13C等磁共振快速回波成像,也可适用于其它需要快速回波成像领域。
背景技术
空气污染(例如,雾霾)加速和加剧肺部疾病问题的恶性发展。医学影像学诊断是肺部疾病最常用也是最直观的一种诊断方式。然而,常用的肺部医学影像诊断技术都存在一定的缺陷,例如胸部X线透视和计算机X射线断层扫描(CT)所放射出的射线都属于高能射线,会对人体造成一定伤害(高能射线会杀死人体内的白细胞),而正电子发射计算机断层(PET)所需的核素也有一定的放射性,所以不宜短时间内多做。超极化129Xe磁共振成像(MRI)是一种先进的无放射性、非侵入的影像学技术,且129Xe能溶入血液和组织,不仅能对肺部通气和肺部表观扩散系数进行成像,还能评价肺部气-气交换和气-血交换功能,从而给临床肺部医学研究引入了一种新的研究手段。
传统MRI是质子(1H)成像,而肺部空腔的质子核自旋浓度比肌肉和脑部组织约小1000倍,因而肺部是其“盲区”。超极化129Xe MRI利用激光光泵和自旋交换技术增强惰性气体(129Xe)的磁共振信号(极化度增强103~105倍),使得对肺部气体成像成为可能。然而,超极化129Xe的气体磁化矢量在采样过程中具有不可再生性,导致其可用磁化矢量强度随每一次射频(RF)脉冲激发的发生而逐步减小。现有技术力图在有限的屏气时间和超极化气体的磁化矢量强度条件下,发展快速有效的成像方法以期最大化MR信号的信噪比(SNR)。
小角度激发方法(例如,定角激发和变角激发)可以在有限的屏气时间和超极化气体的磁化矢量强度下,获得具有较高信噪比的MR图像(见Mugler JP的“Optimization of gradient-echo sequences for hyperpolarized noble gas MRI”,Proceedings of the 6th International Society of Magnetic Resonance in MedicineMeeting,Sydney,Australia,1998,1904页;Lei Zhao等人的“Gradient-echo imagingconsiderations for hyperpolarized 129Xe MRI”,Journal of Magnetic Resonance,SeriesB 113,179-183,1996;以及Yan Xing等人的“Optimal variable flip angle schemesfor dynamic acquisition of exchanging hyperpolarized substrates”,Journal ofMagnetic Resonance,234,75-81,2013)。然而,小角度激发方法存在一定的不足,例如,定角激发方法导致FID幅值线性降低;变角激发方法尽管能保持FID幅值不变,但相对定角激发方法的初始FID幅值,其FID幅值较低。
近年来,超极化129Xe MRI的应用需求,动态成像,以及溶解态成像等方面都对超极化129Xe MRI快速成像提出更高的要求,因此,如何实现快速、高效、适应性强的超极化129Xe MRI成像方法是超极化129Xe MRI临床肺部医学研究的重大挑战。
发明内容
本发明是针对现有基于小角度激发的快速回波成像方法存在的上述技术问题,提供了一种基于定-变角优化激发的超极化磁共振成像方法。
一种基于定-变角优化激发的超极化磁共振成像方法,包括以下步骤:
步骤1,初始化定-变角优化激发方法的射频激发总次数N;
步骤2,初始化定角激发阶段的激发角度,即翻转角α,翻转角α的取值范围为(αmin,αT],其中αmin的取值为arccos(e-1/N),αT的取值为2αmin;
步骤3,根据下式确定定角激发阶段的定角激发次数n*;
其中,方括号[]表示取最接近方括号内数值的正整数,TR为重复时间,T1为从超极化磁化矢量弛豫到热极化平衡态的时间;
步骤4,根据下式确定变角激发阶段的初始翻转角θ*(1)
α≤θ*(1)≤arcsin(exp(TR/T1)·tanα);
步骤5,根据下式确定变角激发阶段的第j次激发角度与第j+1次激发角度之间的关系应满足,其中j=1,…,N-n*:
sinθ*(j+1)=exp(TR/T1)·tanθ*(j);
步骤6,根据步骤2确定的定角激发翻转角α和步骤3确定的定角激发次数n*进行定角激发,采集超极化磁共振(MR)信号;
步骤7、根据步骤4确定的变角激发的初始翻转角θ*(1)、步骤5确定的第j次激发角度与第j+1次激发角度之间的关系及其变角激发的次数进行变角激发,采集超极化MR信号;
步骤8、根据步骤6和步骤7采集的超极化MR信号进行MR图像重建,获得超极化MR图像。
如上所述的射频激发总次数N为64,96,128,192或256。
本发明与现有技术相比,具有以下优点:
1.由于本方法首次提出了一种基于定-变角优化激发的超极化磁共振成像算法,该方法综合利用了定角激发和变角激发方法的优点,能弥补其相应的不足,因此这种快速成像方法可以获得更高信噪比的MR信号,其重建的MR图像质量更佳;
2.由于本方法主要分为三个阶段,首先利用定角激发方法采集超极化MR信号,其次采用变角激发方法采集超极化MR信号,最后根据所述采集的MR信号进行超极化磁共振图像重建,能有效弥补采用单一小角度激发方法的不足。
3.由于本方法在定角激发阶段,通过最优化MR信号的SNR,获得定角激发的最优化步数,避免FID幅值的过度损失,缓解MR图像失真效应。
4.由于本方法在变角激发阶段,通过定角激发和变角激发所满足的横向磁化矢量强度关系式确定变角激发方法的初始角度,校正MR图像伪影效应。
5.有效地弥补小角度激发方法所存在的不足,提高了超极化磁共振的成像速度和MR图像质量。
附图说明
图1是一种基于定-变角优化激发的超极化磁共振成像方法流程图。
图2是一种基于定-变角优化激发的超极化磁共振成像方法的示意图。
图3(a)是采用定角激发(constant flip angle,CFA),变角激发(variable flipangle,VFA)和定-变角优化激发(constant-variable flip angle,CVFA)方法的剩余超极化纵向磁场强度与激发次数之间的关系图(分别用曲线1,曲线2,曲线3表示),而图3(b)是三种激发方法所产生的信号强度与激发次数之间的关系图(分别用曲线1,曲线2,曲线3表示),其中,CFA的激发角度为12°,VFA的初始激发角度为5.07°,CVFA中定角激发阶段的激发角度为8.72°,变角激发阶段的初始激发角度为8.77°,CFA、VFA和CVFA总的激发次数均为128(CVFA中的实验参数选用实施例2中的求取的参数)。图3(a)表明,CFA导致剩余磁场强度衰减最快,VFA衰减最慢。图3(b)表明,CFA导致MR信号强度随激发次数的增加从0.2079快速衰减到0.0126(信号强度经归一化处理),而VFA导致MR信号强度始终保持为0.0884,且当激发次数为40时,CFA和VFA所产生的信号强度近似相等。由此可见,CVFA是CFA和VFA在剩余磁场强度方面和信号强度方面的一种折中处理方法,可以有效避免FID幅值的过度损失,提高MR信号的信噪比和MR图像质量。
图4是在不同的翻转角(激发角)条件下,定角激发的MR信号的信噪比与激发次数之间的关系图,分别用曲线1(翻转角为9度的情况),曲线2(翻转角为12度的情况),曲线3(翻转角为18度的情况)表示,其中,图中实心点分别代表三种翻转角的情况下SNR最大值所在的位置。通过最优化SNR获取定角激发的最优激发次数,可以有效避免FID幅值的过度损失,缓解MR图像失真效应。
具体实施方式
下面通过实施例1~3,并结合附图,对本发明的技术方案作进一步具体的说明。
实施例1:
如图1所示,一种基于定-变角优化激发的超极化磁共振成像方法,该方法分为三个阶段:定角激发阶段,变角激发阶段,以及MR图像重建阶段,包括以下步骤:
步骤1,初始化定-变角优化激发方法的射频(RF)激发次数N,N一般取值为64,96,128,192或256;
步骤2,初始化定角激发阶段的激发角度(翻转角)α,翻转角α的取值范围为(αmin,αT],其中αmin的取值为arccos(e-1/N),具体推导过程参考步骤3,αT的取值为2αmin。
步骤3,确定定角激发阶段的定角激发次数n*。定角激发n次后,MR信号的SNR与下式成正比例关系:
其中,Ts是数据采样周期,α为翻转角,TR为重复时间,T1为从超极化磁化矢量弛豫到热极化平衡态的时间。当翻转角α为定值时,最优化SNR,即且n的取值为正整数,则可得最优化SNR的定角激发次数n*:
其中,[x]表示最接近x的正整数。相对TR而言,超极化气体从超极化平衡态弛豫到热极化平衡态的弛豫时间很长,因此一般可以忽略T1弛豫。若忽略T1弛豫,则TR/T1→0。
定-变角优化激发方法先实施n*次定角激发采集超极化磁共振信号,随后实施N-n*次变角激发采集MR信号,则n*与N的关系应满足(忽略T1弛豫):
步骤4,确定变角激发阶段的初始翻转角θ*(1)。在定角激发阶段,第i次定角激发后,横向磁化矢量强度为:
S(i)=λm0exp[-(i-1)·TR/T1]sinαcosi-1α (4)
其中,λ为超极化增强因子,m0为热极化平衡态的纵向磁化矢量。由于超极化气体热极化平衡态的纵向磁化矢量对超极化纵向磁化矢量的影响可以通过损毁梯度得以消除,因此i次定角激发后,剩余的超极化纵向磁化矢量正比于cosiα。
定-变角优化激发方法从第n*+1次起实施变角激发策略,目的在于避免FID幅值的过度衰减。假若第n*+1次激发依然采用定角激发方法,则第一次变角激发后的横向磁化矢量强度应不小于第n*+1次定角激发后的横向磁化矢量强度,同时也不大于第n*次定角激发后的横向磁化矢量强度。因此,第一次实施变角激发后的横向磁化矢量强度与定角激发的横向磁化矢量强度之间的关系应满足:
S(n*+1)≤S*(1)≤S(n*) (5)
其中,S(n*)表示第n*次定角激发后的横向磁化矢量强度,S*(1)表示第一次变角激发后的横向磁化矢量强度。依据式(5)可得:
即
则定角α与第一次变角θ*(1)之间的关系式为:
α≤θ*(1)≤arcsin(exp(TR/T1)·tanα) (8)
若忽略T1弛豫,则上式可简化为α≤θ*(1)≤arcsin(tanα)。依据式(8)可以确定变角激发阶段的初始翻转角θ*(1),例如,θ*(1)的取值为(α+arcsin(tanα))/2。
步骤5,确定N-n*-1次变角激发的激发角度θ*(i),i=2,…,N-n*。定角激发导致n*个FID幅值逐步降低,而随后N-n*次变角激发的目的在于保持N-n*个FID幅值不变。因此,变角激发阶段的第j次激发角度与第j+1次激发角度之间的关系应满足(j=1,…,N-n*):
若忽略T1弛豫,则上式可简化为sinθ*(j+1)=tanθ*(j)。因此,依据式(9)可以确定N-n*-1次变角激发的激发角度θ*(i),i=2,…,N-n*。
步骤6,根据步骤2确定的定角激发翻转角α和步骤3确定的定角激发次数n*进行定角激发,采集超极化MR信号;
步骤7、根据步骤4确定的变角激发的初始翻转角θ*(1)、步骤5确定的第j次激发角度与第j+1次激发角度之间的关系及其变角激发的次数N-n*进行变角激发,采集超极化MR信号;
步骤8、根据步骤6和步骤7采集的超极化MR信号,并进行MR图像重建,获得超极化磁共振MR图像。
实施例2:
采用实施例1中的技术方案,当N=128时,依据步骤2可得定角激发阶段的激发角度α的取值范围为[7.15°,14.30°],取α=8.72°;依据步骤3可得定角激发次数n*=86;依据步骤4可得变角激发阶段的初始翻转角θ*(1)的范围为[8.72°,8.82°],取θ*(1)=8.77°;依据步骤5可得变角激发的一系列激发角度8.77°,8.87°,8.98°,…,35.15°,44.76°;最后依据步骤6和步骤7采集超极化MR信号,并进行MR图像重建,获得超极化磁共振MR图像。
实施例3:
超极化气体的磁化矢量具有不可恢复性,因此期望第N次激发把剩余纵向磁场矢量全部翻转到x-y平面。满足上述要求的第N次激发的激发角度应为90度,则需要对实施例1中的相关参数—定角激发阶段的激发角度(翻转角)α,定角激发阶段的激发次数n*,和变角激发阶段的初始翻转角θ*(1)进行优化,得到优化后的定-变角激发策略以获取高SNR的超极化磁共振信号。
步骤1,初始化定-变角优化激发方法的射频(RF)激发次数N,N一般取值为64,96,128,192或256;
步骤2,确定优化的定角激发角度(翻转角)αopt,优化的定角激发次数n* opt,和优化的变角激发初始翻转角θ* opt(1)。
充分利用超极化气体不可恢复的磁化矢量需要满足第N-1次激发的激发角度为45度(忽略T1弛豫),则依照式(9)逆推,可得第N-n*-1次激发的激发角度,即变角激发阶段的初始激发角应满足:
为便于定-变角优化激发方法的实施,定角激发阶段优化的激发角度αopt取值为θ* opt(1),则依据式(2)和(10)可以确定优化的定角激发次数n* opt。N与n* opt之间的关系式应满足(忽略T1弛豫):
对于给定的N,利用式(11)可以求得优化的n* opt,
再利用式可以确定定角激发阶段的优化激发角度αopt。
再利用式(10)可以确定变角激发阶段优化的初始激发角θ* opt(1),
再利用式(9)可以确定变角激发阶段的第j次激发角度与第j+1次激发角度之间的关系。
例如,当N=128时,n* opt=85,αopt=8.77°,θ* opt(1)=8.77°;当N=256时,n* opt=170,αopt=6.21°,θ* opt(1)=6.21°。
步骤3,利用步骤2所确定的αopt和n* opt,采用定角激发方法采集超极化MR信号。
步骤4,利用步骤2所确定的θ* opt(1)、式(9)确定的第j次激发角度与第j+1次激发角度之间的关系和变角激发的次数N-n*,采用变角激发方法采集超极化MR信号。
步骤5,利用步骤3和步骤4所采集的超极化MR信号进行MR图像重建,获得超极化129Xe磁共振图像。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (2)
1.一种基于定-变角优化激发的超极化磁共振成像方法,其特征在于,包括以下步骤:
步骤1,初始化定-变角优化激发方法的射频激发总次数N;
步骤2,初始化定角激发阶段的激发角度,即翻转角α,翻转角α的取值范围为(αmin,αT],其中αmin的取值为arccos(e-1/N),αT的取值为2αmin;
步骤3,根据下式确定定角激发阶段的定角激发次数n*;
其中,方括号[]表示取最接近方括号内数值的正整数,TR为重复时间,T1为从超极化磁化矢量弛豫到热极化平衡态的时间;
步骤4,根据下式确定变角激发阶段的初始翻转角θ*(1)
α≤θ*(1)≤arcsin(exp(TR/T1)·tanα);
步骤5,根据下式确定变角激发阶段的第j次激发角度与第j+1次激发角度之间的关系应满足,其中j=1,…,N-n*:
sinθ*(j+1)=exp(TR/T1)·tanθ*(j);
步骤6,根据步骤2确定的定角激发翻转角α和步骤3确定的定角激发次数n*进行定角激发,采集超极化磁共振(MR)信号;
步骤7、根据步骤4确定的变角激发的初始翻转角θ*(1)、步骤5确定的第j次激发角度与第j+1次激发角度之间的关系及其变角激发的次数进行变角激发,采集超极化MR信号;
步骤8、根据步骤6和步骤7采集的超极化MR信号进行MR图像重建,获得超极化MR图像。
2.根据权利要求1所述的一种基于定-变角优化激发的超极化磁共振成像方法,其特征在于,所述的射频激发总次数N为64,96,128,192或256。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410438269.8A CN104181483B (zh) | 2014-08-29 | 2014-08-29 | 一种基于定-变角优化激发的超极化磁共振成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410438269.8A CN104181483B (zh) | 2014-08-29 | 2014-08-29 | 一种基于定-变角优化激发的超极化磁共振成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104181483A true CN104181483A (zh) | 2014-12-03 |
CN104181483B CN104181483B (zh) | 2016-09-07 |
Family
ID=51962685
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410438269.8A Active CN104181483B (zh) | 2014-08-29 | 2014-08-29 | 一种基于定-变角优化激发的超极化磁共振成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104181483B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105301543A (zh) * | 2015-11-16 | 2016-02-03 | 中国科学院武汉物理与数学研究所 | 一种基于单次屏气内多角度激发的超极化角度校准方法 |
CN106355624A (zh) * | 2016-08-26 | 2017-01-25 | 中国科学院武汉物理与数学研究所 | 一种基于变采样率的肺部超极化气体动态成像方法 |
CN107219484A (zh) * | 2016-03-22 | 2017-09-29 | 西门子(深圳)磁共振有限公司 | 多层并行激发弥散成像的激发翻转角度确定方法和装置 |
CN107274461A (zh) * | 2017-06-14 | 2017-10-20 | 合肥工业大学 | 一种超极化气体肺部磁共振图像重建方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5307015A (en) * | 1990-06-29 | 1994-04-26 | The Regents Of The University Of California | NMR relaxometry using variable initial flip angle |
JP6084430B2 (ja) * | 2012-10-30 | 2017-02-22 | 株式会社日立製作所 | 再収束rfパルスのフリップ角制御法及び磁気共鳴イメージング装置 |
-
2014
- 2014-08-29 CN CN201410438269.8A patent/CN104181483B/zh active Active
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105301543A (zh) * | 2015-11-16 | 2016-02-03 | 中国科学院武汉物理与数学研究所 | 一种基于单次屏气内多角度激发的超极化角度校准方法 |
CN105301543B (zh) * | 2015-11-16 | 2017-12-19 | 中国科学院武汉物理与数学研究所 | 一种基于单次屏气内多角度激发的超极化角度校准方法 |
CN107219484A (zh) * | 2016-03-22 | 2017-09-29 | 西门子(深圳)磁共振有限公司 | 多层并行激发弥散成像的激发翻转角度确定方法和装置 |
CN107219484B (zh) * | 2016-03-22 | 2023-09-22 | 西门子(深圳)磁共振有限公司 | 多层并行激发弥散成像的激发翻转角度确定方法和装置 |
CN106355624A (zh) * | 2016-08-26 | 2017-01-25 | 中国科学院武汉物理与数学研究所 | 一种基于变采样率的肺部超极化气体动态成像方法 |
CN106355624B (zh) * | 2016-08-26 | 2018-12-25 | 中国科学院武汉物理与数学研究所 | 一种基于变采样率的肺部超极化气体动态成像方法 |
CN107274461A (zh) * | 2017-06-14 | 2017-10-20 | 合肥工业大学 | 一种超极化气体肺部磁共振图像重建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104181483B (zh) | 2016-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fitts et al. | Arrhythmia insensitive rapid cardiac T1 mapping pulse sequence | |
Möller et al. | MRI of the lungs using hyperpolarized noble gases | |
US7627359B2 (en) | Quantitative contrast enhanced black-blood imaging using quadruple-inversion recovery | |
Zhao et al. | Three-dimensional turbo-spin-echo amide proton transfer MR imaging at 3-Tesla and its application to high-grade human brain tumors | |
Subramanian et al. | Single‐point (constant‐time) imaging in radiofrequency Fourier transform electron paramagnetic resonance | |
Park et al. | Detection of small metastatic brain tumors: comparison of 3D contrast-enhanced whole-brain black-blood imaging and MP-RAGE imaging | |
JP5818631B2 (ja) | 磁気共鳴イメージング装置 | |
CN104181483A (zh) | 一种基于定-变角优化激发的超极化磁共振成像方法 | |
US10151818B2 (en) | Magnetic resonance angiography method and apparatus | |
Zhang et al. | MRI Gibbs‐ringing artifact reduction by means of machine learning using convolutional neural networks | |
CN104507386A (zh) | 使用旋转阵列的永磁体用于便携式磁共振成像的系统和方法 | |
Gallichan | Diffusion MRI of the human brain at ultra-high field (UHF): A review | |
Chen et al. | Comparison of reduced field-of-view diffusion-weighted imaging (DWI) and conventional DWI techniques in the assessment of Cervical carcinoma at 3.0 T: Image quality and FIGO staging | |
Diaz-de-Grenu et al. | MRI detection of tissue pathology beyond atrophy in Alzheimer's disease: Introducing T2-VBM | |
WO2016050938A1 (en) | Magnetic resonance imaging with enhanced bone visualization | |
US20080272781A1 (en) | Q-space sampling method and diffusion spectrum imaging method employing the same | |
JP2014502910A (ja) | インターリーブスピンロッキングイメージング | |
Zhang et al. | Proton/deuterium magnetic resonance imaging of rodents at 9.4 T using birdcage coils | |
CN106646300A (zh) | 一种超极化气体磁共振影像处理方法 | |
Wehrli | The origins and future of nuclear magnetic resonance imaging | |
US20020169372A1 (en) | Magnetic resonance imaging using MT pulse of which duration is shorter | |
US20180059201A1 (en) | Magnetic resonance imaging apparatus | |
CN103354908B (zh) | 用于利用极化转移在高静(b0)磁场下的宽带宽磁共振波谱分析的方法、设备和装置 | |
KR101056451B1 (ko) | 유동성 혈류의 신호를 선택적으로 제거하면서 조직체의 t1-강조 영상을 획득하는 방법 | |
JP5289011B2 (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20180529 Address after: 430000 318 Friendship Avenue, Wuchang District, Wuhan, Hubei. Patentee after: Wuhan Zhongke polarization Medical Technology Co., Ltd. Address before: 430071 Xiaohong Shanxi No. 30, Wuchang District, Wuhan, Hubei Patentee before: Wuhan Inst. of Physics and Mathematics, Chinese Academy of Sciences |