一种CT图像预处理方法技术

本发明专利技术涉及图像预处理技术领域，尤其涉及一种CT图像预处理方法，基于工件CT图像矩阵与空气图像矩阵相减，一定程度上去除伪影及随机噪声，其次基于小波变换对CT图像噪声的类型进行识别并选择适合的滤波器针对性处理噪声。本发明专利技术的有益效果是：针对性去除CT图像中的随机噪声、高斯噪声和椒盐噪声，从而减少去噪中的盲目性；基于空气扫描相减法去噪的同时能一定程度上校正伪影；方法适用性强，当系统参数配置改变，重复方法步骤获得参数。

【技术实现步骤摘要】
一种CT图像预处理方法
本专利技术涉及图像预处理
，尤其涉及一种CT图像预处理方法。
技术介绍
CT技术广泛有效地应用于医学诊断和工业无损检测等领域。但CT图像中噪声，伪影的存在降低了CT图像的质量，进而影响基于图像的无损检测及模式识别的可靠性和有效性，因此对图像进行适当的预处理以减少噪声，伪影的影响是重要的一步。在现有技术中，对于CT图像伪影的去除方法主要有投影正弦图校正法和CT图像校正法两类，这两类方法有较好的伪影抑制效果，但一般都涉及若干参数或阈值的选择,且在CT系统参数改变情况下需重新匹配参数及阈值。对于CT图像的噪声预处理，如果直接盲目地对CT图像噪声进行处理往往达不到去噪效果，可能还会模糊图像细节、边缘。中国专利公开号为CN102800073A，专利技术创造名称为一种锥束CT环形伪影的自动判别与校正方法，该专利技术提供了一种锥束CT环形伪影的自动判别与校正方法，重建所有切片层的空气图像并判定存在环形伪影的空气切片图像层号，将该层的被测物体切片图像与该层的空气切片图像对应像素灰度相减，完成环形伪影校正；而本专利技术首先基于工件CT图像矩阵与空气图像矩阵相减，一定程度上去除伪影及随机噪声，更着重于基于小波变换对CT图像噪声的类型进行识别并选择适合的滤波器针对性处理噪声，方法具有针对性减少去噪中的盲目性。
技术实现思路
本专利技术的目的是为了解决现有技术中存在的缺点，而提出的一种CT图像预处理方法。为了实现上述目的，本专利技术采用了如下技术方案：一种CT图像预处理方法，方法步骤如下：首先基于工件CT图像矩阵与空气图像矩阵相减，一定程度上去除伪影及随机噪声，其次基于小波变换对CT图像噪声的类型进行识别并选择适合的滤波器针对性处理噪声；步骤一，对CT系统配置合适的参数，对模板进行扫描获得CT图像；步骤二，在CT系统同一参数配置下，扫描空气并重建成像；步骤三，在CT系统同一参数配置下，对工件进行切片扫描并重建成像；步骤四，基于空气扫描相减法，且空气扫描相减法不涉及参数配置；步骤五，基于滤波法去噪，提取对角线高频成分直方图信噪比S1,并与参数S_para比较，滤波大于则选用中值滤波器去除椒盐噪声；反之则选用高斯滤波器去除高斯噪声。优选的，所述步骤一中，对CT系统配置合适的参数，对模板进行扫描获得CT图像，步骤一的具体实施如下：在matlab中对CT图像分别加入方差不同的零均值高斯噪声和密度不同的椒盐噪声；基于小波变换，提取对角线高频成分直方图信噪比的均值参数S_para；基于小波变换，提取CT图像高频子带系数的能量分布；采用db4小波分析CT图像，会生成一个低频，三个方向(水平，垂直，对角线)的高频成分，取对角线方向的高频成分，计算其直方图信噪比，高频成分直方图信噪比可表示为：其中，F(i,j)为高频系数矩阵；其中，n为总的高频系数个数；fk是F(i,j)的分布率。优选的，所述步骤二中，在CT系统同一参数配置下，扫描空气并重建成像，从而获得空气图像矩阵dicom0。优选的，所述步骤三中，在CT系统同一参数配置下，对工件进行切片扫描并重建成像，从而获得工件CT图像矩阵dicom_i。优选的，所述步骤四中，基于空气扫描相减法，即工件CT图像矩阵dicom_i与空气图像矩阵dicom0相减，获得去噪图像dicom_i_1，且空气扫描相减法不涉及参数配置。优选的，所述步骤五中，基于滤波法去噪，即对去噪图像dicom_i_1小波变换，提取对角线高频成分直方图信噪比S1,并与参数S_para比较，滤波值大于则选用中值滤波器去除椒盐噪声，反之则选用高斯滤波器去除高斯噪声；从而获得去噪图像dicom_i_2。本专利技术的有益效果是：1、针对性去除CT图像中的随机噪声、高斯噪声和椒盐噪声，从而减少去噪中的盲目性；2、基于空气扫描相减法去噪的同时能一定程度上校正伪影；3、方法适用性强，当系统参数配置改变，重复方法步骤获得参数。附图说明图1为本专利技术提出的一种CT图像预处理方法的方法流程图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。参照图1，一种CT图像预处理方法，方法步骤如下：首先基于工件CT图像矩阵与空气图像矩阵相减，一定程度上去除伪影及随机噪声，其次基于小波变换对CT图像噪声的类型进行识别并选择适合的滤波器针对性处理噪声；步骤一，对CT系统配置合适的参数，对模板进行扫描获得CT图像；步骤二，在CT系统同一参数配置下，扫描空气并重建成像；步骤三，在CT系统同一参数配置下，对工件进行切片扫描并重建成像；步骤四，基于空气扫描相减法，且空气扫描相减法不涉及参数配置；步骤五，基于滤波法去噪，提取对角线高频成分直方图信噪比S1,并与参数S_para比较，滤波大于则选用中值滤波器去除椒盐噪声；反之则选用高斯滤波器去除高斯噪声。所述步骤一中，对CT系统配置合适的参数，对模板进行扫描获得CT图像，步骤一的具体实施如下：在matlab中对CT图像分别加入方差不同的零均值高斯噪声和密度不同的椒盐噪声；基于小波变换，提取对角线高频成分直方图信噪比的均值参数S_para；基于小波变换，提取CT图像高频子带系数的能量分布；采用db4小波分析CT图像，会生成一个低频，三个方向(水平，垂直，对角线)的高频成分，取对角线方向的高频成分，计算其直方图信噪比，高频成分直方图信噪比可表示为：其中，F(i,j)为高频系数矩阵；其中，n为总的高频系数个数；fk是F(i,j)的分布率。所述步骤二中，在CT系统同一参数配置下，扫描空气并重建成像，从而获得空气图像矩阵dicom0；所述步骤三中，在CT系统同一参数配置下，对工件进行切片扫描并重建成像，从而获得工件CT图像矩阵dicom_i；所述步骤四中，基于空气扫描相减法，即工件CT图像矩阵dicom_i与空气图像矩阵dicom0相减，获得去噪图像dicom_i_1，且空气扫描相减法不涉及参数配置；所述步骤五中，基于滤波法去噪，即对去噪图像dicom_i_1小波变换，提取对角线高频成分直方图信噪比S1,并与参数S_para比较，滤波值大于则选用中值滤波器去除椒盐噪声，反之则选用高斯滤波器去除高斯噪声；从而获得去噪图像dicom_i_2。进一步的，本专利技术提供的一种CT图像预处理方法，首先基于工件CT图像矩阵与空气图像矩阵相减，一定程度上去除伪影及随机噪声，方法简单有效；其次基于小波变换对CT图像噪声的类型进行识别并选择适合的滤波器针对性处理噪声，方法具有针对性减少去噪中的盲目性。本专利技术方法可以有效提高CT图像质量，为后本文档来自技高网...

