CN103110421A - 提高磁共振成像信噪比的方法及其*** - Google Patents
提高磁共振成像信噪比的方法及其*** Download PDFInfo
- Publication number
- CN103110421A CN103110421A CN2013100488928A CN201310048892A CN103110421A CN 103110421 A CN103110421 A CN 103110421A CN 2013100488928 A CN2013100488928 A CN 2013100488928A CN 201310048892 A CN201310048892 A CN 201310048892A CN 103110421 A CN103110421 A CN 103110421A
- Authority
- CN
- China
- Prior art keywords
- noise ratio
- signal
- burst length
- repetition
- repetition time
- 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 28
- 238000002595 magnetic resonance imaging Methods 0.000 title claims abstract description 20
- 238000011084 recovery Methods 0.000 claims abstract description 101
- 238000005481 NMR spectroscopy Methods 0.000 claims description 42
- 238000003384 imaging method Methods 0.000 claims description 40
- 238000002059 diagnostic imaging Methods 0.000 abstract description 2
- 229910052698 phosphorus Inorganic materials 0.000 description 8
- 239000011574 phosphorus Substances 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000005415 magnetization Effects 0.000 description 7
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 6
- 238000005457 optimization Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000002950 deficient Effects 0.000 description 2
- -1 phosphorus compound Chemical class 0.000 description 2
- 230000002490 cerebral effect Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 210000004884 grey matter Anatomy 0.000 description 1
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明适用于医学影像技术领域,提供了一种提高磁共振成像信噪比的方法及其***,所述方法包括如下步骤:信噪比信噪比曲线绘制步骤:绘制单位时间内的信噪比-恢复脉冲时间-重复时间曲线;获取参数范围步骤:在所述曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围;确定参数值步骤:根据所述恢复脉冲时间范围及重复时间范围,以及给定的本次扫描时间设置本次扫描的恢复脉冲时间、重复时间以及重复次数。借此,本发明提高了磁共振成像的信噪比。
Description
技术领域
本发明涉及医学影像技术领域,尤其涉及一种提高磁共振成像信噪比的方法及其***。
背景技术
自上世纪80年代初年第一台医用核磁共振成像仪问世以来,核磁共振成像作为一项重要的临床检查手段得到了飞速的发展,在临床诊断工作中发挥着重要作用。同时,核磁共振的成像技术也与日俱进,从成像的对象,设备硬件,脉冲序列,到图像重建,都涌现出很多新的技术。其中磁共振成像脉冲序列的发展处于核心地位,正是由于脉冲序列的发展,在医学磁共振成像仪的应用上才能提高给诊断医生各方面的图像信息。
在磁共振成像脉冲序列中,有一类序列被称为驱动平衡序列,在有些文献中也被称为快速恢复脉冲序列。该类序列在进行数据采集的读出梯度施加完毕后,给予一个180度Y脉冲让剩余的磁化强度矢量重聚,在回波时间在施加一个90度-X脉冲让水平方向上重聚的磁化强度矢量翻转到竖直方向,即被驱使回到磁化强度矢量的初始方向。驱动平衡序列可以与GRE序列、SE序列、FSE序列等序列相结合,以缩短这些序列在获取信号后恢复到平衡态的时间。基于GRE序列的驱动平衡脉冲序列波形示意图如图1所示。
驱动平衡序列的恢复效率与成像参数Tfr,Tfr为第一个90°激励脉冲与最后一个驱动恢复的90°脉冲之间的时间间隔,即恢复脉冲时间、TR(repetition time,重复时间)及检查对象的弛豫常数T1(检查对象的弛豫常数),T2(检查对象的弛豫常数)有关。通常情况下对于T2在100ms一以下,T1在2000ms以下,驱动平衡的恢复效率采用短Tfr与短TR会相对较高。参考文献1:Sylvain Miraux,EricThiaudi_ere,Paul Canioni,et al,Magnetization recovery for signal enhancement:afast imaging DEFT-based technique,J Mag Res,2004(166):28-34(用于信号增强的磁化强度矢量恢复:依照基于驱动平衡的快速成像技术,磁共振杂志,2004(166):28-34.)详细介绍了驱动平衡脉冲序列对信噪比的提升效率与成像参数Tfr、TR,弛豫时间T1与T2之间的关系。但是,当驱动平衡序列应用于信号强度很弱的磷成像时,由于磷化合物T1较长,浓度很低,造成读出带宽与TR对信噪比影响很大,在确定驱动平衡参数时需要考虑读出带宽,T2*衰减,TR以及重复次数之间的关系。
综上可知,现有的磁共振成像技术在实际使用上,显然存在不便与缺陷,所以有必要加以改进。
发明内容
针对上述的缺陷,本发明的目的在于提供一种提高磁共振成像信噪比的方法及其***,以提高磁共振成像的信噪比。
为了实现上述目的,本发明提供一种提高磁共振成像信噪比的方法,包括如下步骤:
信噪比曲线绘制步骤:绘制单位时间内的信噪比-恢复脉冲时间-重复时间曲线;
获取参数范围步骤:在所述曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围;
确定参数值步骤:根据所述恢复脉冲时间范围及重复时间范围,以及给定的本次扫描时间设置本次扫描的恢复脉冲时间、重复时间以及重复次数。
根据所述的方法,在所述确定参数值步骤之后还包括:
在所述给定的本次扫描时间内,根据所述本次扫描的恢复脉冲时间、重复时间以及重复次数执行磁共振扫描。
根据所述的方法,所述信噪比曲线绘制步骤包括:
以解布洛赫方程为基础,获取所述信噪比与恢复脉冲时间、重复时间以及重复次数的关系;
根据所述关系,绘制已修正所述恢复脉冲时间和重复时间因素后的单位时间内的二维信噪比-恢复脉冲时间-重复时间信噪比曲线;
所述获取参数范围步骤包括:
在所述信噪比-恢复脉冲时间-重复时间信噪比曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围。
根据所述的方法,所述确定参数值步骤包括:
在所述极大值区域中的查找所述信噪比的最优区域,并在所述最优区域的任意一点读取对应的恢复脉冲时间和重复时间;或者
在所述极大值区域中的查找所述信噪比最高点,并在所述最高点读取对应的恢复脉冲时间和重复时间;
根据所述对应的恢复脉冲时间和重复时间,以及给定的所述本次扫描时间设置本次扫描的重复次数。
根据所述的方法,所述磁共振成像为驱动平衡脉冲序列成像。
为了实现本发明的另一发明目的本发明还提供了一种提高磁共振成像信噪比的***,包括:
信噪比曲线绘制模块,用于绘制单位时间内的信噪比-恢复脉冲时间-重复时间曲线;
获取参数范围模块,用于在所述曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围;
确定参数值模块,用于根据所述恢复脉冲时间范围及重复时间范围,以及给定的本次扫描时间设置本次扫描的恢复脉冲时间、重复时间以及重复次数。
根据所述的***,所述***还包括:
扫描模块,用于在所述给定的本次扫描时间内,根据所述本次扫描的恢复脉冲时间、重复时间以及重复次数执行磁共振扫描。
根据所述的***,所述信噪比曲线绘制模块包括:
参数关系确定子模块,用于以解布洛赫方程为基础,获取所述信噪比与恢复脉冲时间、重复时间以及重复次数的关系;
信噪比曲线绘制子模块,用于根据所述关系,绘制已修正所述恢复脉冲时间和重复时间因素后的单位时间内的二维信噪比-恢复脉冲时间-重复时间信噪比曲线;
所述获取参数模块在所述信噪比-恢复脉冲时间-重复时间信噪比曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围。
根据所述的***,所述确定参数值模块包括:
第一参数值确定子模块,用于在所述极大值区域中的查找所述信噪比的最优区域,并在所述最优区域的任意一点读取对应的恢复脉冲时间和重复时间;或者
第二参数值确定子模块,用于在所述极大值区域中的查找所述信噪比最高点,并在所述最高点读取对应的恢复脉冲时间和重复时间;
重复次数确定子模块,用于根据所述对应的恢复脉冲时间和重复时间,以及给定的所述本次扫描时间设置本次扫描的重复次数。
根据所述的***,所述磁共振成像为驱动平衡脉冲序列成像。
本发明通过绘制单位时间内的信噪比-恢复脉冲时间-重复时间曲线;然后在所述曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围;并根据所述恢复脉冲时间范围及重复时间范围,以及给定的本次扫描时间设置本次扫描的恢复脉冲时间、重复时间以及重复次数。由此实现了通过获得最佳恢复脉冲施加时间、读出带宽、重复时间和重复次数的组合来达到提高信噪比的目的。此外,在驱动平衡类脉冲序列磁共振成像中,给出了驱动平衡GRE序列单位时间信噪比的最佳分布区域;在该分布指导下可以改善成像质量,缩短扫描时间。具体的,参考恢复脉冲时间与重复时间对成像信噪比的影响,以及二者又和读出梯度与重复次数相互约束的关系,利用布洛赫方程仿真为基础,绘制以修正读出带宽和重复次数因素后信噪比-恢复脉冲时间-重复时间二维单位时间信噪比曲线,在该曲线上获得信噪比最高点的恢复脉冲时间-读出带宽,重复时间-重复次数组合,达到在给定时间内获得最佳的信噪比目的。该方法对于弱成像信号,如磷成像等多核成像有明显改善信噪比作用。
附图说明
图1是现有技术中驱动平衡梯度回波序列的示意图;
图2是本发明第一实施例提供的提高磁共振成像信噪比的***结构示意图;
图3是本发明第二、三、四实施例提供的提高磁共振成像信噪比的***结构示意图;
图4是本发明第五实施例提供的提高磁共振成像信噪比的方法流程图;
图5是本发明一个实施例提供的单位时间信噪比-快速恢复脉冲时间-重复时间二维曲线图
图6是本发明一个实施例提供的无机磷成像单位时间信噪比-快速恢复脉冲时间-重复时间二维曲线图
图7是本发明一个实施例提供的GE3T Signa扫描机执行SPSP_DE_GRE脉冲序列的时序图
图8是本发明一个实施例提供的采用优化参数SPSP_DE_GRE脉冲扫描结果示意图;
图9是现有技术中SPSP_GRE脉冲扫描结果示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
参见图2,在本发明的第一实施例中,提供了一种提高磁共振成像信噪比的***100,包括:
信噪比曲线绘制模块10,用于绘制单位时间内的信噪比-恢复脉冲时间-重复时间曲线;
获取参数范围模块20,用于在所述曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围;
确定参数值模块30,用于根据所述恢复脉冲时间范围及重复时间范围,以及给定的本次扫描时间设置本次扫描的恢复脉冲时间、重复时间以及重复次数。
在该实施例中,首先通过信噪比曲线绘制模块10绘制单位时间内的信噪比-恢复脉冲时间-重复时间曲线;该曲线是归一化单位时间信噪比-恢复脉冲时间-重复时间曲线的绘制。然后,获取参数范围模块20则在所述曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围,该范围包括恢复脉冲时间及重复时间的最优取值区域。最后,由确定参数值模块30根据所述恢复脉冲时间范围及重复时间范围,以及给定的本次扫描时间设置本次扫描的恢复脉冲时间、重复时间以及重复次数。由于所述恢复脉冲时间范围及重复时间范围为恢复脉冲时间及重复时间的最优取值区域,因此最后确定的本次扫描的恢复脉冲时间、重复时间以及重复次数也将是优化获取的,因此通过上述提高磁共振成像信噪比的***100的工作优化了磁共振成像参数,以最终达到提高磁共振成像信噪比的目的。优选的,所述磁共振成像为驱动平衡脉冲序列成像。
参见图3,在本发明的第二实施例中,提高磁共振成像信噪比的***100还包括:
扫描模块40,用于在所述给定的本次扫描时间内,根据所述本次扫描的恢复脉冲时间、重复时间以及重复次数执行磁共振扫描。
在该实施例中,扫描模块40则根据所述本次扫描的恢复脉冲时间、重复时间以及重复次数,在所述给定的本次扫描时间内执行磁共振扫描,由此可以在具体的磁共振成像过程中,实现磁共振成像信噪比的优化。
参见图3,在本发明的第三实施例中,信噪比曲线绘制模块10包括:
参数关系确定子模块11,用于以解Bloch(布洛赫)方程为基础,获取所述信噪比与恢复脉冲时间、重复时间以及重复次数的关系;
信噪比曲线绘制子模块12,用于根据所述关系,绘制已修正所述恢复脉冲时间和重复时间因素后的单位时间内的二维信噪比-恢复脉冲时间-重复时间信噪比曲线;
获取参数模块20在所述信噪比-恢复脉冲时间-重复时间信噪比曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围。
在该实施例中,参数关系确定子模块11是以解Bloch方程为基础获取磁共振成像信噪比的相关参数。Bloch方程是经典力学描述核磁共振现象最为重要的理论基础之一。通过解Bloch方程可以修正所述恢复脉冲时间和重复时间与信噪比的关系。由信噪比曲线绘制子模块12根据参数关系确定子模块11绘制相关的单位时间内的二维信噪比-恢复脉冲时间-重复时间信噪比曲线。
参见图3,在本发明的第四实施例中,确定参数值模块30包括:
第一参数值确定子模块31,用于在所述极大值区域中的查找所述信噪比的最优区域,并在所述最优区域的任意一点读取对应的恢复脉冲时间和重复时间;或者
第二参数值确定子模块32,用于在所述极大值区域中的查找所述信噪比最高点,并在所述最高点读取对应的恢复脉冲时间和重复时间;
重复次数确定子模块33,用于根据所述对应的恢复脉冲时间和重复时间,以及给定的所述本次扫描时间设置本次扫描的重复次数。
在该实施例中,第一参数值确定子模块31可以将所述信噪比曲线中的极大值区域中的部分区域确定为所述信噪比的最优区域,具体的确定可以根据操作者设置的具体取值参数设定。在该最优区域的任意一点对应的恢复脉冲时间和重复时间可以读取为本次磁共振扫描的恢复脉冲时间和重复时间参数。此外,优选的,第二参数值确定子模块32将所述极大值区域中的查找所述信噪比最高点对应的恢复脉冲时间和重复时间读取为本次磁共振扫描的恢复脉冲时间和重复时间参数;最后重复次数确定子模块33可以根据上述两个子模块确定的恢复脉冲时间和重复时间设置本次磁共振扫描的重复次数。
在上述多个实施例中,提高磁共振成像信噪比的***100可以软件单元,硬件单元或软硬件结合单元。并且该***可以设置于磁共振扫描装置中。
参见图4,在本发明的第五实施例中,提供了一种提高磁共振成像信噪比的方法,包括如下步骤:
步骤S401中,信噪比曲线绘制模块10绘制单位时间内的信噪比-恢复脉冲时间-重复时间曲线;该步骤为信噪比曲线绘制步骤;
步骤S402中,获取参数范围模块20在所述曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围;该步骤为获取参数范围步骤:
步骤S403中,确定参数值模块30根据所述恢复脉冲时间范围及重复时间范围,以及给定的本次扫描时间设置本次扫描的恢复脉冲时间、重复时间以及重复次数;该步骤为确定参数值步骤。
在该实施例中,通过获取优化磁共振参数,使得在信噪比较低的非质子成像中能够显著提高成像的时间效率,改善成像质量。
在本发明的第六实施例中,在所述步骤S403之后还包括:
扫描模块40在所述给定的本次扫描时间内,根据所述本次扫描的恢复脉冲时间、重复时间以及重复次数执行磁共振扫描。
在本发明的第六实施例中,在所述步骤S401包括:
参数关系确定子模块11以解Bloch方程为基础,获取所述信噪比与恢复脉冲时间、重复时间以及重复次数的关系;
信噪比曲线绘制子模块12根据所述关系,绘制已修正所述恢复脉冲时间和重复时间因素后的单位时间内的二维信噪比-恢复脉冲时间-重复时间信噪比曲线;
所述步骤S402包括:
获取参数模块20在所述信噪比-恢复脉冲时间-重复时间信噪比曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围。
在该实施例中,磁共振成像序列为驱动平衡序列。基于当驱动平衡序列应用于信号强度很弱的磷成像时,由于磷化合物T1较长,浓度很低,造成读出带宽与TR对信噪比影响很大,在确定驱动平衡参数时需要考虑读出带宽,T2*衰减,TR以及重复次数之间的关系。考虑提出信噪比与Tfr,TR,T1及T2的修正关系,在此关系基础上确定最优磁共振扫描参数。具体的方法实现过程如下:
首先,进行归一化单位时间信噪比-恢复脉冲时间-重复时间曲线的绘制
(1)SNR(Signal/Noise,信噪比)与各成像参数关系;
在存在驱动平衡脉冲情况下,通过解Bloch方程,可得
式中Mxy是稳定状态水平方向磁化强度矢量,M0是初始状态磁化强度矢量,TE是回波时间。由于MRI(Magnetic Resonance Imaging,磁共振成像)的信号强度与Mxy成正比,而噪音水平与上式中各参数无关,所以可以用上式代表驱动平衡磁共振成像的信噪比。
由于在驱动平衡脉冲序列中Tfr受到读出带宽带宽的限制,在其余成像条件相同的情况下,磁共振信噪比与读出带宽的关系为:
带宽与读出时间窗有如下关系:
在预相散梯度,相位编码梯度取硬件***所能提供的最小值情况下,最小回波时间TE可以视作:
TE=0.5×Tread+ΔTmin (公式4)
而在重聚梯度,各平衡梯度均取硬件***所能提供的最小值情况下,读出时间窗与Tfr关系近似为:
Tfr=(Tread+ΔTmin)×2 (公式5)
忽略掉ΔT,TE取最小回波时间,有
TE=Tfr/4 (公式7)
再考虑NEX(number ofexcitation,重复次数)与TR之间的相互约束,这种约束以单位时间信噪比uSNR表示:
综上所述,单位时间信噪比满足以下关系:
公式中的T2*依赖T2、成像时线圈内的对象和匀场的好坏,与T2之间具有如下关系:
可知
在匀场较好的情况下分母中的第二项约在10ms左右。
(2)编程绘制曲线;
根据以上关系,在成像对象为某种具体T1,T2,T2*值时(如在3TMRI扫描机上,大脑灰质T1=1820ms,T2=99ms),以TR与Tfr为自变量,单位时间信噪比为应变量,绘制二维曲线图如图5所示,图中上方的深灰色部分为单位信噪比最优部分。每一色彩灰度级为0.001/ms。
(3)T1,T2对二维曲线图的影响;
SNR/UnitTime-Tfr-TR分布图随T1,T2的不同而改变。也正因为如此,针对不同的T1,T2,分布图可以帮助确定最优的Tfr与TR。
在本发明的第七实施例中,所述步骤S402包括:
第一参数值确定子模块31在所述极大值区域中的查找所述信噪比的最优区域,并在所述最优区域的任意一点读取对应的恢复脉冲时间和重复时间;或者
第二参数值确定子模块32在所述极大值区域中的查找所述信噪比最高点,并在所述最高点读取对应的恢复脉冲时间和重复时间;
重复次数确定子模块33根据所述对应的恢复脉冲时间和重复时间,以及给定的所述本次扫描时间设置本次扫描的重复次数。
在本发明的一个实施例中,通过获得最佳恢复脉冲时间(Tfr)、读出带宽、重复时间TR和重复次数(NEX)的组合来达到提高信噪比的目的。在驱动平衡类脉冲序列磁共振成像中,Tft与TR对成像信噪比(SNR)有很大影响,同时二者又和读出梯度与重复次数相互约束。以Bloch方程仿真为基础,绘制以修正读出带宽和NEX因素后SNR-Tft&TR二维单位时间信噪比曲线,在该曲线上获得SNR最高点的Tfr-读出带宽,TR-NEX组合,达到在给定时间内获得最佳SNR的目的。该技术方案对于弱成像信号,如磷成像等多核成像有明显改善SNR作用。具体的选取极值区域以就近原则确定最优成像参数的操作为:在给定的成像时间和成像矩阵情况下,以图5为参考,在深灰色区域内选择合适的Tfr与TR,结合给定成像时间选择重复次数。优选的,该深灰色区域包括TR的值在200~1800的范围区域。
参见图6~图9,在本发明的一个优选实施例中,在该实施例中是驱动平衡脉冲序列无机磷磁共振成像。成像序列的执行平台是3T GE Signa HDx scanner(GEHealthcare,Waukesha,WI)。扫描对象为半径6cm的100mM Na2HPO4球形水模。执行过程如下:
1、使用质子线圈做3_plane定位相扫描;
2、切换至磷线圈,使用GE2D fidcsi序列中的Spectral Prescan选项进行预扫描,成像模式选择为波谱(Mode=1)。通过参数设置将所用序列设置为一频率选择性射频脉冲。调节匀场梯度,使信号衰减最为缓慢,同时估算出T2*(也可由其它序列测出T2*或由相关经验获得);
3、将无机磷T1,T2和上述T2*代入公式,绘制出曲线图。3T下无机磷Pi的T1=6900ms,T2=153ms。T2*实验初步测量为60ms。绘制为单位时间信噪比-快速恢复脉冲时间-重复时间二维曲线图如图6所示。
4、从图6中获取极大值区域,优选的,极大值区域包括图6中的深灰色区域,即TR的值在500~3500的范围区域。最后确定Tfr与TR的最优区域在TR为1000到1800之间,Tfr为20到40之间;
5、加载具有选择性激励的SPSP_DE_GRE脉冲,根据扫描时间(此处假定不超过6分钟)和最优Tfr与TR范围,设置TR,NEX,Tfr等相关参数。考虑到较大的信号强度易于观察,取TR为1800ms,对应Tfr取为35ms,Matrix32*32,NEX为6。总扫描时间5分46秒;
6、其中加载SPSP_DE_GRE脉冲,脉冲序列如图7所示:
7、使用没有驱动平衡模块的SPSP_GRE脉冲,TR=6000ms,NEX=2,成像时间6分24秒,扫描图像信噪比对比,结果分别图8和图9所示,图8中采用优化参数SPSP_DE_GRE脉冲扫描得到的结果,SNR=12;图9中用于对比的SPSP_GRE脉冲得到的结果,SNR=8。由此可见本发明提供的提高磁共振成像信噪比的方法及其***可以明显的提高磁共振成像信噪比。
在本发明的一个实施例中,提高磁共振成像驱动平衡脉冲序列信噪比的方法,通过下述多个步骤实现:(a)SNR-Tft&TR二维单位时间信噪比曲线的绘制;(b)利用SNR-Tft&TR二维单位时间信噪比曲线获取最佳恢复脉冲施加时间(Tfr)与读出带宽组合和重复时间TR与重复次数(NEX)组合;具体的在(a)中在一定TR情况下读出带宽(或读出时间)修正SNR-Tft曲线;以及在(a)中在一定Tfr情况下读出单位时间信噪比修正(Unit time SNR)-Tft曲线;在(b)一定TR情况下通过读出带宽(或读出时间)修正SNR-Tft曲线提高MRI信噪比的方法;以及在(b)一定Tfr情况下通过单位时间信噪比(Unit time SNR)-Tft曲线提高MRI信噪比的方法。
综上所述,本发明通过绘制单位时间内的信噪比-恢复脉冲时间-重复时间曲线;然后在所述曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围;并根据所述恢复脉冲时间范围及重复时间范围,以及给定的本次扫描时间设置本次扫描的恢复脉冲时间、重复时间以及重复次数。由此实现了通过获得最佳恢复脉冲施加时间、读出带宽、重复时间和重复次数的组合来达到提高信噪比的目的。此外,在驱动平衡类脉冲序列磁共振成像中,给出了驱动平衡GRE序列单位时间信噪比的最佳分布区域;在该分布指导下可以改善成像质量,缩短扫描时间。具体的,参考恢复脉冲时间与重复时间对成像信噪比的影响,以及二者又和读出梯度与重复次数相互约束的关系,利用Bloch方程仿真为基础,绘制以修正读出带宽和重复次数因素后信噪比-恢复脉冲时间-重复时间二维单位时间信噪比曲线,在该曲线上获得信噪比最高点的恢复脉冲时间-读出带宽,重复时间-重复次数组合,达到在给定时间内获得最佳的信噪比目的。该方法对于弱成像信号,如磷成像等多核成像有明显改善信噪比作用。
当然,本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (10)
1.一种提高磁共振成像信噪比的方法,其特征在于,包括如下步骤:
信噪比曲线绘制步骤:绘制单位时间内的信噪比-恢复脉冲时间-重复时间曲线;
获取参数范围步骤:在所述曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围;
确定参数值步骤:根据所述恢复脉冲时间范围及重复时间范围,以及给定的本次扫描时间设置本次扫描的恢复脉冲时间、重复时间以及重复次数。
2.根据权利要求1所述的方法,其特征在于,在所述确定参数值步骤之后还包括:
在所述给定的本次扫描时间内,根据所述本次扫描的恢复脉冲时间、重复时间以及重复次数执行磁共振扫描。
3.根据权利要求1所述的方法,其特征在于,所述信噪比曲线绘制步骤包括:
以解布洛赫方程为基础,获取所述信噪比与恢复脉冲时间、重复时间以及重复次数的关系;
根据所述关系,绘制已修正所述恢复脉冲时间和重复时间因素后的单位时间内的二维信噪比-恢复脉冲时间-重复时间信噪比曲线;
所述获取参数范围步骤包括:
在所述信噪比-恢复脉冲时间-重复时间信噪比曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围。
4.根据权利要求3所述的方法,其特征在于,所述确定参数值步骤包括:
在所述极大值区域中的查找所述信噪比的最优区域,并在所述最优区域的任意一点读取对应的恢复脉冲时间和重复时间;或者
在所述极大值区域中的查找所述信噪比最高点,并在所述最高点读取对应的恢复脉冲时间和重复时间;
根据所述对应的恢复脉冲时间和重复时间,以及给定的所述本次扫描时间设置本次扫描的重复次数。
5.根据权利要求1~4任一项所述的方法,其特征在于,所述磁共振成像为驱动平衡脉冲序列成像。
6.一种提高磁共振成像信噪比的***,其特征在于,包括:
信噪比曲线绘制模块,用于绘制单位时间内的信噪比-恢复脉冲时间-重复时间曲线;
获取参数范围模块,用于在所述曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围;
确定参数值模块,用于根据所述恢复脉冲时间范围及重复时间范围,以及给定的本次扫描时间设置本次扫描的恢复脉冲时间、重复时间以及重复次数。
7.根据权利要求6所述的***,其特征在于,所述***还包括:
扫描模块,用于在所述给定的本次扫描时间内,根据所述本次扫描的恢复脉冲时间、重复时间以及重复次数执行磁共振扫描。
8.根据权利要求6所述的***,其特征在于,所述信噪比曲线绘制模块包括:
参数关系确定子模块,用于以解布洛赫方程为基础,获取所述信噪比与恢复脉冲时间、重复时间以及重复次数的关系;
信噪比曲线绘制子模块,用于根据所述关系,绘制已修正所述恢复脉冲时间和重复时间因素后的单位时间内的二维信噪比-恢复脉冲时间-重复时间信噪比曲线;
所述获取参数模块在所述信噪比-恢复脉冲时间-重复时间信噪比曲线上的极大值区域获取恢复脉冲时间范围及重复时间范围。
9.根据权利要求8所述的***,其特征在于,所述确定参数值模块包括:
第一参数值确定子模块,用于在所述极大值区域中的查找所述信噪比的最优区域,并在所述最优区域的任意一点读取对应的恢复脉冲时间和重复时间;或者
第二参数值确定子模块,用于在所述极大值区域中的查找所述信噪比最高点,并在所述最高点读取对应的恢复脉冲时间和重复时间;
重复次数确定子模块,用于根据所述对应的恢复脉冲时间和重复时间,以及给定的所述本次扫描时间设置本次扫描的重复次数。
10.根据权利要求6~9任一项所述的***,其特征在于,所述磁共振成像为驱动平衡脉冲序列成像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310048892.8A CN103110421B (zh) | 2013-02-07 | 2013-02-07 | 提高磁共振成像信噪比的方法及其*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310048892.8A CN103110421B (zh) | 2013-02-07 | 2013-02-07 | 提高磁共振成像信噪比的方法及其*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103110421A true CN103110421A (zh) | 2013-05-22 |
CN103110421B CN103110421B (zh) | 2016-02-10 |
Family
ID=48408820
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310048892.8A Expired - Fee Related CN103110421B (zh) | 2013-02-07 | 2013-02-07 | 提高磁共振成像信噪比的方法及其*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103110421B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104042216A (zh) * | 2014-07-01 | 2014-09-17 | 中国科学院武汉物理与数学研究所 | 一种基于预扫描和非均匀采样的薄层快速磁共振成像方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1698537A (zh) * | 2004-05-19 | 2005-11-23 | 西北大学 | 冠状动脉的时间分辨的对比度增强的磁共振投影成像 |
CN102116856A (zh) * | 2010-12-30 | 2011-07-06 | 中国科学院深圳先进技术研究院 | 横向弛豫时间测量方法及*** |
US20120268126A1 (en) * | 2011-04-22 | 2012-10-25 | The Regents Of The University Of California | Increasing the perfusion signal to noise ratio (snr) in velocity selective arterial spin labeling (vsasl) in magnetic resonance imaging |
-
2013
- 2013-02-07 CN CN201310048892.8A patent/CN103110421B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1698537A (zh) * | 2004-05-19 | 2005-11-23 | 西北大学 | 冠状动脉的时间分辨的对比度增强的磁共振投影成像 |
CN102116856A (zh) * | 2010-12-30 | 2011-07-06 | 中国科学院深圳先进技术研究院 | 横向弛豫时间测量方法及*** |
US20120268126A1 (en) * | 2011-04-22 | 2012-10-25 | The Regents Of The University Of California | Increasing the perfusion signal to noise ratio (snr) in velocity selective arterial spin labeling (vsasl) in magnetic resonance imaging |
Non-Patent Citations (1)
Title |
---|
彭光明等: "磁共振成像的信噪比:计量方法与影像因素", 《影像技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104042216A (zh) * | 2014-07-01 | 2014-09-17 | 中国科学院武汉物理与数学研究所 | 一种基于预扫描和非均匀采样的薄层快速磁共振成像方法 |
CN104042216B (zh) * | 2014-07-01 | 2015-12-30 | 中国科学院武汉物理与数学研究所 | 一种基于预扫描和非均匀采样的薄层快速磁共振成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103110421B (zh) | 2016-02-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9041394B2 (en) | Magnetic resonance imaging apparatus executing phase-corrected imaging pulse sequence based on data obtained from multiple pulse sub-sequences executed without readout direction gradients but instead using phase or slice encoding direction gradients during readout | |
US20070249929A1 (en) | Systems and methods for magnetic resonance imaging | |
US9046590B2 (en) | Magnetic resonance imaging apparatus phase correction using one or more prescans with identical readout and slice gradients | |
US9541624B2 (en) | Method and apparatus for magnetic resonance imaging | |
CN102914754A (zh) | 用于产生磁共振照片的方法 | |
CN102713658B (zh) | 磁化率梯度绘图 | |
US10001537B2 (en) | Magnetic resonance imaging apparatus and image processing apparatus | |
CN101744616B (zh) | 用于mri中的化学饱和的发射频率校准方法及mri*** | |
US20040227510A1 (en) | Magnetic resonance diffusion imaging with eddy-current compensation | |
JP3866537B2 (ja) | 磁気共鳴撮影装置 | |
EP3078979A2 (de) | Geschwindigkeitskompensierte mr-diffusionsbildgebung | |
CN113017596B (zh) | 一种磁共振多参数定量方法及其应用 | |
CN103083020A (zh) | 磁共振成像设备及其控制方法 | |
US7642778B2 (en) | MRI apparatus | |
JP6762284B2 (ja) | 磁気共鳴イメージング装置およびノイズ除去方法 | |
US10067209B2 (en) | Magnetic resonance imaging apparatus | |
JP3964110B2 (ja) | 磁気共鳴イメージング装置 | |
US10126395B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
DE102017207128B4 (de) | Verfahren zur Aufnahme eines Magnetresonanzdatensatzes, Datenträger sowie Magnetresonanzanlage | |
EP0429295A2 (en) | Magnetic resonance spectroscopy | |
CN103635825A (zh) | 交错黑血和亮血动态对比增强(dce)mri | |
US8148980B2 (en) | Magnetic resonance imaging system and method | |
US7956611B2 (en) | Magnetic resonance imaging apparatus and method | |
CN104155621B (zh) | 一种准确测量静磁场b0分布的方法 | |
US8547101B2 (en) | Magnetic resonance imaging device |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160210 Termination date: 20210207 |
|
CF01 | Termination of patent right due to non-payment of annual fee |