Brigida C. Ferreira

Evaluation of an Epson flatbed scanner to read Gafchromic EBT films for radiation dosimetry.

Gafchromic EBT (EBT) films are becoming increasingly popular due to their advantageous properties. When flatbed colour scanners are used for film dosimetry, a good quality control of the scanning device is a crucial step for accurate results. The proposal of this work was to fully assess the performance of the scanner Epson Expression 10000XL in order to quantify all parameters and needed corrections to minimize dose uncertainties. A standard step tablet, with 32 steps and optical densities from 0.06 to 3.8, was used to check the scanner linearity. The scanner warming-up effect and reproducibility were evaluated by performing 30 consecutive scans plus 20 scans in 15 min intervals. The scanning colour modes: 24 and 48 bits and scanning resolutions from 50 to 300 dpi were tested. A Wiener filter with different pixels regions was applied with the purpose of reducing the film noise. All scans were made in transmission mode with a constant film orientation. The red colour channel was posteriorly extracted from the images to maximize readout sensitivity. Two EBT films were irradiated, perpendicularly and parallel to beam incidence, with a 6 MV photon beam with doses that ranged from 0.2 to 3 Gy. A polynomial expression was used to convert optical density into dose. Dose uncertainty was quantified applying error propagation analysis. A correction for the non-uniform response of the scanner was determined using five films irradiated with a uniform dose. The scanner response was linear until an optical density of approximately 1 which corresponds to doses higher than those of clinical interest for EBT films. The scanner signal stabilized after seven readings. Scanner reproducibility around 0.2% was obtained either with the scanner warm or cold. However, reproducibility was significantly reduced when comparing images digitized with the scanner at different temperatures. Neither the colour depth mode, the scanning resolution, the multiscan option nor the Wiener filter had a significant effect on the shape of the calibration curve. However, a reduction in dose uncertainty was possible by selecting appropriate reading parameters. These are a 48 bit colour depth, a 75 dpi resolution and repeating the scan four times. Finally, the two dimensional Wiener filter applied to a 3 x 3 pixel region to the red component of the image reduced the experimental scan uncertainty to about 0.5% for doses higher than 0.5 Gy. Total scan uncertainty was less than 2% for a perpendicular calibration and reduced to less than 1% for a parallel calibration. A dose over-estimation of around 5% for clinical doses may be made if the image acquired is not corrected for the non-uniform response of the scanner. A protocol to read EBT films using the Epson Expression 10000XL scanner was established for IMRT verification. The contribution for the overall uncertainty in film dosimetry coming from the scanning process was estimated to be around 0.5% for doses higher than 0.5 Gy when reading parameters are optimized. Total scan uncertainty achieved is about 2% when using a perpendicular calibration. It can further be reduced if a parallel calibration is used.

Beam angle optimization for intensity-modulated radiation therapy using a guided pattern search method

Generally, the inverse planning of radiation therapy consists mainly of the fluence optimization. The beam angle optimization (BAO) in intensity-modulated radiation therapy (IMRT) consists of selecting appropriate radiation incidence directions and may influence the quality of the IMRT plans, both to enhance better organ sparing and to improve tumor coverage. However, in clinical practice, most of the time, beam directions continue to be manually selected by the treatment planner without objective and rigorous criteria. The goal of this paper is to introduce a novel approach that uses beams-eye-view dose ray tracing metrics within a pattern search method framework in the optimization of the highly non-convex BAO problem. Pattern search methods are derivative-free optimization methods that require a few function evaluations to progress and converge and have the ability to better avoid local entrapment. The pattern search method framework is composed of a search step and a poll step at each iteration. The poll step performs a local search in a mesh neighborhood and ensures the convergence to a local minimizer or stationary point. The search step provides the flexibility for a global search since it allows searches away from the neighborhood of the current iterate. Beams-eye-view dose metrics assign a score to each radiation beam direction and can be used within the pattern search framework furnishing a priori knowledge of the problem so that directions with larger dosimetric scores are tested first. A set of clinical cases of head-and-neck tumors treated at the Portuguese Institute of Oncology of Coimbra is used to discuss the potential of this approach in the optimization of the BAO problem.