【技术保护点】
1.一种CT图像预处理方法，其特征在于，方法步骤如下：/n首先基于工件CT图像矩阵与空气图像矩阵相减，一定程度上去除伪影及随机噪声，其次基于小波变换对CT图像噪声的类型进行识别并选择适合的滤波器针对性处理噪声；/n步骤一，对CT系统配置合适的参数，对模板进行扫描获得CT图像；/n步骤二，在CT系统同一参数配置下，扫描空气并重建成像；/n步骤三，在CT系统同一参数配置下，对工件进行切片扫描并重建成像；/n步骤四，基于空气扫描相减法，且空气扫描相减法不涉及参数配置；/n步骤五，基于滤波法去噪，提取对角线高频成分直方图信噪比S1,并与参数S_para比较，滤波大于则选用中值滤波器去除椒盐噪声；反之则选用高斯滤波器去除高斯噪声。/n

【技术特征摘要】
1.一种CT图像预处理方法，其特征在于，方法步骤如下：
首先基于工件CT图像矩阵与空气图像矩阵相减，一定程度上去除伪影及随机噪声，其次基于小波变换对CT图像噪声的类型进行识别并选择适合的滤波器针对性处理噪声；
步骤一，对CT系统配置合适的参数，对模板进行扫描获得CT图像；
步骤二，在CT系统同一参数配置下，扫描空气并重建成像；
步骤三，在CT系统同一参数配置下，对工件进行切片扫描并重建成像；
步骤四，基于空气扫描相减法，且空气扫描相减法不涉及参数配置；
步骤五，基于滤波法去噪，提取对角线高频成分直方图信噪比S1,并与参数S_para比较，滤波大于则选用中值滤波器去除椒盐噪声；反之则选用高斯滤波器去除高斯噪声。


2.根据权利要求1所述的一种CT图像预处理方法，其特征在于，所述步骤一中，对CT系统配置合适的参数，对模板进行扫描获得CT图像，步骤一的具体实施如下：
在matlab中对CT图像分别加入方差不同的零均值高斯噪声和密度不同的椒盐噪声；
基于小波变换，提取对角线高频成分直方图信噪比的均值参数S_para；
基于小波变换，提取CT图像高频子带系数的能量分布；
采用db4小波分析CT图像，会生成一个低频，三个方向(水平，垂直，对角线)的高频成分，取对角...

【专利技术属性】
技术研发人员：张广才，何继荣，陈志芸，
申请(专利权)人：扬州哈工博浩智能科技有限公司，
类型：发明
国别省市：江苏;32