CN111685764B - 一种基于磁共振指纹成像的快速定量T1ρ的方法 - Google Patents
一种基于磁共振指纹成像的快速定量T1ρ的方法 Download PDFInfo
- Publication number
- CN111685764B CN111685764B CN202010400589.XA CN202010400589A CN111685764B CN 111685764 B CN111685764 B CN 111685764B CN 202010400589 A CN202010400589 A CN 202010400589A CN 111685764 B CN111685764 B CN 111685764B
- Authority
- CN
- China
- Prior art keywords
- t1rho
- sequence
- magnetic resonance
- spin
- tsl
- 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.)
- Active
Links
Images
Classifications
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Psychiatry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Physiology (AREA)
- High Energy & Nuclear Physics (AREA)
- Signal Processing (AREA)
- Radiology & Medical Imaging (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明属于磁共振成像技术领域，具体为基于磁共振指纹成像的快速定量T1ρ的方法。本发明方法包括在传统MRF序列前增加spin‑lock准备脉冲，通过合理建模，将信号的T1ρ演化与信号的T1、T2演化融合一起，使每个时间点采集的信号都包含T1、T2和T1ρ的成分。设计伪随机变化的序列参数TSL、FA、TR和TE，通过Bloch方程中磁化矢量的演化，建立不同弛豫时间值组合的字典，将实际采集到的图像序列与字典匹配，从而定量出T1ρ、T1和T2值。本发明定量T1ρ的方法与传统指数拟合的方法相比，可以大幅的缩短扫描时间，预期在T1ρ成像中具有较大的应用价值。
Description
技术领域
本发明属于磁共振成像技术领域，具体涉及基于磁共振指纹成像的快速定量T1ρ的方法。
背景技术
磁共振成像 (magnetic resonance imaging, MRI)是一种无创的，具有高空间分辨率、高组织对比度、可任意方向成像的医学影像技术。磁共振成像具有丰富的组织对比度，是临床上重要的诊断工具，被广泛应用于临床诊断、神经科学和精神类疾病等领域。
磁共振指纹成像（magnetic resonance fingerprinting,MRF）最初由Ma等提出(Ma D, Gulani V, Seiberlich N, et al. Magnetic resonance fingerprinting.Nature. 2013;495(7440):187-192)，MRF采用伪随机变化的序列参数进行数据采集,通过布洛赫（Bloch）方程仿真信号演化建立字典，该字典包含大量不同参数（如T1，T2）组合，将采集到的实际磁共振信号用模式识别的技术与字典进行匹配，寻找出最佳的匹配曲线,从而定量出组织的T1、T2值（T1指纵向弛豫时间，T2指横向弛豫时间）。除了可以快速准确定量T1，T2外，MRF还可以应用于水脂分离和波谱中肌酸激酶反应率的定量，以及扩散系数ADC的快速定量，MRF已成为磁共振快速定量技术研发的热门领域。
T1ρ（或T1rho）是指旋转坐标系中的自旋-晶格弛豫时间(spin-latticerelaxation time in the rotating frame), 反映的是大分子的低频运动信息。T1ρ成像中，自旋锁定频率（spin-lock frequency,FSL）一般在几百赫兹到几千赫兹，因此T1ρ对大分子的低频运动很敏感，T1ρ可以反映大分子缓慢运动过程中大分子与氢质子间相互作用的信息。T1ρ的应用主要有T1ρ加权成像，T1ρ mapping和T1ρ分布（T1ρ dispersion）。T1ρ加权成像（或T1ρ mapping）具有不同于传统T1、T2加权成像的软组织对比，被广泛应用于肝脏纤维化分期、关节软骨退变及阿尔兹海默症等疾病的研究。
目前T1ρ的计算一般是通过扫描不同自旋锁定脉冲的时间（spin-lock time,TSL），通过函数：
拟合的方法得到，式中S0是TSL为0时的信号强度，C为常数。指数拟合的方法比较耗时，容易受到运动伪影等的影响，不利于T1ρ成像方法在临床等领域的应用。
发明内容
本发明的目的是为了解决T1ρ成像时间长的问题，提供一种基于磁共振指纹成像的快速定量T1ρ的方法。
为了达到上述发明目的，本发明采用以下技术方案：
通过合理建模，将信号的T1ρ演化与信号的T1、T2演化融合一起，使每个时间点采集的信号都包含T1、T2和T1ρ的成分。设计伪随机变化的序列参数（自旋锁定时间TSL、翻转角FA、重复时间TR和回波时间TE），通过Bloch方程中磁化矢量的演化，建立不同弛豫时间值组合的字典，将实际采集到的图像序列与字典匹配，从而定量出T1ρ、T1和T2值。具体步骤为：
（1）spin-lock准备部分：在传统磁共振指纹成像序列前增加spin-lock准备脉冲；
（2）spin-lock准备部分结束后使用梯度散相残余磁化矢量；
（3）数据采集：使用基于平面回波（EPI）采集方式采集图像数据；也可以是施加180度重聚脉冲后EPI采集方式采集图像数据；也可以是基于磁共振指纹成像技术的采集方式采集图像数据，包括翻转恢复FISP序列、快速自旋回波序列（FSE或RARE）、平衡式稳态自由进动序列（bSSFP）以及螺旋轨迹（Spiral）等采集方式采集图像数据；
（4）伪随机设计序列参数，每次重复是TSL、FA、TR和TE的不同组合，信号按照如下公式演化：
最后将采集的图像序列与字典匹配识别，同时定量出T1ρ、T1和T2。式中S0是TSL为0时的信号强度。
步骤（1）中，spin-lock准备部分，可以是沿着单方向的硬脉冲，为了减少射频场（B1）的不均匀性带来的误差和伪影，可以是由两个方向不同的硬脉冲组合，为了减少静磁场（B0）不均匀性的影响，也可以在两个或者多个自旋锁定脉冲之间增加非选择性180度重聚脉冲。
步骤（2）中，所述使用梯度散相残余磁化矢量，即当spin-lock准备部分结束后在X,Y和Z轴（也可以是某一个或两个轴）施加梯度，以散相XY平面内的残余磁化矢量对后续成像的影响。
步骤（3）中，数据采集使用的是单次激发平面回波序列（EPI），为了减少场不均匀的影响，施加180度重聚脉冲后EPI采集，基于磁共振指纹成像技术的采集方式还可以包括翻转恢复FISP序列、快速自旋回波序列（FSE或RARE）、平衡式稳态自由进动序列（bSSFP）以及螺旋轨迹采集方式（Spiral）等，其中以Spiral采集速度快等优点最为常用。
步骤（4）中，所述伪随机设计序列参数，除传统磁共振指纹成像序列每次重复伪随机改变序列参数（FA、TR和TE）外，融合T1ρ信号演化，伪随机设计TSL参数，使得每次采集的信号都是不同TSL、FA、TR和TE的组合。
考虑到传统磁共振指纹成像方法可以只改变FA、TR和TE中的部分，因此本发明也可以适当将FA、TR和TE中的部分参数固定，伪随机设计参数的个数可以视信号的信噪比和与字典匹配效果来定。
除传统磁共振指纹成像序列每次重复伪随机改变序列参数（FA、TR和TE）外，融合T1ρ信号演化，伪随机设计TSL参数，使得每次采集的信号都是不同TSL、FA、TR和TE的组合。
步骤（4）中，所述图像序列与字典匹配，使用最小二乘法进行匹配；也可以通过比较实际采集图像序列与字典之间的相关性进行匹配；还可以使用点积法进行匹配。
本发明定量T1ρ的方法与传统指数拟合的方法相比，可以大幅的缩短扫描时间，而且本发明快速的定量T1ρ值的同时，能够定量T1和T2值。因此，本发明在T1ρ成像中具有较大的应用价值。
附图说明
图1是本发明的流程图。
图2是本发明的序列时序图。
图3是使用传统饱和恢复法定量出的T1 map， CPMG序列定量出的T2 map，和于RARE采集的T1ρ成像通过指数拟合定量出的T1ρ map。
图4是使用本发明定量出的T1、T2和T1ρ map。
具体实施方式
使用本发明方法对7个T1、T2和T1ρ不同的水模快速定量为例，来说明本发明的具体实施方式。
序列设计：本发明的序列实现是在布鲁克11.7T磁共振系统上实现，软件版本为ParaVision 6.01,接收线圈采用4通道收发一体圆筒线圈。序列每一次激发采用不同的TSL、FA、TR和TE，其中TSL从0ms到50ms,FA从5度到80度，TE从11ms到21ms,TR以Perlin噪声模式变化，范围从25ms到200ms。
传统指纹成像方法一般激发次数设置为1000次时可以获得相对准确的值，本发明也采用1000次的激发次数，对7个水模进行扫描。
使用基于RARE采集的传统T1ρ成像序列对7个水模进行T1ρ map的测量，序列使用11个不同的TSL值进行拟合。
使用饱和恢复法对7个水模进行T1 map的测量，序列使用12个不同的TR值进行拟合。
使用CPMG序列对7个水模进行T2 map的测量，采集20个不同的回波时间信号进行拟合。
根据伪随机设计的序列参数TSL、FA、TR和TE，以及布洛赫方程进行字典建立，本发明中字典建立和字典匹配均使用MathWork公司的MATLAB（2018a）软件上实现。
将上述使用本发明的方法采集到的1000幅不同时间点的图像序列与建立的字典进行匹配，匹配的方法使用最小二乘法，定量出7个水模的T1、T2和T1ρ map。
使用SPSS软件（版本16.0）计算组内相关系数（intraclass correlationcoefficient,ICC）（选择Two-way mixed模型）评价本技术的结果与传统拟合方法的结果的一致性。两种方法测量得到的ICC分别为：T1值ICC=0.842 （p=0.004），T2值ICC=0.834 （p=0.005），T1ρ值ICC=0.919 （p=0.001），表明本发明得到的结果与传统拟合方法得到的结果具有很好的一致性。本发明用很短的数据采集时间就可以同时定量出T1、T2和T1ρ值，大大减少了扫描时间，这对T1ρ成像的广泛应用具有很好的促进作用。
Claims (5)
1.一种基于磁共振指纹成像的快速定量T1ρ的方法，其特征在于，将信号的T1ρ演化与信号的T1、T2演化融合一起，使每个时间点采集的信号都包含T1、T2和T1ρ的成分；设计伪随机变化的序列参数TSL、FA、TR和TE，通过Bloch方程中磁化矢量的演化，建立不同弛豫时间值组合的字典，将实际采集到的图像序列与字典匹配，从而定量出T1ρ值，同时得到T1和T2值；其中，T1指纵向弛豫时间，T2指横向弛豫时间，T1ρ指自旋-晶格弛豫时间；TSL为自旋锁定时间，FA为翻转角，TR为重复时间，TE为回波时间；具体步骤为：
（1）spin-lock准备部分：在传统磁共振指纹成像序列前增加spin-lock准备脉冲；
（2）spin-lock准备部分结束后使用梯度散相残余磁化矢量；
（3）数据采集：使用基于平面回波EPI采集方式采集图像数据；或者施加180度重聚脉冲后EPI采集方式采集图像数据；或者使用基于磁共振指纹成像技术的采集方式采集图像数据，如采用翻转恢复FISP序列、快速自旋回波序列、平衡式稳态自由进动序列及螺旋轨迹采集方式采集图像数据；
（4）伪随机设计序列参数，每次重复是TSL、FA、TR和TE的不同组合，信号按照如下公式演化：
最后将采集的图像序列与字典匹配识别，同时定量出T1ρ、T1和T2，式中S0是TSL为0时的信号强度。
2.根据权利要求1所述的基于磁共振指纹成像的快速定量T1ρ的方法，其特征在于，步骤（1）中所述spin-lock准备部分，是沿着单方向的硬脉冲，或者是由两个方向不同的硬脉冲组合，或者是在多个自旋锁定脉冲之间增加非选择性180度重聚脉冲。
3.根据权利要求2所述的基于磁共振指纹成像的快速定量T1ρ的方法，其特征在于，步骤（2）中所述使用梯度散相残余磁化矢量，是当spin-lock准备部分结束后在X,Y和Z轴中至少一个轴施加梯度，以散相XY平面内的残余磁化矢量对后续成像的影响。
4.根据权利要求2所述的基于磁共振指纹成像的快速定量T1ρ的方法，其特征在于，步骤（4）中所述伪随机设计序列参数，除传统磁共振指纹成像序列每次重复伪随机改变序列参数FA、TR和TE外，融合T1ρ信号演化，伪随机设计TSL参数，使得每次采集的信号都是不同TSL、FA、TR和TE的组合。
5.根据权利要求2所述的基于磁共振指纹成像的快速定量T1ρ的方法，其特征在于，步骤（4）中所述图像序列与字典匹配，使用最小二乘法进行匹配，或者通过比较实际采集图像序列与字典之间的相关性进行匹配，或者使用点积法进行匹配。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010400589.XA CN111685764B (zh) | 2020-05-13 | 2020-05-13 | 一种基于磁共振指纹成像的快速定量T1ρ的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010400589.XA CN111685764B (zh) | 2020-05-13 | 2020-05-13 | 一种基于磁共振指纹成像的快速定量T1ρ的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111685764A CN111685764A (zh) | 2020-09-22 |
| CN111685764B true CN111685764B (zh) | 2021-08-20 |
Family
ID=72477674
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010400589.XA Active CN111685764B (zh) | 2020-05-13 | 2020-05-13 | 一种基于磁共振指纹成像的快速定量T1ρ的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111685764B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112132825A (zh) * | 2020-10-10 | 2020-12-25 | 青岛市市立医院 | 一种基于T2 mapping及T1ρ关节软骨定量成像分析方法 |
| CN112617796B (zh) * | 2020-12-30 | 2024-10-15 | 深圳市联影高端医疗装备创新研究院 | 磁共振功能成像方法、装置、设备、存储介质及系统 |
| CN114217255B (zh) * | 2021-11-29 | 2022-09-20 | 浙江大学 | 一种快速的肝脏多参数定量成像方法 |
| CN114167333B (zh) * | 2021-12-06 | 2022-07-22 | 无锡鸣石峻致医疗科技有限公司 | 一种纵向弛豫时间测定方法、装置、计算机设备及存储介质 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6836114B2 (en) * | 2002-03-15 | 2004-12-28 | The Trustees Of The University Of Pennsylvania | Pulse imaging sequences and methods for T1p-weighted MRI |
| WO2008048641A2 (en) * | 2006-10-17 | 2008-04-24 | The Trustees Of The University Of Pennsylvania | Reducing imaging-scan times for mri systems |
| US20120226136A1 (en) * | 2011-03-04 | 2012-09-06 | The Chinese University Of Hong Kong | Methods and devices of detection, grading, monitoring, and follow-up of fibrosis |
| US10712297B2 (en) * | 2015-04-17 | 2020-07-14 | New York University | Systems and methods for super-resolution surface-layer microscopy using magnetic resonance |
| JP6687383B2 (ja) * | 2015-12-21 | 2020-04-22 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置 |
| CN111090069B (zh) * | 2019-11-21 | 2022-03-29 | 深圳先进技术研究院 | 定量磁共振成像参数确定方法、装置、设备及存储介质 |
-
2020
- 2020-05-13 CN CN202010400589.XA patent/CN111685764B/zh active Active
Non-Patent Citations (3)
| Title |
|---|
| 《Multiparametric or practical quantitative liver MRI: towards millisecond, fat fraction, kilopascal and function era》;Emre Unal 等;《Expert Review of Gastroenterology & Hepatology》;20161209;1-16页 * |
| 《T1rho Fractional-order Relaxation of Human Articular Cartilage*》;Lixian Zou等;《2019 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC)》;20191231;4496-4499页 * |
| RSNA2015腹部影像学;冯翠等;《放射学实践》;20160220(第02期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111685764A (zh) | 2020-09-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111685764B (zh) | 一种基于磁共振指纹成像的快速定量T1ρ的方法 | |
| Bojorquez et al. | What are normal relaxation times of tissues at 3 T? | |
| US9513358B2 (en) | Method and apparatus for magnetic resonance imaging | |
| US8508225B2 (en) | T2-weighted and diffusion-weighted imaging using fast acquisition with double echo (FADE) | |
| US9625547B2 (en) | Magnetic resonance imaging method for the quantification of the T1 and/or T2 relaxation times in a sample | |
| Kobayashi et al. | Diffusion-weighting caused by spoiler gradients in the fast imaging with steady-state precession sequence may lead to inaccurate T2 measurements in MR fingerprinting | |
| US8797031B2 (en) | MR imaging system for discriminating between imaged tissue types | |
| CN108663644B (zh) | 基于残差网络重建的单扫描定量磁共振t2*成像方法 | |
| Foxley et al. | Improving diffusion-weighted imaging of post-mortem human brains: SSFP at 7 T | |
| CN111505553B (zh) | 磁共振成像系统和方法 | |
| US8558546B2 (en) | Relaxometry | |
| Heule et al. | Variable flip angle T1 mapping in the human brain with reduced t2 sensitivity using fast radiofrequency‐spoiled gradient echo imaging | |
| Landman et al. | Resolution of crossing fibers with constrained compressed sensing using traditional diffusion tensor MRI | |
| WO2021196865A1 (zh) | 一种振荡梯度准备的3d梯度自旋回波成像方法及设备 | |
| US20110025327A1 (en) | Method for radiofrequency mapping in magnetic resonance imaging | |
| US7495440B2 (en) | Q-space sampling method and diffusion spectrum imaging method employing the same | |
| US7319328B1 (en) | System, method, and computer-readable medium for magnetic resonance diffusion anisotropy image processing | |
| Leroi et al. | Simultaneous proton density, T1, T2, and flip‐angle mapping of the brain at 7 T using multiparametric 3D SSFP imaging and parallel‐transmission universal pulses | |
| Cao et al. | DTI‐MR fingerprinting for rapid high‐resolution whole‐brain T1, T2, proton density, ADC, and fractional anisotropy mapping | |
| Assländer et al. | On multi-path longitudinal spin relaxation in brain tissue | |
| US8717023B2 (en) | Relaxometry quantification self-justification fitting | |
| Dong et al. | SNR-efficient distortion-free diffusion relaxometry imaging using ACcelerated Echo-train shifted EPTI (ACE-EPTI) | |
| Makhijani et al. | Improved blood suppression in three‐dimensional (3D) fast spin‐echo (FSE) vessel wall imaging using a combination of double inversion‐recovery (DIR) and diffusion sensitizing gradient (DSG) preparations | |
| Iyyakkunnel et al. | Complex B1+ mapping with Carr‐Purcell spin echoes and its application to electrical properties tomography | |
| Dong et al. | Single-shot Echo Planar Time-resolved Imaging for multi-echo functional MRI and distortion-free diffusion imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |