Computed tomography using level set method and algebraic reconstruction technique

Abstract

In this paper, a novel reconstruction technique based on level set method and algebraic reconstruction technique is proposed for multiphase flow computed tomography (CT) system. The curvature-driven noise reduction method is inserted into the conventional iteration procedure of algebraic reconstruction technique to improve the image quality and convergence speed with limited projection data. By evolving the image as a set of iso-intensity contours after each updation, the sufficient number of iterations for acceptable results is reduced by 80%–90%, while the image quality is enhanced obviously. Quantitative evaluation of image quality is given by using both relative image error and correlation coefficient. The resultant images can be utilized to detect flow regimes for monitoring industrial multiphase flow. Laboratory results demonstrate the feasibility of the proposed method. Phantoms of four typical flow regimes can be reconstructed from few-view projection data efficiently, and the corresponding image errors and correlation coefficients are acceptable for the cases tested in this paper.

Keywords

computed tomography / multiphase flow / image reconstruction

References

[1] International Atomic Energy Agency. Industrial Process Gamma Tomography[R]. 2008, Austria: IAEA.

[2] Pan X. C., Sidky E. Y., Vannier M. Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction[J]. Inverse Problems, 2009, 25(12): 123009

[3] Zhang Lifeng. Research on Exciting-Measuring Modes and Image Reconstruction Algorithms for Electrical Tomography[D]. 2010, Tianjin: School of Electrical Engineering and Automation, Tianjin University.

[4] Wu C., Cheng Y., Ding Y., et al. A novel Xray computed tomography method for fast measurement of multiphase flow [J]. Chemical Engineering Science, 2007, 62(12): 4325-4335.

[5] Maad R., Johansen G. A. Experimental analysis of highspeed gamma-ray tomography performance[J]. Measurement Science and Technology, 2008, 19(8): 1-10.

[6] Wang H., Xue Qian. Algebraic reconstruction algorithm based on expanding projection matrix and blind deconvolution[J]. Journal of Tianjin University, 2011, 44(5): 385-390.

[7] Zhuang Tiange. Principles and Algorithms of CT[M]. 1992, Shanghai: Shanghai Jiaotong University Press.

[8] Li Z., Cong P., Wu Haifeng. Research of ART method in CT image reconstruction[J]. Nuclear Electronics and Detection Technology, 2005, 25(2): 184-186.

[9] Atkinson C., Soria J. An efficient simultaneous reconstruction technique for tomographic particle image velocimetry[J]. Experiments in Fluids, 2009, 47(4/5): 553-568.

[10] Osher S., Sethian J. A. Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton-Jacobi formulations[J]. Journal of Computational Physics, 1988, 79, 12-49.

[11] Sethian J. A. Level Set Methods and Fast Marching Methods[M]. 1999, New York: Cambridge University Press.

[12] Malladi R., Sethian J. A. A unified approach to noise removal, image enhancement, and shape recovery[J]. IEEE Transactions on Image Processing, 1996, 5(11): 1554-1568.

[13] Johansen G. A. Nuclear tomography methods in industry [J]. Nuclear Physics A, 2005, 752, 696c-705c.