论文部分内容阅读
计算机断层成像(Computed Tomography,简称CT)是通过无损方式获取物体内部结构信息的一种技术手段,已被广泛地应用到医疗诊断、工业无损探伤、航空航天等领域。在数据采集的速度、空间分辨率以及射线的利用率等方面,锥束CT明显优于二维平行束和扇束CT,是现代CT成像领域内一个活跃的研究方向。在医学成像领域内,数字减影血管照影(Digital Subtraction Angiography,简称DSA)设备能够采集锥束投影数据并可对高对比度物体进行重建。近年来,由于平板检测器的出现,利用DSA设备生成低对比度物体断层图像的技术引起了广泛的关注。锥束投影的精确重建面临算法复杂、数值计算稳定性差及计算量剧增的难题,是锥束CT发展的一个瓶颈,也是制约DSA设备锥束断层重建质量的一个主要因素。目前,针对适合DSA设备扫描轨迹的锥束重建算法的研究工作较少,本文基于滤波反投影锥束重建算法对此进行一些探索性的研究,主要工作包括以下几个方面:DSA生成断层图像首选的扫描轨迹是满足短扫描重建条件的圆弧轨迹,目前阶段主要采用短扫描加权的FDK算法。FDK算法属于近似锥束重建,在较大锥角情况下其重建结果会产生锥束伪像。为此,本文首先研究重排类型的FDK衍生算法,并将重建质量最好的FDK-SLANT算法应用到无重排步骤的锥束投影重建中,改善了FDK算法存在的低对比度衰减伪像。该类型算法的特点是通过一些固定斜平面与检测器平面的交线确定不同投影视角上投影数据的滤波线。受该算法的启发,在斜平面绕转轴旋转过程中形成近似平面空间中,重新分析了FDK算法的推导过程,进而提出了一种利用转角及其增量关系改进经典FDK算法的新方法,并且将改进方法进一步应用到重排类型的T-FDK算法中。改进算法能够在保持FDK类型算法执行效率的同时,有效地解决了在锥角较大情况下FDK算法存在的低对比度衰减伪像,显著扩大了低对比度物体的锥向重建范围。改进后的T-FDK算法能够获得较好的重建结果,然而该算法需要将锥束投影数据重排为三维平行束投影。超短扫描轨迹FDK类型算法通过建立锥束投影与三维平行束投影之间的重排计算关系,避免了T-FDK算法的重排步骤。依据重建公式是否存在轨迹参数偏导数的计算项,超短扫描圆弧轨迹算法可分为两类,其中不含有轨迹参数偏导数计算项的重建算法数值计算精度高,然而其重建结果具有较小的锥向重建范围。在三维Radon逆变换分解的基础上,本文推导出一个新的偏导数计算与轨迹参数无关的锥束重建公式,使两类算法获得一致的重建结果,并且进一步将转角增量关系的改进措施推广引入到超短扫描轨迹的两类重建算法中,扩大了低对比度物体的锥向重建范围。当重建密度差较大的物体时,圆弧轨迹的FDK类型算法在较小锥角的切片图像上也会产生严重的涂污伪像。针对该问题,本文研究了Katsevich重建理论在圆弧轨迹上的应用,其重建计算的关键在于构造因子的确定。主要研究了基于几何图形法确定构造因子大小的方法,通过对检测器平面的有效划分,探讨了构造因子在检测器平面上的分布规律,并由此推导了满足投影数据采样精度的Katsevich类型重建公式,从而有效避免了FDK类型算法所存在的涂污伪像。不同于FDK类型算法,在同一个投影视角上,Katsevich类型算法需要沿三个方向对投影数据进行滤波,其中两个方向为非水平方向,而重建公式中基于Hilbert变换的滤波属于非局部运算,同时DSA检测器采集到的投影数据也存在轴向及横向截断。为此,论文进一步研究了利用保形的分段三次Hermite插值函数外插值截断的投影数据,从而允许Katsevich类型算法重建截断的投影数据。单圆弧轨迹所采集的锥束投影,理论上不能提供完备的Radon数据,目前已有的重建算法尚不能完善解决该轨迹上的伪像问题。为了确保锥束投影数据的完备性,本文研圆弧加直线扫描轨迹上Katsevich类型的锥束重建算法。Katsevich本人已讨论了Katsevich理论应用到该轨迹上的两个基本问题:构造因子的取值和PI线决定的有效反投影区域。本文对此作了进一步分析,研究了圆弧加直线轨迹锥束精确重建时轨迹扫描区间所需要满足的条件;补充了圆弧与直线轨迹交点位置处构造因子的分析,进而修正Katsevich所获得的重建公式;利用边界过渡的权重函数改善反投影区域边界投影数据的阶跃变化;推导了平板检测器下Katsevich类型的重建公式。由于在重建计算中含对轨迹参数的偏导数计算,当重建物体小且密布分布时,容易产生计算误差,影响重建图像的质量。本文利用超短扫描路径上重建算法的研究结果以及直线轨迹的特点,进一步推导出不含有轨迹偏导数的重建公式,改善了重建图像的质量。