Fast parallel algorithm for three-dimensional distance-driven model in iterative computed tomography reconstruction

doi:10.1088/1674-1056/24/2/028703

›› 2015, Vol. 24 ›› Issue (2): 28703-028703.doi: 10.1088/1674-1056/24/2/028703

• INTERDISCIPLINARY PHYSICS AND RELATED AREAS OF SCIENCE AND TECHNOLOGY • 上一篇下一篇

Fast parallel algorithm for three-dimensional distance-driven model in iterative computed tomography reconstruction

陈建林, 李磊, 王林元, 蔡爱龙, 席晓琦, 张瀚铭, 李建新, 闫镔

National Digital Switching System Engineering & Technological R & D Center, Zhengzhou 450002, China

收稿日期:2014-06-25 修回日期:2014-09-11 出版日期:2015-02-05 发布日期:2015-02-05
基金资助:
Projected supported by the National High Technology Research and Development Program of China (Grant No. 2012AA011603) and the National Natural Science Foundation of China (Grant No. 61372172).

Fast parallel algorithm for three-dimensional distance-driven model in iterative computed tomography reconstruction

Chen Jian-Lin (陈建林), Li Lei (李磊), Wang Lin-Yuan (王林元), Cai Ai-Long (蔡爱龙), Xi Xiao-Qi (席晓琦), Zhang Han-Ming (张瀚铭), Li Jian-Xin (李建新), Yan Bin (闫镔)

National Digital Switching System Engineering & Technological R & D Center, Zhengzhou 450002, China

Received:2014-06-25 Revised:2014-09-11 Online:2015-02-05 Published:2015-02-05
Contact: Yan Bin E-mail:tom.yan@gmail.com
Supported by:
Projected supported by the National High Technology Research and Development Program of China (Grant No. 2012AA011603) and the National Natural Science Foundation of China (Grant No. 61372172).

摘要/Abstract

摘要： The projection matrix model is used to describe the physical relationship between reconstructed object and projection. Such a model has a strong influence on projection and backprojection, two vital operations in iterative computed tomographic reconstruction. The distance-driven model (DDM) is a state-of-the-art technology that simulates forward and back projections. This model has a low computational complexity and a relatively high spatial resolution; however, it includes only a few methods in a parallel operation with a matched model scheme. This study introduces a fast and parallelizable algorithm to improve the traditional DDM for computing the parallel projection and backprojection operations. Our proposed model has been implemented on a GPU (graphic processing unit) platform and has achieved satisfactory computational efficiency with no approximation. The runtime for the projection and backprojection operations with our model is approximately 4.5 s and 10.5 s per loop, respectively, with an image size of 256×256×256 and 360 projections with a size of 512×512. We compare several general algorithms that have been proposed for maximizing GPU efficiency by using the unmatched projection/backprojection models in a parallel computation. The imaging resolution is not sacrificed and remains accurate during computed tomographic reconstruction.

关键词: computed tomography, iterative reconstruction, parallelizable algorithm, distance-driven model

Abstract: The projection matrix model is used to describe the physical relationship between reconstructed object and projection. Such a model has a strong influence on projection and backprojection, two vital operations in iterative computed tomographic reconstruction. The distance-driven model (DDM) is a state-of-the-art technology that simulates forward and back projections. This model has a low computational complexity and a relatively high spatial resolution; however, it includes only a few methods in a parallel operation with a matched model scheme. This study introduces a fast and parallelizable algorithm to improve the traditional DDM for computing the parallel projection and backprojection operations. Our proposed model has been implemented on a GPU (graphic processing unit) platform and has achieved satisfactory computational efficiency with no approximation. The runtime for the projection and backprojection operations with our model is approximately 4.5 s and 10.5 s per loop, respectively, with an image size of 256×256×256 and 360 projections with a size of 512×512. We compare several general algorithms that have been proposed for maximizing GPU efficiency by using the unmatched projection/backprojection models in a parallel computation. The imaging resolution is not sacrificed and remains accurate during computed tomographic reconstruction.

Key words: computed tomography, iterative reconstruction, parallelizable algorithm, distance-driven model

中图分类号: (X-ray imaging)

87.59.-e

07.85.-m (X- and γ-ray instruments) 87.57.Q- (Computed tomography)

陈建林, 李磊, 王林元, 蔡爱龙, 席晓琦, 张瀚铭, 李建新, 闫镔. Fast parallel algorithm for three-dimensional distance-driven model in iterative computed tomography reconstruction[J]. , 2015, 24(2): 28703-028703.

Chen Jian-Lin (陈建林), Li Lei (李磊), Wang Lin-Yuan (王林元), Cai Ai-Long (蔡爱龙), Xi Xiao-Qi (席晓琦), Zhang Han-Ming (张瀚铭), Li Jian-Xin (李建新), Yan Bin (闫镔). Fast parallel algorithm for three-dimensional distance-driven model in iterative computed tomography reconstruction[J]. Chin. Phys. B, 2015, 24(2): 28703-028703.

参考文献 25

