GPU-based high-performance computing for radiation therapy

Physics in Medicine and Biology

Volumn 59, Issue 4, 2014, Pages

  Jia, Xun ^a   Ziegenhein, Peter ^b   Jiang, Steve B ^a  

^a UNIVERSITY OF TEXAS SOUTHWESTERN MEDICAL CENTER   (United States)

^b INSTITUTE OF CANCER RESEARCH   (United Kingdom)

Author keywords

graphics processing units; high performance computing; radiation therapy

Indexed keywords

ACCELERATION FACTORS; CLINICAL ENVIRONMENTS; COMPUTATIONAL POWER; GRAPHICS PROCESSING UNIT; HARDWARE STRUCTURES; HIGH-PERFORMANCE COMPUTING; LOW COSTS; PROGRAMMING MODELS;

COMPUTER GRAPHICS; PROGRAM PROCESSORS; RADIOTHERAPY;

COMPUTER GRAPHICS EQUIPMENT;

COMPUTER ANALYSIS; COMPUTER ASSISTED RADIOTHERAPY; COMPUTER GRAPHICS; DIAGNOSTIC IMAGING; EQUIPMENT; HUMAN; METHODOLOGY; REVIEW;

COMPUTER GRAPHICS; COMPUTING METHODOLOGIES; DIAGNOSTIC IMAGING; HUMANS; RADIOTHERAPY, COMPUTER-ASSISTED;

EID: 84893845759     PISSN: 00319155     EISSN: 13616560     Source Type: Journal    
DOI: 10.1088/0031-9155/59/4/R151     Document Type: Review

References (123)

