利索能及
我要发布
收藏
专利号: 2014107582906
申请人: 中国科学院深圳先进技术研究院
专利类型:发明专利
专利状态:已下证
更新日期:2024-12-09
缴费截止日期: 暂无
联系人

摘要:

权利要求书:

1.一种磁共振弥散张量去噪方法,其特征在于,包括:采集K空间数据;

根据采集的K空间数据,求取弥散张量的最大似然估计值;

利用弥散张量和弥散参数的稀疏性,将求取的最大似然估计值作为初始值,计算弥散张量的最大后验概率估计值;

根据求取的所述最大后验概率估计值,计算弥散系数;

其中,利用弥散张量和弥散参数的稀疏性,将求取的最大似然估计值作为初始值,计算弥散张量的最大后验概率估计值,包括:按照以下公式计算弥散张量的最大后验概率估计值:其中, R(·)表示作用于弥散张量D的惩罚函数,G(·)表示作用于弥散参数FA的惩罚函数,γ1、γ2表示正则化参数、D1,D2,D3,D4,D5,D6代表不同方向的6个弥散张量。

2.如权利要求1所述的方法,其特征在于,根据采集的K空间数据,求取弥散张量的最大似然估计值,包括:按照以下公式计算弥散张量的最大似然估计值:其中,表示弥散张量的最大似然估计值,dm表示第m个弥散加权图像所对应的K空间数据,F表示傅立叶编码矩阵,I0表示无弥散加权的参考图像, 表示第m个弥散加权图像的相位,b表示弥散加权因子,gm=(gxm,gym,gzm)T表示第m个梯度向量。

3.如权利要求2所述的方法,其特征在于,R(·)约束每个弥散张量的稀疏性,或者R(·)约束弥散张量的联合稀疏性。

4.如权利要求3所述的方法,其特征在于:

在R(·)约束每个弥散张量的稀疏性的情况下,在R(·)约束弥散张量的联合稀疏性的情况下:其中,Ψ为运算符,表示稀疏变换,ψi表示Ψ的第i列,||·||1表示求L1范数,||·||2表示求L2范数。

5.一种磁共振弥散张量去噪装置,其特征在于,包括:K空间数据采集模块,用于采集K空间数据;

最大似然估计求取模块,用于根据采集的K空间数据,求取弥散张量的最大似然估计值;

最大后验概率估计模块,用于利用弥散张量和弥散参数的稀疏性,将求取的最大似然估计值作为初始值,计算弥散张量的最大后验概率估计值;

弥散系数计算模块,用于根据求取的所述最大后验概率估计值,计算弥散系数;

其中,所述最大后验概率估计模块,具体用于按照以下公式计算弥散张量的最大后验概率估计值:其中, R(·)表示作用于弥散张量D的惩罚函数,G(·)表示作用于弥散参数FA的惩罚函数,γ1、γ2表示正则化参数、D1,D2,D3,D4,D5,D6代表不同方向的6个弥散张量。

6.如权利要求5所述的装置,其特征在于,所述最大似然估计求取模块,具体用于按照以下公式计算弥散张量的最大似然估计值:其中, 表示弥散张量的最大似然估计值,dm表示第m个弥散加权图像所对应的K空间数据,F表示傅立叶编码矩阵,I0表示无弥散加权的参考图像, 表示第m个弥散加权图像的相位,b表示弥散加权因子,gm=(gxm,gym,gzm)T表示第m个梯度向量。

7.如权利要求6所述的装置,其特征在于,R(·)约束每个弥散张量的稀疏性,或者R(·)约束弥散张量的联合稀疏性。

8.如权利要求7所述的装置,其特征在于:

在R(·)约束每个弥散张量的稀疏性的情况下,在R(·)约束弥散张量的联合稀疏性的情况下:其中,Ψ为运算符,表示稀疏变换,ψi表示Ψ的第i列,||·||1表示求L1范数,||·||2表示求L2范数。