In circular full-cover cone-beam computed tomography (CT), the field-of-view (FOV) is limited by the width of planar detector, resulting in low imaging efficiency for large object. The FOV can be doubled by half-cover scanning, in which the back-projection filtration (BPF) algorithm based on the concept of PI-line is the best choice for image reconstruction. However, the integral intervals of different PI-lines are unequal in the BPF algorithm, leading to heavy communication consuming and calculation. As a result, the reconstruction efficiency by use of the BPF algorithm is low. In this paper, an efficient image reconstruction strategy based on the BPF algorithm for flat object is proposed. With the method, we demonstrate that the inequality of integral interval of PI-line can be ignored in the discrete implementation of the BPF algorithm when the thickness of flat object is less than 2R sin(2π/Np) (R is the scanning radius and Np is the number of uniform sampled projections in a full circle). Compared with the original BPF algorithm for half-cover scanning, our method has two major advantages: the first one is that the outer loop is the sample angle while the inner loop is the PI-line, which reduces the communication consuming for computer significantly;the second one is that the derivative of projection, back-projection and inverse Hilbert transform along the PI-line can be computed using parallel computing techniques readily. The results of numerical simulation and real data experiment indicate that the computational efficiency of the proposed method is 5.6 times that for original BPF algorithm and the reconstruction errors of the two methods are comparable.% 全覆盖圆轨迹扫描的成像视野受探测器宽度限制,对于大物体的成像效率较低.半覆盖扫描可以将成像视野扩展近1倍,图像重建首推使用反投影滤波型算法.反投影滤波型算法按PI线重建,各PI线积分区间的不一致性导致通信和计算消耗大,影响重建效率.针对半覆盖成像中扁平形状物体的重建问题,提出了一种改进的反投影滤波型算法,且证明了当扁平物体的厚度小于2Rsin(2π/Np)(R为扫描半径, Np为圆扫描一周均匀采集的投影数量)时, PI线积分区间的不一致性在数值计算过程中的误差是可以忽略的.改进后的算法相比原半覆盖反投影滤波算法具有两个明显的优势:一是数值计算过程中角度循环移至PI线循环之外,算法的通信需求显著降低;二是投影数据求导、反投影和沿PI线滤波三个步骤均能够并行计算,算法的并行性得到增强.数值仿真与实际数据的实验结果表明,本文算法与原半覆盖反投影滤波算法的重建精度相当,但计算效率提高了4.6倍.
展开▼