To compensate for artifacts in reconstructed images of industrial computed tomography (ICT), Grangeat type exact cone beam reconstruction algorithm was studied. Grangeat algorithm based on flat detector was implemented by the steps of computing first partial derivatives of cone-beam projection data along rows and lines respectively, first partial derivatives of Radon data, back-projection and so on. Due to the missed data of a circular scan orbit, the incomplete Radon space was analyzed. The nearest-neighbor interpolation, linear interpolation, 3-points quadratic interpolation and one-dimensional linear regression interpolation methods were used to estimate missed data respectively, and the effects of different methods on the image quality were examined. The compute simulation results prove the validity of proposed methods.
Feldkamp, Davis L A, Kress J. Practical cone-beam algorithm[J]. Journal of the Optical Society of America, 1984, 1(1):1612-619
Grangeat P. Mathematical framework of cone beam 3D reconstruction via the first derivative of the Radon transform Louis A K,Natterer F. Mathematical Methods in Tomography. Herman GT: Springer Verlag, 1991:66-97
Lee Seung Wook, Cho Gyuseong, Wang Ge. Artifacts associated with implementation of the Grangeat formula[J].Med Phys.2002, 29(12):2871-2880
����,��־ǿ,����,��.Բ������ؽ��������ǹ�ϵ���о�[J].CT������Ӧ���о�, 2003,12(3):9-16 Zeng Kai, Chen Zhiqiang, Zhang Li, et al. The study of relationship between circular orbit cone beam reconstruction error and cone angle[J]. CT Theory and Applications, 2003,12(3):9-16(in Chinese)
Chen Zhiqiang, Li Liang, Kang Kejun, et al. New developments of exact cone-beam CT reconstruction algorithms[J]. CT Theory and Applications, 2005,14(3):65-70
Marr R B, Chen C N, lauterbur P C��On two approaches to 3D reconstruction in NMR zeugmatography Natterer F. Mathematical Aspects of Computerized Tomography. Herman GT: Springer Verlag, 1981:225-240