一种基于流体力学微分方程的脑纤维非对称重构方法技术

一种基于流体力学微分方程的脑纤维非对称重构方法，将纤维束看作流束，通过引入流体力学中的相关概念对纤维重构进行建模；以流体力学中的散度概念来描述纤维在体素内的连续性；通过延伸散度在体素间的概念来描述纤维束的空间一致性；通过这两个约束条件最优化计算非对称型脑纤维微结构重构。实验分别在模拟磁共振数据和实际临床数据上与目前流行的纤维重构方法进行对比，实验结果证明本发明专利技术提出的方法提高了纤维重构的精确性。

An asymmetric reconstruction method of brain fibers based on differential equations of fluid mechanics

An asymmetric reconstruction method of brain fiber based on the differential equation of fluid mechanics, which regards the fiber bundle as a flow bundle, and models the fiber reconstruction by introducing the related concepts in fluid mechanics; describes the continuity of fiber in voxel by using the concept of divergence in fluid mechanics; describes the spatial consistency of fiber bundle by using the concept of divergence between voxels; and describes the spatial consistency of fiber bundle by using these two constraint strips The results show that the reconstruction of the asymmetric brain fiber microstructure can be calculated by the optimization method. The experimental results show that the method proposed in the invention improves the accuracy of fiber reconstruction.

【技术实现步骤摘要】
一种基于流体力学微分方程的脑纤维非对称重构方法技术邻域本专利技术涉及计算机图形学下的医学成像、神经解剖学邻域，尤其是一种脑纤维微结构重构方法。
技术介绍
大脑是控制人类进行逻辑思维、学习记忆、运动及情感等复杂功能活动的综合器官，对人脑工作机理的探究是当前科学研究的前沿和热点。脑白质纤维重构技术通过将具有解剖学意义的大脑纤维空间微结构信息进行成像，是目前在活体大脑中纤维组织微结构成像唯一的非侵入无创方法，已成为大脑认知机理探究、神经类疾病病理分析以及脑手术导航等脑科学研究的重要技术手段。脑纤维微结构重构是脑纤维成像的基础步骤，为纤维束跟踪提供精确的纤维方向估计。传统重构方法往往仅依靠单体素信息进行体素的微结构重构，且模型的类型基本为中心对称型模型，限制了纤维方向重构的准确性。提出新的更为准确的脑纤维微结构重构模型是研究的热点。
技术实现思路
为了克服已有脑纤维微结构重构方法对于单体素信息的依赖性、精度较低以及模型为中心对称的问题，本专利技术提出一种结合体素的邻域信息，基于流体力学微分方程的非对称型脑纤维微结构重构方法。本专利技术解决其技术问题所采用的技术方案是：一种基于流体力学微分方程的脑纤维非对称重构方法，包括如下步骤：步骤一、非对称性纤维轨迹分布(fibertrajectorydistribution,FTD)函数：将纤维束看作流束，流束由一系列的流线集合构成S＝{si,i＝1,···,n}，空间任意一点(x,y,z)的流线切线为该点的流场，以流场来表示体素内任意一点的纤维方向：以整个体素内的流场分布表示FTD，使用一组三元二次多项式来近似体素内的FTD：υ(x,y,z)＝AC(x,y,z)(2)其中，系数矩阵A定义如下：以C表示C(x,y,z)＝[x2,y2,z2,xy,xz,yz,x,y,z,1]T；步骤二、结合纤维束在邻域上的空间连续性，过程如下：2.1、体素内FTD的连续性约束假设水分子在同一纤维束中的扩散位移保持连续性，使用连续不可压缩流体理论，通过引入扩散张量矢量场上纤维流的发散概念来描述纤维轨迹的空间连续性：当FTD不属于神经纤维束的起始或终止区域时，divΩ满足下式：divΩ＝0(5)联立公式(2)、(3)、(4)、(5)得：2.2、FTD体素间的空间连续性约束虽然公式(6)能保证体素内的FTD满足连续性，但是不能表示体素间的纤维束连续性，同一纤维束在相邻体素间是一致的，即对应的FTD应满足体素间的连续性，提出了一种FTD体素间的一致性函数来表征纤维束在相邻体素间的连续性，假设体素是单位立方体，以Nc＝(c0,c1,…,c5)表示中心体素c相邻的六个体素，以Ac＝(A0,A1,…,A5)表示邻域体素的FTD系数，对于流束经过的相邻体素c与ci，其交面任意点的连续性函数divΓ(x,y,z)定义如下：其中Γ为流束经过的相邻体素连接面，S为流束，υ与υi分别为c与ci的FTD，将ci与c映射到同一坐标系得：其中ajk和分别为A和Ai中的元素，整个面Γ上的连续性函数divΨ即为divΓ(x,y,z)的面积分：每一体素有六个相接邻域体素，体素的体素间连续性函数为六个连接面的连续性函数之和：联立式(2)-式(9)得：其中ajk表示中心体素的FTD系数矩阵A中的元素，表示邻域第i个体素的FTD系数矩阵Ai中的元素；步骤三、FTD的计算通过最小化体素内与体素间连续性函数计算FTD旨在使得纤维轨迹分布最贴切于纤维方向分布函数(fiberorientationdistribution,FOD)，FTD中的系数矩阵A可通过最优化如下代价函数来计算：其中Φ(υ(x,y,z))为FOD在点(x,y,z)的概率，为简化计算，我们取中心体素的26邻域C＝[c1,c2,…,c26]中FOD的峰值P＝[p1,p2,…,p26]作为Φ(υ(x,y,z))的近似，这样公式(11)简化如下：流场系数A求得后，即求得体素的FTD。本专利技术的有益效果为：提高了纤维重构的精确性。附图说明图1为中心体素与邻域体素的示意图。具体实施方式下面结合附图对本专利技术做进一步说明。参照图1，一种基于流体力学微分方程的脑纤维非对称重构方法，包括如下步骤：步骤一、非对称性纤维轨迹分布(fibertrajectorydistribution,FTD)函数：将纤维束看作流束，流束由一系列的流线集合构成S＝{si,i＝1,···,n}，空间任意一点(x,y,z)的流线切线为该点的流场，以流场来表示体素内任意一点的纤维方向：以整个体素内的流场分布表示FTD，使用一组三元二次多项式来近似体素内的FTD：υ(x,y,z)＝AC(x,y,z)(2)其中，系数矩阵A定义如下：以C表示C(x,y,z)＝[x2,y2,z2,xy,xz,yz,x,y,z,1]T；步骤二、结合纤维束在邻域上的空间连续性，由于在dMRI数据获取期间，大脑神经纤维可看成是不变形式的，对应于流体力学中的相关概念，即将流束看作不可压缩恒定流；过程如下：2.1、体素内FTD的连续性约束假设水分子在同一纤维束中的扩散位移保持连续性，使用连续不可压缩流体理论，通过引入扩散张量矢量场上纤维流的发散概念来描述纤维轨迹的空间连续性：当FTD不属于神经纤维束的起始或终止区域时，divΩ满足下式：divΩ＝0(5)联立公式(2)、(3)、(4)、(5)得：2.2、FTD体素间的空间连续性约束虽然公式(6)能保证体素内的FTD满足连续性，但是不能表示体素间的纤维束连续性，同一纤维束在相邻体素间是一致的，即对应的FTD应满足体素间的连续性，提出了一种FTD体素间的一致性函数来表征纤维束在相邻体素间的连续性，假设体素是单位立方体，以Nc＝(c0,c1,…,c5)表示中心体素c相邻的六个体素，以Ac＝(A0,A1,…,A5)表示邻域体素的FTD系数，对于流束经过的相邻体素c与ci，其交面任意点的连续性函数divΓ(x,y,z)定义如下：其中Γ为流束经过的相邻体素连接面，S为流束，υ与υi分别为c与ci的FTD，将ci与c映射到同一坐标系得：其中ajk和分别为A和Ai中的元素，整个面Γ上的连续性函数divΨ即为divΓ(x,y,z)的面积分：每一体素有六个相接邻域体素，体素的体素间连续性函数为六个连接面的连续性函数之和：联立式(2)-式(9)得：其中ajk表示中心体素的FTD系数矩阵A中的元素，表示邻域第i个体素的FTD系数矩阵Ai中的元素；步骤三、FTD的计算通过最小化体素内与体素间连续性函本文档来自技高网...