[1]	Yang F Q, Zhang D H, Huang K D, Wang K and Xu Z 2014 Acta Phys. Sin. 63 058701 (in Chinese)
[2]	Wang L Y, Zhang H M, Cai A L, Yan B, Li L and Hu G E 2013 Acta Phys. Sin. 62 198701 (in Chinese)
[3]	Bian J, Siewerdsen J, Han X, Sidky E Y, Prince J L, Pelizzari C A and Pan X C 2010 Phys. Med. Biol. 55 6575
[4]	Zhang H M, Wang L Y, Yan B, Li L, Xi X Q and Lu L Z 2013 Chin. Phys. B 22 078701
[5]	Wang L Y, Li L, Yan B, Jiang C S, Wang H Y and Bao S L 2010 Chin. Phys. B 19 088106
[6]	Jia X and Lou Y F 2011 Phys. Med. Biol. 47 77
[7]	De Man B and Fessler J A 2010 Biomedical Mathematics: Promising Directions in Imaging, Therapy Planning, and Inverse Problems (Medical Physics Publishing) pp. 113-140
[8]	Siddon R L 1985 Med. Phys. 12 252
[9]	Joseph P M 1983 IEEE Trans. Med. Imaging 1 192
[10]	Wang G, Vannier M W and Cheng P C 1999 Microsc. Microanal. 5 58
[11]	De Man B and Basu S 2002 IEEE Nuclear Science Symposium Medical Imaging Conference, Norfolk, IEEE, p. 1477
[12]	Lewitt R M 1990 J. Opt. Soc. Am. 7 1834
[13]	Yu H and Wang G 2012 J. X-Ray Sci. Technol. 20 229
[14]	Zhang O, Connor Y and Fessler J A 2006 IEEE Trans. Med. Imaging 25 582
[15]	Vazquez F, Garzon E M and Fernandez 2011 The Comput. J. 54 1861
[16]	Herman G T 1980 Image Reconstruction from Projections: the Fundamentals of Computerized Tomography (New York: Academic Press) p. 20
[17]	Zeng G L and Gullberg G T 2000 IEEE Trans. Med. Imaging 19 548
[18]	Nguyen V G and Lee S J 2012 Opt. Eng. 51 093203
[19]	Long Y, Fessler J A and Balter J M 2010 IEEE Trans. Med. Imaging 29 1839
[20]	De Man B and Basu S 2004 Phys. Med. Biol. 49 2463
[21]	Andersen A H and Kak A C 1984 Ulwason. Imag. 6 81
[22]	Gao H 2012 Med. Phys. 39 7110
[23]	Gordon R, Bender R and Herman G T 1970 J. Theor. Biol. 29 471
[24]	Wang G and Jiang M 2004 J. X-Ray Sci. Technol. 12 169
[25]	Daniel G, Ulrich H, Ralf S, Elmar S and Ulrich S 2009 SIAM J. Imag. Sci. 31 4204

Fast parallel algorithm for three-dimensional distance-driven model in iterative computed tomography reconstruction

Fast parallel algorithm for three-dimensional distance-driven model in iterative computed tomography reconstruction

RichHTML

PDF (PC)

赞

可视化

摘要/Abstract

引用本文

使用本文

参考文献 25

相关文章 13

Metrics

本文评价

推荐阅读 0

[1]	卫晨希, 吴朝, 魏文彬, 鲍园, 骆荣辉, 王磊, 刘刚, 田扬超. Single-shot grating-based x-ray differential phase contrast imaging with a modified analyzer grating[J]. 中国物理B, 2017, 26(10): 108701-108701.
[2]	张瀚铭, 王林元, 李磊, 闫镔, 蔡爱龙, 胡国恩. Optimization-based image reconstruction in x-ray computed tomography by sparsity exploitation of local continuity and nonlocal spatial self-similarity[J]. 中国物理B, 2016, 25(7): 78701-078701.
[3]	金朝, 张瀚铭, 闫镔, 李磊, 王林元, 蔡爱龙. Novel Fourier-based iterative reconstruction for sparse fan projection using alternating direction total variation minimization[J]. 中国物理B, 2016, 25(3): 38701-038701.
[4]	李盼云, 张凯, 黄万霞, 袁清习, 王研, 鞠在强, 吴自玉, 朱佩平. Cosine fitting radiography and computed tomography[J]. 中国物理B, 2015, 24(6): 68704-068704.
[5]	程任翔, 陶超, 刘晓峻. Enhancement of photoacoustic tomography in the tissue with speed-of-sound variance using ultrasound computed tomography[J]. 中国物理B, 2015, 24(11): 114301-114301.
[6]	鲍园, 王研, 高昆, 王志立, 朱佩平, 吴自玉. Investigation of noise properties in grating-based x-ray phase tomography with reverse projection method[J]. 中国物理B, 2015, 24(10): 108702-108702.
[7]	黄魁东, 张华, 史仪凯, 张亮, 徐哲. Scatter correction method for cone-beam CT based on interlacing-slit scan[J]. , 2014, 23(9): 98106-098106.
[8]	胡勇, 李静超, 沈明武, 史向阳. Formation of multifunctional Fe₃O₄/Au composite nanoparticles for dual-mode MR/CT imaging applications[J]. 中国物理B, 2014, 23(7): 78704-078704.
[9]	张凯, 李德娥, 洪友丽, 朱佩平, 袁清习, 黄万霞, 高昆, 周红章, 吴自玉. Penetrating view of nano-structures in Aleochara verna spermatheca and flagellum by hard X-ray microscopy[J]. 中国物理B, 2013, 22(7): 76801-076801.
[10]	韩玉, 闫镔, 李磊, 宇超群, 李建新, 包尚联. Multiple helical scans and the reconstruction of over FOV-sized objects in cone-beam CT[J]. 中国物理B, 2012, 21(6): 68701-068701.
[11]	汪先超, 闫镔, 李磊, 胡国恩 . Cone-beam local reconstruction based on Radon inversion transform[J]. 中国物理B, 2012, 21(11): 118702-118702.
[12]	贾鹏翔, 张峰, 闫镔, 包尚联. Scattering correction method for panel detector based cone beam computed tomography system[J]. 中国物理B, 2010, 19(8): 87802-087802.
[13]	王林元, 李磊, 闫镔, 江成顺, 王浩宇, 包尚联. An algorithm for computed tomography image reconstruction from limited-view projections[J]. 中国物理B, 2010, 19(8): 88106-088106.