[发明专利]磁共振弥散张量去噪方法和系统

权利要求书

1.一种磁共振弥散张量去噪方法，其包括：

图像数据获取步骤：获取磁共振弥散加权图像所对应的K空间数据；

去噪步骤：基于磁共振弥散加权成像模型和采样噪声的高斯分布性质，利用弥散相对各向异性的稀疏性，采用最大后验概率估计的方法由所述K空间数据获得每一个空间位置所对应的去噪后的弥散张量矩阵；

弥散参数计算步骤：基于所述去噪后的弥散张量矩阵，获得弥散参数图。

2.根据权利要求1所述的磁共振弥散张量去噪方法，其特征在于，所述去噪步骤包括：

基于磁共振弥散加权成像模型和采样噪声的高斯分布性质，利用弥散相对各向异性的稀疏性，采用最大后验概率估计的方法构建去噪函数模型，所述去噪函数模型参见如下述公式(1)：

    公式(1)

其中，表示弥散张量矩阵的估计值，RA表示弥散相对各向异性，即RA=D12+D22+D32+2D42+2D52+2D62D1+D2+D3-D1+D2+D33,]]>D_i(i＝1,…,6)分别为所述弥散张量矩阵D中的元素，D=D1D4D5D4D2D6D5D6D3;]]>G(·)是作用于弥散相对各向异性的稀疏约束函数，λ为相应的正则化参数；d_m为第m个弥散加权图像所对应的K空间数据；F表示傅立叶编码矩阵；表示第m个弥散加权图像，其中，I₀表示无弥散加权的参考图像，为第m个弥散加权图像的相位，b是弥散加权因子，g_m是第m个弥散加权图像所对应的弥散梯度向量g_m＝(g_xm，g_ym，g_zm)^T；

利用所述K空间数据，求解所述公式(1)，获得每一个空间位置所对应的去噪后的弥散张量矩阵。