GPU-based ultra fast IMRT plan optimization

Abstract

The widespread adoption of on-board volumetric imaging in cancer radiotherapy has stimulated research efforts to develop online adaptive radiotherapy techniques to handle the inter-fraction variation of the patient's geometry. Such efforts face major technical challenges to perform treatment planning in real-time. To overcome this challenge, we are developing a supercomputing online re-planning environment (SCORE) at the University of California San Diego (UCSD). As part of the SCORE project, this paper presents our work on the implementation of an intensity modulated radiation therapy (IMRT) optimization algorithm on graphics processing units (GPUs). We adopt a penalty-based quadratic optimization model, which is solved by using a gradient projection method with Armijo's line search rule. Our optimization algorithm has been implemented in CUDA for parallel GPU computing as well as in C for serial CPU computing for comparison purpose. A prostate IMRT case with various beamlet and voxel sizes was used to evaluate our implementation. On an NVIDIA Tesla C1060 GPU card, we have achieved speedup factors of 20-40 without losing accuracy, compared to the results from an Intel Xeon 2.27 GHz CPU. For a specific 9-field prostate IMRT case with 5x5 mm^2 beamlet size and 2.5x2.5x2.5 mm^3 voxel size, our GPU implementation takes only 2.8 seconds to generate an optimal IMRT plan. Our work has therefore solved a major problem in developing online re-planning technologies for adaptive radiotherapy.

submitted to Physics in Medicine and Biology

References in corpus (1)

• GPU-based ultra fast dose calculation using a finite pencil beam model

Cited by in corpus (20)

• GPU-based Low Dose CT Reconstruction via Edge-preserving Total Variation Regularization
• Implementation and evaluation of various demons deformable image registration algorithms on GPU
• GPU-based fast Monte Carlo simulation for radiotherapy dose calculation
• GPU-based Iterative Cone Beam CT Reconstruction Using Tight Frame Regularization
• A GPU Tool for Efficient, Accurate, and Realistic Simulation of Cone Beam CT Projections
• Development of a GPU-based Monte Carlo dose calculation code for coupled electron-photon transport
• Multi-GPU Accelerated Multi-Spin Monte Carlo Simulations of the 2D Ising Model
• Multicriteria VMAT optimization
• Four-dimensional Cone Beam CT Reconstruction and Enhancement using a Temporal Non-Local Means Method
• GPU-based ultra-fast direct aperture optimization for online adaptive radiation therapy
• Ultra-fast treatment plan optimization for volumetric modulated arc therapy (VMAT)
• Fast Monte Carlo Simulation for Patient-specific CT/CBCT Imaging Dose Calculation
• GPU-based fast gamma index calcuation
• Beam Orientation Optimization for Intensity Modulated Radiation Therapy using Adaptive l1 Minimization
• A GPU-based finite-size pencil beam algorithm with 3D-density correction for radiotherapy dose calculation
• A GPU-based multi-criteria optimization algorithm for HDR brachytherapy
• High-dimensional Automated Radiation Therapy Treatment Planning via Bayesian Optimization
• Massively parallelized replica-exchange simulations of polymers on GPUs
• The fixed-point iteration method for IMRT optimization with truncated dose deposition coefficient matrix
• Conditional Gradient Methods for Convex Optimization with General Affine and Nonlinear Constraints