RESUMEN
In the past decade, several different radiotherapy treatment plan evaluation and optimization schemes have been proposed as viable approaches, aiming for dose escalation or an increase of healthy tissue sparing. In particular, it has been argued that dose-mass plan evaluation and treatment plan optimization might be viable alternatives to the standard of care, which is realized through dose-volume evaluation and optimization. The purpose of this investigation is to apply dose-mass optimization to a cohort of lung cancer patients and compare the achievable healthy tissue sparing to that one achievable through dose-volume optimization. Fourteen non-small cell lung cancer (NSCLC) patient plans were studied retrospectively. The range of tumor motion was less than 0.5 cm and motion management in the treatment planning process was not considered. For each case, dose-volume (DV)-based and dose-mass (DM)-based optimization was performed. Nine-field step-and-shoot IMRT was used, with all of the optimization parameters kept the same between DV and DM optimizations. Commonly used dosimetric indices (DIs) such as dose to 1% the spinal cord volume, dose to 50% of the esophageal volume, and doses to 20 and 30% of healthy lung volumes were used for cross-comparison. Similarly, mass-based indices (MIs), such as doses to 20 and 30% of healthy lung masses, 1% of spinal cord mass, and 33% of heart mass, were also tallied. Statistical equivalence tests were performed to quantify the findings for the entire patient cohort. Both DV and DM plans for each case were normalized such that 95% of the planning target volume received the prescribed dose. DM optimization resulted in more organs at risk (OAR) sparing than DV optimization. The average sparing of cord, heart, and esophagus was 23, 4, and 6%, respectively. For the majority of the DIs, DM optimization resulted in lower lung doses. On average, the doses to 20 and 30% of healthy lung were lower by approximately 3 and 4%, whereas lung volumes receiving 2000 and 3000 cGy were lower by 3 and 2%, respectively. The behavior of MIs was very similar. The statistical analyses of the results again indicated better healthy anatomical structure sparing with DM optimization. The presented findings indicate that dose-mass-based optimization results in statistically significant OAR sparing as compared to dose-volume-based optimization for NSCLC. However, the sparing is case-dependent and it is not observed for all tallied dosimetric endpoints.
Asunto(s)
Carcinoma de Pulmón de Células no Pequeñas/radioterapia , Planificación de la Radioterapia Asistida por Computador/métodos , Humanos , Movimiento , Órganos en Riesgo/efectos de la radiación , Dosificación RadioterapéuticaRESUMEN
The purpose of this work is to evaluate the modeling of carbon fiber couch attenuation properties with a commercial treatment planning system (TPS, Pinnacle3, v8.0d). A carbon fiber couch (Brain-Lab) was incorporated into the TPS by automatic contouring of all transverse CT slices. The couch shape and dimensions were set according to the vendor specifications. The couch composition was realized by assigning appropriate densities to the delineated contours. The couch modeling by the TPS was validated by absolute dosimetric measurements. A phantom consisting of several solid water slabs was CT scanned, the CT data set was imported into the TPS, and the carbon fiber couch was auto-contoured. Open (unblocked) field plans for different gantry angles and field sizes were generated. The doses to a point at 3 cm depth, placed at the linac isocenter, were computed. The phantom was irradiated according to the dose calculation setup and doses were measured with an ion chamber. In addition, percent depth dose (PDD) curves were computed as well as measured with radiographic film. The calculated and measured doses, transmissions, and PDDs were cross-compared. Doses for several posterior fields (0 degree, 30 degrees, 50 degrees, 75 degrees, 83 degrees) were calculated for 6 and 18 MV photon beams. For model validation a nominal field size of 10 x 10 cm2 was chosen and 100 MU were delivered for each portal. The largest difference between computed and measured doses for those posterior fields was within 1.7%. A comparison between computed and measured transmissions for the aforementioned fields was performed and the results were found to agree within 1.1%. The differences between computed and measured doses for different field sizes, ranging from 5 x 5 cm2 to 25 x 25 cm2 in 5 cm increments, were within 2%. Measured and computed PDD curves with and without the couch agree from the surface up to 30 cm depth. The PDDs indicate a surface dose increase resulting from the carbon fiber couch field modification. The carbon fiber couch attenuation for individual posterior oblique fields (75 degrees) can be in excess of 8% depending on the beam energy and field size. When the couch is contoured in Pinnacle3 its attenuation properties are modeled to within 1.7% with respect to measurements. These results demonstrate that appropriate contouring together with relevant density information for the contours is sufficient for adequate modeling of carbon fiber supporting devices by modern commercial treatment planning systems.
Asunto(s)
Carbono/efectos de la radiación , Modelos Químicos , Planificación de la Radioterapia Asistida por Computador/métodos , Radioterapia de Intensidad Modulada/instrumentación , Fibra de Carbono , Dosificación Radioterapéutica , Sensibilidad y Especificidad , Tomografía Computarizada por Rayos XRESUMEN
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.
Asunto(s)
Algoritmos , Neoplasias de Cabeza y Cuello/radioterapia , Método de Montecarlo , Planificación de la Radioterapia Asistida por Computador/métodos , Radioterapia de Intensidad Modulada/métodos , Neoplasias de Cabeza y Cuello/patología , Humanos , Valor Predictivo de las Pruebas , Dosificación RadioterapéuticaRESUMEN
This study quantifies the dose prediction errors (DPEs) in dynamic IMRT dose calculations resulting from (a) use of an intensity matrix to estimate the multi-leaf collimator (MLC) modulated photon fluence (DPE(IGfluence) instead of an explicit MLC particle transport, and (b) handling of tissue heterogeneities (DPE(hetero)) by superposition/convolution (SC) and pencil beam (PB) dose calculation algorithms. Monte Carlo (MC) computed doses are used as reference standards. Eighteen head-and-neck dynamic MLC IMRT treatment plans are investigated. DPEs are evaluated via comparing the dose received by 98% of the GTV (GTV D 98%), the CTV D 95%, the nodal D 90%, the cord and the brainstem D 02%, the parotid D 50%, the parotid mean dose (D (Mean)), and generalized equivalent uniform doses (gEUDs) for the above structures. For the MC-generated intensity grids, DPE(IGfluence) is within +/- 2.1% for all targets and critical structures. The SC algorithm DPE(hetero) is within +/- 3% for 98.3% of the indices tallied, and within +/- 3.4% for all of the tallied indices. The PB algorithm DPE(hetero) is within +/- 3% for 92% of the tallied indices. Statistical equivalence tests indicate that PB DPE(hetero) requires a +/- 3.6% interval to state equivalence with the MC standard, while the intervals are < 1.5% for SC DPE(hetero) and DPE(IGfluence). Overall, these results indicate that SC and MC IMRT dose calculations which use MC-derived intensity matrices for fluence prediction do not introduce significant dose errors compared with full Monte Carlo dose computations; however, PB algorithms may result in clinically significant dose deviations.
Asunto(s)
Algoritmos , Artefactos , Modelos Biológicos , Radiometría/métodos , Planificación de la Radioterapia Asistida por Computador/métodos , Anisotropía , Carga Corporal (Radioterapia) , Simulación por Computador , Modelos Estadísticos , Método de Montecarlo , Especificidad de Órganos , Dosificación Radioterapéutica , Radioterapia Conformacional , Efectividad Biológica Relativa , Reproducibilidad de los Resultados , Dispersión de Radiación , Sensibilidad y Especificidad , Programas InformáticosRESUMEN
A hybrid dose-computation method is designed which accurately accounts for multileaf collimator (MLC)-induced intensity modulation in intensity modulated radiation therapy (IMRT) dose calculations. The method employs Monte Carlo (MC) modeling to determine the fluence modulation caused by the delivery of dynamic or multisegmental (step-and-shoot) MLC fields, and a conventional dose-computation algorithm to estimate the delivered dose to a phantom or a patient. Thus, it determines the IMRT fluence prediction accuracy achievable by analytic methods in the limit that the analytic method includes all details of the MLC leaf transport and scatter. The hybrid method is validated and benchmarked by comparison with in-phantom film dose measurements, as well as dose calculations from two in-house, and two commercial treatment planning system analytic fluence estimation methods. All computation methods utilize the same dose algorithm to calculate dose to a phantom, varying only in the estimation of the MLC modulation of the incident photon energy fluence. Gamma analysis, with respect to measured two-dimensional (2D) dose planes, is used to benchmark each algorithm's performance. The analyzed fields include static and dynamic test patterns, as well as fields from ten DMLC IMRT treatment plans (79 fields) and five SMLC treatment plans (29 fields). The test fields (fully closed MLC, picket fence, sliding windows of different size, and leaf-tip profiles) cover the extremes of MLC usage during IMRT, while the patient fields represent realistic clinical conditions. Of the methods tested, the hybrid method most accurately reproduces measurements. For the hybrid method, 79 of 79 DMLC field calculations have gamma < 1 (3%/3 mm) for more than 95% of the points (per field) while for SMLC fields, 27 of 29 pass the same criteria. The analytic energy fluence estimation methods show inferior pass rates, with 76 of 79 DMLC and 24 of 29 SMLC fields having more than 95% of the test points with gamma < or = 1 (3%/3 mm). Paired one-way ANOVA tests of the gamma analysis results found that the hybrid method better predicts measurements in terms of both the fraction of points with gamma < or = 1 and the average gamma for both 2%/2 mm and 3%/3 mm criteria. These results quantify the enhancement in accuracy in IMRT dose calculations when MC is used to model the MLC field modulation.