Feasibility-guided automated planning for stereotactic treatments of prostate cancer.

Significant improvements in plan quality using automated planning have been previously demonstrated. The aim of this study was to develop an optimal automated class solution for stereotactic radiotherapy (SBRT) planning of prostate cancer using the new Feasibility module implemented in the pinnacle evolution. Twelve patients were retrospectively enrolled in this planning study. Five plans were designed for each patient. Four plans were automatically generated using the 4 proposed templates for SBRT optimization implemented in the new pinnacle evolution treatment planning systems, differing for different settings of dose-fallout (low, medium, high and veryhigh). Based on the obtained results, the fifth plan (feas) was generated customizing the template with the optimal criteria obtained from the previous step and integrating in the template the "a-priori" knowledge of OARs sparing based on the Feasibility module, able to estimate the best possible dose-volume histograms of OARs before starting optimization. Prescribed dose was 35 Gy to the prostate in 5 fractions. All plans were generated with a full volumetric-modulated arc therapy arc and 6MV flattening filter-free beams, and optimized to ensure the same target coverage (95% of the prescription dose to 98% of the target). Plans were assessed according to dosimetric parameters and planning and delivery efficiency. Differences among the plans were evaluated using a Kruskal-Wallis 1-way analysis of variance. The requests for more aggressive objectives for dose falloff parameters (from low to veryhigh) translated in a statistically significant improvement of dose conformity, but at the expense of a dose homogeneity. The best automated plans in terms of best trade-off between target coverage and OARs sparing among the 4 plans automatically generated by the SBRT module were the high plans. The veryhigh plans reported a significant increase of high-doses to prostate, rectum, and bladder that was considered dosimetrically and clinically unacceptable. The feas plans were optimized on the basis on high plans, reporting significant reduction of rectum irradiation; Dmean, and V18 decreased by 19% to 23% (p = 0.031) and 4% to 7% (p = 0.059), respectively. No statistically significant differences were found in femoral heads and penile bulb irradiation for all dosimetric metrics. feas plans showed a significant increase of MU/Gy (mean: 368; p = 0.004), reflecting an increased level of fluence modulation. Thanks to the new efficient optimization engines implemented in pinnacle evolution (L-BFGS and layered graph), mean planning time was decreased to less than 10 minutes for all plans and all techniques. The integration of dose-volume histograms a-priori knowledge provided by the feasibility module in the automated planning process for SBRT planning has shown to significantly improve plan quality compared to generic protocol values as inputs.

Reproducibility and stability of spirometer-guided deep inspiration breath-hold in left-breast treatments using an optical surface monitoring system.

The aim of this study was to evaluate the reproducibility and stability of left breast positioning during spirometer-guided deep-inspiration breath-hold (DIBH) radiotherapy using an optical surface imaging system (AlignRT). The AlignRT optical tracking system was used to monitor five left-sided breast cancer patients treated using the Active Breathing Coordinator spirometer with DIBH technique. Treatment plans were created using an automated hybrid-VMAT technique on DIBH CTs. A prescribed dose of 60 Gy to the tumor bed and 50 Gy to the breast in 25 fractions was planned. During each treatment session, the antero-posterior (VRT), superior-inferior (LNG), and lateral (LAT) motion of patients was continuously recorded by AlignRT. The intra-breath-hold stability and the intra- and inter-fraction reproducibility were analyzed for all breath-holds and treatment fractions. The dosimetric impact of the residual motion during DIBH was evaluated from the isocenter shifts amplitudes obtained from the 50%, 90%, and 100% cumulative distribution functions of intra-fractional reproducibility. The positional variations of 590 breath-holds as measured by AlignRT were evaluated. The mean intra-breath-hold stability during DIBH was 1.0 ± 0.4 mm, 2.1 ± 1.9 mm, and 0.7 ± 0.5 mm in the VRT, LNG, and LAT directions, with a maximal value of 8.8 mm in LNG direction. Similarly, the mean intra-breath-hold reproducibility was 1.4 ± 0.8 mm, 1.7 ± 1.0 mm, and 0.8 ± 0.5 mm in the VRT, LNG, and LAT directions, with a maximal value of 4.1 mm in LNG direction. Inter-fractional reproducibility showed better reliability, with difference in breathing levels in all fractions of 0.3 mm on average. Based on tolerance limits corresponding to the 90% cumulative distribution level, gating window widths of 1 mm, 2 mm, and 5 mm in the LAT, VRT, and LNG directions were considered an appropriate choice. In conclusion, despite the use of a dedicated spirometer at constant tidal volume, a non-negligible variability of the breast surface position has been reported during breath-holds. The real-time monitoring of breast surface using surface-guided optical technology is strongly recommended.

