[发明专利]一种CT图像预处理方法

权利要求书

1.一种CT图像预处理方法，其特征在于，方法步骤如下：

首先基于工件CT图像矩阵与空气图像矩阵相减，一定程度上去除伪影及随机噪声，其次基于小波变换对CT图像噪声的类型进行识别并选择适合的滤波器针对性处理噪声；

步骤一，对CT系统配置合适的参数，对模板进行扫描获得CT图像；

步骤二，在CT系统同一参数配置下，扫描空气并重建成像；

步骤三，在CT系统同一参数配置下，对工件进行切片扫描并重建成像；

步骤四，基于空气扫描相减法，且空气扫描相减法不涉及参数配置；

步骤五，基于滤波法去噪，提取对角线高频成分直方图信噪比S1,并与参数S_para比较，滤波大于则选用中值滤波器去除椒盐噪声；反之则选用高斯滤波器去除高斯噪声。

2.根据权利要求1所述的一种CT图像预处理方法，其特征在于，所述步骤一中，对CT系统配置合适的参数，对模板进行扫描获得CT图像，步骤一的具体实施如下：

在matlab中对CT图像分别加入方差不同的零均值高斯噪声和密度不同的椒盐噪声；

基于小波变换，提取对角线高频成分直方图信噪比的均值参数S_para；

基于小波变换，提取CT图像高频子带系数的能量分布；

采用db4小波分析CT图像，会生成一个低频，三个方向(水平，垂直，对角线)的高频成分，取对角线方向的高频成分，计算其直方图信噪比，高频成分直方图信噪比可表示为：

其中，F(i,j)为高频系数矩阵；

其中，n为总的高频系数个数；

f_k是F(i,j)的分布率。

3.根据权利要求1所述的一种CT图像预处理方法，其特征在于，所述步骤二中，在CT系统同一参数配置下，扫描空气并重建成像，从而获得空气图像矩阵dicom0。

4.根据权利要求1所述的一种CT图像预处理方法，其特征在于，所述步骤三中，在CT系统同一参数配置下，对工件进行切片扫描并重建成像，从而获得工件CT图像矩阵dicom_i。

5.根据权利要求1所述的一种CT图像预处理方法，其特征在于，所述步骤四中，基于空气扫描相减法，即工件CT图像矩阵dicom_i与空气图像矩阵dicom0相减，获得去噪图像dicom_i_1，且空气扫描相减法不涉及参数配置。

6.根据权利要求1所述的一种CT图像预处理方法，其特征在于，所述步骤五中，基于滤波法去噪，即对去噪图像dicom_i_1小波变换，提取对角线高频成分直方图信噪比S1,并与参数S_para比较，滤波值大于则选用中值滤波器去除椒盐噪声，反之则选用高斯滤波器去除高斯噪声；从而获得去噪图像dicom_i_2。