Developing a phenomenological model of the proton trajectory within a heterogeneous medium required for proton imaging.

To develop an accurate phenomenological model of the cubic spline path estimate of the proton path, accounting for the initial proton energy and water equivalent thickness (WET) traversed. Monte Carlo (MC) simulations were used to calculate the path of protons crossing various WET (10-30 cm) of different material (LN300, water and CB2-50% CaCO3) for a range of initial energies (180-330 MeV). For each MC trajectory, cubic spline trajectories (CST) were constructed based on the entrance and exit information of the protons and compared with the MC using the root mean square (RMS) metric. The CST path is dependent on the direction vector magnitudes (|P0,1|). First, |P0,1| is set to the proton path length (with factor Λ(Norm)(0,1) = 1.0). Then, two optimal factor Λ(0,1) are introduced in |P0,1|. The factors are varied to minimize the RMS difference with MC paths for every configuration. A set of Λ(opt)(0,1) factors, function of WET/water equivalent path length (WEPL), that minimizes the RMS are presented. MTF analysis is then performed on proton radiographs of a line-pair phantom reconstructed using the CST trajectories. Λ(opt)(0,1) was fitted to the WET/WEPL ratio using a quadratic function (Y = A + BX(2) where A = 1.01,0.99, B = 0.43,- 0.46 respectively for Λ(opt)(0), Λ(opt)(1)). The RMS deviation calculated along the path, between the CST and the MC, increases with the WET. The increase is larger when using Λ(Norm)(0,1) than Λ(opt)(0,1) (difference of 5.0% with WET/WEPL = 0.66). For 230/330 MeV protons, the MTF10% was found to increase by 40/16% respectively for a thin phantom (15 cm) when using the Λ(opt)(0,1) model compared to the Λ(Norm)(0,1) model. Calculation times for Λ(opt)(0,1) are scaled down compared to MLP and RMS deviation are similar within standard deviation.B ased on the results of this study, using CST with the Λ(opt)(0,1) factors reduces the RMS deviation and increases the spatial resolution when reconstructing proton trajectories.

Adaptation of the CVT algorithm for catheter optimization in high dose rate brachytherapy.

PURPOSE: An innovative, simple, and fast method to optimize the number and position of catheters is presented for prostate and breast high dose rate (HDR) brachytherapy, both for arbitrary templates or template-free implants (such as robotic templates). METHODS: Eight clinical cases were chosen randomly from a bank of patients, previously treated in our clinic to test our method. The 2D Centroidal Voronoi Tessellations (CVT) algorithm was adapted to distribute catheters uniformly in space, within the maximum external contour of the planning target volume. The catheters optimization procedure includes the inverse planning simulated annealing algorithm (IPSA). Complete treatment plans can then be generated from the algorithm for different number of catheters. The best plan is chosen from different dosimetry criteria and will automatically provide the number of catheters and their positions. After the CVT algorithm parameters were optimized for speed and dosimetric results, it was validated against prostate clinical cases, using clinically relevant dose parameters. The robustness to implantation error was also evaluated. Finally, the efficiency of the method was tested in breast interstitial HDR brachytherapy cases. RESULTS: The effect of the number and locations of the catheters on prostate cancer patients was studied. Treatment plans with a better or equivalent dose distributions could be obtained with fewer catheters. A better or equal prostate V100 was obtained down to 12 catheters. Plans with nine or less catheters would not be clinically acceptable in terms of prostate V100 and D90. Implantation errors up to 3 mm were acceptable since no statistical difference was found when compared to 0 mm error (p > 0.05). No significant difference in dosimetric indices was observed for the different combination of parameters within the CVT algorithm. A linear relation was found between the number of random points and the optimization time of the CVT algorithm. Because the computation time decrease with the number of points and that no effects were observed on the dosimetric indices when varying the number of sampling points and the number of iterations, they were respectively fixed to 2500 and to 100. The computation time to obtain ten complete treatments plans ranging from 9 to 18 catheters, with the corresponding dosimetric indices, was 90 s. However, 93% of the computation time is used by a research version of IPSA. For the breast, on average, the Radiation Therapy Oncology Group recommendations would be satisfied down to 12 catheters. Plans with nine or less catheters would not be clinically acceptable in terms of V100, dose homogeneity index, and D90. CONCLUSIONS: The authors have devised a simple, fast and efficient method to optimize the number and position of catheters in interstitial HDR brachytherapy. The method was shown to be robust for both prostate and breast HDR brachytherapy. More importantly, the computation time of the algorithm is acceptable for clinical use. Ultimately, this catheter optimization algorithm could be coupled with a 3D ultrasound system to allow real-time guidance and planning in HDR brachytherapy.