Prediction and classification of VMAT dosimetric accuracy using plan complexity and log-files analysis.

PURPOSE: We presented different machine learning models based on log files analysis and complexity indexes to predict and classify the dosimetric accuracy of VMAT plans. METHODS: A total of 1302 VMAT arcs from 651 treatment plans were analyzed using the modulation complexity score (MCS) and the dynamic log-files generated by the linac. Predicted and measured fluences were compared using γ-analysis in terms of mean γ-values (γmean) and γ-pass rate (γ%). A kernel regression model was developed aiming to predict individual γ% and γmean values. Multinomial logistic regression (LR), Naïve-Bayes (NB) and support vector machine (SVM) models were developed based on MCS values to classify QA results as "pass" (γ%greater than90 % and γmean < 0.5), "control" (80 % < Î³% < 90 % and 0.50 < Î³mean < 0.75) and "fail" (γ% < 80 % and γmean > 0.75). Training, validation and testing groups were used to evaluate the model reliability. A complexity-based traffic light protocol was implemented to flag pass (green light), control (orange light) and failed plans (red light). RESULTS: Prediction accuracy of residuals for γ% was 2.1 % and 2.2 % in the training and testing cohorts, respectively. For 2 %(local)/2mm, both γ% and γmean classification performances reported weighted precision, recall and F1-values greater than 90 % for all machine learning models. The optimal MCS threshold value for the identification of failed plans was 0.130, with a sensibility and specificity of 0.994 and 0.952, respectively. The optimal MCS threshold for reliable plans was 0.270, with a sensitivity and specificity of 0.925 and 0.922, respectively. CONCLUSIONS: Machine learning can accurately predict the dosimetric accuracy of VMAT treatments, representing an efficient tool to assist patient-specific QA.

Prediction of VMAT delivery accuracy using plan modulation complexity score and log-files analysis.

The purpose of this study was to develop a predictive model based on plan complexity metrics and linac log-files analysis to classify the dosimetric accuracy of VMAT plans. A total of 612 VMAT plans, corresponding to 1224 arcs, were analyzed. All VMAT arcs underwent pre-treatment verification that was performed by means of the dynamic log-files generated by the linac. The comparison of predicted (by TPS) and measured (by log-files) integral fluences was performed usingγ-analysis in terms of the percentage of points withγ-value smaller than one (γ%) and using a stringent 2%(local)/2 mm criteria. Thisγ-analysis was performed by a commercial software LinacWatch. The action limits (AL) were derived from the mean values, standard deviations and the confidence limit (CL) of theγ% distribution. A plan complexity metric, the modulation complexity score (MCS), based on the aperture beam area variability and leaf sequence variability was used as input variable of the model. A binary logistic regression (LR) model was developed to classify QA results as 'pass' (γ% ≥ AL) or 'fail' (γ% < AL). Receiver operator characteristics (ROC) curves were used to determine the optimal MCS threshold to flag 'failed' plans that need to be re-optimized. The model reliability was evaluated stratifying the plans in training, validation and testing groups. The confidence and action limits forγ% were found 20.1% and 79.9%, respectively. The accuracy of the model for the training and testing dataset was 97.4% and 98.0%, respectively. The optimal MCS threshold value for the identification of failed plans was 0.142, providing a true positive rate able to flag the plans failing QA of 91%. In clinical routine, the use of this MCS threshold may allow the prompt identification of overly modulated plans, then reducing the number of QA failures and improving the quality of VMAT plans used for treatment.

Epid-based in vivo dose verification for lung stereotactic treatments delivered with multiple breath-hold segmented volumetric modulated arc therapy.

