首页 / 专利库 / 数学与统计 / 扩散张量 / 一种高阶张量的纤维取向分布估计稀疏去卷积方法

一种高阶张量的纤维取向分布估计稀疏去卷积方法

阅读:747发布:2020-05-22

专利汇可以提供一种高阶张量的纤维取向分布估计稀疏去卷积方法专利检索,专利查询,专利分析的服务。并且一种高阶张量的 纤维 取向分布估计稀疏去卷积方法,包括如下步骤:1)使用高阶张量?纤维取向分布HOT?ODF卷积表示扩散 信号 衰减轮廓信号的纤维响应函数;2)l2惩罚项;3)l0和l1惩罚项;4)纤维取向分布域的稀疏,实现纤维取向分布估计。本 发明 提供一种分布域充分表征实际稀疏、准确性较高的高阶张量的纤维取向分布估计稀疏去卷积方法。,下面是一种高阶张量的纤维取向分布估计稀疏去卷积方法专利的具体信息内容。

1.一种高阶张量的纤维取向分布估计稀疏去卷积方法,其特征在于:所述方法包括如下步骤:
1)使用高阶张量-纤维取向分布HOT-ODF卷积表示扩散信号衰减轮廓信号的纤维响应函数:
其中,D(υ)为HOT-ODF基函数,S(g)为扩散信号强度,So为未加扩散梯度的脉冲,其中g是磁场梯度方向,υ是在球体S2上的三维单位向量;R(υ,g)代表扩散信号衰减,是一个单一相干的定向纤维组测定的轴向对称响应函数:
其中,μ表示对扩散敏感性系数b和各向异性相互作用的影响程度信号衰减;R(υ,g)能够从一个z轴对齐旋转对称的纤维组区域的信号衰减来估计;HOT-ODF函数D(υ)通过笛卡尔张量近似,对于l阶张量:
υ=(υ1,υ2,υ3)T         (4)
υi(i=1,2,3)代表三维单位向量υ的第i列向量;r,s都代表阶次,drs代表扩散系数,D(υ)简化成:
xj是与drs相关的系数, 是第j个张量单项式,j是系数:
其中是公式(5)中的项数m是j的上限,是有界的:
m≤(l+1)(l+2)/2        (7)
给定一个磁场梯度gi,i=1,2,...和扩散加权的b值b相关联的DW-MRI信号衰减Si/S0的一个数据集,第l阶张量的系数可通过最小化下列能量函数E相对于所述未知多项式系数估计:
估计所述HOT系数在公式(8)构成形式的离散逆问题:
y=Ax+ξ         (9)
其中,A是n×m的矩阵且 x是包含未知系数xj的m维向量,y是包含
yi=Si/S0的n维向量,ξ代表噪声和测量误差;
2)l2惩罚项
球谐函数的Tikhonov正则化方式,即非负约束球形去卷积,涉及最小化两个项数的加权和:
其中,是正则化参数,L是包含在精确解的平滑先验信息约束矩阵;Tikhonov正则化导致了整体上减少了在偏差-方差权衡上的均方误差;Tikhonov正则化方法的问题是,x估计的解决方案中每个元素通常是非零和非负的,而真实的解决方案涉及非零元素子集;因此公式(10)的解决方法是非稀疏的;
3)l0和l1惩罚项
纤维取向的真实分布被认为是稀疏,假定只有少数x的元素为非零,它们产生的观测信号y的最大稀疏通过最小化l0表示:
其中,ξ代表阈值,||x||0代表x中非零元素的数目,考虑到一个体素内只存在几个纤维群体,然后承担一个体素内小于k个纤维束的群体;公式(11)中的元素能够被整理为寻找一个给定稀疏值k可能的最小误差:
而最小化l0问题则是一个非确定性多项式问题,它放宽为l1问题;通过||ωx||1近似||x||0,就变成了l1问题:
其中,||ωx||1代表权重先验信息的凸约束;x的稀疏元素趋向于零和权重ω一起迭代更新;权重迭代 ω(t+1)代表权重ω的第t+1次迭代, 代表第i个稀疏元素的第t次迭代,通过参数λ≥0平衡以最大限度地减少错误和稀疏这两项双重目标:
公式(14)的问题中,纤维取向分布的稀疏约束被施加在x中;
4)纤维取向分布域的稀疏
稀疏表示是通过最大限度地减少l0问题:
T
其中,D(υ,x)=xF,x=[x1,x2,...,xm],F=[f1(υ),f2(υ),...,fm(υ)]是从HOT-ODF在公式(5)计算出纤维取向分布的值;公式(15)同样放宽为l1问题:
D(υ,x)的稀疏元素趋向于零同ω一起迭代更新;权重ω的修正方案定义为:
其中, 是第i个权重的第t+1次迭代,β是一个给定的向量
正元素与Di大小相同;δ是一个正常数;公式(14)和(15)是拉索型最优化问题,应对该拉索型最优化问题的方法是利用弹性网模型:
其中,参数λ用于平衡稀疏约束项和噪声项;假设λ是标准化D的高斯分布,其中是标准化D的最小值;max(D)-min(D)在单位空间中通常是恒定不变
的,所以设置λ=e-[D-min(D)],参数α∈[0,1]控制稀疏和平滑度的量,参数ψ是用来约束向量元素:
其中,Γ是可控的阈值,低于该阈值相应的纤维取向密度幅值定为零,实现纤维取向分布估计。

