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

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

Multiple helical scans and the reconstruction of over FOV-sized objects in cone-beam CT

In cone-beam computed tomography(CBCT),there are often cases where the size of the specimen is larger than the field of view(FOV)(referred to as over FOV-sized(OFS)).To acquire the complete projection data for OFS objects,some scan modes have been developed for long objects and short but over-wide objects.However,these modes still cannot meet the requirements for both longitudinally long and transversely wide objects.In this paper,we propose a multiple helical scan mode and a corresponding reconstruction algorithm for both longitudinally long and transversely wide objects.The simulation results show that our model can deal with the problem and that the results are acceptable,while the OFS object is twice as long compared with the FOV in the same latitude.     

Novel Fourier-based iterative reconstruction for sparse fan projection using alternating direction total variation minimization

Sparse-view x-ray computed tomography(CT) imaging is an interesting topic in CT field and can efficiently decrease radiation dose. Compared with spatial reconstruction, a Fourier-based algorithm has advantages in reconstruction speed and memory usage. A novel Fourier-based iterative reconstruction technique that utilizes non-uniform fast Fourier transform(NUFFT) is presented in this work along with advanced total variation(TV) regularization for a fan sparse-view CT. The proposition of a selective matrix contributes to improve reconstruction quality. The new method employs the NUFFT and its adjoin to iterate back and forth between the Fourier and image space. The performance of the proposed algorithm is demonstrated through a series of digital simulations and experimental phantom studies. Results of the proposed algorithm are compared with those of existing TV-regularized techniques based on compressed sensing method, as well as basic algebraic reconstruction technique. Compared with the existing TV-regularized techniques, the proposed Fourier-based technique significantly improves convergence rate and reduces memory allocation, respectively.     

一种迭代的锥束计算机层析成像系统几何全参数标定方法EI北大核心CSCD

计算机层析成像(CL)系统对于扁平状物体的检测具有独特的优势,几何参数标定是CL系统获得高质量图像的重要步骤。现有CL系统几何参数标定方法仅能标定部分几何参数,难以一次完成所有几何参数的标定。借鉴计算机断层成像(CT)系统几何参数标定方法,提出了一种高精度的迭代标定方法,利用简单体模快速标定CL系统的所有几何参数。对CT系统几何参数标定方法用于CL系统的有效性进行分析,确定对CT系统几何参数标定方法敏感的CL系统几何参数;根据CL系统投影几何关系,构建新的反映CL实际系统与理想系统间误差的非线性最小二乘目标函数,利用迭代法不断优化敏感参数和其他关键参数。实验结果表明,本文方法能够有效标定CL系统所有几何参数,敏感参数以及受敏感参数影响的其他几何参数标定精度均有明显提升。利用已标定的几何参数重建三维Shepp-Logan体模和印刷电路板体模,重建图像中没有几何伪影,验证了算法的有效性。     

Scattering correction method for panel detector based cone beam computed tomography system

<正>A scattering correction method for a panel detector based cone beam computed tomography system is presented. First,the x-ray spectrum of the system is acquired by using the Monte Carlo simulation method.Secondly,scattered photon distribution is calculated and stored as correction matrixes by using the Monte Carlo simulation method according to scanned objects and computed tomography system specialties.Thirdly,scattered photons are removed from projection data by correction matrixes.A comparison of reconstruction image between before and after scattering correction demonstrates that the scattering correction method is effective for the panel detector based cone beam computed tomography system.     

基于加权TV正则化的X射线CT系统能谱估计方法

在CT硬化伪影校正、双能CT图像重建以及CT辐射剂量计算等实际应用中,X射线的能谱信息具有重要的作用。然而,由于透射测量方程组系数矩阵的病态性、X光子的统计涨落和噪声的干扰,使得EM等能谱估计方法难以获得较精细的能谱刻画。针对该问题,提出一种基于加权TV正则化的X射线CT系统能谱估计方法。首先采用能谱能量范围内带有不同K-edge的材料作为体模以降低投影测量方程之间的相关性。然后,利用CT成像系统的几何参数来获取投影测量数据所对应的准确透射长度信息来减小投影方程的测量误差。最后,综合利用透射衰减测量数据的保真性、轫致辐射能谱部分的连续性、特征辐射能谱部分的离散性、能谱的非负性和归一性以及平均有效衰减系数等信息,使用加权TV正则化方法建立目标函数,利用正则化理论中的L曲线准则通过黄金分割参数搜索策略求得最优正则化参数及对应的能谱估计结果。分别通过不同统计波动模型下的仿真能谱和实际测量数据对算法性能进行了验证,结果表明该方法与EM等方法相比,不受初始能谱的影响,有效提高了能谱估计结果的稳定性和准确性。     

