论文摘要
计算机断层成像技术(简称CT)在医学、航空、安检等领域中有着广泛的应用。在成像原理上,工业CT和医学CT都是基于Radon变换。经过多年的发展,CT的扫描方式已经由最初的二维平行束扫描发展到三维锥束扫描。三维锥束CT与二维CT相比,具有射线利用率高、纵向分辨率高、扫描时间短等优点。图像重建方法是CT技术中较为关键的部分。图像重建有两类基本方法:解析法和迭代法。解析法的优点是能够在较短的时间内重建出高质量的图像,但在投影数据不完全或高噪声的情况下,该算法不再适用。而迭代法在投影旋转分度数较少或高噪声的情况下仍能重建出较高质量的图像。迭代重建法的缺点是速度慢,重建图像的时间与投影的数目、重建图像的大小、迭代的次数相关,但随着计算机性能的不断提高,迭代法所需重建时间较长的缺点将降为次要矛盾。本文对工业CT(ICT)锥束扫描方式下的迭代重建算法进行了研究。在工业CT应用中,有时会遇到缺失投影数据的情况,这种情况下应用迭代重建算法能够重建出高质量的图像,本文比较研究了不同迭代算法在不完全投影数据情况下的重建结果。此外由于锥束CT采用高密度面阵探测器,相邻探测单元之间存在射线散射引起的串扰,所以投影数据中存在大量噪声。因此我们对不同迭代算法的抗噪性也进行了对比,并针对代数重建(ART)算法在同样情况下容易受噪声影响,通过选择松弛因子、投影次序和迭代的初值,改善了ART算法的抗噪性。在代数迭代算法中,联合代数迭代算法(SART)是一种重要的图像重建方法。该算法的主要思想是通过一个角度下的所有投影测量值和计算值的偏差来修正像素值,通过迭代计算使测量值和计算值总的误差变小,但该算法对重建图像域来说是无约束的,因此迭代到一定次数后,待重建图像中的噪声会出现逐渐增加的情况。为了解决这一问题,针对三维锥束CT情况,研究了一种增加惩罚项联合代数迭代方法。计算机仿真实验表明:与无惩罚项的联合代数迭代算法相比,该方法在锥束CT图像重建中,能够在抑制噪声的同时提高重建图像的收敛速度。在统计迭代重建算法中,有序子集期望值最大重建算法(OSEM)能够在较短的时间内重建出高质量的图像,当投影数据中含有噪声时,子集大小的选取将会影响到重建图像的质量和收敛速度。针对三维锥束CT,论文中研究了一种基于图像分块的变子集OSEM重建方法,该方法将图像空间分割成等大小的图像块,然后在迭代过程中,对于不同的图像块用变化的子集进行图像重建。仿真实验表明:在锥束CT图像重建中,该方法能够在抑制噪声的同时提高重建图像的收敛速度。