一种非负高阶张量拟牛顿搜索的纤维方向分布估计方法
专利汇可以提供一种非负高阶张量拟牛顿搜索的纤维方向分布估计方法专利检索，专利查询，专利分析的服务。并且一种非负高阶张量拟 牛 顿搜索的 纤维 方向分布估计方法，包括以下步骤：读取脑部磁共振数据，获取施加梯度方向g的磁共振 信号 S(g)和未施加梯度方向的磁共振信号S0，以及梯度方向数据，选取所需的感兴趣区域，并计算该区域的扩散衰减信号S(g)/S0；将感兴趣区域内的每个 体素 的扩散衰减信号S(g)/S0逐个建模为具有扩散形态的椭球分布模型；通过计算张量系数向量c得到扩散函数D(v)，再计算每个 采样 点的扩散函数值，最后将扩散函数值拟合成扩散模型，搜索极值并计算纤维方向。本 发明 采用了用高阶笛卡尔张量拟合纤维方向分布函数平方根的方式，保证了所得纤维方向分布函数的非负性，且 角 度 分辨率 较高，实验效果好。，下面是一种非负高阶张量拟牛顿搜索的纤维方向分布估计方法专利的具体信息内容。
1.一种非负高阶张量拟牛顿搜索的纤维方向分布估计方法，其特征在于：所述的纤维方向分布估计方法包括以下步骤：
(1)数据预处理：读取脑部扩散加权磁共振数据，获取施加梯度方向g时的磁共振信号S(g)和未施加梯度方向时的磁共振信号S0，以及相应的梯度方向数据，选取感兴趣区域，并计算该区域的扩散衰减信号S(g)/S0；
(2)将感兴趣区域内每个体素中的扩散衰减信号逐个建模为具有扩散形态的椭球分布模型，建模过程如下：
2.1)体素微结构建模：将扩散衰减信号假设为沿重建向量v的单条纤维信号响应函数R(v,g)与扩散函数D(v)在球面上的卷积：
其中， 近似看作一个高斯分布函数，g＝{gi∈R1×3|i＝1,...,n}为梯度方向，v＝{vp∈R1×3|p＝1,...,K}为在单位球面上采样的重建向量，R1×3表示维度为1×3的实数域矩阵空间，n和K分别表示梯度方向和重建向量的个数，μ＝εb是表征扩散效率ε与扩散敏感系数b共同影响的一个参数，扩散函数D(v)的表达式如下：
drs表示单项式 的系数，l为高阶张量的阶数，r,s分别表示重建向量v＝(vx,vy,vz)的基方向vx和vy的指数，λj表示第j个张量的张量系数，j＝1,...,m，m表示张量的个数表示第j个张量单项式，并满足 c是由m个张
T
量系数组成的系数向量，表达式为c＝[λ1,λ2,..λm]，F(v)＝[f1(v),f2(v),...fm(v)]由m个张量单项式构成；
2.2)数学模型：
扩散加权磁共振信号有n个扩散梯度方向gi，i＝1,...,n，并且沿重建向量v进行重建，那么系数向量c通过最小化下面的代价函数J(c)求得：
其中，Ei＝S(gi)/S0是第i个扩散梯度方向gi上的衰减信号； 是一个
m×m维的矩阵，其值只与扩散梯度方向gi、重建向量v以及参数μ有关，对于每一个扩散梯度方向gi，都有一个Qi矩阵与之对应；
(3)计算张量系数向量c，得到扩散函数D(v)，再计算每个采样点处的扩散函数值，最后将扩散函数值拟合成扩散模型，搜索极值并计算纤维方向。
2.如权利要求1所述的一种非负高阶张量拟牛顿搜索的纤维方向分布估计方法，其特征在于：所述步骤(3)中，所述张量系数向量c的计算包括以下步骤：
3.1)在单位半球面上均匀采样321个离散的点，以球心为原点获取这321个重建向量v，计算单条纤维响应函数R(v,g)的值，设定高阶张量模型的阶数l，计算单项式矩阵F(v)，进而计算出步骤2.2)中的矩阵Qi；
3.2)使用BFGS拟牛顿搜索算法迭代求解2.2)中的最小化问题，步骤如下：
步骤3.2.1已知代价函数J(c)，选取一个初始点c1作为第一次迭代的搜索起始点，迭代次数计为k＝1，并设置最大迭代次数k_max，计算J(c)的梯度向量表达式：
步骤3.2.2对于第k次迭代，令 表示第k次迭代的梯度向量，计算拟牛顿方向dk＝-Hktk作为第k次迭代的搜索方向，以ck为起点，沿方向dk进行一维搜索，求得本次搜索的可接受步长αk，矩阵Hk表示的是第k次迭代时ck点处Hesse矩阵的逆的近似矩阵，其更新方法如下：
其中，I是单位阵，δk＝ck+1-ck是相邻两次迭代的解的差，βk＝tk+1-tk是相邻两次迭代的梯度向量之差；
步骤3.2.3更新张量系数向量ck+1＝ck+αkdk，计算对应的代价函数值J(ck+1)，当达到最大迭代次数或相邻两次迭代的代价函数值满足 时，则终止迭代，否则，返回步骤3.2.2，进入下一次迭代；参数σ为一个较小的常数，用于判断算法是否收敛到一个局部极小值。
3.如权利要求2所述的一种非负高阶张量拟牛顿搜索的纤维方向分布估计方法，其特征在于：所述步骤(3)还包括3.3)将得到的张量系数用于拟合扩散函数，获取纤维方向分布函数模型，搜索极值并计算纤维方向，步骤如下：
步骤3.3.1对十二面体进行5次细分，得到10242个球面上的相邻等距点，以球心为原点得到相应个数的重建向量V，通过3.2)得到的张量系数向量可以求得扩散函数步骤3.3.2由前面得到的扩散函数求得在10242个重建向量上的扩散函数值，即纤维方向分布函数的值，通过搜索纤维方向分布函数值中的极值点来获取纤维的主方向，极值点的搜索方法如下：
对每一个重建向量Vq，q＝1,…,10242，在10242个重建向量中搜索出与Vq的夹角小于θ的所有向量，比较Vq与这些向量所对应的纤维方向分布函数值的大小，若Vq所对应的值最大，则判断Vq为该体素的一个极值方向；依次遍历所有重建向量，最后得到N个极值方向，这些方向即为纤维的主方向。
一种非负高阶张量拟牛顿搜索的纤维方向分布估计方法
技术领域
背景技术
发明内容
具体实施方式
度方向，v＝{vp∈R |p＝1,...,K}为在单位球面上采样的重建向量，R 表示维度为1×3的实数域矩阵空间，n和K分别表示梯度方向和重建向量的个数，μ＝εb是表征扩散效率ε与扩散敏感系数b共同影响的一个参数，扩散函数D(v)的表达式如下：
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种基于扩散峰度成像的径向峰度各项异性定量方法 | 2020-05-22 | 838 |
| 基于结构先验信息的人体心脏左心室扩散张量估计方法 | 2020-05-17 | 945 |
| 一种扩散张量成像大脑白质纤维聚类方法 | 2020-05-14 | 328 |
| 基于绝对平方梯度的结构张量迹约束下的扩散滤波方法 | 2020-05-14 | 135 |
| 一种非负高阶张量拟牛顿搜索的纤维方向分布估计方法 | 2020-05-26 | 802 |
| 一种扩散张量磁共振成像方法 | 2020-05-11 | 395 |
| 扩散程度可控的各向异性扩散图像去噪增强方法 | 2020-05-19 | 819 |
| 用于DTI纤维跟踪的自动种子点选择 | 2020-05-25 | 674 |
| 一种利用双能减影数据分离胸片图像中成分的方法及装置 | 2020-05-25 | 105 |
| 用于脑白质纤维跟踪的高阶扩散张量混合稀疏成像方法 | 2020-05-18 | 1033 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
