Performance of a hybrid MC dose algorithm for IMRT optimization dose evaluation

Adaptive proton therapy

Radiation therapy treatments are typically planned based on a single image set, assuming that the patient's anatomy and its position relative to the delivery system remains constant during the course of treatment. Similarly, the prescription dose assumes constant biological dose-response over the treatment course. However, variations can and do occur on multiple time scales. For treatment sites with significant intra-fractional motion, geometric changes happen over seconds or minutes, while biological considerations change over days or weeks. At an intermediate timescale, geometric changes occur between daily treatment fractions. Adaptive radiation therapy is applied to consider changes in patient anatomy during the course of fractionated treatment delivery. While traditionally adaptation has been done off-line with replanning based on new CT images, online treatment adaptation based on on-board imaging has gained momentum in recent years due to advanced imaging techniques combined with treatment delivery systems. Adaptation is particularly important in proton therapy where small changes in patient anatomy can lead to significant dose perturbations due to the dose conformality and finite range of proton beams. This review summarizes the current state-of-the-art of on-line adaptive proton therapy and identifies areas requiring further research.

NeuralDAO: Incorporating neural network generated dose into direct aperture optimization for end-to-end IMRT planning

PURPOSE

Thecurrent practice in intensity-modulated radiation therapy (IMRT) planning almost always includes different dose calculation strategies for plan optimization and final dose verification. The accurate Monte Carlo (MC) dose algorithm is considered to be time-consuming for the optimization. Thus a fast, simplified dose algorithm is used in the optimization. The significant differences between the optimized dose and the delivered dose lead to tediously planning loops and potentially suboptimal solutions. This work aims to develop an IMRT optimization algorithm to minimize the dose discrepancy so that the delivered dose can be optimized in a holistic, end-to-end manner.

METHODS

The proposed algorithm, namely NeuralDAO, integrates a neural dose network into the column generation (CG) direct aperture optimization (DAO) formulation for step-and-shoot IMRT planning. The neural dose network is designed and trained to produce doses of MC-level accuracy within few milliseconds. Its differentiability is fully exploited to compute gradients for identifying potential aperture shapes. A prototype of NeuralDAO was developed in PyTorch and available to the public. Five lung patient cases have been studied. Dosimetric accuracy was compared with the MC dose. Plan quality and time were compared with a state-of-the-art (SoA) dose-correct algorithm. Statistical analysis was performed by Wilcoxon signed-rank test.

RESULTS

The average gamma passing rate at 2 mm/2% is 99.7% between the optimized and delivered doses. The convergence process produced by NeuralDAO is virtually identical to that produced by an MC-based DAO. The average dose calculation time is 12.1 ms for an aperture on GPU. One session of optimization took 10-36 min. Compared with the SoA, better conformity index and homogeneity index were observed for the target. The esophagus was significantly spared. Significant reductions were observed for the replanning number and the planning time.

CONCLUSIONS

A new DAO algorithm based on the neural dose network has been developed. The results suggest that this algorithm minimizes the discrepancy between the optimized and delivered doses, which offers a promising approach to reduce the time and effort required in IMRT planning. This work demonstrates the possibility of applying the neural network in IMRT optimization. It is of great potential to extend this algorithm to other treatment modalities.

Performance of a hybrid Monte Carlo-Pencil Beam dose algorithm for proton therapy inverse planning

PURPOSE

Analytical algorithms have a limited accuracy when modeling very heterogeneous tumor sites. This work addresses the performance of a hybrid dose optimizer that combines both Monte Carlo (MC) and pencil beam (PB) dose engines to get the best trade-off between speed and accuracy for proton therapy plans.

METHODS

The hybrid algorithm calculates the optimal spot weights (w) by means of an iterative optimization process where the dose at each iteration is computed by using a precomputed dose influence matrix based on the conventional PB plus a correction term c obtained from a MC simulation. Updates of c can be triggered as often as necessary by calling the MC dose engine with the last corrected values of w as input. In order to analyze the performance of the hybrid algorithm against dose calculation errors, it was applied to a simplistic water phantom for which several test cases with different errors were simulated, including proton range uncertainties. Afterwards, the algorithm was used in three clinical cases (prostate, lung, and brain) and benchmarked against full MC-based optimization. The influence of different stopping criteria in the final results was also investigated.

