针对小动物的锥束CT重建的研究及验证

基于经典3D Shepp-Logan模型,应用射线追踪和射线衰减定律得到锥束CT的模拟仿真数据.利用FDK三维重建算法实现锥束CT图像重建,采用仿真图像对重建算法进行了验证,通过搭建的锥束CT系统和自制仿体进行了层析成像和图像重建.结果表明:算法可较为准确地重建出原物体的空间分布和密度关系,针对仿真数据的中心断层重建图像的空间分辨率可达到8 lp/mm.     

一种利用小体积偏置探测器的新型锥束X射线CT系统

在传统的CT系统中, 系统的硬件成本和计算量都是非常巨大的. 提出了一种利用偏移放置的面阵探测器的锥束CT系统, 在这种系统中X射线束仅仅覆盖被扫描物体的一半体积, 投影数据在探测器方向上是截断的, 探测器尺寸和投影数据量都减少为传统CT系统的一半. 在这种新型扫描方式下, 现有的CT重建算法都不能处理. 因此, 提出了一种BPF形式的直接反投影重建方法. 该方法不需要对投影数据重排, 直接反投影滤波重建出最终的图像. 因此, 该算法在数学上更简洁, 计算速度更快. 最后, 数值模拟实验和真实CT系统实验结果均验证了该系统和重建算法能够获得高质量的CT图像.     

一种圆轨迹锥束CT中截断投影数据的高效重建算法

本文基于数据重排方法, 提出了T-BPF (Tent-BPF)算法, 该算法先将锥束投影数据重排成平行投影数据, 然后使用一种推导的BPF型算法重建重排后的平行投影数据. T-BPF算法将原BPF算法反投影中变化的角度积分限变成固定的, 反投影中各层循环之间没有了相关性, 这意味着T-BPF算法较原BPF算法具有更好的可并行性. 实验结果显示: 使用GPU对256³的Shepp-Logan体模的图像重建进行并行加速, T-BPF算法在保证重建质量的前提下, 加速比达到了1036, 较原BPF算法有很大提升. T-BPF算法为截断投影数据的3D图像快速重建提供了方法. 关键词： X射线光学 CT 图像重建 GPU     

光栅成像的锥束相位CT重建算法

X射线相位衬度成像对弱吸收物质有着独有的高密度分辨率,在医学、生物学、材料学等领域显示出良好的应用前景.但是其中的干涉成像法、衍射增强成像法和同轴成像法具有一定局限性,难以被广泛应用.光栅成像克服了以往相位衬度成像的缺陷,也使相位衬度成像向锥束成像发展成为可能.本文致力于锥束相位X射线计算机断层摄影术(CT)重建算法的研究,根据基于光栅的锥束相位成像的特点,利用吸收CT中FDK (Feldkamp-Davis-Kress)重建算法的思想,发展出适用于锥束相位成像的CT重建方法.该方法为滤波反投影类型,以相位一阶导数像为投影值,直接重建物体的相位项.通过仿真验证了算法的正确性.     

螺旋锥束计算机断层成像倾斜扇束反投影滤波局部重建算法

螺旋锥束计算机断层成像(CT)作为常用的临床诊断工具,如何尽可能地减少其辐射剂量是热点研究领域之一.局部成像利用准直器减小射线直照区域,能够有效降低CT辐射剂量.然而,局部成像会造成投影数据横向截断,产生局部重建问题.现有螺旋反投影滤波(BPF)算法只能实现局部曲面重建,难以实现局部体区域重建.在圆轨迹扇束BPF算法的基础上,通过加权修正和坐标扩展,提出了螺旋锥束CT倾斜扇束反投影滤波(TFB-BPF)重建算法.该算法把重建区域按层划分,对各层构建倾斜的扇形束几何,并使用经过加权修正的TFB-BPF算法逐层进行重建.该算法最大的优势是滤波线(与原始螺旋锥束BPF算法中PI线等价)在二维平面内选择,算法更加简洁高效,并且能够应用于局部体区域的重建.实验结果表明,算法能够有效重建物体局部体区域,并且重建图像质量较好,没有明显的截断伪影.     

