CN103976735A - 基于磁共振的黑血电影成像方法 - Google Patents

基于磁共振的黑血电影成像方法 Download PDF

Info

Publication number
CN103976735A
CN103976735A CN201410184612.0A CN201410184612A CN103976735A CN 103976735 A CN103976735 A CN 103976735A CN 201410184612 A CN201410184612 A CN 201410184612A CN 103976735 A CN103976735 A CN 103976735A
Authority
CN
China
Prior art keywords
proton
target area
gradient
motion
imaging
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
Application number
CN201410184612.0A
Other languages
English (en)
Other versions
CN103976735B (zh
Inventor
郭华
董莉
苑纯
张喆
戴二鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN201410184612.0A priority Critical patent/CN103976735B/zh
Publication of CN103976735A publication Critical patent/CN103976735A/zh
Application granted granted Critical
Publication of CN103976735B publication Critical patent/CN103976735B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明提出一种基于磁共振的黑血电影成像方法,包括:在磁共振的磁化准备阶段通过速度选择准备VESP序列抑制目标区域内的血流信号;通过心电门控同步心电信号;以及基于心电信号,对其中血流信号被抑制的目标区域进行图像采集并根据所采集到的第一图像生成目标区域的黑血电影图像。本发明实施例的基于磁共振的黑血电影成像方法,抑制了血流信号,消除了血液流动所诱发的流动伪影,从而能够清晰确定血管壁和血管腔之间的边界,实现对血管壁边界的精确描绘,最终实现对血管动态特性和功能特性的准确评估。

Description

基于磁共振的黑血电影成像方法
技术领域
本发明涉及磁共振技术领域,特别涉及一种基于磁共振的黑血电影成像方法。
背景技术
相关的研究表明,高血压、动脉粥样硬化、中风等心血管疾病的发生与动脉僵硬存在重要的联系。原因是随着动脉管壁的弹性降低,动脉存储血液的顺应性减弱,从而造成血压增加,增加心血管疾病发生的风险。因此,评价动脉管壁的动力学特性,对于相关的心血管疾病的预防、监测以及治疗,都具有重要的临床意义。
目前主要使用亮血电影成像技术研究动脉动力学特性。借助电影成像技术,可以获取心动周期内不同时间点的若干帧图像,通过测量管腔面积在心动周期内的变化,可以量化分析动脉动力学特性。但是,亮血电影成像受到两个方面的限制:容易受到流动伪影的干扰,影响定量分析的准确性;不能区分血管内外管壁,因此不能提供对应的血管壁信息。
相对于亮血成像技术,黑血成像技术具有更加优越的特性。黑血成像技术,是指通过特定的磁化准备序列,抑制成像区域的血流信号,使其最终在图像上表现为“暗信号”的成像方式。黑血成像技术主要有如下两个方面的优势:首先,由于血流信号被抑制,流动伪影大幅度减少;同时,血流与管壁之间的对比度增加,分界更加清晰,因而可以更加准确地评估血管的动力学特性。但是目前由于血流抑制效率、信噪比以及成像速度等因素的限制,难以同时满足抑制血流信号和快速获取动力学图像两个方面的要求,黑血技术仍局限应用于静态成像,提供结构学信息,而在电影成像中的应用受到了限制。
发明内容
本发明旨在至少解决现有技术中存在的技术问题之一。
有鉴于此,本发明需要提供一种基于磁共振的黑血电影成像方法,通过该方法抑制了血流信号,消除了血液流动所诱发的流动伪影,从而能够清晰确定血管壁和血管腔之间的边界,实现对血管壁边界的精确描绘,最终实现对血管动态特性和功能特性的准确评估。
根据本发明的一个的实施例,提出一种基于磁共振的黑血电影成像方法,包括:在所述磁共振的磁化准备阶段通过速度选择准备VESP序列抑制目标区域内的血流信号;通过心电门控同步心电信号;以及基于所述心电信号,对其中所述血流信号被抑制的目标区域进行图像采集并根据所采集到的第一图像生成所述目标区域的黑血电影图像。
根据本发明的一个的实施例,所述VESP序列包括T2准备脉冲序列、散相梯度以及在所述T2准备脉冲序列后的破坏梯度,其中,所述T2准备脉冲序列依次包括一个90°x脉冲,一个或多个180°y脉冲,以及一个90°-x脉冲;所述散相梯度在读出方向、相位编码方向和选层方向中的至少一个方向施加。
根据在本发明的一个实施例,所述通过速度选择准备VESP序列抑制目标区域内的血流信号,包括:将所述90°x脉冲作用于所述目标区域,以使所述目标区域的血流中运动的质子和所述目标区域的血管中静止的质子翻转至横截面;将所述一个或多个180°y脉冲和所述散相梯度共同作用于所述目标区域,以将翻转后的运动的质子散相;将所述90°-x脉冲作用于翻转后的静止的质子和散相后剩余的运动的质子,以将所述翻转后的静止的质子和所述散相后剩余的运动的质子翻转至垂直方向;将所述破坏梯度作用于所述目标区域以清除所述横截面残余的磁矩。
根据本发明的一个实施例,所述散相梯度为满足以下条件的任意流速敏感梯度:
m 0 = ∫ G ( t ) dt = 0 ,
以使所述散相梯度不引起所述静止质子的相位变化;以及
m 1 = ∫ G ( t ) · tdt ≠ 0 ,
以使所述散相梯度引起所述运动质子散相;其中,m0,m1分别为所述散相梯度的零阶梯度矩和一阶梯度矩。
根据本发明的一个实施例,所述目标区域内的运动的质子分别具有不同的流动速度,在所述将所述一个或多个180°y脉冲和所述散相梯度共同作用于所述目标区域的过程中,所述散相梯度作用于所述具有不同流动速度的运动的质子时,所述目标区域内单个像素内具有不同流动速度的运动的质子分别产生不同的相位变化,发生散相;当所述单个像素内的运动的质子发生散相时,所述单个像素内的运动的质子的磁矩相互抵消,所述单个像素对应的血流信号被抑制。
根据本发明的一个实施例,所述目标区域内单个运动质子的相位变化为:
φ = γ · v ( r → ) · ∫ G ( t ) · tdt = γ · v ( r → ) · m 1 , 其中, m 1 = ∫ G ( t ) · tdt 为所述散相梯度的一阶梯度矩,γ为所述单个运动的质子的旋磁比,为所述单个运动质子的流动速度。
所述心电门控同步心电信号的方式主要包括前瞻式(Prospective)和/或回顾式(Retrospective)。
根据本发明的一个实施例,在所述通过心电门控同步心电信号之前,还包括:通过压脂序列抑制所述目标区域的脂肪信号。
根据本发明的一个实施例,通过快速图像采集方法对所述其中所述血流信号被抑制的目标区域进行图像采集,其中,所述快速图像采集方法主要包括平面回波成像、快速自旋回波成像、螺旋桨式成像、螺旋式成像、变密度螺旋式成像、自旋回波成像、稳态自由进动成像、梯度回波成像、平行成像。
本发明实施例的黑血电影成像方法,是一种基于速度选择性准备(VElocity SelectivePreparation,VESP)黑血技术和快速图像采集技术的黑血电影成像方法,具有以下优点:
(1)通过在磁共振的磁化准备阶段抑制目标区域内的血流信号,消除了血流所诱发的流动伪影,从而能够清晰确定血管壁和血管腔之间的边界,实现对血管壁边界的精确描绘,最终实现对血管动态特性和功能特性的更加准确的评估;
(2)通过在磁共振的磁化准备阶段抑制目标区域内的血流信号,可以同时区分出血管壁内外管壁,进而计算血管壁面积或者血管壁厚度,并分析其在心动周期内的变化;
(3)通过使用心电门控同步心电信号,避免了在一个心动周期内连续采集一幅图像,相对于静态成像方式而言,可以减少血管搏动引起的图像模糊;
(4)使用VESP黑血技术,并可以与多种快速成像技术结合,成像速度快,信噪比高,具有更加广泛的适用性。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明一个实施例的基于磁共振的黑血电影成像方法的流程图;
图2为根据本本发明一个实施例的VESP序列的示意图;
图3为根据本发明一个实施例的通过速度选择准备VESP序列抑制目标区域内的血流信号的方法流程图;
图4为根据本发明一个实施例的通过前瞻式和回顾式两种同步心电信号的示意图;
图5a和图5b为根据本发明一个实施例的通过回顾式同步心电信号过程中采集图像信号的示意图。
具体实施方式
下面参考附图描述根据本发明实施例的基于磁共振的黑血电影成像方方法,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
为了研究动脉管壁的动力学系统及其与相关心血管疾病的关系,本发明中,发明人首次将黑血成像技术应用于电影成像中,即首次提出了基于磁共振的黑血电影成像技术,相对于亮血电影成像技术而言,黑血电影成像技术在研究动脉壁管动力学特性方面更加有优势。黑血电影成像技术能够抑制血流信号,从而消除血流引起的流动伪影影响,进而实现对血管壁边界的精确描绘,最终实现对血管壁功能参数的准确评估。同时黑血电影技术可以区分血管壁内外边界,从而能够提供详细的血管壁信息。
基于上述研究,本发明的实施例提出了一种基于磁共振的黑血电影成像方法,图1为根据本发明一个实施例的基于磁共振的黑血电影成像方法的流程图。
如图1所示,根据本发明实施例的基于磁共振的黑血电影成像方法,包括步骤S101:在磁共振的磁化准备阶段通过速度选择准备VESP序列抑制目标区域内的血流信号。
下面将对步骤S101进行详细说明。在步骤S101中,目标区域为黑血电影图像的采集区域,可根据不同研究和应用,将目标区域设置为不同的位置,如可为颈总动脉、腹主动脉、心脏等。本发明的实施例对目标区域的位置不做限定。
在本发明的一个实施例中,VESP序列包括T2准备脉冲序列、散相梯度以及在T2准备脉冲序列后的破坏梯度,其中,T2准备脉冲序列依次包括一个90°x脉冲,一个或多个180°y脉冲,以及一个90°-x脉冲。其中,90°x脉冲是沿预设坐标系的x轴方向入射的90°脉冲,180°y脉冲是沿预设坐标系的y轴方向入射的180°脉冲,90°-x脉冲是沿预设坐标系的-x轴方向入射的90°脉冲。应当理解,本发明的实施例对脉冲的形状和实现方式不做限定,如矩形脉冲、SINC脉冲或者绝热(adiabatic)脉冲等。
在本发明的实施例中,T2准备脉冲序列中的脉冲也可依次分别为上述脉冲的等效脉冲,例如,如果90°x脉冲为顺时针方向上的,其等效脉冲为逆时针方向上的270°x脉冲。本发明实施例中使用的脉冲仅为示例性的,并非仅能通过上述示例来实现,本领域技术人员可根据上述示例进行相应的扩展或变化,这些均应包含在本发明的保护范围之内。
在本发明的实施例中,散相梯度在K空间的读出方向、相位编码方向和选层方向中的至少一个方向施加,在本发明的一个优选的实施例中,可在读出方向、相位编码方向和选层方向三个方向均施加散相梯度。
具体地,图2为根据本本发明一个实施例的VESP序列的示意图。如图2所示,该VESP序列包括由一个90°x脉冲,两个180°y脉冲,以及一个90°-x脉冲组成的T2准备脉冲序列、从读出方向(Gx)、相位编码方向(Gy)和选层方向(Gz)三个方向施加的散相梯度以及破坏梯度S。
图3为根据本发明一个实施例的通过速度选择准备VESP序列抑制目标区域内的血流信号的方法流程图。如图3本发明实施例的通过速度选择准备VESP序列抑制目标区域内的血流信号包括以下步骤:
S201,将90°x脉冲作用于目标区域,以使目标区域的血流中运动的质子和目标区域的血管中静止的质子翻转至横截面。
在本发明的实施例中,当90°x脉冲作用于目标区域时,目标区域的血流中运动的质子和目标区域的血管中静止的质子受到90°x脉冲作用被翻转至横截面,即上述预设坐标系的xy平面。
S202,将一个或多个180°y脉冲和散相梯度共同作用于目标区域,以将翻转后的运动的质子散相。
在本发明的实施例中,散相梯度为满足以下条件的任意流速敏感梯度:
m 0 = ∫ G ( t ) dt = 0 ,
m 1 = ∫ G ( t ) · tdt ≠ 0 ,
其中,m0,m1分别为散相梯度的零阶梯度矩和一阶梯度矩。散相梯度的零阶梯度矩为0,则该散相梯度不会引起静止质子的相位变化,而一阶梯度矩不为0,则该散相梯度可以使运动质子散相。
具体地,通过一个180°y脉冲使翻转后的运动的和静止的质子重聚,消除B0场不均匀可能带来的影响;通过将散相梯度作用于目标区域内均匀分布的运动的质子,以使目标区域内的运动的质子发生相位变化,其中,目标区域内单个质子的相位变化为:
φ = γ · v ( r → ) · ∫ G ( t ) · tdt = γ · v ( r → ) · m 1 m 0 = ∫ G ( t ) dt = 0 ,
其中,为散相梯度的一阶梯度矩,γ为单个运动的质子的旋磁比,为单个运动质子的流动速度。
从上述单个质子的相位变化公式可以看出不同的流动速度会导致不同的相位变化,由于目标区域内的运动的质子具有不同的流动速度,因此,当散相梯度作用于目标区域时,目标区域内的运动的质子分别产生不同的相位变化,即发生散相,此时目标区域内单个像素内的运动质子的磁矩相互抵消,相应地,该像素内的血流信号被抑制。
S203,将90°-x脉冲作用于翻转后的静止的质子和散相后剩余的运动的质子,以将翻转后的静止的质子和散相后剩余的运动的质子翻转至垂直方向。
在本发明的实施例中,当90°-x脉冲作用于目标区域时,目标区域内的翻转后的静止的质子和散相后剩余的运动的质子受到90°-x脉冲作用被翻转至垂直方向,即上述预设坐标系的z方向。
S204,将破坏梯度作用于目标区域以清除横截面残余的磁矩。
在本发明的实施例中,在梯度磁场作用于目标区域后,由于T1弛豫的原因,在垂直方向会有一定量的磁矩恢复。因此,施加完S203的90°-x脉冲后,横截面静止的质子和散相后剩余的运动的质子翻转至垂直方向,同时该部分的磁矩会相应地被翻转到横截面,为了减少对后面信号采集可能造成的影响,可通过将破坏梯度作用于目标区域,清除该部分磁矩。
应当理解,在本发明的实施例中,为了加强血流信号的抑制效果,可在磁共振的磁化准备阶段通过多个VESP序列多次对血流信号进行抑制。
S102,通过心电门控同步心电信号。
在本发明的实施例中,心电门控同步心电信号的方式主要包括前瞻式和/或回顾式,也可为其他可同步心电信号的方式,本发明对此不做限定。
在本发明的实施例中,前瞻式同步心电信号过程中,在检测到R波后经过一个触发延时之后,开始采集图像信号,并在采集过程中将采集到的图像信号分配到相应的图像帧;回顾式同步心电信号过程中,在整个同步心电信号的过程中一直采集图像信号,并在每个心动周期结束时将在本轮心动周期内采集到的图像信号分别分配到各个图像帧。具体地,通过前瞻式和回顾式两种同步心电信号的方式可如图4所示,如图4所示,前瞻式同步心电信号过程中在每次检测到R波之后都需要经过触发延时后才进入图像采集,并在图像采集过程中经采集到的图像信号分配到1,2,3,4,5,6,……;而回顾式同步心电信号的过程中一直持续采集图像信号,并将在每两个R波之间的图像信号分别分配到特定的图像帧(1,2,……,13)。
S103,基于心电信号,对其中血流信号被抑制的目标区域进行图像采集并根据所采集到的第一图像生成目标区域的黑血电影图像。
具体地,在本发明的一个实施例中,当目标区域的血流信号被抑制后,此时对目标区域进行图像采集即可根据采集到的图像生成目标区域的黑血电影图像。
在本发明的实施例中,主要通过快速采集方法对其中血流信号被抑制的目标区域进行图像采集,从而成像速度快,信噪比高,具有更加广泛的适用性。其中,快速图像采集方法主要包括平面回波成像(EPI)、快速自旋回波成像(FSE)、螺旋桨式成像(PROPELLER)、螺旋式成像(Spiral)、变密度螺旋式成像(VDS)、自旋回波成像(SE)、稳态自由进动成像(bSSFP)、梯度回波成像(SPGR)、平行成像等。本发明的实施例对快速图像采集方法的种类不限。
本发明实施例的基于磁共振的黑血电影成像方法,通过在磁共振的磁化准备阶段抑制目标区域内的血流信号,并基于心电同步信号对该目标区域进行图像采集生成黑血电影图像。在磁共振的磁化准备阶段抑制了血液信号,消除了血液流动所诱发的流动伪影,从而能够清晰确定血管壁和血管腔之间的边界,实现对血管壁边界的精确描绘,最终实现对血管动态特性和功能特性的准确评估,并且可以同时区分出血管壁内外管壁,进而计算血管壁面积或者血管壁厚度,并分析其在心动周期内的变化。此外,通过使用心电门控同步心电信号,避免了在一个心动周期内连续采集一幅图像,相对于静态成像方式而言,可以减少血管搏动引起的图像模糊。
由于本发明实施例的方法可以监测动脉动力学特性的改变,因此本发明实施例的方法可以用于一些心血管相关疾病的预防、诊断和治疗,如高血压、动脉粥样硬化、动脉瘤以及其他大动脉壁畸形等;同时也可以用于心脏成像,比如,心功能评价等。
在本发明的一个实施例中,本发明实施例的基于磁共振的黑血电影成像方法在通过心电门控同步心电信号之前还包括:
通过压脂序列抑制目标区域的脂肪信号。从而在去除了血流造成的流动伪影之外,还去除了脂肪信号的干扰,从而使黑血电影图像中血管壁分界更加清晰,有利于对血管壁动态特性和功能特性的准确评估。
具体地,图5a和图5b为根据本发明一个实施例的通过回顾式同步心电信号过程中采集图像信号的示意图。如图5a和图5b所示,在同步心电信号过程中在两个R波之间按照一定的频率多次采集图像信号,并在每次采集图像信号(图5a和图5b中用白色块表示)之前均有抑制血液信号的过程(图5a和图5b中灰色块VESP表示)和抑制脂肪信号的过程(图5a和图5b中黑色块FS表示)。其中,在图5a中,在抑制血液信号和抑制脂肪信号后,通过梯度回波成像(SPGR)方法进行快速图像采集;在图5b中,在抑制血液信号和抑制脂肪信号后,通过稳态自由进动成像(bSSFP)方法进行快速图像采集。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为,表示包括一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分,并且本发明的优选实施方式的范围包括另外的实现,其中可以不按所示出或讨论的顺序,包括根据所涉及的功能按基本同时的方式或按相反的顺序,来执行功能,这应被本发明的实施例所属技术领域的技术人员所理解。
在流程图中表示或在此以其他方式描述的逻辑和/或步骤,例如,可以被认为是用于实现逻辑功能的可执行指令的定序列表,可以具体实现在任何计算机可读介质中,以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用,或结合这些指令执行系统、装置或设备而使用。就本说明书而言,"计算机可读介质"可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下:具有一个或多个布线的电连接部(电子装置),便携式计算机盘盒(磁装置),随机存取存储器(RAM),只读存储器(ROM),可擦除可编辑只读存储器(EPROM或闪速存储器),光纤装置,以及便携式光盘只读存储器(CDROM)。另外,计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质,因为可以例如通过对纸或其他介质进行光学扫描,接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序,然后将其存储在计算机存储器中。
应当理解,本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中,多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如,如果用硬件来实现,和在另一实施方式中一样,可用本领域公知的下列技术中的任一项或他们的组合来实现:具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路,具有合适的组合逻辑门电路的专用集成电路,可编程门阵列(PGA),现场可编程门阵列(FPGA)等。
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,该程序在执行时,包括方法实施例的步骤之一或其组合。
此外,在本发明各个实施例中的各功能单元可以集成在一个处理模块中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。
上述提到的存储介质可以是只读存储器,磁盘或光盘等。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同限定。

Claims (9)

1.一种基于磁共振的黑血电影成像方法,其特征在于,包括:
在所述磁共振的磁化准备阶段通过速度选择准备VESP序列抑制目标区域内的血流信号;
通过心电门控同步心电信号;以及
基于所述心电信号,对其中所述血流信号被抑制的目标区域进行图像采集并根据所采集到的第一图像生成所述目标区域的黑血电影图像。
2.如权利要求1所述的方法,其特征在于,所述VESP序列包括T2准备脉冲序列、散相梯度以及在所述T2准备脉冲序列后的破坏梯度,其中,
所述T2准备脉冲序列依次包括一个90°x脉冲,一个或多个180°y脉冲,以及一个90°-x脉冲;
所述散相梯度在读出方向、相位编码方向和选层方向中的至少一个方向施加。
3.如权利要求2所述的方法,其特征在于,所述通过速度选择准备VESP序列抑制目标区域内的血流信号,包括:
将所述90°x脉冲作用于所述目标区域,以使所述目标区域的血流中运动的质子和所述目标区域的血管中静止的质子翻转至横截面;
将所述一个或多个180°y脉冲和所述散相梯度共同作用于所述目标区域,以将翻转后的运动的质子散相;
将所述90°-x脉冲作用于翻转后的静止的质子和散相后剩余的运动的质子,以将所述翻转后的静止的质子和所述散相后剩余的运动的质子翻转至垂直方向;
将所述破坏梯度作用于所述目标区域以清除所述横截面残余的磁矩。
4.如权利要求2所述的方法,其特征在于,所述散相梯度为满足以下条件的任意流速敏感梯度:
m 0 = ∫ G ( t ) dt = 0 ,
以使所述散相梯度不引起所述静止质子的相位变化;以及
m 1 = ∫ G ( t ) · tdt ≠ 0 ,
以使所述散相梯度引起所述运动质子散相;其中,
m0,m1分别为所述散相梯度的零阶梯度矩和一阶梯度矩。
5.如权利要求3所述的方法,其特征在于,所述目标区域内的运动的质子分别具有不同的流动速度,在所述将所述一个或多个180°y脉冲和所述散相梯度共同作用于所述目标区域的过程中,
所述散相梯度作用于所述具有不同流动速度的运动的质子时,所述目标区域内单个像素内具有不同流动速度的运动的质子分别产生不同的相位变化,发生散相;
当所述单个像素内的运动的质子发生散相时,所述单个像素内的运动的质子的磁矩相互抵消,所述单个像素对应的血流信号被抑制。
6.如权利要求5所述的方法,其特征在于,其中,所述目标区域内单个运动质子的相位变化为:
φ = γ · v ( r → ) · ∫ G ( t ) · tdt = γ · v ( r → ) · m 1 m 0 = ∫ G ( t ) dt = 0 ,
其中,为所述散相梯度的一阶梯度矩,γ为所述单个运动的质子的旋磁比,为所述单个运动质子的流动速度。
7.如权利要求1所述的方法,其特征在于,所述心电门控同步心电信号的方式主要包括前瞻式和/或回顾式。
8.如权利要求1所述的方法,其特征在于,在所述通过心电门控同步心电信号之前,还包括:
通过压脂序列抑制所述目标区域的脂肪信号。
9.如权利要求1所述的方法,其特征在于,通过快速图像采集方法对所述其中所述血流信号被抑制的目标区域进行图像采集,其中,
所述快速图像采集方法主要包括平面回波成像、快速自旋回波成像、螺旋桨式成像、螺旋式成像、变密度螺旋式成像、自旋回波成像、稳态自由进动成像、梯度回波成像、平行成像。
CN201410184612.0A 2014-05-04 2014-05-04 基于磁共振的黑血电影成像方法 Expired - Fee Related CN103976735B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410184612.0A CN103976735B (zh) 2014-05-04 2014-05-04 基于磁共振的黑血电影成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410184612.0A CN103976735B (zh) 2014-05-04 2014-05-04 基于磁共振的黑血电影成像方法

Publications (2)

Publication Number Publication Date
CN103976735A true CN103976735A (zh) 2014-08-13
CN103976735B CN103976735B (zh) 2016-03-23

Family

ID=51269087

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410184612.0A Expired - Fee Related CN103976735B (zh) 2014-05-04 2014-05-04 基于磁共振的黑血电影成像方法

Country Status (1)

Country Link
CN (1) CN103976735B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104305959A (zh) * 2014-10-13 2015-01-28 清华大学 基于速度选择射频脉冲的磁共振黑血成像方法及系统
CN105596003A (zh) * 2015-12-24 2016-05-25 深圳先进技术研究院 一种下肢深静脉血栓磁共振成像方法和装置
WO2019169671A1 (zh) * 2018-03-05 2019-09-12 奥泰医疗系统有限责任公司 一种脂肪压制黑血磁共振成像方法
CN111202519A (zh) * 2020-01-17 2020-05-29 首都医科大学宣武医院 一种在体血栓软硬度检测的方法及其系统
CN112014782A (zh) * 2020-08-06 2020-12-01 深圳先进技术研究院 磁共振成像方法、装置及计算机存储介质
CN112986878A (zh) * 2019-12-18 2021-06-18 上海联影医疗科技股份有限公司 冠脉成像方法以及磁共振成像系统
US11918335B2 (en) 2020-08-06 2024-03-05 Shenzhen Institutes Of Advanced Technology Magnetic resonance imaging method, apparatus, and computer storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1543326A (zh) * 2001-11-15 2004-11-03 ��ʽ���綫֥ 在预脉冲和成象脉冲列之间需要长等待时间的核磁共振机
US20090261825A1 (en) * 2008-04-18 2009-10-22 Case Western Reserve University Dark blood balanced steady state free precession imaging
CN102551721A (zh) * 2011-12-12 2012-07-11 中国科学院深圳先进技术研究院 血管成像方法和系统
CN103110420A (zh) * 2013-01-31 2013-05-22 深圳先进技术研究院 磁共振血管成像方法和系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1543326A (zh) * 2001-11-15 2004-11-03 ��ʽ���綫֥ 在预脉冲和成象脉冲列之间需要长等待时间的核磁共振机
US20090261825A1 (en) * 2008-04-18 2009-10-22 Case Western Reserve University Dark blood balanced steady state free precession imaging
CN102551721A (zh) * 2011-12-12 2012-07-11 中国科学院深圳先进技术研究院 血管成像方法和系统
CN103110420A (zh) * 2013-01-31 2013-05-22 深圳先进技术研究院 磁共振血管成像方法和系统

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104305959A (zh) * 2014-10-13 2015-01-28 清华大学 基于速度选择射频脉冲的磁共振黑血成像方法及系统
CN105596003A (zh) * 2015-12-24 2016-05-25 深圳先进技术研究院 一种下肢深静脉血栓磁共振成像方法和装置
CN105596003B (zh) * 2015-12-24 2019-08-02 深圳先进技术研究院 一种下肢深静脉血栓磁共振成像方法和装置
WO2019169671A1 (zh) * 2018-03-05 2019-09-12 奥泰医疗系统有限责任公司 一种脂肪压制黑血磁共振成像方法
CN112986878A (zh) * 2019-12-18 2021-06-18 上海联影医疗科技股份有限公司 冠脉成像方法以及磁共振成像系统
CN111202519A (zh) * 2020-01-17 2020-05-29 首都医科大学宣武医院 一种在体血栓软硬度检测的方法及其系统
CN112014782A (zh) * 2020-08-06 2020-12-01 深圳先进技术研究院 磁共振成像方法、装置及计算机存储介质
US11918335B2 (en) 2020-08-06 2024-03-05 Shenzhen Institutes Of Advanced Technology Magnetic resonance imaging method, apparatus, and computer storage medium

Also Published As

Publication number Publication date
CN103976735B (zh) 2016-03-23

Similar Documents

Publication Publication Date Title
CN103976735B (zh) 基于磁共振的黑血电影成像方法
US8577443B2 (en) Magnetic resonance imaging apparatus and image processing apparatus
US8577442B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
US8417007B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP6381879B2 (ja) 磁気共鳴イメージング装置
US8301224B2 (en) System and method for automatic, non-invasive diagnosis of pulmonary hypertension and measurement of mean pulmonary arterial pressure
JP6521910B2 (ja) 磁気共鳴イメージング装置
KR100646914B1 (ko) 자기공명 이미징장치
WO2014085288A1 (en) Noninvasive 4-d time-resolved dynamic magnetic resonance angiography
JP6109978B2 (ja) 磁気共鳴イメージング装置
US8581582B2 (en) MRI non-contrast time-slip angiography using variably positioned cine sub-sequence
JP2009160378A (ja) 磁気共鳴イメージング装置
JP6713988B2 (ja) 血流アーチファクトを抑制させたDixonMR撮像
JP5901905B2 (ja) 磁気共鳴イメージング装置
Stankovic et al. Kt GRAPPA-accelerated 4D flow MRI of liver hemodynamics: influence of different acceleration factors on qualitative and quantitative assessment of blood flow
JP2000201903A5 (zh)
JP6991728B2 (ja) 画像処理装置、磁気共鳴イメージング装置及び画像処理方法
JP6513413B2 (ja) 医用画像処理装置及び磁気共鳴イメージング装置
Corrias et al. Heart applications of 4D flow
JP2018503453A (ja) 磁化率重み付け磁気共鳴イメージングのための方法及びシステム
Richards et al. The role of 4-dimensional flow in the assessment of bicuspid aortic valve and its valvulo-aortopathies
Haraldsson et al. Improved estimation and visualization of two‐dimensional myocardial strain rate using MR velocity mapping
CN105686829B (zh) 在检查对象的周期运动的情况下的形变计算
KR100732790B1 (ko) 프리펄스와 촬상펄스열 사이에 긴 대기시간을 요하는자기공명 이미징
CN103083021B (zh) 一种最佳m1值的获取方法、装置及医疗设备

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: 20160323

Termination date: 20170504

CF01 Termination of patent right due to non-payment of annual fee