首页 / 专利库 / 数学与统计 / 扩散张量 / 一种基于扩散张量成像的脑部特征提取方法

一种基于扩散张量成像的脑部特征提取方法

阅读:802发布:2020-05-13

专利汇可以提供一种基于扩散张量成像的脑部特征提取方法专利检索,专利查询,专利分析的服务。并且一种基于 扩散张量 成像的脑部特征提取方法,包括以下步骤:1),采用基于三维谐振子的重构与估计模型:该模型采用多壳q空间中 水 分子扩散方向数据的映射的方法,简称SHORE模型;2),计算两个基于SHORE模型的标量映射:计算基于SHORE模型下的均方位移和零位移概率,均方位移表示在t时刻水分子距离初始 位置 的距离,即水分子随机运动在空间范围内的度量,标量映射是通过水分子在一定时间内最低限度扩散的概率 密度 ;3),数据的 可视化 显示:通过调用Matplotlib二维图像数据显示工具来呈现SHORE模型通过映射后的图像。本 发明 结合扩散张量成像技术,更清晰、更有效、更能体现脑部特征结构。,下面是一种基于扩散张量成像的脑部特征提取方法专利的具体信息内容。

1.一种基于扩散张量成像的脑部特征提取方法,其特征在于:所述方法包括以下步骤:
1),采用基于三维谐振子的重构与估计模型:
该模型采用多壳q空间中分子扩散方向数据的映射的方法,简称SHORE模型;首先q空间磁共振衰减信号在三维量子学谐振子问题中扩展为正交基,也就是通过三维扩展方案来表示一维的q空间数据;而衰减信号的傅里叶反变换得到能够提供连接信息的平均传播算子;
2),计算两个基于SHORE模型的标量映射:
计算基于SHORE模型下的均方位移和零位移概率,均方位移简称MSD,指的是表示在t时刻水分子距离初始位置的距离,即水分子随机运动在空间范围内的度量,零位移概率简称rtop,所述标量映射是通过水分子在一定时间内最低限度扩散的概率密度
3),数据的可视化显示:
通过调用Matplotlib二维图像数据显示工具来呈现SHORE模型通过映射后的图像;所述Matplotlib是一个二维绘图库并且能够应用于python脚本。
2.如权利要求1所述的一种基于扩散张量成像的脑部特征提取方法,其特征在于:所述步骤3)中,调用Matplotlib中对应功能的类。
3.如权利要求1或2所述的一种基于扩散张量成像的脑部特征提取方法,其特征在于:
所述步骤1)中,建立基于SHORE模型包括以下步骤:
首先读取数据,建立模型:
E(q)=S(q)/S0   (1)
E(q)为q空间磁共振衰减信号,S(q)为扩散信号强度,而S0未加扩散梯度的脉冲,q为时一个实验控制参数:
q=(2π)-1γδG   (2)
γ为磁旋比,δ为扩散梯度脉冲的持续时间,G代表扩散梯度向量;q空间磁共振衰减信号正交基基础下在三维量子力学谐振子问题中扩展:
其中,cnlm为变换系数,Φnlm(u,q)为基函数,n为Φnlm(u,q)的径向阶,l球面谐波函数的阶次,L为l的最大值,N为n的最大值:
u表示每个体素常数的估计数据,ζ是一个比例系数,Γ为Gamma函数, 是广义拉盖尔多项式, 是球面谐波函数,m为度。
4.如权利要求3所述的一种基于扩散张量成像的脑部特征提取方法,其特征在于:所述步骤2)中,计算基于SHORE的标量映射,包括以下步骤:
2.1计算系踪传播算子,衰减信号的傅里叶反变换能够得到系踪传播算子:
R代表净位移向量,是傅里叶共轭,P(R)为系踪传播算子;
2.2计算基于SHORE模型下的均方位移:
首先计算水分子的平均扩散位移矢量,它与体素中平均扩散系数相关,并通过爱因斯坦扩散方程计算:
其中 是扩散位移矢量,Δ表示扩散时间,MD是平均扩散系数,均方位移可以通过概率密度函数的加权平方位移求得:
其中概率密度函数的求法如下:
其中D为扩散张量;
2.3计算SHORE模型下的零位移概率
零位移概率标量映射是通过水分子在一定时间内最低限度扩散的概率密度;零位移概率是给定的概率密度函数中水分子零位移的概率,即:
其中Po零位移概率。
5.如权利要求1或2所述的一种基于扩散张量成像的脑部特征提取方法,其特征在于:
所述步骤3)中,将步骤2)中计算得到的数据放入Matplotlib库中进行绘图,最后通过Matplotlib工具显示出来。

说明书全文

一种基于扩散张量成像的脑部特征提取方法

技术领域