RESULTS

The hybrid algorithm achieved excellent results provided that the estimated range in a homogeneous material is the same for the two dose engines involved, i.e., PB and MC. For the three patient cases, the hybrid plans were clinically equivalent to those obtained with full MC-based optimization. Only a single update of c was needed in the hybrid algorithm to fulfill the clinical dose constraints, which represents an extra computation time to obtain c that ranged from 1 (brain) to 4 min (lung) with respect to the conventional PB-based optimization, and an estimated average gain factor of 14 with respect to full MC-based optimization.

CONCLUSION

The hybrid algorithm provides an improved trade-off between accuracy and speed. This algorithm can be immediately considered as an option for improving dose calculation accuracy of commercial analytical treatment planning systems, without a significant increase in the computation time (≪5 min) with respect to current PB-based optimization.

Patient-specific CT dosimetry calculation: a feasibility study

Current estimation of radiation dose from computed tomography (CT) scans on patients has relied on the measurement of Computed Tomography Dose Index (CTDI) in standard cylindrical phantoms, and calculations based on mathematical representations of “standard man”. Radiation dose to both adult and pediatric patients from a CT scan has been a concern, as noted in recent reports. The purpose of this study was to investigate the feasibility of adapting a radiation treatment planning system (RTPS) to provide patient‐specific CT dosimetry. A radiation treatment planning system was modified to calculate patient‐specific CT dose distributions, which can be represented by dose at specific points within an organ of interest, as well as organ dose‐volumes (after image segmentation) for a GE Light Speed Ultra Plus CT scanner. The RTPS calculation algorithm is based on a semi‐empirical, measured correction‐based algorithm, which has been well established in the radiotherapy community. Digital representations of the physical phantoms (virtual phantom) were acquired with the GE CT scanner in axial mode. Thermoluminescent dosimeter (TLDs) measurements in pediatric anthropomorphic phantoms were utilized to validate the dose at specific points within organs of interest relative to RTPS calculations and Monte Carlo simulations of the same virtual phantoms (digital representation). Congruence of the calculated and measured point doses for the same physical anthropomorphic phantom geometry was used to verify the feasibility of the method. The RTPS algorithm can be extended to calculate the organ dose by calculating a dose distribution point‐by‐point for a designated volume. Electron Gamma Shower (EGSnrc) codes for radiation transport calculations developed by National Research Council of Canada (NRCC) were utilized to perform the Monte Carlo (MC) simulation. In general, the RTPS and MC dose calculations are within 10% of the TLD measurements for the infant and child chest scans. With respect to the dose comparisons for the head, the RTPS dose calculations are slightly higher (10%–20%) than the TLD measurements, while the MC results were within 10% of the TLD measurements. The advantage of the algebraic dose calculation engine of the RTPS is a substantially reduced computation time (minutes vs. days) relative to Monte Carlo calculations, as well as providing patient‐specific dose estimation. It also provides the basis for a more elaborate reporting of dosimetric results, such as patient specific organ dose volumes after image segmentation.

PACS numbers: 87.55.D‐, 87.57.Q‐, 87.53.Bn, 87.55.K‐

Evaluation of dose prediction errors and optimization convergence errors of deliverable-based head-and-neck IMRT plans computed with a superposition/convolution dose algorithm

