CN109814056B - 一种获取磁共振精准定量图像的方法 - Google Patents
一种获取磁共振精准定量图像的方法 Download PDFInfo
- Publication number
- CN109814056B CN109814056B CN201910105084.8A CN201910105084A CN109814056B CN 109814056 B CN109814056 B CN 109814056B CN 201910105084 A CN201910105084 A CN 201910105084A CN 109814056 B CN109814056 B CN 109814056B
- Authority
- CN
- China
- Prior art keywords
- maps
- magnetic resonance
- acquiring
- images
- data
- 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
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明提供一种磁共振精准T1maps的重建方法,其主要目的在于解决传统射频发射场校正方法存在的一些缺陷,并且本发明可以利用采集的一组数据直接获得精准T1maps而不需要计算射频场的B1+maps。该方法利用多个翻转角的梯度回波序列采集数据,在考虑射频发射场不均性的情况下,联合采集到的数据结合梯度回波信号表达式计算T 1 maps。该方法不像传统方法那样,需要先获得射频发射场的B1+maps然后再带入到信号表达式中计算T1maps。传统的DAM校正方法是使用长TR的条件,要求TR>>T1来达到图像信号正比于sin(α)和sin(2α)的目的,而本方法则是直接利用多个图像信号表达式直接计算组织的T1值,因而,本方法没有对长TR条件的依赖,并且相较于传统的一些校正方法很容易实现。
Description
技术领域:本发明属于磁共振成像技术领域,是一种用于获取磁共振精准定量图像的重建方法。
背景技术:定量磁共振成像是近些年发展起来一种新的磁共振成像技术,它是利用采集的多幅磁共振加权像结合磁共振信号表达式来计算某种参数的定量图谱的方法。由于实际磁共振成像设备中信号的影响因素有很多,使用理论表达式计算的定量图像往往存在一定误差,其中在高场磁共振设备中射频发射场(B1+场)的不均匀性就是最重要的误差源之一。为了获得准确的定量图像,通常需要对射频发射场进行校正。目前,已有的射频场不均匀性校正的方法有很多,如传统的双角法(DAM)、双角法的改进方法(如饱和双角法SDAM)以及利用组织参考区域来获得全局拟合B1+maps的方法等。这些方法有着诸多的限制,如DAM需要长TR的条件,这限制了磁共振扫描的速度;而利用参考组织的定量值获得全局拟合B1+maps的方法则需要事先知道参考组织的定量值或者有相应的金标准方法测量参考组织的实际定量值,然而目前国际上并没有一个统一的标准的各种组织各种参数的定量值体系,也没有适合任意组织不同参数的金标准测量方法;此外,上述的这些方法一般只能对单一参数的定量图像进行校正。本发明提出了一种,利用多个翻转角的梯度回波序列结合信号表达式计算精准T1maps(T1定量图像)的方法。该方法克服了传统DAM方法在B1+场校正时对长TR要求的缺陷,并且不需要进行B1+maps的计算便可直接通过解方程形式获取精准的T1maps。
发明内容:
发明目的:本发明提供一种磁共振精准T1maps的重建方法,其主要目的在于解决传统射频发射场校正方法存在的一些缺陷,并且本发明可以利用采集的一组数据直接获得精准T1maps而不需要计算射频场的B1+maps。
技术方案:
一种获取磁共振精准定量图像的方法,其特征在于:
包括:
数据采集:
分别利用三个翻转角的梯度回波序列采集数据,三个翻转角满足倍数关系即α,2α和3α;
重建:
建立方程:
式中:
I1、I2和I3为数据采集步骤中的各个不同的梯度回波序列采集到的图像信号表达式;
t=cosB1αB1为翻转角与名义角之间的比例系数;
对上式求解得到t的图像信号表达式,进行平滑滤波;
选用I1和I2作除法运算,得T1的计算表达式,式如下:
所述的获取磁共振精准定量图像的方法,其特征在于:在选用I1和I2作除法运算的基础上,选用I1和I3或者I2和I3进行相同的运算;将各组T1maps求平均值以进一步降低噪声的影响。
所述的获取磁共振精准定量图像的方法:选用双回波梯度回波序列采集数据,得到R2star-maps。
所述的获取磁共振精准定量图像的方法:选用双回波梯度回波序列采集数据,结合所述R2star-maps得到PDmaps。
所述的获取磁共振精准定量图像的方法:利用计算T1maps的前一步的t,t=cosB1α经过余弦函数反变换得到B1maps。
所述的获取磁共振精准定量图像的方法:选用多回波梯度回波序列采集数据,使用任意回波的数据再结合求均值的降噪方法,获得多幅的T1maps图像,对得到的多幅T1maps求均值可以得到更加准确且噪声更低的定量图像。
优点及效果:
为了获得精准T1maps,本发明提供了一种用于获取磁共振精准T1maps的重建方法。该方法利用多个翻转角的梯度回波序列采集数据,在考虑射频发射场不均性的情况下,联合采集到的数据结合梯度回波信号表达式计算T1maps。该方法不像传统方法那样,需要先获得射频发射场的B1+maps然后再带入到信号表达式中计算T1maps。传统的DAM校正方法是使用长TR的条件,要求TR>>T1来达到图像信号正比于sin(α)和sin(2α)的目的,而本方法则是直接利用多个图像信号表达式直接计算组织的T1值,因而,本方法没有对长TR条件的依赖,并且相较于传统的一些校正方法很容易实现,对于满足VFA(variable flip angle)方法进行定量图像重建需求的TR,本方法均适用。此外,本方法是依据表达式利用“盛金公式”计算的t的解析解,因此,本方法得到的定量值的精确度很高。
附图说明:
图1校正前获得的T1maps。
图2为本方法获得的T1maps。
图3为选取的对比区域示意图。
具体实施方式:
本发明是通过以下技术方案来实现的:
数据采集:分别使用三个不同翻转单回波的梯度回波序列扫描数据,规定三个翻转满足倍数关系,即α,2α和3α。
重建:在考虑射频发射不均匀情况下,引入实际翻转角与名义角之间的比例系数B1,联合多组方程求解T1,算法公式如下:
(1)利用不同角相同回波的信号表达式,假设三幅图像分别为I1、I2和I3,联立三幅图像的信号表达式可得到t的三次方程,这里t=cosB1α,表达式如下:
(2)对上式进行计算,这里可以使用“盛金公式”求解,由于图像存在噪声,且经过多步数学运算,噪声会被放大,因此需要对得到的t的图像进行平滑滤波;
(3)利用三个翻转角中任意两个角的信号表达式作除法运算,这里选用I1和I2作除法运算,可得到T1的计算表达式,表达式如下:
如图所示,图1为利用传统VFA方法在不考虑射频场不均匀性时计算得到的T1maps,图2是使用本方法获得的T1maps,图3为选取的脑部对比区域。
本方法与文献中T1maps校正结果的对比见下表:
上表是本方法获得的T1maps与文献T1maps校正结果的对比,从表中的对比结果可以看出,使用本发明获得的T1maps与一些文献中的通过不同校正方法获得的T1maps的结果基本一致。相比较于未校正的T1maps,本方法结合VFA方法不但能够获得精准的T1maps,而且不像传统校正方法那样需要先通过一些方法获得B1+maps然后再对T1maps校正。因此本方法不但简便而且具有实际的应用价值和意义。
在考虑射频发射场不均匀情况下,本发明通过分析梯度回波序列采集信号时信号的表达式,提出了一种求解精准T1maps的方法,属于磁共振成像技术领域。在引入实际翻转角与名义翻转角之间的比例关系后,通过解方程的形式计算T1值,因此本方法不但可以解决传统方法进行定量成像时需要对射频场不均匀性校正的问题,而且可以获得精准的T1maps。
上表中引述的文献:
[1]Deichmann R.Fast high-resolution T1mapping of the human brain[J].Magn Reson Med,2005,54(1):20-27.
[2]Rooney WD,Johnson G,Li X,et al.Magnetic field andtissuedependencies of human brain longitudinal 1H2O relaxation in vivo[J].Magnetic Resonance in Medicine,2010,57(2):308–318.
[3]Orospeusquens AM,Laurilar M,Shah NJ.Magnetic field dependence ofthe distribution of NMR relaxation times in the living human brain[J].Magnetic Resonance Materials in Physics Biology and Medicine.2008,21:131–147.
Claims (6)
2.根据权利要求1所述的获取磁共振精准定量图像的方法,其特征在于:在选用I1和I2作除法运算的基础上,选用I1和I3或者I2和I3进行相同的运算;将各组T1maps求平均值以进一步降低噪声的影响。
3.根据权利要求1所述的获取磁共振精准定量图像的方法:选用双回波梯度回波序列采集数据,得到R2star-maps。
4.根据权利要求3所述的获取磁共振精准定量图像的方法:选用双回波梯度回波序列采集数据,结合所述R2star-maps得到PDmaps。
5.根据权利要求1所述的获取磁共振精准定量图像的方法:利用计算T1maps的前一步的t,t=cosB1α经过余弦函数反变换得到B1maps。
6.根据权利要求2所述的获取磁共振精准定量图像的方法:选用多回波梯度回波序列采集数据,使用任意回波的数据再结合求均值的降噪方法,获得多幅的T1maps图像,对得到的多幅T1maps求均值可以得到更加准确且噪声更低的定量图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910105084.8A CN109814056B (zh) | 2019-02-01 | 2019-02-01 | 一种获取磁共振精准定量图像的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910105084.8A CN109814056B (zh) | 2019-02-01 | 2019-02-01 | 一种获取磁共振精准定量图像的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109814056A CN109814056A (zh) | 2019-05-28 |
CN109814056B true CN109814056B (zh) | 2020-11-10 |
Family
ID=66606310
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910105084.8A Active CN109814056B (zh) | 2019-02-01 | 2019-02-01 | 一种获取磁共振精准定量图像的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109814056B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113017596B (zh) * | 2021-03-09 | 2022-11-11 | 深圳高性能医疗器械国家研究院有限公司 | 一种磁共振多参数定量方法及其应用 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102759723A (zh) * | 2012-04-11 | 2012-10-31 | 华东师范大学 | 一种磁共振t2*图像生成方法 |
JP2015188635A (ja) * | 2014-03-28 | 2015-11-02 | 株式会社日立製作所 | 磁気共鳴イメージング装置 |
CN105699923A (zh) * | 2015-09-12 | 2016-06-22 | 北京大学 | 一种无创动态测量组织r2,r2*和r2`参数图像的磁共振成像方法 |
CN108294753A (zh) * | 2018-01-03 | 2018-07-20 | 上海东软医疗科技有限公司 | 磁共振定量信息图的获取方法和装置 |
CN108802648A (zh) * | 2018-04-03 | 2018-11-13 | 上海东软医疗科技有限公司 | 一种基于梯度回波的磁共振定量成像方法和装置 |
CN108896943A (zh) * | 2018-05-10 | 2018-11-27 | 上海东软医疗科技有限公司 | 一种磁共振定量成像方法和装置 |
CN109100669A (zh) * | 2018-07-12 | 2018-12-28 | 厦门大学 | 基于重叠回波的单扫描同步磁共振扩散及t2成像方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10054653B2 (en) * | 2016-05-03 | 2018-08-21 | Siemens Healthcare Gmbh | Magnetic resonance method and apparatus for quantitative simultaneous multi-slice assessment of tissue displacement, deformation, and related biomarker parameters |
JP7039194B2 (ja) * | 2017-06-26 | 2022-03-22 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置および医用画像処理装置 |
-
2019
- 2019-02-01 CN CN201910105084.8A patent/CN109814056B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102759723A (zh) * | 2012-04-11 | 2012-10-31 | 华东师范大学 | 一种磁共振t2*图像生成方法 |
JP2015188635A (ja) * | 2014-03-28 | 2015-11-02 | 株式会社日立製作所 | 磁気共鳴イメージング装置 |
CN105699923A (zh) * | 2015-09-12 | 2016-06-22 | 北京大学 | 一种无创动态测量组织r2,r2*和r2`参数图像的磁共振成像方法 |
CN108294753A (zh) * | 2018-01-03 | 2018-07-20 | 上海东软医疗科技有限公司 | 磁共振定量信息图的获取方法和装置 |
CN108802648A (zh) * | 2018-04-03 | 2018-11-13 | 上海东软医疗科技有限公司 | 一种基于梯度回波的磁共振定量成像方法和装置 |
CN108896943A (zh) * | 2018-05-10 | 2018-11-27 | 上海东软医疗科技有限公司 | 一种磁共振定量成像方法和装置 |
CN109100669A (zh) * | 2018-07-12 | 2018-12-28 | 厦门大学 | 基于重叠回波的单扫描同步磁共振扩散及t2成像方法 |
Non-Patent Citations (3)
Title |
---|
Single-Shot T2 Mapping Through OverLapping-Echo Detachment (OLED) Planar Imaging;Congbo Cai等;《IEEE Transactions on Biomedical Engineering》;20171031;第64卷(第10期);第2450-2461页 * |
定量磁化率成像的基本原理及方法概述;郑志伟等;《磁共振成像》;20160630;第7卷(第6期);第454-460页 * |
心肌纵向弛豫时间和细胞外容积定量MR成像;陈铟铟等;《放射学时间》;20160831;第31卷(第8期);第794-798页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109814056A (zh) | 2019-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107024670B (zh) | 磁共振系统的校正方法及装置 | |
CN103238082B (zh) | 使用多点Dixon技术和低分辨率校准的MR成像 | |
US7888935B1 (en) | K-space trajectory estimation in spiral MRI system and related method thereof | |
US20170016972A1 (en) | Fast Prospective Motion Correction For MR Imaging | |
CN102914754B (zh) | 用于产生磁共振照片的方法 | |
US20130187649A1 (en) | System for Motion Corrected MR Diffusion Imaging | |
JPH0213434A (ja) | 磁気共鳴イメージング方法およびシステム | |
CA2964945A1 (en) | System and method for producing distortion free magnetic resonance images using dual-echo echo-planar imaging | |
US20110025327A1 (en) | Method for radiofrequency mapping in magnetic resonance imaging | |
US20190324102A1 (en) | System and method for stimulated echo based mapping (stem) using magnetic resonance imaging (mri) | |
US7395108B2 (en) | System for concurrent MRI imaging and magnetic field homogeneity measurement | |
WO2017132182A1 (en) | Systems and methods for joint reconstruction of quantitative t2 and adc maps | |
US9316711B2 (en) | System and method for accelerated magnetic resonance imaging using spectral sensitivity | |
CN109814056B (zh) | 一种获取磁共振精准定量图像的方法 | |
US8928317B2 (en) | System and method for controlling apparent timing dependencies for T2-weighted MRI imaging | |
US10228434B2 (en) | Multi-shot echo planar imaging using reordered segments and RF excitation pulse phase and slice profiles matched across interleaves | |
US11402454B2 (en) | Correction of distorted diffusion-weighted magnetic resonance image data | |
US10254367B2 (en) | Magnetic resonance imaging method and apparatus with motion-corrected model-based acceleration of parameter mapping | |
US11874353B2 (en) | Multi-shot echo planar imaging using reordered segments and recursive radio frequency pulse design giving matched slice profiles across segments | |
Latta et al. | Influence of k-space trajectory corrections on proton density mapping with ultrashort echo time imaging: Application for imaging of short T2 components in white matter | |
Milford et al. | A novel method for T2 quantification in presence of B1 inhomogeneities | |
JPH0549610A (ja) | 磁気共鳴診断装置 | |
WO2020223567A1 (en) | Simultaneous multi-slice mrsi using density weighted concentric ring acquisition | |
US20240331225A1 (en) | Methods, apparatuses, systems, and computer-readable mediums for biomarker quantification using free-breathing stack-of-radial imaging | |
Kowalik et al. | Rapid breath‐hold assessment of myocardial velocities using spiral UNFOLD‐ed SENSE tissue phase mapping |
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 |