说明书全文

一种高阶张量的纤维取向分布估计稀疏去卷积方法

技术领域

[0001] 本发明涉及计算机图形学下的医学成像、神经解剖学领域,尤其是一种纤维取向分布估计方法。

背景技术

[0002] 磁共振弥散加权成像(DW-MRI)是一种非侵入性的成像技术能够揭示人脑的微观结构信息;扩散张量成像(DTI)通常可用来近似估计扩散函数但至少需要7组DW-MR图像;然而,体素经常受部分体积效应的影响,同时扩散张量模型未能精确地表示神经束几何形状中多个纤维群体的体素;建立准确的估计方法需要对纤维进行取向,解决神经结构地区复合纤维模型就必须使用高分辨率扩散成像(HARDI)这一新方法;基于模型的方法依赖于一个更复杂的模型来表征扩散加权信号衰减;这样的方法包括多张量模型,定向函数球谐去卷积,以及高阶张量(HOT)。

发明内容

[0003] 为了克服已有纤维取向分布方法的分布域无法充分表征实际稀疏、准确性较低的不足,本发明提供一种分布域充分表征实际稀疏、准确性较高的高阶张量的纤维取向分布估计稀疏去卷积方法。
[0004] 本发明解决其技术问题所采用的技术方案是:
[0005] 一种高阶张量的纤维取向分布估计稀疏去卷积方法,包括如下步骤:
[0006] 1)使用高阶张量-纤维取向分布HOT-ODF卷积表示扩散信号衰减轮廓信号的纤维响应函数:
[0007]
[0008] 其中,D(υ)为HOT-ODF基函数,S(g)为扩散信号强度,So为未加扩散梯度的脉冲,其中g是磁场梯度方向,υ是在球体S2上的三维单位向量;R(υ,g)代表扩散信号衰减,是一个单一相干的定向纤维组测定的轴向对称响应函数:
[0009]
[0010] 其中,μ表示对扩散敏感性系数b和各向异性相互作用的影响程度信号衰减;R(υ,g)能够从一个z轴对齐旋转对称的纤维组区域的信号衰减来估计;HOT-ODF函数D(υ)通过笛卡尔张量近似,对于l阶张量:
[0011]
[0012] υ=(υ1,υ2,υ3)T  (4)
[0013] υi(i=1,2,3)代表三维单位向量υ的第i列向量;r,s都代表阶次,drs代表扩散系数,D(υ)简化成
[0014]
[0015] xj是与drs相关的系数, 是第j个张量单项式,j是系数:
[0016]
[0017] 其中是公式(5)中的项数m是j的上限,是有界的:
[0018] m≤(l+1)(l+2)/2  (7)
[0019] 给定一个磁场梯度gi,i=1,2,...和扩散加权的b值b相关联的DW-MRI信号衰减Si/S0的一个数据集,第l阶张量的系数可通过最小化下列能量函数E相对于所述未知多项式系数估计:
[0020]
[0021] 估计所述HOT系数在公式(8)构成形式的离散逆问题:
[0022] y=Ax+ξ  (9)
[0023] 其中,A是n×m的矩阵且 x是包含未知系数xj的m维向量,y是包含yi=Si/S0的n维向量,ξ代表噪声和测量误差;
[0024] 2)l2惩罚项
[0025] 球谐函数的Tikhonov正则化方式,即非负约束球形去卷积,涉及最小化两个项数的加权和:
[0026]
[0027] 其中,是正则化参数,L是包含在精确解的平滑先验信息约束矩阵;Tikhonov正则化导致了整体上减少了在偏差-方差权衡上的均方误差;Tikhonov正则化方法的问题是,x估计的解决方案中每个元素通常是非零和非负的,而真实的解决方案涉及非零元素子集;因此公式(10)的解决方法是非稀疏的;
[0028] 3)l0和l1惩罚项
[0029] 纤维取向的真实分布被认为是稀疏,假定只有少数x的元素为非零,它们产生的观测信号y的最大稀疏通过最小化l0表示:
[0030]
[0031] 其中,ξ代表阈值,||x||0代表x中非零元素的数目,考虑到一个体素内只存在几个纤维群体,然后承担一个体素内小于k个纤维束的群体;公式(11)中的元素能够被整理为寻找一个给定稀疏值k可能的最小误差:
[0032]
[0033] 而最小化l0问题则是一个非确定性多项式问题,它放宽为l1问题;通过||ωx||1近似||x||0,就变成了l1问题:
[0034]
[0035] 其中,||ωx||1代表权重先验信息的凸约束;x的稀疏元素趋向于零和权重ω一起迭代更新;权重迭代 代表权重ω的第t+1次迭代, 代表第i个稀疏元素的第t次迭代,通过参数λ≥0平衡以最大限度地减少错误和稀疏这两项双重目标:
[0036]
[0037] 公式(14)的问题中,纤维取向分布的稀疏约束被施加在x中;
[0038] 4)纤维取向分布域的稀疏
[0039] 稀疏表示是通过最大限度地减少l0问题:
[0040]
[0041] 其中,D(υ,x)=xF,x=[x1,x2,...,xm]T,F=[f1(υ),f2(υ),...,fm(υ)]是从HOT-ODF在公式(5)计算出纤维取向分布的值;公式(15)同样放宽为l1问题:
[0042]
[0043] D(υ,x)的稀疏元素趋向于零同ω一起迭代更新;权重ω的修正方案定义为:
[0044]
[0045] 其中, 是第i个权重的第t+1次迭代,β是一个给定的向量正元素与Di大小相同;δ是一个正常数;公式(14)和(15)是拉索型最优化问题,应对该拉索型最优化问题的方法是利用弹性网模型:
[0046]
[0047] 其中,参数λ用于平衡稀疏约束项和噪声项;假设λ是标准化D的高斯分布,其中是标准化D的最小值;max(D)-min(D)在单位空间中通常是恒定不变的,所以设置λ=e-[D-min(D)],参数a∈[0,1]控制稀疏和平滑度的量,参数ψ是用来约束向量元素:
[0048]
[0049] 其中,Γ是可控的阈值,低于该阈值相应的纤维取向密度幅值定为零。
[0050] 本发明的技术构思为:高阶张量拟合模型的扩散估计,提出了一个稀疏高阶张量正规化模型,它通过直接对纤维取向分布的表示实施稀疏约束,而不是将其强加给基函数系数的方法。
[0051] 结合稀疏模型l2惩罚项的稳定作用和l1惩罚项的稀疏激励效应以充分表征在纤维取向分布领域实际稀疏,另通过迭代加权正则化方法来解决去卷积问题;最后去卷积技术与现有l2或l1正则化的方法对比测试合成数据和真实人的大脑。
[0052] 本发明的有益效果主要表现在:分布域充分表征实际稀疏、准确性较高。

