[发明专利]一种弥散张量磁共振图像张量域非局部均值去噪方法

摘要

本发明公开一种弥散张量磁共振图像张量域非局部均值去噪方法，属于数字图像处理和应用数学交叉学科技术领域，解决弥散张量磁共振图像易于被噪声影响的问题。首先，依次遍历该弥散张量磁共振图像的体素，并以遍历到的每个体素为中心设定相应搜索区域；然后将搜索区域内的所有体素逐个与中心体素进行张量矩阵相似性比较；最后，按照张量矩阵相似性高低赋予搜索区域内的体素不同权值，计算加权均值张量矩阵，获得该中心体素的去噪结果。解决了弥散张量磁共振图像易于被噪声影响的问题。

主权项

一种弥散张量磁共振图像张量域非局部均值去噪方法，其特征在于，如下步骤：(1)在背景区域，计算背景区域的灰度直方图，并使用高斯函数拟合灰度直方图，根据高斯函数的方差确定噪声标准差，并设定衰减系数h和搜索区域半径Ω；(2)根据噪声标准差、衰减系数h和搜索区域半径Ω依次遍历弥散张量磁共振图像中的所有体素，并以遍历到的每个体素为中心体素，设定半径为Ω的正方搜索区域Q；(3)将正方搜索区域Q中的所有体素依次与中心体素进行张量矩阵相似性比较，获得描述张量矩阵相似性的相似性距离d；(4)根据衰减系数h、正方搜索区域Q和张量矩阵相似性的相似性距离d获得归一化参数，计算公式如下：Z(p)=Σq∈Qexp(-d2(p,q)h2),]]>其中，p为中心体素，q为比较体素，d(p,q)为中心体素与比较体素的张量矩阵相似性的相似性距离；(5)根据归一化参数Z(p)、衰减系数h和张量矩阵相似性的相似性距离d计算各体素的权值，计算公式如下：w(p,q)=1Z(p)exp(-d2(p,q)h2);]]>(6)根据各体素权值w(p,q)，计算被遍历到体素的加权均值张量矩阵I，获得该中心体素的去噪结果，计算公式如下：I(p)=exp(Σq∈Qw(p,q)log(V(q))),]]>其中，I(p)为中心体素去噪后张量矩阵，V(q)为比较体素的张量矩阵；(7)重复上述步骤，直到遍历完该弥散张量磁共振图像的所有体素。