论文部分内容阅读
CT (Computed Tomography),即计算机断层成像技术,是医学和工业中的一种优质无损诊断技术,其本质是利用投影数据重建物体断层图像。高分辨率重建、大物体重建和三维重建的工程需求使得投影数据量激增,使重建出的CT图像的数据量也激增,从而重建速度就成为一个瓶颈。精密部件的检测对图像重建提出了更高的精度要求,如何直接重建出满足工程目标的CT图像也是CT应用中亟待解决的关键问题。本文致力于重建精度的提高、重建速度的提高以及目标重建等三个方面的研究。论文首先研究了滤波反投影(FBP)算法的离散实现公式,利用信号处理理论及定积分的定义,推导出了它应该包含投影采样间隔和角度采样间隔两个参数,并作为整个公式的因子。利用定积分数值计算方法,并根据反投影的特点,提出了基于辛普森公式的精确反投影方法。将插值方法应用于反投影的过程,仿真结果表明线性插值比样条插值效果更好,并给出了定性解释。在分析斜变滤波器特点的基础上,提出了一种斜变滤波器的设计原则,分析了频域加窗和时域加窗对重建的影响,确定了斜变滤波器单位脉冲响应的长度选择方法和一种理想斜变滤波器单位脉冲响应的中心值计算方法。仿真结果表明,频域中指数窗性能最优而时域中矩形窗性能最优,单位脉冲响应的长度为(?)-1时,CT图像无空气伪像。根据投影信号和滤波器单位脉冲响应都为长序列的特点,提出了基于算术傅里叶变换(AFT)的滤波过程的加速方法,以及基于快速哈达玛变换(FHT)的滤波过程的加速方法。仿真结果表明,AFT滤波算法比FFT滤波算法快30%,而FHT滤波算法比FFT滤波算法快近50%。本文研究了反投影过程的三种加速方法:对每个投影预插值到更精细的间隔,反投影时采用最快的零阶插值方法;采用被动式射线驱动法反投影;以空间换时间的策略来预存储三角函数值而避免反投影时复杂的三角运算。在此基础上,提出了一种组合反投影加速方法。仿真实验表明,比传统方法快4倍多。在“特征重建”概念的基础上,提出了“目标重建”的概念,并设计了目标重建的三种实现方法:近似法、空间域精确法、频域精确法。实际实验表明:近似法重建效果一般,而精确法重建效果良好,且可以有效避免非白噪声对CT图像后处理的影响。论文从成像理论的角度,以解析式图像重建算法为研究对象,研究了加速方法、精度提高方法以及目标重建的方法,研究结果对解析式图像重建算法具有很好的推动作用。