工业ct三维图重建与分割算法研究.pdf

万方数据   Study on Industrial CT 3D Image Reconstruction and Segmentation Algorithm A Thesis Submitted to Chongqing University in Partial Fulfillment of the Requirement for the Doctor’s Degree of Engineering By Zou Yongning Supervised by Prof. Li Jianwei Specialty: Instrument Science and Technology College of Optoelectronic Engineering of Chongqing University, Chongqing, China December, 2014 万方数据 万方数据 中文摘要 摘 要 工业 CT 是使用X 射线来生成被测物体外部和内部二维或三维表示的技术。 工业CT 在许多领域已被用于零部件的内部检测，主要应用有缺陷检测、故障分析、 测量、装配分析等。锥束CT(CBCT)通过平板探测器采集投影数据，能够实现快速 体积成像，并且具有各向同性的高分辨率和较高的X 射线使用效率。锥束 CT 主 要有圆形轨迹锥束 CT 和螺旋轨迹锥束 CT 两种形式。本文主要研究了工业锥束 CT 重建算法，并提出了改进方法，同时研究了三维CT 图像缺陷分割方法。 由于制作工艺问题，实际平板探测器获取的图像不一定理想，可能存在非均 匀性，即使通过采集求平均也不能消去，需要对数据进行校正才可以消除。平板 探测器的非均匀性存在于整个平面上，平板探测器越大，非均匀性问题就越突出。 围绕这个问题，本文提出了非均匀性校正的锥束CT 重建算法。首先，在重建前对 投影数据进行非均匀性校正，然后采用校正后的数据进行重建。采用仿真数据和 实际采集的数据进行验证，结果表明提出的算法能显著改善重建图像的质量。 因为存在X 射线散射和电子噪声干扰，平板探测器采集的投影数据可能被许 多噪声污染。如果不去噪就进行图像重建，那么重建图像的信噪比就会降低。而 且螺旋轨迹锥束CT 的Katsevich 重建算法需要计算投影数据的导数，如果直接对 含噪的投影数据进行求导，会放大噪声在投影数据中的影响，重建图像的质量会 严重下降，所以含噪投影数据需要提前进行滤波处理。中值滤波能够较好地去除 信号或图像中的椒盐噪声，但是该滤波方法可能会破坏图像的一些重要细节。本 文提出采用冲击滤波器改进的算法来去噪，然后使用Katsevich 算法对滤波后的投 影数据进行图像重建。实验证实，采用该方法重建图像的质量明显好于中值滤波 的重建图像。 在锥束CT 发展过程中，采集数据的尺寸越来越大，因此锥束CT 的图像重建 成为一项非常耗时的计算任务。GPU 集群是高性能并行计算设备，利用多个计算 机资源去同时处理计算任务，使计算时间极大地降低。本文采用GPU 集群技术来 加速图像重建同时维护图像质量