CT解析重建方法进展：从圆轨迹扫描到多源直线扫描成像.docxVIP

CT解析重建方法进展：从圆轨迹扫描

到多源直线扫描成像

摘要：从1971年Hounsfield完成世界上首次X射线计算机断层成像（CT）扫描以来，CT技术在扫描方式上经历了5代发展，并在医疗、安检等领域深刻地影响着人们的生活。在扫描方式不断进步的同时，重建方法也在不断发展。与迭代算法相比，解析算法具有重建速度快、误差分析简单以及占用运算资源小等特点，成为CT系统中常用的成像算法。本文以CT扫描方式的发展为脉络介绍解析重建方法的发展，从经典的圆轨迹和螺旋轨迹扫描出发，到非标准轨迹扫描，并介绍当前有关多源直线CT的研究进展。最后，我们展望CT系统的发展在未来所面临的挑战。

关键词：X射线计算机断层成像；CT解析重建；非标准轨迹扫描；直线轨迹扫描；分布式光源CT

X射线计算机断层成像（computedtomography，CT）利用不同角度下X射线穿透被检物体获得的投影数据，通过算法重建出物体断层图像或三维图像，解决了X射线透视成像过程中存在的混叠和遮挡问题，可以实现空间各点数值的精确求解，是X射线成像技术的一次重要飞跃。从1971年Hounsfield完成世界上首次CT扫描至今[1]，经过50年的发展，CT已经成为一种主流的非入侵式成像手段，给医学诊断、工业检测和安全检查等领域带来了革命性的影响。伴随成像需求的变化，CT技术经历了5代发展（图1），从最初的“笔束平移-旋转”扫描，发展到现在的多层螺旋锥束扫描、电子束扫描等。随着CT技术在三维成像、医学导航、快速安检等领域的深入应用，低剂量成像、定量成像、快速成像等新需求也给CT成像技术提出了更高的要求，但以圆轨迹扫描和螺旋轨迹扫描为代表的标准扫描轨迹CT系统在更快的成像速度、更大的通道尺寸等方面的发展遇到瓶颈。近年来，非标准轨迹扫描及分布式多光源成像大大拓展了CT系统的应用场景，在解决大通道成像、高速成像等问题上表现出了巨大的潜力。

CT成像模式的发展带动了CT图像重建算法的发展。传统CT重建算法可以分为解析重建算法与迭代重建算法两大类。与迭代算法相比，解析算法具有重建速度快、误差分析简单以及占用运算资源小等特点，是当前CT系统中应用的主流算法。本文依托CT扫描模式

的发展，介绍解析重建方法的发展，从经典的圆轨迹和螺旋轨迹扫描出发，到非标准轨迹扫描成像，并介绍了当前分布式多源直线CT的研究进展。最后，展望了CT成像系统在未来面临的挑战。

图1各代CT系统扫描示意图

Fig.1TheillustrationsofCTsystemsacrossthegenerations

1标准轨迹CT扫描模式与重建方法

在CT扫描投影和数据重建过程中，主要涉及投影数据和图像数据在物体空间、拉东空

间和傅里叶空间之间的相互变换（图2）。

拉东变换（Radontransform）和中心切片定理（centralslicetheorem）是解析重建的理论基础。1917年，奥地利数学家拉东[2]首次提出，由物体某断面各方向的物理特性参数线积分值，可计算该断面的物理特性参数分布。在二维成像中，CT重建就是物体空间与拉东空间相互变换的过程。二维图像的中心切片定理揭示了傅里叶

图2三种空间之间的关系（二维成像）

Fig.2Therelationshipbetweenthethree

domains(2DImaging)

空间是拉东空间与物体空间之间的桥梁，定理指出：二维函数f(x,y)在角度θ下的平行束

投影p(t,θ)关于探测器方向t的一维傅里叶变换(ξ,θ)给定了f(x,y)的二维傅里叶空间

f?(ξx,ξy)中与探测器方向平行的过原点的一个切片，其数学表达为：

(ξ,θ)=f?(ξcosθ,ξsinθ)。（1）

本节主要针对圆轨迹和螺旋轨迹两种标准轨迹扫描模式及其解析重建算法进行介绍。

1.1圆轨迹CT扫描与重建方法

在圆轨迹CT扫描中，光源与探测器绕扫描中心做圆轨迹运动，这是最为经典的扫描模式之一。常用的圆轨迹CT解析重建算法有滤波反投影算法（平行束/扇形束扫描）与FDK（Feldkamp-Davis-Kress）算法（锥束扫描）等。

滤波反投影算法[3]（filteredbackprojection，FBP）因形式简洁、易于实现、重建速度较快而得到广泛应用。平行束扫描投影的FBP重建过程可表示为：

其中，内层积分对应投影的滤波步骤（h