Treatment plan comparison between helical tomotherapy and MLC-based IMRT using radiobiological measures

The rapid implementation of advanced treatment planning and delivery technologies for radiation therapy has brought new challenges in evaluating the most effective treatment modality. Intensity-modulated radiotherapy (IMRT) using multi-leaf collimators (MLC) and helical tomotherapy (HT) are becoming popular modes of treatment delivery and their application and effectiveness continues to be investigated. Presently, there are several treatment planning systems (TPS) that can generate and optimize IMRT plans based on user-defined objective functions for the internal target volume (ITV) and organs at risk (OAR). However, the radiobiological parameters of the different tumours and normal tissues are typically not taken into account during dose prescription and optimization of a treatment plan or during plan evaluation. The suitability of a treatment plan is typically decided based on dosimetric criteria such as dose-volume histograms (DVH), maximum, minimum, mean and standard deviation of the dose distribution. For a more comprehensive treatment plan evaluation, the biologically effective uniform dose (D) is applied together with the complication-free tumour control probability (P(+)). Its utilization is demonstrated using three clinical cases that were planned with two different forms of IMRT. In this study, three different cancer types at different anatomical sites were investigated: head and neck, lung and prostate cancers. For each cancer type, a linac MLC-based step-and-shoot IMRT plan and a HT plan were developed. The MLC-based IMRT treatment plans were developed on the Philips treatment-planning platform, using the Pinnacle 7.6 software release. For the tomotherapy HiArt plans, the dedicated tomotherapy treatment planning station was used, running version 2.1.2. By using D as the common prescription point of the treatment plans and plotting the tissue response probabilities versus D for a range of prescription doses, a number of plan trials can be compared based on radiobiological measures. The applied plan evaluation method shows that in the head and neck cancer case the HT treatment gives better results than MLC-based IMRT in terms of expected clinical outcome P(+) of 62.2% and 46.0%, D to the ITV of 72.3 Gy and 70.7 Gy, respectively). In the lung cancer and prostate cancer cases, the MLC-based IMRT plans are better over the clinically useful dose prescription range. For the lung cancer case, the HT and MLC-based IMRT plans give a P(+) of 66.9% and 72.9%, D to the ITV of 64.0 Gy and 66.9 Gy, respectively. Similarly, for the prostate cancer case, the two radiation modalities give a P(+) of 68.7% and 72.2%, D to the ITV of 86.0 Gy and 85.9 Gy, respectively. If a higher risk of complications (higher than 5%) could be allowed, the complication-free tumour control could increase by over 40%, 2% and 30% compared to the initial dose prescription for the three cancer cases, respectively. Both MLC-based IMRT and HT can encompass the often-large ITV required while they minimize the volume of the organs at risk receiving high doses. Radiobiological evaluation of treatment plans may provide an improved correlation of the delivered treatment with the clinical outcome by taking into account the dose-response characteristics of the irradiated targets and normal tissues. There may exist clinical cases, which may look dosimetrically similar but in radiobiological terms may be quite different. In such situations, traditional dose-based evaluation tools can be complemented by the use of P(+)--D diagrams to effectively evaluate and compare treatment plans.

The impact of different dose–response parameters on biologically optimized IMRT in breast cancer