The purpose of this study is to evaluate dose prediction errors (DPEs) and optimization convergence errors (OCEs) resulting from use of a superposition/convolution dose calculation algorithm in deliverable intensity-modulated radiation therapy (IMRT) optimization for head-and-neck (HN) patients. Thirteen HN IMRT patient plans were retrospectively reoptimized. The IMRT optimization was performed in three sequential steps: (1) fast optimization in which an initial nondeliverable IMRT solution was achieved and then converted to multileaf collimator (MLC) leaf sequences; (2) mixed deliverable optimization that used a Monte Carlo (MC) algorithm to account for the incident photon fluence modulation by the MLC, whereas a superposition/convolution (SC) dose calculation algorithm was utilized for the patient dose calculations; and (3) MC deliverable-based optimization in which both fluence and patient dose calculations were performed with a MC algorithm. DPEs of the mixed method were quantified by evaluating the differences between the mixed optimization SC dose result and a MC dose recalculation of the mixed optimization solution. OCEs of the mixed method were quantified by evaluating the differences between the MC recalculation of the mixed optimization solution and the final MC optimization solution. The results were analyzed through dose volume indices derived from the cumulative dose-volume histograms for selected anatomic structures. Statistical equivalence tests were used to determine the significance of the DPEs and the OCEs. Furthermore, a correlation analysis between DPEs and OCEs was performed. The evaluated DPEs were within +/- 2.8% while the OCEs were within 5.5%, indicating that OCEs can be clinically significant even when DPEs are clinically insignificant. The full MC-dose-based optimization reduced normal tissue dose by as much as 8.5% compared with the mixed-method optimization results. The DPEs and the OCEs in the targets had correlation coefficients greater than 0.71, and there was no correlation for the organs at risk. Because full MC-based optimization results in lower normal tissue doses, this method proves advantageous for HN IMRT optimization.

The effect of statistical noise on IMRT plan quality and convergence for MC-based and MC-correction-based optimized treatment plans

Monte Carlo (MC) is rarely used for IMRT plan optimization outside of research centres due to the extensive computational resources or long computation times required to complete the process. Time can be reduced by degrading the statistical precision of the MC dose calculation used within the optimization loop. However, this eventually introduces optimization convergence errors (OCEs). This study determines the statistical noise levels tolerated during MC-IMRT optimization under the condition that the optimized plan has OCEs <100 cGy (1.5% of the prescription dose) for MC-optimized IMRT treatment plans.Seven-field prostate IMRT treatment plans for 10 prostate patients are used in this study. Pre-optimization is performed for deliverable beams with a pencil-beam (PB) dose algorithm. Further deliverable-based optimization proceeds using: (1) MC-based optimization, where dose is recomputed with MC after each intensity update or (2) a once-corrected (OC) MC-hybrid optimization, where a MC dose computation defines beam-by-beam dose correction matrices that are used during a PB-based optimization. Optimizations are performed with nominal per beam MC statistical precisions of 2, 5, 8, 10, 15, and 20%. Following optimizer convergence, beams are re-computed with MC using 2% per beam nominal statistical precision and the 2 PTV and 10 OAR dose indices used in the optimization objective function are tallied. For both the MC-optimization and OC-optimization methods, statistical equivalence tests found that OCEs are less than 1.5% of the prescription dose for plans optimized with nominal statistical uncertainties of up to 10% per beam. The achieved statistical uncertainty in the patient for the 10% per beam simulations from the combination of the 7 beams is ~3% with respect to maximum dose for voxels with D>0.5D(max). The MC dose computation time for the OC-optimization is only 6.2 minutes on a single 3 Ghz processor with results clinically equivalent to high precision MC computations.

Computational Challenges for Image-Guided Radiation Therapy: Framework and Current Research

It is arguable that the imaging and delivery hardware necessary for delivering real-time adaptive image-guided radiotherapy is available on high-end linear accelerators. Robust and computationally efficient software is the limiting factor in achieving highly accurate and precise radiotherapy to the constantly changing anatomy of a cancer patient. The limitations are not caused by the availability of algorithms but rather issues of reliability, integration, and calculation time. However, each of the software components is an active area of research and development at academic and commercial centers. This article describes the software solutions in 4 broad areas: deformable image registration, adaptive replanning, real-time image guidance, and dose calculation and accumulation. Given the pace of technological advancement, the integration of these software solutions to develop real-time adaptive image-guided radiotherapy and the associated challenges they bring will be implemented to varying degrees by all major manufacturers over the coming years.