[发明专利]一种用于脑部微组织结构测量的扩散弛豫谱成像方法

说明书

本发明公开了一种用于脑部微组织结构测量的扩散弛豫谱成像方法，包括：获取目标对象在多回波时间和多b值下的扩散磁共振数据；建立多组织脑微结构的扩散弛豫谱成像模型；通过基于球面平均技术的多参数联合优化方法获得多组织脑微结构的扩散弛豫时间和扩散敏感系数谱；将得到的弛豫时间、脑组织体积分数以及分布函数进行表征得到脑组织微观结构参数。本发明可提供组织成分特定的非T2加权的神经突内分数和自由水分数，解决了传统成分分数的TE依赖性。

技术领域

本发明属于磁共振成像和医学图像处理技术领域，具体涉及一种用于脑部微组织结构测量的扩散弛豫谱成像方法。

背景技术

脑部微组织结构测量的目的是为了实现对神经元微结构的特异性，以产生无创成像的独特临床生物标记物。由于生物组织十分复杂，微观结构建模通常需要基于特定的假设对组织结构进行模拟，以得到其感兴趣的目标指标。对于微观参数特异性成像，实现此目的的一种有效方法是增加数据的测量维度，以更好的区分组织微成分，如利用扩散和T1、T2弛豫和波谱成像之间的协同作用。相关进展包括弛豫扩散相关成像技术，已被应用于动物和离体大脑中。最近，这种多维度成像方法与微观结构建模相结合，以改善对人脑成像的组织微结构参数的估计。但是，多维度成像在提供特异性参数的同时也进一步增加了数据采集以及模型拟合的困难。

另外，在现有的扩散磁共振数据采集方案下，微观结构参数对扩散信号变化不敏感，限制了模型参数的特异性。其中，常见的多成分模型只考虑了不同成分间的扩散特性差异，而未考虑弛豫时间等其他特性的差异，而扩散磁共振数据也经常以单个TE进行采集。这导致从模型得出的定量指标会受到来自T1和T2加权的干扰。当成份之间存在T1或T2加权差异时，模型不能将成份之间的弛豫和扩散特性区分开来。在组织中神经突内和神经突外成份之间，T2的差异相比来说要小得多，但也被证明会导致在白质中一些扩散模型的导出参数具有回波时间依赖性这将减少模型衍生的定量指标的特异性，及研究结果的可解释性，特别是对于研究伴随着T2弛豫时间大幅改变的大脑发育和成熟的过程。

综上，脑组织微结构的特征在于扩散率和横向弛豫率之间的异质性。扩散磁共振成像及其高阶模型能够无创探测脑组织微结构。其中，受限、受阻和自由成分的水分子之间的信号差异可以通过高阶扩散模型来表征，从而估计体素内的细胞体、轴突纤维和自由成分的相对比例，从而用于检测大脑的发育、退化和疾病驱动的脑组织微结构改变。然而，组织成分的准确表征不仅受组织依赖的扩散系数的影响，也受到横向弛豫率的影响。明确考虑弛豫-扩散的耦合关系可能会显著改善脑组织微结构的表征。具有单一回波时间的标准扩散磁共振成像技术主要提供有关扩散率的信息，而弛豫扩散磁共振成像涉及多个回波时间和多个扩散加权强度，能够探测弛豫和扩散率之间的组织特征的耦合关系，因此提出一种用于脑部微组织结构测量的扩散弛豫谱成像方法具有重要意义。

发明内容

为了突破现有技术的挑战，本发明提供一种扩散弛豫谱成像方法，探测弛豫-扩散的耦合关系，并定量测量脑组织微结构参数，解决了传统成分分数的TE依赖性。

实现本发明目的的技术解决方案为：

一种用于脑部微组织结构测量的扩散弛豫谱成像方法，包括以下步骤：

步骤1、通过磁共振设备对脑部进行扩散加权成像扫描，根据不同b值和不同回波时间的设定，进行采样，得到多b值多回波时间的扩散磁共振成像图像；

步骤2、基于球面平均扩散弛豫谱成像模型对多b值多回波时间的扩散磁共振成像图像进行拟合得到具有脑组织特异性和b值依赖的弛豫时间，以及脑组织体积分数；

步骤3、基于扩散弛豫谱成像模型和步骤2中得到的弛豫时间对多b值多回波时间的扩散磁共振图像进行拟合得到无T2加权的具有脑组织特异性的纤维方向分布函数；

步骤4、将步骤2得到的弛豫时间、脑组织体积分数以及步骤3得到的分布函数进行表征得到脑组织微观结构参数。