The full potential of biologically optimized radiation therapy can only be maximized with the prediction of individual patient radiosensitivity prior to treatment. Unfortunately, the available biological parameters, derived from clinical trials, reflect an average radiosensitivity of the examined populations. In the present study, a breast cancer patient of stage I-II with positive lymph nodes was chosen in order to analyse the effect of the variation of individual radiosensitivity on the optimal dose distribution. Thus, deviations from the average biological parameters, describing tumour, heart and lung response, were introduced covering the range of patient radiosensitivity reported in the literature. Two treatment configurations of three and seven biologically optimized intensity-modulated beams were employed. The different dose distributions were analysed using biological and physical parameters such as the complication-free tumour control probability (P(+)), the biologically effective uniform dose (D), dose volume histograms, mean doses, standard deviations, maximum and minimum doses. In the three-beam plan, the difference in P(+) between the optimal dose distribution (when the individual patient radiosensitivity is known) and the reference dose distribution, which is optimal for the average patient biology, ranges up to 13.9% when varying the radiosensitivity of the target volume, up to 0.9% when varying the radiosensitivity of the heart and up to 1.3% when varying the radiosensitivity of the lung. Similarly, in the seven-beam plan, the differences in P(+) are up to 13.1% for the target, up to 1.6% for the heart and up to 0.9% for the left lung. When the radiosensitivity of the most important tissues in breast cancer radiation therapy was simultaneously changed, the maximum gain in outcome was as high as 7.7%. The impact of the dose-response uncertainties on the treatment outcome was clinically insignificant for the majority of the simulated patients. However, the jump from generalized to individualized radiation therapy may significantly increase the therapeutic window for patients with extreme radio sensitivity or radioresistance, provided that these are identified. Even for radiosensitive patients a simple treatment technique is sufficient to maximize the outcome, since no significant benefits were obtained with a more complex technique using seven intensity-modulated beams portals.

Comparison of the helical tomotherapy and MLC-based IMRT radiation modalities in treating brain and cranio-spinal tumors

The investigation of the clinical efficacy and effectiveness of Intensity Modulated Radiotherapy (IMRT) using Multileaf Collimators (MLC) and Helical Tomotherapy (HT) has been an issue of increasing interest over the past few years. In order to assess the suitability of a treatment plan, dosimetric criteria such as dose-volume histograms (DVH), maximum, minimum, mean, and standard deviation of the dose distribution are typically used. Nevertheless, the radiobiological parameters of the different tumors and normal tissues are often not taken into account. The use of the biologically effective uniform dose (D̿) together with the complication-free tumor control probability (P+) were applied to evaluate the two radiation modalities. Two different clinical cases of brain and cranio-spinal axis cancers have been investigated by developing a linac MLC-based step-and-shoot IMRT plan and a Helical Tomotherapy plan. The treatment plans of the MLC-based IMRT were developed on the Philips treatment planning station using the Pinnacle 7.6 software release while the dedicated Tomotherapy treatment planning station was used for the HT plan. With the use of the P+ index and the D̿ concept as the common prescription point, the different treatment plans were compared based on radiobiological measures. The tissue response probabilities were plotted against D̿ for a range of prescription doses. The applied plan evaluation method shows that in the brain cancer, the HT treatment gives slightly better results than the MLC-based IMRT in terms of optimum expected clinical outcome (P+ of 66.1% and 63.5% for a D̿ to the PTV of 63.0 Gy and 62.0 Gy, respectively). In the cranio-spinal axis cancer, the HT plan is significantly better compared to the MLC-based IMRT plan over the clinically useful dose prescription range (P+ of 84.1% and 28.3% for a D̿ to the PTV of 50.6 Gy and 44.0 Gy, respectively). If a higher than 5% risk for complications could be allowed, the complication-free tumor control could be increased by almost 30% compared to the initial dose prescription. In comparison to MLC based-IMRT, HT can better encompass the often large PTV while minimizing the volume of the OARs receiving high dose. A radiobiological treatment plan evaluation can provide a closer association of the delivered treatment with the clinical outcome by taking into account the dose-response relations of the irradiated tumors and normal tissues. The use of P — D̿ diagrams can complement the traditional tools of evaluation such as DVHs, in order to compare and effectively evaluate different treatment plans.

Response-probability volume histograms and iso-probability of response charts in treatment plan evaluation

