一种基于扩散张量成像的脑部特征提取方法
【技术领域】
[0001] 本发明涉及计算机图形学下的医学成像、神经解剖学领域，尤其一种基于扩散张 量成像的脑部特征提取方法。
【背景技术】
[0002] 随着时代的发展，医学影像技术的进步，扩散张量成像技术在神经科学的研究中 占越来越大的影响力，拥有先进的神经影像技术是这个时代不可或缺的；扩散张量成像作 为新兴的一种描述大脑结构的方法，同时也是唯一一种活体检测人脑结构的方法，在神经 医学领域主要是对脑组织结构特征的研究；目前，扩散张量成像技术正被广泛的应用于精 神科疾病和诊断的辅助手段，甚至可以用于术前手术方案的制定，可以说其在医学领域的 贡献有着无可替代的优势;所以对基于扩散张量的算法研究对于脑科学有着重大的意义。
[0003] 在扩散张量成像中的每个体素都对应于一个二阶张量，所以张量数据的表示和有 用信息的提取是扩散张量成像技术的关键因素;到目前为止，研究者提出了越来越多中不 同的关于扩散张量成像技术下的脑部特征提取算法，而对扩散张量数据的处理，以及脑部 特征结构的可视化表达，是扩散张量数据研究的热点。
【发明内容】
[0004] 为了实现更好的脑部特征结构的可视化，本发明结合扩散张量成像技术，提出了 一种更清晰、更有效、更能体现脑部特征结构的基于扩散张量成像的脑部特征提取方法。
[0005] 为了解决上述技术问题提供的技术方案为：
[0006] -种基于扩散张量成像的脑部特征提取方法，所述方法包括以下步骤：
[0007] 1 )，采用基于三维谐振子的重构与估计模型：
[0008] 该模型采用多壳q空间中水分子扩散方向数据的映射的方法，简称SHORE模型；首 先q空间磁共振衰减信号在三维量子力学谐振子问题中扩展为正交基，也就是通过三维扩 展方案来表示一维的q空间数据;而衰减信号的傅里叶反变换得到能够提供连接信息的平 均传播算子；
[0009] 2)，计算两个基于SHORE模型的标量映射：
[0010]计算基于SHORE模型下的均方位移和零位移概率，均方位移简称MSD，指的是表示 在t时刻水分子距离初始位置的距离，即水分子随机运动在空间范围内的度量，零位移概率 简称rtop，所述标量映射是通过水分子在一定时间内最低限度扩散的概率密度；
[0011] 3)，数据的可视化显示：
[0012] 通过调用Matplotlib二维图像数据显示工具来呈现SHORE模型通过映射后的图 像;所述Matplotlib是一个二维绘图库并且能够应用于python脚本。
[0013] 进一步，所述步骤3)中，调用Matplotlib中对应功能的类。
[0014] 再进一步，所述步骤1)中，建立基于SHORE模型包括以下步骤：
[0015]首先读取数据，建立模型：
[0016] E(q)=S(q)/So (1)
[0017] E (q)为q空间磁共振衰减信号，S (q)为扩散信号强度，而So未加扩散梯度的脉冲，q 为时一个实验控制参数：
[0018] q=(2ji)-1T5G (2)
[0019] γ为磁旋比，δ为扩散梯度脉冲的持续时间，G代表扩散梯度向量;q空间磁共振衰 减信号正交基基础下在三维量子力学谐振子问题中扩展：
[0021 ] 其中，Cnlm为变换系数，Φη1ω(ιι，q)为基函数，η为〇 nlm(u，q)的径向阶，1球面谐波函 数的阶次，L为1的最大值，N为η的最大值：
[0023] u表示每个体素常数的估计数据，ζ是一个比例系数，Γ为Ga_a函数，
[0024] Ztif是广义拉盖尔多项式，Γ(?)是球面谐波函数，m为角度。
[0025] 所述步骤2)中，计算基于SHORE的标量映射，包括以下步骤：
[0026] 2.1计算系踪传播算子，衰减信号的傅里叶反变换能够得到系踪传播算子：
[0028] R代表净位移向量，是傅里叶共辄，P(R)为系踪传播算子；
[0029] 2.2计算基于SHORE模型下的均方位移：
[0030]首先计算水分子的平均扩散位移矢量，它与体素中平均扩散系数相关，并通过爱 因斯坦扩散方程计算：
[0032]其中妄是扩散位移矢量，Δ表示扩散时间，MD是平均扩散系数，均方位移可以通过 概率密度函数的加权平方位移求得：
[0034]其中概率密度函数的求法如下：
[0036]其中D为扩散张量；
[0037] 2.3计算SHORE模型下的零位移概率
[0038]零位移概率标量映射是通过水分子在一定时间内最低限度扩散的概率密度;零位 移概率是给定的概率密度函数中水分子零位移的概率，gp:
[0040] 其中P。零位移概率。
[0041]所述步骤3)中，将步骤2)中计算得到的数据放入Matplotlib库中进行绘图，最后 通过Matplotlib工具显示出来。
[0042]本发明的有益效果为:结合扩散张量成像技术，更清晰、更有效、更能体现脑部特 征结构。
【附图说明】
[0043]图1为本发明中的流程图；
【具体实施方式】
[0044]以下将结合附图对本发明做进一步详细说明：
[0045] 参照图1，一种基于扩散张量成像的脑部特征提取方法，所述方法包括以下步骤：
[0046] 1 )，采用基于三维谐振子的重构与估计模型：
[0047]该模型采用多壳q空间中水分子扩散方向数据的映射的方法，简称SHORE模型；首 先q空间磁共振衰减信号在三维量子力学谐振子问题中扩展为正交基，也就是通过三维扩 展方案来表示一维的q空间数据;而衰减信号的傅里叶反变换得到能够提供连接信息的平 均传播算子;过程如下：
[0048]首先读取数据，建立模型：
[0049] E(q)=S(q)/So (1)
[0050] E (q)为q空间磁共振衰减信号，S (q)为扩散信号强度，而So未加扩散梯度的脉冲，q 为时一个实验控制参数：
[0051] q=(2ji)-1T5G (2)
[0052] γ为磁旋比，δ为扩散梯度脉冲的持续时间，G代表扩散梯度向量;q空间磁共振衰 减信号正交基基础下在三维量子力学谐振子问题中扩展：
[0054] 其中，Cnlm为变换系数，Φη1ω(ιι，q)为基函数，η为〇