Optimization of the regularization parameter in the Dual Annealing method used for the reconstruction of energy spectrum of electron beam generated by the AQURE mobile accelerator.

INTRODUCTION: The shape of the energy spectrum is an essential component of any electron beam Monte Carlo model. Due to specialized equipment and the long measurement time for the direct methods for determining the energy spectrum, attractive alternatives are backward spectrum reconstructions from the measured data. One such approach is solving the first-degree Fredholm integral equation with appropriate regularization. It makes it possible to calculate the depth distribution as the sum of the distributions from monoenergetic beams. This study aims to determine the optimal value of the regularization parameter for the problem of determining the spectrum of the electron beam produced by a mobile accelerator used during intraoperative radiotherapy. MATERIAL AND METHODS: The Geant4 package was used to generate the distributions of deep doses for monoenergetic beams for two models with different degrees of complexity, i.e. simple (theoretical) and full (for the mobile accelerator). The dose distributions for four different shapes of energy spectrum (for each model) were obtained similarly. They were established as the reference data for further calculations. The Dual Annealing optimization method was used to obtain the reconstructed spectrum. The multiple optimizations that differ by the regularization parameter (ranging from 0 to 1) were performed. For each reconstruction, similarity indicators of the energy spectrum and the dose distribution to the referenced data were calculated to determine the optimal regularization parameters. RESULTS: Optimal regularization parameters determined by similarity indicators for the spectrum and the dose distribution differ for geometry models considered in the study. The regularization parameter for the simple geometry ranged from 0.03 to 0.05, while for full geometry, they were from 0.05 to 0.06. The results for conventional linear accelerators found in the literature range from 0.5 to 1.1. CONCLUSION: The Dual Annealing optimization method can be effectively used to solve the Fredholm equation with Tikhonov regularization to reconstruct an electron beam's energy spectrum. The regularization parameter value depends on the beam-forming system. Its value for the mobile accelerator considered in the study ranges from 0.05 to 0.06, depending on the nominal beam energy value.

Impact of the Intra- and Inter-observer Variability in the Delineation of Parotid Glands on the Dose Calculation During Head and Neck Helical Tomotherapy.

The intra- and inter-observer variability in delineation of the parotids on the kilo-voltage computed tomography (kVCT) and mega-voltage computed tomography (MVCT) were examined to establish their impact on the dose calculation during adaptive head and neck helical tomotherapy (HT). Three observers delineated left and right parotids for ten randomly selected patients with oropharynx cancer treated on HT. The pre-treatment kVCT and the MVCT from the first fraction of irradiation were selected to delineation. The delineation procedure was repeated three times by each observer. The parotids were delineated according to the institutional protocol. The analyses included intra-observer reproducibility and inter-structure, -observer and -modality variability of the volume and dose. The differences between the left and right parotid outlines were not statistically significant (p > 0.3). The reproducibility of the delineation was confirmed for each observer on the kVCT (p > 0.2) and on the MVCT (p > 0.1). The inter-observer variability of the outlines was significant (p < 0.001) as well as the inter-modality variability (p < 0.006). The parotids delineated on the MVCT were 10% smaller than on the kVCT. The inter-observer variability of the parotids delineation did not affect the average dose (p = 0.096 on the kVCT and p = 0.176 on the MVCT). The dose calculated on the MVCT was higher by 3.3% than dose from the kVCT (p = 0.009). Usage of the institutional protocols for the parotids delineation reduces intra-observer variability and increases reproducibility of the outlines. These protocols do not eliminate delineation differences between the observers, but these differences are not clinically significant and do not affect average doses in the parotids. The volumes of the parotids delineated on the MVCT are smaller than on the kVCT, which affects the differences in the calculated doses.

Impact of the Intra- and Inter-observer Variability in the Delineation of Parotid Glands on the Dose Calculation During Head and Neck Helical Tomotherapy.

The intra- and inter-observer variability in delineation of the parotids on the kilo-voltage computed tomography (kVCT) and mega-voltage computed tomography (MVCT) were examined to establish their impact on the dose calculation during adaptive head and neck helical tomotherapy (HT). Three observers delineated left and right parotids for ten randomly selected patients with oropharynx cancer treated on HT. The pre-treatment kVCT and the MVCT from the first fraction of irradiation were selected to delineation. The delineation procedure was repeated three times by each observer. The parotids were delineated according to the institutional protocol. The analyses included intra-observer reproducibility and inter-structure, -observer and -modality variability of the volume and dose. The differences between the left and right parotid outlines were not statistically significant (p > 0.3). The reproducibility of the delineation was confirmed for each observer on the kVCT (p > 0.2) and on the MVCT (p > 0.1). The inter-observer variability of the outlines was significant (p < 0.001) as well as the inter-modality variability (p < 0.006). The parotids delineated on the MVCT were 10% smaller than on the kVCT. The inter-observer variability of the parotids delineation did not affect the average dose (p = 0.096 on the kVCT and p = 0.176 on the MVCT). The dose calculated on the MVCT was higher by 3.3% than dose from the kVCT (p = 0.009). Usage of the institutional protocols for the parotids delineation reduces intra-observer variability and increases reproducibility of the outlines. These protocols do not eliminate delineation differences between the observers, but these differences are not clinically significant and do not affect average doses in the parotids. The volumes of the parotids delineated on the MVCT are smaller than on the kVCT, which affects the differences in the calculated doses.