[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)/S0  (1)
[0017] E(q)为q空间磁共振衰减信号,S(q)为扩散信号强度,而S0未加扩散梯度的脉冲,q为时一个实验控制参数:
[0018] q=(2π)-1γδG  (2)
[0019] γ为磁旋比,δ为扩散梯度脉冲的持续时间,G代表扩散梯度向量;q空间磁共振衰减信号正交基基础下在三维量子力学谐振子问题中扩展:
[0020]
[0021] 其中,cnlm为变换系数,Φnlm(u,q)为基函数,n为Φnlm(u,q)的径向阶,l球面谐波函数的阶次,L为l的最大值,N为n的最大值:
[0022]
[0023] u表示每个体素常数的估计数据,ζ是一个比例系数,Γ为Gamma函数,[0024] 是广义拉盖尔多项式, 是球面谐波函数,m为度。
[0025] 所述步骤2)中,计算基于SHORE的标量映射,包括以下步骤:
[0026] 2.1计算系踪传播算子,衰减信号的傅里叶反变换能够得到系踪传播算子:
[0027]
[0028] R代表净位移向量,是傅里叶共轭,P(R)为系踪传播算子;
[0029] 2.2计算基于SHORE模型下的均方位移:
[0030] 首先计算水分子的平均扩散位移矢量,它与体素中平均扩散系数相关,并通过爱因斯坦扩散方程计算:
[0031]
[0032] 其中 是扩散位移矢量,Δ表示扩散时间,MD是平均扩散系数,均方位移可以通过概率密度函数的加权平方位移求得:
[0033]
[0034] 其中概率密度函数的求法如下:
[0035]
[0036] 其中D为扩散张量;
[0037] 2.3计算SHORE模型下的零位移概率
[0038] 零位移概率标量映射是通过水分子在一定时间内最低限度扩散的概率密度;零位移概率是给定的概率密度函数中水分子零位移的概率,即:
[0039]
[0040] 其中Po零位移概率。
[0041] 所述步骤3)中,将步骤2)中计算得到的数据放入Matplotlib库中进行绘图,最后通过Matplotlib工具显示出来。
[0042] 本发明的有益效果为:结合扩散张量成像技术,更清晰、更有效、更能体现脑部特征结构。附图说明
[0043] 图1为本发明中的流程图

具体实施方式

[0044] 以下将结合附图对本发明做进一步详细说明:
[0045] 参照图1,一种基于扩散张量成像的脑部特征提取方法,所述方法包括以下步骤:
[0046] 1),采用基于三维谐振子的重构与估计模型:
[0047] 该模型采用多壳q空间中水分子扩散方向数据的映射的方法,简称SHORE模型;首先q空间磁共振衰减信号在三维量子力学谐振子问题中扩展为正交基,也就是通过三维扩展方案来表示一维的q空间数据;而衰减信号的傅里叶反变换得到能够提供连接信息的平均传播算子;过程如下:
[0048] 首先读取数据,建立模型:
[0049] E(q)=S(q)/S0  (1)
[0050] E(q)为q空间磁共振衰减信号,S(q)为扩散信号强度,而S0未加扩散梯度的脉冲,q为时一个实验控制参数:
[0051] q=(2π)-1γδG  (2)
[0052] γ为磁旋比,δ为扩散梯度脉冲的持续时间,G代表扩散梯度向量;q空间磁共振衰减信号正交基基础下在三维量子力学谐振子问题中扩展:
[0053]
[0054] 其中,cnlm为变换系数,Φnlm(u,q)为基函数,n为Φnlm(u,q)的径向阶,l球面谐波函数的阶次,m为角度,L为l的最大值,N为n的最大值:
[0055]
[0056] u表示每个体素常数的估计数据,ζ是一个比例系数,Γ为Gamma函数, 是广义拉盖尔多项式, 是球面谐波函数;
[0057] 2),计算两个基于SHORE模型的标量映射:
[0058] 计算基于SHORE模型下的均方位移和零位移概率,均方位移简称MSD,指的是表示在t时刻水分子距离初始位置的距离,即水分子随机运动在空间范围内的度量,零位移概率简称rtop,所述标量映射是通过水分子在一定时间内最低限度扩散的概率密度;过程如下:
[0059] 2.1计算系踪传播算子,衰减信号的傅里叶反变换能够得到系踪传播算子:
[0060]
[0061] R代表净位移向量,是傅里叶共轭,P(R)为系踪传播算子;
[0062] 2.2计算两个基于SHORE模型的标量映射:
[0063] 计算基于SHORE模型下的均方位移和零位移概率,均方位移首先计算水分子的平均扩散位移矢量,它与体素中平均扩散系数相关,并通过爱因斯坦扩散方程计算:
[0064]
[0065] 其中 是扩散位移矢量,Δ表示扩散时间,MD是平均扩散系数,均方位移MSD可以通过概率密度函数的加权平方位移求得:
[0066]
[0067] 其中概率密度函数的求法如下:
[0068]
[0069] 其中D为扩散张量;
[0070] 2.3计算SHORE模型下的零位移概率
[0071] 零位移概率这种标量映射是通过水分子在一定时间内最低限度扩散的概率密度;零位移概率是给定的概率密度函数中水分子零位移的概率,即:
[0072]
[0073] 其中Po为零位移概率;
[0074] 3),数据的可视化显示:
[0075] 通过调用Matplotlib二维图像数据显示工具来呈现SHORE模型通过映射后的图像;Matplotlib是一个二维绘图库并且能够应用于python脚本,这里将步骤二中计算得到的数据放入Matplotlib库中进行绘图,最后通过Matplotlib工具显示出来。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