论文部分内容阅读
CT (Computed Tomography)是一种先进的无损检测技术,在其很多应用中,经常会遇到感兴趣区域重建、大物体重建、高分辨率重建等实际应用需求,这些往往都会产生局部重建问题。截断投影数据的局部重建是CT重建领域的一个难点,传统的经典重建算法用于局部重建会在重建图像中产生严重的截断伪影。对于实际CT系统的大规模数据量的重建,现有的局部重建算法在重建效率、并行性质、可实现性方面存在诸多的缺陷和不足,很难满足实际CT系统快速准确重建的需求。本文针对不同的具体应用,在保证重建质量的前提下,研究并提出了重建效率高、并行性质好且可方便地用于实际CT系统的局部重建算法。主要研究内容和创新点如下:1.对于PI线段上没有发生数据截断的局部重建问题,提出了基于数据重排方法的反投影滤波(Backprojection Filtration, BPF)型局部重建算法(T-BPF算法)。该算法先将锥束投影数据重排成帐篷状平行投影数据,然后推导了一种BPF型算法对重排后的帐篷状平行投影数据进行重建。T-BPF算法在保证重建质量的前提下,重建效率比原BPF算法有所提高,而且反投影的运算中各层循环之间没有了相关性,适合使用并行计算进行加速。数字仿真和真实数据重建的实验结果证实了该算法的正确性及其优势。2.对于投影数据任意截断的一般局部重建问题,提出了基于Radon逆变换的抑制截断伪影的近似局部重建算法。该方法是将2D Radon逆变换推广用于3D锥束CT的局部重建,其滤波过程分为两步:第一步是对投影数据求导,该步骤是局部的,投影数据截断情况下的结果仍然是准确的;第二步是对求导后的投影数据进行Hilbert滤波。该方法具有与ATRACT (Approximate Truncation Resistant Algorithm for Computed Tomography)相当的抑制截断伪影的能力,且其重建效率是ATRACT的两倍以上。数字仿真和真实数据的重建实验证实了该方法的有效性。新提出的基于Radon逆变换的局部重建算法为圆轨迹锥束CT对重建图像精度要求不高的一般局部重建问题提供了一种简单且有效的方法。3.针对扁平状特殊结构物体的局部重建问题,当PI线段上没有发生数据截断时,提出了一种改进的BPF算法,用于超短扫描的扁平状物体的局部重建。该方法可以大大减少数据采集和成像的时间,提高实际CT系统的使用寿命。该方法是在原BPF算法基础上的一种改进,即固定反投影中的投影角度积分限。通过推导发现如果扁平物体的厚度小于0.0349R(R是圆轨迹的扫描半径),改进算法的误差是可以忽略的。当被重建物体的厚度超出可重建的范围,可以适当增大扫描半径来进行重建。数字仿真和真实数据的实验结果显示该方法是扁平状物体局部重建的一个很好的选择,因为它不仅具有与其它重建算法相当的重建质量,而且可以大大减少成像时间,提高CT系统的使用寿命。4.针对扁平状特殊结构物体的局部重建问题,当投影数据任意截断时,提出了一种基于FDK算法的扁平状物体的局部重建方法。通过研究发现如果扁平物体的厚度和局部区域在水平方向上的长度满足一定条件,简洁、高效的FDK算法可以处理该局部重建问题。通过数字仿真实验得到了一个可以比较准确重建扁平状物体局部区域的经验性条件,即:局部区域在水平方向上的长度大于该物体在水平方向上的长度的1/8,物体厚度小于该物体在水平方向上的长度的1/13。数字仿真和真实数据的重建结果证实了在这个条件下FDK算法可以很好地实现扁平状物体的局部重建。5.针对重建物体超出扫描视野的CT图像重建问题,考虑一种视野只覆盖一半物体的扫描结构,对于该扫描方式下任意投影角度都发生截断的投影数据,结合BPF算法和反投影图像可以通过不同投影角度下的投影数据得到这个发现,提出了一种重建半覆盖扫描下截断投影数据的算法。该算法不需要数据重排这个操作,且对于每个点的反投影图像只需要180度范围内投影数据,大大提高了算法的重建效率,重建图像的空间分辨率也有所提高。数字仿真和真实数据的实验验证了算法的正确性及其优点。最后,对本文工作进行了总结,并对CT图像局部重建算法的发展进行了展望,指出了未来的研究方向和构想。