[发明专利]一种局部参数自适应的工业CT图像分割方法

说明书

本发明公开了一种局部参数自适应的工业CT图像分割方法，提供一种能较好地对灰度不均匀图像进行分割的算法，即基于Phansalkar算法，利用遗传算法优化找到使目标像素方差与背景像素方差的绝对差值达到最小时的Phansalkar算法参数，将优化所得参数作为Phansalkar算法参数对当前窗口分割。此算法能较好地分割灰度不均匀图像，且此方法无参数，可移植性好。

技术领域

本发明属于工业CT图像处理领域，涉及一种局部参数自适应的工业CT图像分割方法。

背景技术

工业CT图像分割技术是工业CT无损检测技术的重要环节，其精度直接影响工业CT检测精度，因此，提高工业CT图像分割精度十分重要。工业CT图像分割的难点在于其分割精度不足，不能满足工业测量要求。主要原因在于，工业CT图像中包含大量干扰信息，例如噪声、伪影和低对比度等，导致图像分割时易产生误分割，导致最终检测精度下降。

工业CT图像分割算法均是由自然图像分割算法演变而来，因为自然图像分割算法不强调精度，因此，需要针对工业CT图像的特征对自然图像分割算法进行改进。为提高工业CT图像分割结果的精度，需要开发可同时处理图像噪声、伪影、并对低对比度图像处理结果优的算法。现有的算法主要分为模板法、阈值分割法、模糊聚类法、活动轮廓模型法等。

其中，模板法将等检测件的标准件进行扫描，得到其CT图像作为标准图像，而后通过找出标准件与被测试件不同的方法，分割工业CT图像，此算法只适用于工业CT图像中的缺陷分割。而阈值法因其对伪影、噪声适应性较差，只适用于无伪影噪声的图像。模糊聚类法对初始聚类中心敏感易陷入极值且没有考虑空间信息对噪声敏感。活动轮廓模型法抗噪性较好，对伪影处理结果不理想。

目前大多数分割方法无法完成具有伪影和噪声的图像的高精度分割，并且易产生误分割。阮健等人在《CT理论与应用研究》(2010,19(1):56-61)的文章“一种工业CT图像的分割算法”提出了一种适用于工业CT图像的分割算法，首先利用最大类间方差法和图像处理方法处理了外层伪影，然后利用聚类迭代的方法处理中心空气，得到感兴趣的区域。实验结果表明，此算法能够对具有先验知识的工业断层图像准确地提取感兴趣的区域。邱钊等人在《计算机工程》(2004,30(8):159-161)中的文章“基于边缘信息的工业CT图像分割法”提出了采用多尺度二维小波变换模极大值作为分裂合并图像分割的一致性度量,采用四叉树结构作为图像块的初始划分,设置合适的门限值抑制噪声的干扰并提取出图像边缘信息。向新程等人在《原子能科学技术》(2007,41(3))中的文章“二维直方图阈值法分割工业CT图像的研究”提出了对CT图像先后采用了比特窗算法进行预处理和二维直方图阈值法进行分割，实现一类工件缺陷的分割。然而上述方法对零件的伪影处理结果仍不能满足高精度分割要求，且对零件形状有一定封闭性要求，可移植性差。

发明内容

针对传统工业CT图像处理方法对图像中的伪影和噪声的分割结果精度不足的问题，本发明提供一种局部参数自适应的工业CT图像分割方法，对图像中的每一个局部窗口，改变Phansalkar算法的参数，同时计算当前参数组对应的分割结果，并计算与分割结果中目标背景对应的原图像中的目标像素方差σ_obj和背景像素方差σ_bac，而后计算σ_obj和σ_bac的差的绝对值σ_dif,令σ_dif为目标函数，利用遗传算法优化找到使σ_dif达到最小值的参数组，相应分割结果为当前窗口分割结果。此算法能较好地分割灰度不均匀图像，且此方法无参数，可移植性好。

本发明解决其技术问题所采用的技术方案包括以下步骤：