[发明专利]一种口腔CBCT图像去噪方法及系统

权利要求书

1.一种口腔CBCT图像去噪方法，其特征在于，包括如下步骤：

S01.从一组CBCT序列图像中选择一幅待去噪图像g_i，其中，i表示一组CBCT序列图像的第i幅图像，且i1；

S02.对图像g_i进行低通滤波算法去噪，得到去噪后的图像g′_i；

S03.对图像g_i和图像g′_i做差值处理，得到差值图像矩阵diff1，对差值图像矩阵diff1进行绝对值处理abs1＝|diff1|，对绝对值处理后的差值图像矩阵abs1进行归一化处理，得到归一化图像矩阵Norm1：

Norm1_(j,k)＝[abs1_(j,k)-Min(abs1)]/[Max(abs1)-Min(abs1)]

其中，Norm1_(j,k)表示归一化处理后的图像矩阵Norm1中第j行第k列的值，abs1_(j,k)表示绝对值处理后的差值图像矩阵abs1中第j行第k列的值，Max(abs1)表示矩阵abs1中的最大值，Min(abs1)表示矩阵abs1中的最小值；

S04.对第i-1幅图像g_i-1和图像g_i做差值处理，得到差值图像矩阵diff2，对差值图像矩阵diff2进行绝对值处理abs2＝|diff2|，对绝对值处理后的差值图像矩阵abs2进行归一化处理，得到归一化图像矩阵Norm2：

Norm2_(j,k)＝[abs2_(j,k)-Min(abs2)]/[Max(abs2)-Min(abs2)]

其中，Norm2_(j,k)表示归一化处理后的图像矩阵Norm2中第j行第k列的值，abs2_(j,k)表示绝对值处理后的差值图像矩阵abs2中第j行第k列的值，Max(abs2)表示矩阵abs2中的最大值，Min(abs2)表示矩阵abs2中的最小值；

S05.对第i+1幅图像g_i+1和图像g_i做差值处理，得到差值图像矩阵diff3，对差值图像矩阵diff3进行绝对值处理abs3＝|diff3|，对绝对值处理后的差值图像矩阵abs3进行归一化处理，得到归一化图像矩阵Norm3：

其中，Norm3_(j,k)表示归一化处理后的图像矩阵Norm3中第j行第k列的值，abs3_(j,k)表示绝对值处理后的差值图像矩阵abs3中第j行第k列的值，Max(abs3)表示矩阵abs3中的最大值，Min(abs3)表示矩阵abs3中的最小值；

S06.对归一化图像矩阵Norm1、Norm2、Norm3进行加权处理和中值滤波处理、归一化处理，得到权重图像矩阵W；

S07.利用权重图像矩阵W，对图像g_i和图像g′_i采用如下公式进行加权处理，得到最终去噪后的图像：

最终去噪后的图像＝W_(j,k)×g_i(j,k)+(1-W_(j,k))×g′_i(j,k)

其中，W_(j,k)表示权重图像矩阵W中第j行第k列的值，g_i(j,k)表示待去噪图像g_i第j行第k列的像素值，g′_i(j,k)表示去噪后图像g′_i第j行第k列的像素值。

2.根据权利要求1所述的一种口腔CBCT图像去噪方法，其特征在于，所述步骤S02中的低通滤波算法包括高斯滤波算法、均值滤波算法、中值滤波算法。

3.根据权利要求1所述的一种口腔CBCT图像去噪方法，其特征在于，所述步骤S06中的加权处理方法中图像矩阵Norm1、Norm2、Norm3各自的权重之和为1。