【技术保护点】
1.一种基于流体力学微分方程的脑纤维非对称重构方法，其特征在于，所述方法包括如下步骤：/n步骤一、非对称性纤维轨迹分布FTD函数：/n将纤维束看作流束，流束由一系列的流线集合构成S＝{s

【技术特征摘要】
1.一种基于流体力学微分方程的脑纤维非对称重构方法，其特征在于，所述方法包括如下步骤：
步骤一、非对称性纤维轨迹分布FTD函数：
将纤维束看作流束，流束由一系列的流线集合构成S＝{si,i＝1,…,n}，空间任意一点(x,y,z)的流线切线为该点的流场，以流场来表示体素内任意一点的纤维方向：



以整个体素内的流场分布表示FTD，使用一组三元二次多项式来近似体素内的FTD：
υ(x,y,z)＝AC(x,y,z)(2)
其中，系数矩阵A定义如下：



以C表示C(x,y,z)＝[x2,y2,z2,xy,xz,yz,x,y,z,1]T；
步骤二、结合纤维束在邻域上的空间连续性，过程如下：
2.1、体素内FTD的连续性约束
假设水分子在同一纤维束中的扩散位移保持连续性，使用连续不可压缩流体理论，通过引入扩散张量矢量场上纤维流的发散概念来描述纤维轨迹的空间连续性：



当FTD不属于神经纤维束的起始或终止区域时，divΩ满足下式：
divΩ＝0(5)
联立公式(2)、(3)、(4)、(5)得：



2.2、FTD体素间的空间连续性约束
虽然公式(6)能保证体素内的FTD满足连续性，但是不能表示体素间的纤维束连续性，同一纤维束在相邻体素间是一致的，即对应的FTD应满足体素间的连续性，提出了一种FTD体素间的一致性函数来表征纤维束在相邻体素...

【专利技术属性】
技术研发人员：冯远静，李茂，何建忠，陈余凯，
申请(专利权)人：浙江工业大学，
类型：发明
国别省市：浙江;33