We evaluated an EPID-based in-vivo dosimetry (IVD) method for the dose verification and the treatment reproducibility of lung SBRT-VMAT treatments in clinical routine. Ten patients with lung metastases treated with Elekta VMAT technique were enrolled. All patients were irradiated in five consecutive fractions, with total doses of 50 Gy. Set-up was carried out with the Elekta stereotactic body frame. Eight patients were simulated and treated using the Active Breath Control (ABC) system, a spirometer enabling patients to maintain a breath-hold at a predetermined lung volume. Two patients were simulated and treated in free-breathing using an abdominal compressor. IVD was performed using the SOFTDISO software. IVD tests were evaluated by means of (a) ratio R between daily in-vivo isocenter dose and planned dose and (b) γ-analysis between EPID integral portal images in terms of percentage of points with γ-value smaller than one (γ% ) and mean γ-values (γmean ) using a 3%(global)/3 mm criteria. Alert criteria of ±5% for R ratio, γ%  < 90%, and γmean  > 0.67 were chosen. 50 transit EPID images were acquired. For the patients treated with ABC spirometer, the results reported a high level of accuracy in dose delivery with 100% of tests within ±5%. The γ-analysis showed a mean value of γmean equal to 0.21 (range: 0.04-0.56) and a mean γ% equal to 96.9 (range: 78-100). Relevant discrepancies were observed only for the two patients treated without ABC, mainly due to a blurring dose effect due to residual respiratory motion. Our method provided a fast and accurate procedure in clinical routine for verifying delivered dose as well as for detecting errors.

Linac-based extracranial radiosurgery with Elekta volumetric modulated arc therapy and an anatomy-based treatment planning system: Feasibility and initial experience.

We reported our initial experience in using Elekta volumetric modulated arc therapy (VMAT) and an anatomy-based treatment planning system (TPS) for single high-dose radiosurgery (SRS-VMAT) of liver metastases. This study included a cohort of 12 patients treated with a 26-Gy single fraction. Single-arc VMAT plans were generated with Ergo++ TPS. The prescription isodose surface (IDS) was selected to fulfill the 2 following criteria: 95% of planning target volume (PTV) reached 100% of the prescription dose and 99% of PTV reached a minimum of 90% of prescription dose. A 1-mm multileaf collimator (MLC) block margin was added around the PTV. For a comparison of dose distributions with literature data, several conformity indexes (conformity index [CI], conformation number [CN], and gradient index [GI]) were calculated. Treatment efficiency and pretreatment dosimetric verification were assessed. Early clinical data were also reported. Our results reported that target and organ-at-risk objectives were met for all patients. Mean and maximum doses to PTVs were on average 112.9% and 121.5% of prescribed dose, respectively. A very high degree of dose conformity was obtained, with CI, CN, and GI average values equal to 1.29, 0.80, and 3.63, respectively. The beam-on-time was on average 9.3 minutes, i.e., 0.36min/Gy. The mean number of monitor units was 3162, i.e., 121.6MU/Gy. Pretreatment verification (3%-3mm) showed an optimal agreement with calculated values; mean γ value was 0.27 and 98.2% of measured points resulted with γ < 1. With a median follow-up of 16 months complete response was observed in 12/14 (86%) lesions; partial response was observed in 2/14 (14%) lesions. No radiation-induced liver disease (RILD) was observed in any patients as well no duodenal ulceration or esophagitis or gastric hemorrhage. In conclusion, this analysis demonstrated the feasibility and the appropriateness of high-dose single-fraction SRS-VMAT in liver metastases performed with Elekta VMAT and Ergo++ TPS. Preliminary clinical outcomes showed a high rate of local control and minimum incidence of acute toxicity.

Real time transit dosimetry for the breath-hold radiotherapy technique: an initial experience.

INTRODUCTION: The breath-hold is one of the techniques to obtain the dose escalation for lung tumors. However, the change of the patient's breath pattern can influence the stability of the inhaled air volume, IAV, used in this work as a surrogate parameter to assure the tumor position reproducibility during dose delivery. MATERIALS AND METHOD: In this paper, an Elekta active breathing coordinator has been used for lung tumor irradiation. This device is not an absolute spirometer and the feasibility study here presented developed (i) the possibility to select a specific range epsilon of IAV values comfortable for the patient and (ii) the ability of a transit signal rate S(t), obtained by a small ion-chamber positioned on the portal image device, to supply in real time the in vivo isocenter dose reproducibility. Indeed, while the selection of the IAV range depends on the patient's ability to follow instructions for breath-hold, the S(t) monitoring can supply to the radiation therapist a surrogate of the tumor irradiation reproducibility. RESULTS: The detection of the S(t) in real time during breath-hold was used to determine the interfraction isocenter dose variations due to the reproducibility of the patient's breathing pattern. The agreement between the reconstructed and planned isocenter dose in breath-hold at the interfraction level was well within 1.5%, while in free breathing a disagreement up to 8% was observed. The standard deviation of the S(t) in breath-hold observed at the intrafraction level is a bit higher than the one obtained without the patient and this can be justified by the presence of a small residual tumor motion as heartbeat. CONCLUSION: The technique is simple and can be implemented for routine use in a busy clinic.