Evaluation of image-guidance strategies for prostate cancer.

In this study, set-up accuracy and time consumption of different image-guidance protocols used for prostate cancer patients were compared. Set-up corrections from 60 prostate cancer patients treated on helical tomotherapy (HT) were used to simulate four types of image-guidance protocols which were based on: (i) a limited number of imaging sessions (IG-1), (ii) reduced registration tasks during daily imaging (IG-2), or (iii) and (iv) mixed methods of imaging (IG-3, IG-4). Each protocol was evaluated for three referencing scenarios based on the first fraction, first three fractions and first five fractions. Residual set-up error, the difference between the average set-up correction and the actual correction required, was used to evaluate the accuracy of each protocol. The first five fractions referencing scenario provides the highest reduction of the margins for each image-guidance protocol evaluated in this study. The first type of protocol is the shortest way to the effective correction of the systematic component of set-up error. For the second type of the protocol, the control of the residual errors is better and, as a result, the reduction of the margins is more significant than that obtained for the first one. Moreover, the second type of the protocol provides the highest accuracy of delivered dose. The result obtained for the fourth type of protocol does not decrease the calculated margins or increase their accuracy in correspondence to the no image guidance scheme. The fourth type of the protocol is not recommended as a protocol to be used to increase the conformity of the dose. The choice of the rest protocols should be validated in the context of (i) institutional practice regarding patient set-up procedure and its time consumption, (ii) acceptable balance between the amount of the dose delivered to the organ at risk and the additional imaging dose and (iii) patient anatomical conditions.

Tomotherapy: implications on daily workload and scheduling patients based on three years' institutional experience.

Helical tomotherapy (HT) was introduced at the Greater Poland Cancer Centre (GPCC) in April 2009. Retrospective analysis included data from the treatments performed for the first 656 patients treated with HT between May 2009 and May 2012 at the GPCC. In order to evaluate the implications on daily workload and scheduling of patients, stepwise regression and time analysis for each component of the overall treatment time, such as positioning, imaging, registration, and irradiation were performed. A detailed analysis included: (1) learning curves and optimized time needed for positioning and registration; (2) relation between irradiation time and parameters used for plan creation; and (3) average time of daily imaging. The irradiation component has the highest influence on the overall treatment time (R = 0.911). The lowest influence was observed for the imaging (R = 0.670). The learning curve for positioning was 7 months while the reduction of the average daily time needed for registration was observed even after two years. The irradiation time strongly depends on the planning parameters. Changing the pitch from 0.215 to 0.287 for pelvic cancer cases decreased the average daily beam-on time per patient by about 2 minutes. Similar changes for head and neck reduced this time by 1.3 minutes. The limitation in the usage of 1 cm field width only for complex cases, lower than 10 cm in the cranio-caudal direction, reduced the beam-on time per patient by 2 minutes. The average overall treatment time decreased from 21.5 minutes per patient in the first year of the HT usage to 13.8 minutes per patient in current practice. Our current practice shows that for a group of patients including mainly those with pelvis and head and neck cancers, the HT treatment takes approximately 15 minutes per patient allowing 40 patients to be treated within 10 hours.

B-Spline registration based on new concept of an intelligent masking procedure and GPU computations for the head and neck adaptive tomotherapy.

The deformable image registration (DIR) procedure has been optimized for helical tomotherapy. The data on registration shifts obtained on matching planning image with pre-treatment megavoltage CT are used in our software for acceleration of the first step (rigid registration) of the DIR procedure and for implementation of the B-Spline algorithm with intelligent masking. Priorities of the masks were automatically calculated based on disagreement detected during rigid registration. Evaluation tasks included: (a) comparison of accuracy and rate for schemes of pre-registered and non-registered images; (b) qualification of the effectiveness of the intelligent masking process, and (c) determination of acceleration of achievable with GPU computing. A specially designed head and neck phantom used for evaluation included structures with controlled changes of position, volume, density, and shape. Re-contouring procedures were performed with an Adaptive Planning software (Tomotherapy Inc.). No statistical difference was observed in accuracy of DIR based on structure position match on the tomotherapy unit and non pre-registered images (p > 0.7). Using pre-registered data reduces the total time required for execution of the elastic registration procedure by 5%. These data are also necessary for intelligent masking procedure during B-Spine registration. Intelligent masking procedure increases accuracy of the registration for a masked structure (p < 0.04) without decreasing the accuracy in non-masked tissues and additionally reduces the total time by 13%. GPU computations speed up procedure 30 times. GPU computing of the DIR in current status of our investigation could be realized in a relatively short time after pre-treatment imaging. The proposed approach can be used in the routine assessment of anatomic changes occurring in healthy tissue during the course of radiotherapy. Further developments will be concentrated on the full integration of DIR computations in the imaging and treatment process of helical tomotherapy.