具体实施方式

[0053] 下面对本发明作进一步描述。
[0054] 一种高阶张量的纤维取向分布估计稀疏去卷积方法,包括如下步骤:
[0055] 1)扩散加权成像的图像与分子有关,其测量结果可以呈现水分子扩散在球体沿不同空间方向的敏感度;使用高阶张量-纤维取向分布(HOT-ODF)卷积表示扩散信号衰减轮廓信号的纤维响应函数:
[0056]
[0057] D(υ)为HOT-ODF基函数,S(g)为扩散信号强度,So为未加扩散梯度的脉冲,其中g是磁场梯度方向,υ是在球体S2上的三维单位向量;R(υ,g)代表扩散信号衰减,是一个单一相干的定向纤维组测定的轴向对称响应函数:
[0058]
[0059] 其中μ表示对扩散敏感性系数b和各向异性相互作用的影响程度信号衰减;R(υ,g)能够从一个z轴对齐旋转对称的纤维组区域的信号衰减来估计;HOT-ODF函数D(υ)通过笛卡尔张量近似,对于l阶张量:
[0060]
[0061] υ=(υ1,υ2,υ3)T  (4)
[0062] υi(i=1,2,3)代表三维单位向量υ的第i列向量;r,s都代表阶次,drs代表扩散系数,D(υ)又可以简化成:
[0063]
[0064] xj是与drs相关的系数, 是第j个张量单项式,j是系数:
[0065]
[0066] 其中是公式(5)中的项数m是j的上限,是有界的:
[0067] m≤(l+1)(l+2)/2  (7)
[0068] 给定一个磁场梯度gi,i=1,2,...和扩散加权的b值b相关联的DW-MRI信号衰减Si/S0的一个数据集,第l阶张量的系数可通过最小化下列能量函数E相对于所述未知多项式系数估计:
[0069]
[0070] 估计所述HOT系数在公式(8)构成形式的离散逆问题:
[0071] y=Ax+ξ  (9)
[0072] A是n×m的矩阵且 x是包含未知系数xj的m维向量,y是包含yi=Si/S0的n维向量,ξ代表噪声和测量误差,在这里代表着Rician噪声;这问题关系到张量阶次l和磁场梯度方向的数量n。
[0073] 2)l2惩罚项
[0074] 一个常见的正则化方法是最小二乘法的惩罚;球谐函数的Tikhonov正则化方式,即非负约束球形去卷积,涉及最小化两个项数的加权和:
[0075]
[0076] 其中 是正则化参数,L是包含在精确解的平滑先验信息约束矩阵;Tikhonov正则化导致了整体上减少了在偏差-方差权衡上的均方误差;Tikhonov正则化方法的问题是,x估计的解决方案中每个元素通常是非零和非负的,而真实的解决方案涉及非零元素子集;因此公式(10)的解决方法通常是非稀疏的。
[0077] 3)l0和l1惩罚项
[0078] 这步骤是解决观测数据只占用x元素一个小子集的方法;纤维取向的真实分布可以被认为是稀疏假定只有少数x的元素为非零,它们产生的观测信号y的最大稀疏通过最小化l0表示:
[0079]
[0080] 其中ξ代表阈值,||x||0代表x中非零元素的数目,考虑到一个体素内只存在几个纤维群体,然后承担一个体素内小于k个纤维束的群体;公式(11)中的元素能够被整理为寻找一个给定稀疏值k可能的最小误差:
[0081]
[0082] 而最小化l0问题则是一个非确定性多项式问题,它通常放宽为l1问题;通过||ωx||1近似||x||0,这就变成了l1问题:
[0083]
[0084] 其中||ωx||1代表权重先验信息的凸约束;x的稀疏元素趋向于零和权重ω一起迭代更新;权重迭代 ω(t+1)代表权重ω的第t+1次迭代, 代表第i个稀疏元素的第t次迭代它也可以通过参数λ≥0平衡以最大限度地减少错误和稀疏这两项双重目标:
[0085]
[0086] 公式(14)的问题中,纤维取向分布的稀疏约束被施加在x中。
高效检索全球专利

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

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

申请试用

分析报告

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

申请试用

QQ群二维码
意见反馈