Comparison of measured and computed portal dose for IMRT treatment.

A new 2D array Seven 29T model (PTW, Freiburg), equipped with 729 vented plane-parallel ion chambers, projected for pretreatment verification of radiotherapy plans, was used as a detector for the transmitted or portal dose measurements below a Rando phantom. The dosimetric qualities of the 2D array make it attractive for measuring transmitted dose maps from step-and-shoot intensity-modulated radiotherapy (IMRT). It is well known that for step-and-shoot IMRT beams that use a small number of monitor units (MUs) per sequence, the early and recent electronic portal imaging devices (EPIDs) present a different response at X-ray start-up that affects the accuracy of the measured transmitted dose. The comparison of portal doses measured to those calculated by a commercial treatment-planning system (TPS) can verify correct dose delivery during treatment. This direct validation was tested by irradiating a simulated head tumor in a Rando anthropomorphic phantom by step-and-shoot IMRT beams. The absolute transmitted doses on a plane orthogonal to the beam central axis below the phantom were measured by the 2D array calibrated in terms of dose to water and compared with the computed portal dose extracted by custom software. In a previous paper, the comparison between the IMRT portal doses, computed by a commercial TPS and measured by a linear array that supplied a 1 mm spatial dose resolution, was carried out. The gamma-index analysis supplied an agreement of more than 95% of the dose point with acceptance criteria, in terms of dose difference, DeltaDmax, and distance agreement, deltadmax, equal to 4% and 4 mm, respectively. In this paper, we verify the possible use of the PTW 2D array for measurements of the transmitted doses during several fractions of head and neck tumor radiotherapy. There are two advantages in the use of this 2D array as a portal dose device for the IMRT quality assurance program: first is the ability to perform absolute dose comparisons for hundreds of measurement positions to verify the correct dose delivery in several fractions of the therapy; second is the efficiency in time to detect these kinds of dose distributions within the field of view area of the CT scanner.

In-vivo portal dosimetry by an ionization chamber.

As all methods for in-vivo dosimetry require special efforts many physicists are often discouraged in verifying the middle dose in a patient along the beam central axis. This work reports a practical method for the determination of the middle dose value, D(m), on the central beam axis, using a signal S(t), obtained by a small thimble ion-chamber positioned at the center of the electronic portal imaging device, and irradiated by the X-ray beam transmitted through the patient. The use of a stable ion-chamber reduces many of the disadvantages associated to the use of diodes as their periodic recalibration and time consuming positioning. The method makes use of a set of correlation functions obtained by the S(t) and D(m) ratios, determined by irradiating a water-equivalent phantom with 6 MV, 10 MV and 5 MV X-ray beams. Several tests were carried out in phantoms with asymmetric inhomogeneities. The method here proposed is based on the determination of the water-equivalent thickness of the patient, along the beam central axis, by the treatment planning system that makes use of the electron densities obtained by a computer tomography scanner, that works with calibrated Hounsfield numbers. This way, it is therefore possible to compare the dose, D(m, TPS), obtained by a treatment planning system, with the in-vivo dose D(m) value, both defined at density middle point (identified along the beam central axis, where the thick material, in terms of g cm(-2), above and below, is the same). The method has been applied for the in-vivo dosimetry of 30 patients, treated with conformed beams for pelvic tumor, checking: anterior-posterior or posterior-anterior irradiations and lateral-lateral irradiations. For every checked field at least five measurements were carried out. Applying a correct quality assurance program based on the tests of the patient set-up, machine settings and calculations, results showed that the method is able to verify agreements between the dose D(m,TPS) and the in-vivo dose value D(m), within 4% for 95% of the 240 measurements carried out in-vivo.