1. Agostinelli S et al 2003 GEANT4 - a simulation toolkit Nucl. Instrum. Methods Phys. Res. A 506 250-303
2. Armbrust M et al 2010 A view of cloud computing Commun. ACM 53 50-8
3. Badal A and Badano A 2009 Accelerating Monte Carlo simulations of photon transport in a voxelized geometry using a massively parallel graphics processing unit Med. Phys. 36 4878-80
4. Bangert M, Ziegenhein P and Oelfke U 2012 Characterizing the combinatorial beam angle selection problem Phys. Med. Biol. 57 6707-23
5. Bell N and Garland M 2008 Efficient sparse matrix-vector multiplication on CUDA NVIDIA Technical Report NVR-2008-004
6. Bian J G et al 2010 Evaluation of sparse-view reconstruction from flat-panel-detector cone-beam CT Phys. Med. Biol. 55 6575-99
7. Cabello J, Gillam J E and Rafecas M 2012 High performance 3D PET reconstruction using spherical basis functions on a polar grid Int. J. Biomed. Imaging 2012 452910
8. Candes E J, Romberg J and Tao T 2006 Robust uncertainty principles: exact signa (Pubitemid 43174093)
9. Chen G T Y, Singh R P, Castro J R, Lyman J T and Quivey J M 1979 Treatment planning for heavy ion radiotherapy Int. J. Radiat. Oncol. Biol. Phys. 5 1809-19 (Pubitemid 10187528)
10. Chou C-Y, Chuo Y-Y, Hung Y and Wang W 2011 A fast forward projection using multithreads for multirays on GPUs in medical image reconstruction Med. Phys. 38 4052-65
11. de Greef M, Crezee J, van Eijk J C, Pool R and Bel A 2009 Accelerated ray tracing for radiotherapy dose calculations on a GPU Med. Phys. 36 4095-102
12. Dongarra J et al 1994 MPI - a message passing interface standard Int. J. Supercomput. Appl. High Perform. Comput. 8 165
13. Donoho D L 2006 Compressed sensing IEEE Trans. Inf. Theory 52 1289-306
14. Fang J, Varbanescu A L and Sips H 2011 A comprehensive performance comparison of CUDA and OpenCL 2011 Int. Conf. on Parallel Processing pp 216-25
15. Feldkamp L A, Davis L C and Kress J W 1984 Practical cone beam algorithm J. Opt. Soc. Am. A 1 612-9
16. Folkerts M, Jia X and Jiang S B 2012 A fast GPU-optimized DRR calculation algorithm for iterative CBCT reconstruction in preparation
17. Folkerts M et al 2010 MO-FF-A4-05: implementation and evaluation of various DRR algorithms on GPU Med. Phys. 37 3367
18. Fujimoto R, Kurihara T and Nagamine Y 2011 GPU-based fast pencil beam algorithm for proton therapy Phys. Med. Biol. 56 1319-28
19. Gao H, Li R, Lin Y and Xing L 2012 4D cone beam CT via spatiotemporal tensor framelet Med. Phys. 39 6943-6
20. Gautier Q et al 2013 Development of a GPU research platform for automatic treatment planning and adaptive radiotherapy re-planning Med. Phys. 40 534
21. Gifford K A, Horton J L, Wareing T A, Failla G and Mourtada F 2006 Comparison of a finite-element multigroup discrete-ordinates code with Monte Carlo for radiotherapy calculations Phys. Med. Biol. 51 2253-65
22. Golub G H and van Loan C F 1996 Matrix Computation (Baltimore, MD: John Hopkins University Press)
23. Gu X, Jelen U, Li J, Jia X and Jiang S B 2011a A GPU-based finite-size pencil beam algorithm with 3D-density correction for radiotherapy dose calculation Phys. Med. Biol. 56 3337-50
24. Gu X, Jia X and Jiang S B 2011b GPU-based fast gamma index calculation Phys. Med. Biol. 56 1431-41
25. Gu X et al 2009 GPU-based ultra-fast dose calculation using a finite size pencil beam model Phys. Med. Biol. 54 6287-97
26. Gu X et al 2010 Implementation and evaluation of various demons deformable image registration algorithms on a GPU Phys. Med. Biol. 55 207-19
27. Gu X et al 2013 A contour-guided deformable image registration algorithm for adaptive radiotherapy Phys. Med. Biol. 58 1889-901
28. Hissoiny S, Ozell B, Bouchard H and Despres P 2011a GPUMCD: a new GPU-oriented Monte Carlo dose calculation platform Med. Phys. 38 754-64
29. Hissoiny S, Ozell B and Despres P 2010 A convolution-superposition dose calculation engine for GPUs Med. Phys. 37 1029-37
30. Hissoiny S, Ozell B, Despres P and Carrier J F 2011b Validation of GPUMCD for low-energy brachytherapy seed dosimetry Med. Phys. 38 4101-7
31. Hissoiny S, Raaijmakers A J, Ozell B, Despres P and Raaymakers B W 2011c Fast dose calculation in magnetic fields with GPUMCD Phys. Med. Biol. 56 5119-29
32. Hofmann H G, Keck B, Rohkohl C and Hornegger J 2011 Comparing performance of many-core CPUs and GPUs for static and motion compensated reconstruction of C-arm CT data Med. Phys. 38 468-73
33. Jacobs F, Sundermann E, de Sutter B, Christiaens M and Lemahieu I 1998 A fast algorithm to calculate the exact radiological path through a pixel or voxel space J. Comput. Inf. Technol. CIT 6 89-94
34. Jacques R, Taylor R, Wong J and McNutt T 2010 Towards real-time radiation therapy: GPU accelerated superposition/convolution Comput. Methods Programs Biomed. 98 285-92
35. Jacques R, Wong J, Taylor R and McNutt T 2011 Real-time dose computation: GPU-accelerated source modeling and superposition/convolution Med. Phys. 38 294-305
36. Jacques R A, Taylor R H, Wong J W and McNutt T R 2008 Towards real-time radiation therapy: superposition/convolution at interactive rates Int. J. Radiat. Oncol. Biol. Phys. 72 S667
37. Jaffray D A, Drake D G, Moreau M, Martinez A A and Wong J W 1999 A radiographic and tomographic imaging system integrated into a medical linear accelerator for localization of bone and soft-tissue targets Int. J. Radiat. Oncol. Biol. Phys. 45 773-89 (Pubitemid 29462661)
38. Jaffray D A, Siewerdsen J H, Wong J W and Martinez A A 2002 Flat-panel cone-beam computed tomography for image-guided radiation therapy Int. J. Radiat. Oncol. Biol. Phys. 53 1337-49 (Pubitemid 34787872)
39. Jahnke L, Fleckenstein J, Wenz F and Hesser J 2012 GMC: a GPU implementation of a Monte Carlo dose calculation based on Geant4 Phys. Med. Biol. 57 1217-29
40. Jia X, Dong B, Lou Y and Jiang S B 2011a GPU-based iterative cone-beam CT reconstruction using tight frame regularization Phys. Med. Biol. 56 3787-807
41. Jia X, Gu X, Graves Y J, Folkerts M and Jiang S B 2011b GPU-based fast Monte Carlo simulation for radiotherapy dose calculation Phys. Med. Biol. 56 7017-31
42. Jia X, Lou Y, Li R, Song W Y and Jiang S B 2010a GPU-based fast cone beam CT reconstruction from undersampled and noisy projection data via total variation Med. Phys. 37 1757-60
43. Jia X, Schuemann J, Paganetti H and Jiang S B 2012a GPU-based fast Monte Carlo dose calculation for proton therapy Phys. Med. Biol. 57 7783-97
44. Jia X, Schumann J, Paganetti H and Jiang S B 2013 Development of gPMC v2.0, a GPU-based Monte Carlo dose calculation package for proton radiotherapy Med. Phys. 40 498
45. Jia X, Tian Z, Lou Y, Sonke J-J and Jiang S B 2012b Four-dimensional cone beam CT reconstruction and enhancement using a temporal nonlocal means method Med. Phys. 39 5592-602
46. Jia X, Yan H, Cervino L, Folkerts M and Jiang S B 2012c A GPU tool for efficient, accurate, and realistic simulation of cone beam CT projections Med. Phys. 39 7368-78
47. Jia X, Yan H, Gu X and Jiang S B 2012d Fast Monte Carlo simulation for patient-specific CT/CBCT imaging dose calculation Phys. Med. Biol. 57 577-90
48. Jia X et al 2010b Development of a GPU-based Monte Carlo dose calculation code for coupled electron-photon transport Phys. Med. Biol. 55 3077
49. Jia X et al 2011c GPU-based cone beam CT reconstruction via total variation regularization J. X-Ray Sci. Technol. 19 139
50. Ju T, Simpson T, Deasy J O and Low D A 2008 Geometric interpretation of the gamma dose distribution comparison technique: interpolation-free calculation Med. Phys. 35 879-87 (Pubitemid 351322669)
51. Kakimoto T, Dohi K, Shibata Y and Oguri K 2012 Performance comparison of GPU programming frameworks with the striped Smith-Waterman algorithm Comput. Archit. News 40 70-5
52. Kawrakow I 2000 Accurate condensed history Monte Carlo simulation of electron transport: I. EGSnrc, the new EGS4 version Med. Phys. 27 485-98 (Pubitemid 30171797)
53. Keyes RW, Romano C, Arnold D and Luan S 2010 Radiation therapy calculations using an on-demand virtual cluster via cloud computing http://arxiv.org/abs/1009.5282
54. Kim J, Li S, Zhao Y and Movsas B 2007 Real-time intensity-based deformable fusion on PC graphics hardware Int. Conf. on the Use of Computers in Radiotherapy
55. Kohno R et al 2003 Experimental evaluation of validity of simplified Monte Carlo method in proton dose calculations Phys. Med. Biol. 48 1277-88 (Pubitemid 36700464)
56. Kohno R et al 2011 Clinical implementation of a GPU-based simplified Monte Carlo method for a treatment planning system of proton beam therapy Phys. Med. Biol. 56 N287-94
57. Li J S, Shahine B, Fourkal E and Ma C M 2005 A particle track-repeating algorithm for proton beam dose calculation Phys. Med. Biol. 50 1001-10 (Pubitemid 40477320)
58. Li R et al 2010 Real-time volumetric image reconstruction and 3D tumor localization based on a single x-ray projection image for lung cancer radiotherapy Med. Phys. 37 2822-6
59. Li R et al 2011 3D tumor localization through real-time volumetric x-ray imaging for lung cancer radiotherapy Med. Phys. 38 2783-94
60. Low D A and Dempsey J F 2003 Evaluation of the gamma dose distribution comparison method Med. Phys. 30 2455-64 (Pubitemid 37129751)
61. Low D A, Harms W B, Mutic S and Purdy J A 1998 A technique for the quantitative evaluation of dose distributions Med. Phys. 25 656-61 (Pubitemid 28234885)
62. Lu W 2010 A non-voxel-based broad-beam (NVBB) framework for IMRT treatment planning Phys. Med. Biol. 55 7175
63. Lu W and Chen M 2010 Fluence-convolution broad-beam (FCBB) dose calculation Phys. Med. Biol. 55 7211-29
64. Luu J et al 2009 FPGA-based Monte Carlo computation of light absorption for photodynamic cancer therapy FCCM'09: Proc. 17th IEEE Symp. on Field Programmable Custom Computing Machines pp 157-64
65. Maier A et al 2011 Three-dimensional anisotropic adaptive filtering of projection data for noise reduction in cone beam CT Med. Phys. 38 5896-909
66. Men C, Jia X and Jiang S B 2010a GPU-based ultra-fast direct aperture optimization for online adaptive radiation therapy Phys. Med. Biol. 55 4309
67. Men C, Romeijn H E, Jia X and Jiang S B 2010b Ultrafast treatment plan optimization for volumetric modulated arc therapy (VMAT) Med. Phys. 37 5787-91
68. Men C, Romeijn H E, Tas Z C and Dempsey J F 2007 An exact approach to direct aperture optimization in IMRT treatment planning Phys. Med. Biol. 52 7333-52 (Pubitemid 350242843)
69. Men C et al 2009 GPU-based ultrafast IMRT plan optimization Phys. Med. Biol. 54 6565
70. Moore G E 1965 Cramming more components onto integrated circuits Electronics 38 114
71. Mueller K and Yagel R 2000 Rapid 3-D cone-beam reconstruction with the simultaneous algebraic reconstruction technique (SART) using 2-D texture mapping hardware IEEE Trans. Med. Imaging 19 1227-37 (Pubitemid 32870330)
72. Mueller K, Yagel R and Wheller J J 1999 Anti-aliased three-dimensional cone-beam reconstruction of low-contrast objects with algebraic methods IEEE Trans. Med. Imaging 18 519-37 (Pubitemid 30507530)
73. Noel P B et al 2010 GPU-based cone beam computed tomography Comput. Methods Programs Biomed. 98 271-7
74. Okitsu Y, Ino F and Hagihara K 2010 High-performance cone beam reconstruction using CUDA compatible GPUs Parallel Comput. 36 129-41
75. Pallipuram V K, Bhuiyan M and Smith M C 2012 A comparative study of GPU programming models and architectures using neural networks J. Supercomput. 61 673-718
76. Park J C et al 2011 Ultra-fast digital tomosynthesis reconstruction using general-purpose GPU programming for image-guided radiation therapy Technol. Cancer Res. Treat. 10 295-306
77. Park J C et al 2012 Fast compressed sensing-based CBCT reconstruction using Barzilai-Borwein formulation for application to on-line IGRT Med. Phys. 39 1207-17
78. Pasciak A S and Ford J R 2008 High-speed evaluation of track-structure Monte Carlo electron transport simulations Phys. Med. Biol. 53 5539-53
79. Peng F et al 2012 A new column-generation-based algorithm for VMAT treatment plan optimization Phys. Med. Biol. 57 4569-88
80. Perl J, Shin J, Schumann J, Faddegon B and Paganetti H 2012 TOPAS: an innovative proton Monte Carlo platform for research and clinical applications Med. Phys. at press
81. Persoon L C G G, Podesta M, van Elmpt W J C, Nijsten S M J J G and Verhaegen F 2011 A fast three-dimensional gamma evaluation using a GPU utilizing texture memory for on-the-fly interpolations Med. Phys. 38 4032-5
82. Podlozhnyuk V 2007 Image convolution with CUDA http://developer.download. nvidia.com/compute/DevZone/C/html/C/src/convolutionSeparable/doc/ convolutionSeparable.pdf
83. Poole C M, Cornelius I, Trapp J V and Langton C M 2012 Radiotherapy Monte Carlo simulation using cloud computing technology Australas. Phys. Eng. Sci. Med. 35 497-502
84. Pratx G and Xing L 2011a GPU computing in medical physics: a review Med. Phys. 38 2685-97
85. Pratx G and Xing L 2011b Monte Carlo simulation of photon migration in a cloud computing environment with MapReduce J. Biomed. Opt. 16 125003
86. Rangan K 2008 The cloud wars: $100+ billion at stake Technical Report (In Merrill Lynch 7 May 2008)
87. Salvat F, Fernández-Varea J M and Sempau J 2009 PENELOPE-2008: a code system for Monte Carlo simulation of electron and photon transport (Issy-les-Moulineaux: OECD-NEA)
88. Samant S S, Xia J Y, Muyan-Ozcelilk P and Owens J D 2008 High performance computing for deformable image registration: towards a new paradigm in adaptive radiotherapy Med. Phys. 35 3546-53
89. Shackleford J A, Kandasamy N and Sharp G C 2010 On developing B-spline registration algorithms for multi-core processors Phys. Med. Biol. 55 6329-51
90. Sharp G C, Kandasamy N, Singh H and Folkert M 2007 GPU-based streaming architectures for fast cone-beam CT image reconstruction and demons deformable registration Phys. Med. Biol. 52 5771-83 (Pubitemid 47485412)
91. Siddon R L 1985 Fast calculation of the exact radiological path for a 3-dimensional CT array Med. Phys. 12 252-5 (Pubitemid 15089015)
92. Sidky E Y and Pan X C 2008 Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization Phys. Med. Biol. 53 4777-807
93. Siegele L 2008 Let it rise: a special report on corporate IT Economist Magazine
94. Sonke J J, Zijp L, Remeijer P and van Herk M 2005 Respiratory correlated cone beam CT Med. Phys. 32 1176-86 (Pubitemid 40590074)
95. Spoerk J, Bergmann H, Wanschitz F, Dong S and Birkfellner W 2007 Fast DRR splat rendering using common consumer graphics hardware Med. Phys. 34 4302-8 (Pubitemid 350029652)
96. Su C-L, Chen P-Y, Lan C-C, Huang L-S and Wu K-H 2012 Overview and comparison of OpenCL and CUDA technology for GPGPU APCCAS'12: IEEE Asia Pacific Conf. on Circuits and Systems pp 448-51
97. Thirion J P 1998 Image matching as a diffusion process: an analogy with Maxwell's demons Med. Image Anal. 2 243-60
98. Tian Z, Jia X, Dong B, Lou Y and Jiang S B 2011a Low-dose 4DCT reconstruction via temporal nonlocal means Med. Phys. 38 1359-65
99. Tian Z, Jia X, Yuan K, Pan T and Jiang S B 2011b Low-dose CT reconstruction via edge-preserving total variation regularization Phys. Med. Biol. 56 5949-67
100. Townson R W et al 2013 GPU-based Monte Carlo radiotherapy dose calculation using phase-space sources Phys. Med. Biol. 58 4341-56
101. Tyagi N, Bose A and Chetty I J 2004 Implementation of the DPM Monte Carlo code on a parallel architecture for treatment planning applications Med. Phys. 31 2721-5 (Pubitemid 39287890)
102. Wang X, Yan H, Cervino L, Jiang S B and Jia X 2014 Towards the clinical implementation of iterative conebeam CT reconstruction for radiation therapy using a multi-GPU system Phys. Med. Biol. in preparation
103. Watt A and Watt M 1992 Advanced Animation and Rendering Techniques: Theory and Practice (Reading, MA: Addison-Wesley)
104. Webb S 1989 Optimization of conformal radiotherapy dose distributions by simulated annealing Phys. Med. Biol. 34 1349-70 (Pubitemid 19245988)
105. Weber R, Gothandaraman A, Hinde R J and Peterson G D 2011 Comparing hardware accelerators in scientific applications: a case study IEEE Trans. Parallel and Distrib. Syst. 22 58-68
106. Weng X et al 2003 A vectorized Monte Carlo code for radiotherapy treatment planning dose calculation Phys. Med. Biol. 48 N111-20 (Pubitemid 36503202)
107. Wu J, Kim M, Peters J, Chung H and Samant S S 2009 Evaluation of similarity measures for use in the intensity-based rigid 2D-3D registration for patient positioning in radiotherapy Med. Phys. 36 5391-403
108. Xiao K, Chen D Z, Hu X S and Zhou B 2012 Efficient implementation of the 3D-DDA ray traversal algorithm on GPU and its application in radiation dose calculation Med. Phys. 39 7619-25
109. Xu F and Mueller K 2004 Towards a unified framework for rapid 3D computed tomography on commodity GPUs IEEE Nucl. Sci. Symp. Conf. Rec. vol 4 pp 2757-9 (Pubitemid 40088805)
110. Xu F and Mueller K 2005 Accelerating popular tomographic reconstruction algorithms on commodity PC graphics hardware IEEE Trans. Nucl. Sci. 52 654-63 (Pubitemid 41172556)
111. Xu F and Mueller K 2007 Real-time 3D computed tomographic reconstruction using commodity graphics hardware Phys. Med. Biol. 52 3405-19 (Pubitemid 46830633)
112. Xu F et al 2010 On the efficiency of iterative ordered subset reconstruction algorithms for acceleration on GPUs Comput. Methods Programs Biomed. 98 261-70
113. Xu W and Mueller K 2009 Accelerating regularized iterative CT reconstruction on commodity graphics hardware (GPU) ISBI'09: IEEE Int. Symp. on Biomedical Imaging: from Nano to Macro pp 1287-90
114. Xu W and Mueller K 2010 Evaluating popular non-linear image processing filters for their use in regularized iterative CT NSS/MIC'10: IEEE Nuclear Science Symp. and Medical Imaging Conf. pp 2864-5
115. Yan D, Vicini F, Wong J and Martinez A 1997 Adaptive radiation therapy Phys. Med. Biol. 42 123-32 (Pubitemid 27055957)
116. Yan H, Cervino L, Jia X and Jiang S B 2012 A comprehensive study on the relationship between the image quality and imaging dose in low-dose cone beam CT Phys. Med. Biol. 57 2063-80
117. Yan H, Godfrey D J and Yin F-F 2008 Fast reconstruction of digital tomosynthesis using on-board images Med. Phys. 35 2162-9 (Pubitemid 351614830)
118. Yan H, Ren L, Godfrey D J and Yin F-F 2007 Accelerating reconstruction of reference digital tomosynthesis using graphics hardware Med. Phys. 34 3768-76 (Pubitemid 47482878)
119. Yepes P, Randeniya S, Taddei P J and Newhauser W D 2009 Monte Carlo fast dose calculator for proton radiotherapy: application to a voxelized geometry representing a patient with prostate cancer Phys. Med. Biol. 54 N21-8
120. Yepes P P, Mirkovic D and Taddei P J 2010 A GPU implementation of a track-repeating algorithm for proton radiotherapy dose calculations Phys. Med. Biol. 55 7107-20
121. Zhen X et al 2012 CT to cone-beam CT deformable registration with simultaneous intensity correction Phys. Med. Biol. 57 6807-26
122. Zhuge Y, Cao Y, Udupa J K and Miller R W 2011 Parallel fuzzy connected image segmentation on GPU Med. Phys. 38 4365-71
123. Ziegenhein P, Kamerling C P, Bangert M, Kunkel J and Oelfke U 2013 Performance-optimized clinical IMRT planning on modern CPUs Phys. Med. Biol. 58 3705-15