一种基于重叠回波的单扫描定量磁共振T2成像方法，涉及磁共振成像的方法。通过在单次扫描中加入两个有相同偏转角的激发脉冲来产生两个不同演化时间的回波信号，尽管两个回波信号的演化时间不同，导致两个回波信号的T2加权不同，这两个回波信号来自同一个成像切片，可通过两个回波信号之间的先验知识：两者结构类似、联合边缘的稀疏性来分离这两个回波信号，利用稀疏变换配合相应的分离算法对这两个回波信号进行分离；最后对分离得到的两个信号进行T2计算得到定量T2图像。获得了单次扫描的定量T2成像，将定量T2成像的时间由秒级甚至分钟级，减少到ms级，并且得到的T2图像质量能够与常规的单次扫描EPI序列得到的图像质量相当。
【技术实现步骤摘要】
本专利技术涉及磁共振成像的方法，尤其是涉及一种基于重叠回波的单扫描定量磁共 振Τ2成像方法。
技术介绍
磁共振参数成像（Τ2成像、12$成像和扩散成像）因其能够提供丰富 的组织特征属性的定量信息而在临床诊断上有着广泛的应用（ B. Zhao, F. Lam,and Z.P.Liang, "Model-based MR Parameter mapping with sparsity constraints:parameter estimation and performance bounds,'，IEEE Trans. Med. Imag ?，vol. 33, no. 9, pp. 1832-1844, 2014)，比如：心肌梗死的诊断、测量肝脏中铁含量是否过 量等。特别地，T2弛豫时间的定量分析在精神病学以及神经病科等临床医学磁共振成像 中引起了越来越多的关注。然而，磁共振参数成像在其成像过程中往往需要获取一系列对 比加权的图像，因而一般其获取数据的时间很长。虽然现在有很多不同的成像方法的提出 来克服上述问题，比如：下采样自旋回波磁共振成像（Spin-Echo MRI)、梯度自旋回波磁共 振成像（Gradient Spin Echo MRI)、基于压缩感知的平行成像（Parallel Imaging with CS)等。同时，一些基于模型的重建方法和基于布洛赫仿真的重建方法相继被提出，以此 来进一步加快成像的速度。但是多次激发的磁共振参数成像方法在获取阶段仍然需要耗 费数秒的时间，这样对不可重复的神经活动进行实时参数成像就变成了几乎不可能实现的 任务。因此，单扫描的多回波的平面回波成像（echo-planar imaging，EPI)的成像方法 被提出（ S. Posse, S. Wiese, D. Gembris, K. Mathiak, C. Kessler, M. L. Grosse-Ruyken，B. Elghawaghi, T. Richards, S. R. Dager, and V. G. Kiselev, ^Enhancement of BOLD-Contrast Sensitivity by Single-Shot Multi-Echo Functional MR Imaging,'，Magn. Reson. Med.，vol. 42, pp. 87 - 97, 1999)，此方法通过将一系列对比加权图像的获取包含在一次扫 描中所获得的多个回波中。然而这种方法存在局限性，一方面是这种方法需要延长回波链， 必然导致增加获取的时间与信号的衰减；另一方面是这种方法的实现与常规EPI方法相比 是以延长重复时间（TR)为代价的，这就可能需要牺牲所得回波图像的空间分辨率；而最重 要的是这种方法目前只能用于T2*定量成像，尚无法用于T2定量成像。此外，尽管有不同 的快速T2成像方法相继被提出，包括梯度自旋回波序列，但是这些方法都是用多次激发序 列来进行T2成像，这样不仅效果差强人意，成像效率也亟待提高。
技术实现思路
本专利技术的目的在于提供一种基于重叠回波的单扫描定量磁共振T2成像方法。 本专利技术包括如下步骤： (1)在磁共振成像仪操作台上，打开磁共振成像仪中的操作软件，首先对成像物体 进行感兴趣区域定位，然后进行调谐、匀场、功率校正和频率校正； (2)导入事先编译好的0LED成像序列；根据具体的实验情况，设置脉冲序列的各 个参数； 所述0LED成像序列的结构依次为：翻转角为α的片选脉冲、，翻转角 为α的片选脉冲、ΤΕ/2、180°重聚脉冲、采样回波链； 将两个小角度激发脉冲结合两个回波移位梯度匕和G 2，使得两个回波在Κ空间的 中心产生偏移，所述180°重聚脉冲以及两个小角度激发脉冲都与层选梯度Gss相结合进行 层选；第二个小角度激发脉冲前后分别施加回波延时和，所述180°重聚 脉冲前后有X，y，z三个方向的破坏梯度作用； 所述采样回波链是由分别作用在X，y方向的梯度链组成；X方向的梯度链由一系 列正负梯度构成，且每个梯度的面积是第一个回波移位梯度匕的三倍；y方向的梯度链是 由一系列大小相等的"blips"梯度构成，且所述"blips"梯度的总面积等于移位梯度面积 的四倍； 在所述采样回波链之前，X和y方向分别施加了重聚梯度GroJP G "，所述的面 积是X方向第一个梯度面积的一半，方向与X方向第一个梯度方向相反；所述G&的面积是 所有所述"blips"梯度的总面积的一半，方向与所述"blips"梯度方向相反； (3)执行步骤（2)设置好的所述0LED成像序列，进行数据采样；数据采样完成后 得到两个回波信号的K空间数据； (4)对步骤（3)得到的K空间数据进行分析以及对回波信号磁化矢量M+的演化进 行理论推导，在第二个回波移位梯度G2之后，可以得到以下式子： 式中是自旋密度分布，α是激发脉冲的翻转角，通过实验发现，当α = 45° 时，两个回波信号的强度都相对较高，，其 中h分别对应的第一个和第二个回波移位梯度的持续时间，γ是磁旋比；由上式可 知，实际上有三个被不同相位调制的回波信号，其中第一项是由第二个激发脉冲产生的，而 后两项是由第一个激发脉冲产生的；然而要分离出这三个信号相对单扫描获取的信号来说 是非常复杂的，通过分析后两项，可知后两项的回波中心位置是不一样的，且最后一项信号 强度相对第二项来说比较小，因此最后一项的回波可以通过简单的处理而被忽略； (5)对步骤（4)得到的回波信号用下述分离算法进行分离，根据傅里叶变换理论， 两个回波信号在图像域的线性相位是不一样的，此外，虽然两个回波信号由于演化时间不 同导致Τ2加权不同，但是它们是来自同一个图像层；因此，利用两者的图像结构相似的先 验信息可以对两个回波信号进行联合重建，重建算法如下： 其中，Xl，χ2分别是从第一个和第二个回波信号中重建出来的图像；是尺度因子，x1(]，x2(^ v别是第一个和第二个回波信号的初始图像ρ λ2和λ3分别是拉格朗日乘数法可调整约束权重；▽是梯度算子；第一项是保真项，第二项和 第三项是对第一幅和第二幅图像的稀疏性约束，最后一项是两幅图像轮廓相似性约束；这 两幅图像有如下关系： 其中：分别是第一幅和第二幅图像的线性相位位移；X。是原始信 号，是由包含第一个和第二个回波信号的原始信号进行傅里叶逆变换得到的，通过迭代算 法求解上述式子就可以得到分离后的第一个和第二个回波信号产生的图像； (6)步骤（5)分离出来的图像进行T2成像计算；对于单扫描的T2成像方法来说， 只需要两幅不同的回波时间图像就可行了，T2的值直接通过T2弛豫方程求得： 其中是校正因子，；S^S2分别是第一个回波信号和 第二个回波信号的图像强度，加入全变分（Total Variation)外推法来增强图像的分辨率， 而且设定一个阈值，当得到的数值低于阈值时，会被认为是噪声而被忽略，同样当计算出的 T2值过大时也是不合理的，也会被省略；最后通过T2成像计算得到了具有较好分辨率的高 品质的T2图像。 本专利技术提供一种在一次扫描的情况下，获得重叠的回波信号，然后利用分离算法 对重叠的信号进行分离，最后进行T2计算，就能获得与本文档来自技高网...
【技术保护点】
一种基于重叠回波的单扫描定量磁共振T2成像方法，其特征在于包括如下步骤：(1)在磁共振成像仪操作台上，打开磁共振成像仪中的操作软件，首先对成像物体进行感兴趣区域定位，然后进行调谐、匀场、功率校正和频率校正；(2)导入事先编译好的OLED成像序列；根据具体的实验情况，设置脉冲序列的各个参数；所述OLED成像序列的结构依次为：翻转角为α的片选脉冲、(TE2‑TE1)/2、翻转角为α的片选脉冲、TE1/2、180°重聚脉冲、采样回波链；将两个小角度激发脉冲结合两个回波移位梯度G1和G2，使得两个回波在K空间的中心产生偏移，所述180°重聚脉冲以及两个小角度激发脉冲都与层选梯度Gss相结合进行层选；第二个小角度激发脉冲前后分别施加回波延时(TE2‑TE1)/2和TE1/2，所述180°重聚脉冲前后有x，y，z三个方向的破坏梯度作用；所述采样回波链是由分别作用在x，y方向的梯度链组成；x方向的梯度链由一系列正负梯度构成，且每个梯度的面积是第一个回波移位梯度G1的三倍；y方向的梯度链是由一系列大小相等的“blips”梯度构成，且所述“blips”梯度的总面积等于移位梯度面积的四倍；在所述采样回波链之前，x和y方向分别施加了重聚梯度Gror和Gar，所述Gror的面积是x方向第一个梯度面积的一半，方向与x方向第一个梯度方向相反；所述Gar的面积是所有所述“blips”梯度的总面积的一半，方向与所述“blips”梯度方向相反；(3)执行步骤(2)设置好的所述OLED成像序列，进行数据采样；数据采样完成后得到两个回波信号的K空间数据；(4)对步骤(3)得到的K空间数据进行分析以及对回波信号磁化矢量M+的演化进行理论推导，在第二个回波移位梯度G2之后，可以得到以下式子：M+=∫r→ρ(r→)sinα{-i cosαeiθ2+12e-δTE/T2(r→)[(1+cosα)ei(θ2-θ1)+(1-cosα)ei(θ1+θ2)]}dr→]]>式中是自旋密度分布，α是激发脉冲的翻转角，通过实验发现，当α＝45°时，两个回波信号的强度都相对较高，δTE＝(TE2‑TE1)/2，其中δ1，δ2分别对应的第一个和第二个回波移位梯度的持续时间，γ是磁旋比；由上式可知，实际上有三个被不同相位调制的回波信号，其中第一项是由第二个激发脉冲产生的，而后两项是由第一个激发脉冲产生的；然而要分离出这三个信号相对单扫描获取的信号来说是非常复杂的，通过分析后两项，可知后两项的回波中心位置是不一样的，且最后一项信号强度相对第二项来说比较小，因此最后一项的回波可以通过简单的处理而被忽略；(5)对步骤(4)得到的回波信号用下述分离算法进行分离，根据傅里叶变换理论，两个回波信号在图像域的线性相位是不一样的，此外，虽然两个回波信号由于演化时间不同导致T2加权不同，但是它们是来自同一个图像层；因此，利用两者的图像结构相似的先验信息可以对两个回波信号进行联合重建，重建算法如下：{x1,x2}=argminx1,x2[||x1-x10||22+λ1||▿x1||1+λ2||▿x2||1+λ3||▿(x1-βx2)||1]]]>其中，x1，x2分别是从第一个和第二个回波信号中重建出来的图像；是尺度因子，x10，x20分别是第一个和第二个回波信号的初始图像；λ1，λ2和λ3分别是拉格朗日乘数法可调整约束权重；是梯度算子；第一项是保真项，第二项和第三项是对第一幅和第二幅图像的稀疏性约束，最后一项是两幅图像轮廓相似性约束；这两幅图像有如下关系：其中，分别是第一幅和第二幅图像的线性相位位移；x0是原始信号，是由包含第一个和第二个回波信号的原始信号进行傅里叶逆变换得到的，通过迭代算法求解上述式子就可以得到分离后的第一个和第二个回波信号产生的图像；(6)步骤(5)分离出来的图像进行T2成像计算；对于单扫描的T2成像方法来说，只需要两幅不同的回波时间图像就可行了，T2的值直接通过T2弛豫方程求得：T2(r)=-ΔTEln(μS2(r→)S1(r→)),]]>其中是校正因子，ΔTE＝TE2‑TE1；S1和S2分别是第一个回波信号和第二个回波信号的图像强度，加入全变分(Total Variation)外推法来增强图像的分辨率，而且设定一个阈值，当得到的数值低于阈值时，会被认为是噪声而被忽略，同样当计算...
【技术特征摘要】
【专利技术属性】
技术研发人员：蔡聪波，曾益清，陈人楷，丁兴号，蔡淑惠，陈忠，
申请(专利权)人：厦门大学，
类型：发明
国别省市：福建;35
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。