PURPOSE This study aims at demonstrating a new method for treatment plan evaluation and comparison based on the radiobiological response of individual voxels. This is performed by applying them on three different cancer types and treatment plans of different conformalities. Furthermore, their usefulness is examined in conjunction with traditionally applied radiobiological and dosimetric treatment plan evaluation criteria. METHODS Three different cancer types (head and neck, breast and prostate) were selected to quantify the benefits of the proposed treatment plan evaluation method. In each case, conventional conformal radiotherapy (CRT) and intensity modulated radiotherapy (IMRT) treatment configurations were planned. Iso-probability of response charts was produced by calculating the response probability in every voxel using the linear-quadratic-Poisson model and the dose-response parameters of the corresponding structure to which this voxel belongs. The overall probabilities of target and normal tissue responses were calculated using the Poisson and the relative seriality models, respectively. The 3D dose distribution converted to a 2 Gy fractionation, D2(GY) and iso-BED distributions are also shown and compared with the proposed methodology. Response-probability volume histograms (RVH) were derived and compared with common dose volume histograms (DVH). The different dose distributions were also compared using the complication-free tumor control probability, P+, the biologically effective uniform dose, D, and common dosimetric criteria. RESULTS 3D Iso-probability of response distributions is very useful for plan evaluation since their visual information focuses on the doses that are likely to have a larger clinical effect in that particular organ. The graphical display becomes independent of the prescription dose highlighting the local radiation therapy effect in each voxel without the loss of important spatial information. For example, due to the exponential nature of the Poisson distribution, cold spots in the target volumes or hot spots in the normal tissues are much easier to be identified. Response-volume histograms, as DVH, can also be derived and used for plan comparison. RVH are advantageous since by incorporating the radiobiological properties of each voxel they summarize the 3D distribution into 2D without the loss of relevant information. Thus, more clinically relevant radiobiological objectives and constraints could be defined and used in treatment planning optimization. These measures become increasingly important when dose distributions need to be designed according to the microscopic biological properties of tumor and normal tissues. CONCLUSIONS The proposed methods do not aim to replace quantifiers like the probabilities of total tissue response, which ultimately are the quantities of interest to evaluate treatment success. However, iso-probability of response charts and response-probability volume histograms illustrates more clearly the difference in effectiveness between different treatment plans than the information provided by alternative dosimetric data. The use of 3D iso-probability of response distributions could serve as a good descriptor of the effectiveness of a dose distribution indicating primarily the regions in a tissue that dominate its response.

Radiobiological evaluation of forward and inverse IMRT using different fractionations for head and neck tumours

PurposeTo quantify the radiobiological advantages obtained by an Improved Forward Planning technique (IFP) and two IMRT techniques using different fractionation schemes for the irradiation of head and neck tumours. The conventional radiation therapy technique (CONVT) was used here as a benchmark.MethodsSeven patients with head and neck tumours were selected for this retrospective planning study. The PTV1 included the primary tumour, PTV2 the high risk lymph nodes and PTV3 the low risk lymph nodes. Except for the conventional technique where a maximum dose of 64.8 Gy was prescribed to the PTV1, 70.2 Gy, 59.4 Gy and 50.4 Gy were prescribed respectively to PTV1, PTV2 and PTV3. Except for IMRT2, all techniques were delivered by three sequential phases. The IFP technique used five to seven directions with a total of 15 to 21 beams. The IMRT techniques used five to nine directions and around 80 segments. The first, IMRT1, was prescribed with the conventional fractionation scheme of 1.8 Gy per fraction delivered in 39 fractions by three treatment phases. The second, IMRT2, simultaneously irradiated the PTV2 and PTV3 with 59.4 Gy and 50.4 Gy in 28 fractions, respectively, while the PTV1 was boosted with six subsequent fractions of 1.8 Gy. Tissue response was calculated using the relative seriality model and the Poisson Linear-Quadratic-Time model to simulate repopulation in the primary tumour.ResultsThe average probability of total tumour control increased from 38% with CONVT to 80% with IFP, to 85% with IMRT1 and 89% with IMRT2. The shorter treatment time and larger dose per fraction obtained with IMRT2 resulted in an 11% increase in the probability of control in the PTV1 with respect to IFP and 7% relatively to IMRT1 (p < 0.05). The average probability of total patient complications was reduced from 80% with CONVT to 61% with IFP and 31% with IMRT. The corresponding probability of complications in the ipsilateral parotid was 63%, 42% and 20%; in the contralateral parotid it was 50%, 20% and 9%; in the oral cavity it was 2%, 15% and 4% and in the mandible it was 1%, 5% and 3%, respectively.ConclusionsA significant improvement in treatment outcome was obtained with IMRT compared to conventional radiation therapy. The practical and biological advantages of IMRT2, employing a shorter treatment time, may outweigh the small differences obtained in the organs at risk between the two IMRT techniques. This technique is therefore presently being used in the clinic for selected patients with head and neck tumours. A significant improvement in the quality of the dose distribution was obtained with IFP compared to CONVT. Thus, this beam arrangement is used in the clinical routine as an alternative to IMRT.