结合透射CT的康普顿背散射图像重建

在康普顿散射成像(CST)技术中可以结合透射成像重建出衰减系数来消除散射重建的非线性,但这样得到的投影矩阵带有误差。而CST重建问题的不适定性对噪声和投影矩阵的误差非常敏感,重建结果会有较大误差。针对此问题,基于压缩感知理论提出了一种新的CST重建算法。新方法将图像重建问题归结为一个图像的全变分(TV)最小化问题,并使用收敛速度较快的基于交替方向法的Split-Bregman方法进行求解。在仿真实验中,通过与代数重建技术(ART)进行比较,在测量数据充足和测量数据不足两种情况下,本文算法都具有更好的重建质量,证明了所提算法在重建精度和抗噪性能方面的优势。     

Fast local reconstruction by selective backprojection for low dose in dental computed tomography

The high radiation dose in computed tomography (CT) scans increases the lifetime risk of cancer, which becomes a major clinical concern. The backprojection-filtration (BPF) algorithm could reduce the radiation dose by reconstructing the images from truncated data in a short scan. In a dental CT, it could reduce the radiation dose for the teeth by using the projection acquired in a short scan, and could avoid irradiation to the other part by using truncated projection. However, the limit of integration for backprojection varies per PI-line, resulting in low calculation efficiency and poor parallel performance. Recently, a tent BPF has been proposed to improve the calculation efficiency by rearranging the projection. However, the memory-consuming data rebinning process is included. Accordingly, the selective BPF (S-BPF) algorithm is proposed in this paper. In this algorithm, the derivative of the projection is backprojected to the points whose x coordinate is less than that of the source focal spot to obtain the differentiated backprojection. The finite Hilbert inverse is then applied to each PI-line segment. S-BPF avoids the influence of the variable limit of integration by selective backprojection without additional time cost or memory cost. The simulation experiment and the real experiment demonstrated the higher reconstruction efficiency of S-BPF.     

A modified interval subdividing based geometric calibration method for interior tomography

The interior tomography is commonly met in practice, whereas the self-calibration method for geometric parameters remains far from explored. To determine the geometry of interior tomography, a modified interval subdividing based method, which was originally developed by Tan et al.,[11]was presented in this paper. For the self-calibration method, it is necessary to obtain the reconstructed image with only geometric artifacts. Therefore, truncation artifacts reduction is a key problem for the self-calibration method of an interior tomography. In the method, an interior reconstruction algorithm instead of the Feldkamp–Davis–Kress(FDK) algorithm was employed for truncation artifact reduction. Moreover, the concept of a minimum interval was defined as the stop criterion of subdividing to ensure the geometric parameters are determined nicely. The results of numerical simulation demonstrated that our method could provide a solution to the selfcalibration for interior tomography while the original interval subdividing based method could not. Furthermore, real data experiment results showed that our method could significantly suppress geometric artifacts and obtain high quality images for interior tomography with less imaging cost and faster speed compared with the traditional geometric calibration method with a dedicated calibration phantom.     

Cone-beam local reconstruction based on a Radon inversion transformation

The local reconstruction from truncated projection data is one area of interest in image reconstruction for computed tomography(CT),which creates the possibility for dose reduction.In this paper,a filtered-backprojection(FBP) algorithm based on the Radon inversion transform is presented to deal with the three-dimensional(3D) local reconstruction in the circular geometry.The algorithm achieves the data filtering in two steps.The first step is the derivative of projections,which acts locally on the data and can thus be carried out accurately even in the presence of data truncation.The second step is the nonlocal Hilbert filtering.The numerical simulations and the real data reconstructions have been conducted to validate the new reconstruction algorithm.Compared with the approximate truncation resistant algorithm for computed tomography(ATRACT),not only it has a comparable ability to restrain truncation artifacts,but also its reconstruction efficiency is improved.It is about twice as fast as that of the ATRACT.Therefore,this work provides a simple and efficient approach for the approximate reconstruction from truncated projections in the circular cone-beam CT.