锥束射线RT扫描大视场三维CT成像方法研究

受探测器尺寸和重建锥角的限制,基于单圆轨道扫描和FDK算法的三维CT成像方法难以对大物体进行三维CT成像。提出了一种基于平板探测器,采用锥束RT扫描方式的成像方法,推导了其三维CT成像算法。首先把对回转物体进行三次锥束平移透照的投影数据重新排成为更大范围的倾斜平行投影数据,然后经过加权滤波和反投影重建,获得物体的三维层析结果。仿真结果表明,该方法在重构大物体图像时,获得了满意的图像质量。水平方向视场是单圆轨道FDK算法视场的2.88倍,在垂直方向减少了锥角过大形成的伪影,提高了物体的重构精度。     

基于最小区域的锥束ART算法快速图像重建

针对锥束ART算法重建速度慢的问题,考虑到工业CT重建目标尺寸差异较大的特点,提出了一种基于最小重建区域的快速三维图像重建方法。由不同视角下的锥束投影构建最小区域包络,对最小区域包络进行膨胀处理,以消除重建目标边界附近所出现的伪影,在此基础上,给出了射线穿过最小重建区域的体素索引及权因子的计算。该方法能够根据重建目标的尺寸自适应地确定最小重建区域,从而最大限度地减少了计算量。实验结果表明,该方法不仅大幅度提高了锥束ART算法的重建速度,而且有效地提高了重建质量。     

一种基于Radon逆变换的半覆盖螺旋锥束CT重建算法

半覆盖螺旋锥束计算机断层成像能够扩展传统螺旋锥束计算机断层成像的成像视野, 实现小面板探测器成像超视野物体. 但是, 各角度下的投影都存在数据截断, 会使重建结果中产生截断伪影, 降低图像质量. 本文提出了一种基于Radon逆变换的半覆盖螺旋锥束重建算法, 该算法在滤波时先使用局部算子, 再使用全局算子. 局部算子不受数据截断的影响, 并且降低了运算后的值由于数据截断所造成的不连续性, 因此, 减少了随后全局运算产生的截断误差. 仿真和实际实验结果均验证了本文算法的有效性. 和现有算法的对比也表明, 本文算法针对半覆盖螺旋锥束投影具有更强的截断伪影抑制能力, 能够有效提高重建图像的质量.     

X射线CT转台单侧两次螺旋扫描单层重排重建算法

锥束计算机断层成像(CT)在临床医学及工业无损检测等诸多领域获得了越来越广泛的应用。但是,受面阵探测器尺寸等硬件条件限制,锥束CT的成像视野有限,难以满足实际应用中对大尺寸物体成像的需求。为了扩展锥束CT的视野,提出了锥束CT转台单侧两次螺旋扫描单层重排重建算法。该算法在转台的同一侧进行两次螺旋扫描,通过统一成像几何条件,将两组螺旋锥形束投影重排为多层平行束投影,利用平行束投影的对称性去除数据的横向截断,得到一组完全覆盖物体横截面的平行束投影数据,再通过滤波反投影方式得到重建图像。实验结果表明,该方法的成像视野扩展率可达2.56,且成像质量与传统全覆盖算法相当。     

改进的大锥角锥束CT圆轨迹反投影滤波算法北大核心CSCD

锥束CT圆轨迹反投影滤波(BPF)算法能够有效处理数据截断问题,在感兴趣区域成像研究中具有重要的地位。但是当锥角增大时,重建图像的顶端和底端会出现严重的伪影,从而影响图像质量。为解决该问题,从扇束圆轨迹BPF算法出发,重新推导了锥束圆轨迹BPF算法,通过对锥束伪影成因分析,提出了一种基于锥角大小和轴向距离大小的权重函数,修正了重建过程中差分反投影图像求解部分,达到抑制伪影的目的。实验结果表明,改进的大锥角锥束CT圆轨迹BPF算法能够在锥角增大时有效地抑制锥角效应,提高重建图像质量。