Assessment and topographic characterization of locoregional recurrences in head and neck tumours.

PurposeTo evaluate the differences between three methods of classification of recurrences in patients with head and neck tumours treated with Radiation Therapy (RT).Materials and methods367 patients with head and neck tumours were included in the study. Tumour recurrences were delineated in the CT images taken during patient follow-up and deformable registration was used to transfer this volume into the planning CT. The methods used to classify recurrences were: method CTV quantified the intersection volume between the recurrence and the Clinical Target Volume (CTV); method TV quantified the intersection between the Treated Volume and the recurrence (for method CTV and TV, recurrences were classified in-field if more than 95% of their volume were inside the volume of interest, marginal if the intersection was between 20-95% and outfield otherwise); and method COM was based on the position of the Centre Of Mass of the recurrence. A dose assessment in the recurrence volume was also made.ResultsThe 2-year Kaplan-Meier locoregional recurrence incidence was 10%. Tumour recurrences occurred in 22 patients in a mean time of 16.5 ± 9.4 months resulting in 28 recurrence volumes. The percentage of in-field recurrences for methods CTV, TV and COM was 7%, 43% and 50%, respectively. Agreement between the three methods in characterizing individually in-field and marginal recurrences was found only in six cases. Methods CTV and COM agreed in 14. The percentage of outfield recurrences was 29% using all methods. For local recurrences (in-field or marginal to gross disease) the average difference between the prescribed dose and D98% in the recurrence volume was -5.2 ± 3.5% (range: -10.1%-0.9%).ConclusionsThe classification of in-field and marginal recurrences is very dependent on the method used to characterize recurrences. Using methods TV and COM the largest percentage of tumour recurrences occurred in-field in tissues irradiated with high doses.

Simulated annealing applied to IMRT beam angle optimization: A computational study

Electing irradiation directions to use in IMRT treatments is one of the first decisions to make in treatment planning. Beam angle optimization (BAO) is a difficult problem to tackle from the mathematical optimization point of view. It is highly non-convex, and optimization approaches based on gradient descent methods will probably get trapped in one of the many local minima. Simulated Annealing (SA) is a local search probabilistic procedure that is known to be able to deal with multimodal problems. SA for BAO was retrospectively applied to ten clinical examples of treated cases of head-and neck tumors signalized as complex cases where proper target coverage and organ sparing proved difficult to achieve. The number of directions to use was considered fixed and equal to 5 or 7. It is shown that SA can lead to solutions that significantly improve organ sparing, even considering a reduced number of angles, without jeopardizing tumor coverage.

A graphic user interface toolkit for specification, report and comparison of dose-response relations and treatment plans using the biologically effective uniform dose

A toolkit (BEUDcal) has been developed for evaluating the effectiveness and for predicting the outcome of treatment plans by calculating the biologically effective uniform dose (BEUD) and complication-free tumor control probability. The input for the BEUDcal is the differential dose-volume histograms of organs exported from the treatment planning system. A clinical database is built for the dose-response parameters of different tumors and normal tissues. Dose-response probabilities of all the examined organs are illustrated together with the corresponding BEUDs and the P(+) values. Furthermore, BEUDcal is able to generate a report that simultaneously presents the radiobiological evaluation together with the physical dose indices, showing the complementary relation between the physical and radiobiological treatment plan analysis performed by BEUDcal. Comparisons between treatment plans for helical tomotherapy and multileaf collimator-based intensity modulated radiotherapy of a lung patient were demonstrated to show the versatility of BEUDcal in the assessment and report of dose-response relations.