Study of lung density corrections in a clinical trial (RTOG 88-08). Radiation Therapy Oncology Group

Accuracy Evaluation of EPL and ETAR Algorithms in the Treatment Planning Systems using CIRS Thorax Phantom

Background:

It is recommended for each set of radiation data and algorithm that subtle deliberation is done regarding dose calculation accuracy. Knowing the errors in dose calculation for each treatment plan will result in an accurate estimate of the actual dose achieved by the tumor.

Objective:

This study aims to evaluate the equivalent path length (EPL) and equivalent tissue air ratio (ETAR) algorithms in radiation dose calculation.

Material and Methods:

In this experimental study, the TEC-DOC 1583 guideline was used. Measurements and calculations were obtained for each algorithm at specific points in thorax CIRS phantom for 6 and 18 MVs and results were compared.

Results:

In the EPL, calculations were in agreement with measurements for 27 points and differences between them ranged from 0.1% to 10.4% at 6 MV. The calculations were in agreement with measurements for 21 points and differences between them ranged from 0.4% to 13% at 18 MV. In ETAR, calculations were also in consistent with measurements for 21 points, and differences between them ranged from 0.1% to 9% at 6 MV. Moreover, for 18 MV, the calculations were in agreement with measurements for 17 points and differences between them ranged from 0% to 11%.

Conclusion:

For the EPL algorithm, more dose points were in consistent with acceptance criteria. The errors in the ETAR were 1% to 2% less than the EPL. The greatest calculation error occurs in low-density lung tissue with inhomogeneities or in high-density bone. Errors were larger in shallow depths. The error in higher energy was more than low energy beam.

The dependence of inhomogeneity correction factors on photon beam quality index performed with the Anisotropic Analytical Algorithm

Purpose: The purpose of the study was to investigate the dependence of tissue inhomogeneity correction factors (ICFs) on the photon beam quality index (QI).

Materials and Methods: Heterogeneous phantoms, comprising semi-infinite slabs of the lung (0.10, 0.20, 0.26 and 0.30 g/cm3), adipose tissue (0.92 g/cm3) and bone (1.85 g/cm3) in water, were constructed in the Eclipse treatment planning system. Several calculation models of 6 MV and 15 MV photon beams for quality index (TPR20,10) = 0.670±k*0.01 and TPR20,10 = 0.760±k*0.01, k = -3, -2, -1, 0, 1, 2, 3 respectively were built in the Eclipse. The ICFs were calculated with the anisotropic analytical algorithm (AAA) for several beam sizes and points lying at several depths inside of and below inhomogeneities of different thicknesses.

Results: The ICFs increased for lung and adipose tissues with increasing beam quality (TPR20,10), while decreased for bone. Calculations with AAA predict that the maximum difference in ICFs of 1.0% and 2.5% for adipose and bone tissues, respectively. For lung tissue, changes of ICFs of a maximum of 9.2% (6 MV) and 13.8% (15 MV). For points where charged particle equilibrium exists, a linear dependence of ICFs on TPR20,10 was observed. If CPE doesn’t exist, the dependence became more complex. For points inside of the low-density inhomogeneity, the dependence of the ICFs on energy was not linear but the changes of ICFs were smaller than 3.0%. Measurements results carried out with the CIRS phantom were consistent with the calculation results.

Conclusions: A negligible dependence of the ICFs on energy was found for adipose and bone tissue. For lung tissue, in the CPE region, the dependence of ICFs on different beam quality indexes with the same nominal energy may not be neglected, however, this dependence was linear. Where there is no CPE, the dependence of the ICFs on energy was more complicated.

Evaluation of Lung Density and Its Dosimetric Impact on Lung Cancer Radiotherapy: A Simulation Study

BACKGROUND

The dosimetric parameters required in lung cancer radiation therapy are taken from a homogeneous water phantom; however, during treatment, the expected results are being affected because of its inhomogeneity. Therefore, it becomes necessary to quantify these deviations.

OBJECTIVE

The present study has been undertaken to find out inter- and intra- lung density variations and its dosimetric impact on lung cancer radiotherapy using Monte Carlo code FLUKA and PBC algorithms.

MATERIAL AND METHODS

Density of 100 lungs was recorded from their CT images along with age. Then, after PDD calculated by FLUKA MC Code and PBC algorithm for virtual phantom having density 0.2 gm/cm3 and 0.4 gm/cm3 (density range obtained from CT images of 100 lungs) using Co-60 10 x10 cm2 beams were compared.

RESULTS

Average left and right lung densities were 0.275±0.387 and 0.270±0.383 respectively. The deviation in PBC calculated PDD were (+)216%, (+91%), (+)45%, (+)26.88%, (+)14%, (-)1%, (+)2%, (-)0.4%, (-)1%, (+)1%, (+)4%, (+)4.5% for 0.4 gm/cm3 and (+)311%, (+)177%, (+)118%, (+)90.95%, (+)72.23%, (+)55.83% ,(+)38.85%, (+)28.80%, (+)21.79%, (+)15.95%, (+)1.67%, (-) 2.13%, (+)1.27%, (+)0.35%, (-)1.79%, (-)2.75% for 0.2 gm/cm3 density mediums at depths of 1mm, 2mm, 3mm, 4mm, 5mm, 6 mm, 7 mm, 8mm, 9mm,10mm, 15mm, 30mm, 40mm, 50mm, 80mm and 100 mm, respectively.

CONCLUSION

Large variations in inter- and intra- lung density were recorded. PBC overestimated the dose at air/lung interface as well as inside lung. The results of Monte Carlo simulation can be used to assess the performance of other treatment planning systems used in lung cancer radiotherapy.

Lung Cancer Radiation Therapy: Monte Carlo Investigation of “Under Dose” by High Energy Photons

Loss of electronic equilibrium in lung tissue causes a build-up region in the tumor. Increasing the photon energy increases the depth at which electronic equilibrium is reestablished within the lung tumor. This study uses the Monte Carlo code PENELOPE for simulations of radiation treatment of tumor surrounded by lung. Six MV photons were compared to 15 MV photons using four beam arrangements in both homogeneous and heterogeneous media. The experimental results demonstrate that for every beam arrangement in heterogeneous media 15 MV photons delivered 5% to 10% lower dose to the tumor periphery than 6 MV photons. The simulations also show that in axial coplanar treatment plans, the loss of electronic equilibrium was greatest in the coronal plane. In conclusion there is a tumor sparing effect at the tumor-lung interface that is a function of beam energy. As an alternative to increasing beam energy, the addition of multiple beam angles with lower energy photons improved target coverage. If higher energy beams are required for patients with large separation, then adding multiple beam angles does offer some improved target coverage. The non-coplanar technique with the lower energy photons covered the tumor with a greatest isodose at the tumor periphery without tangential sparing in the coronal plane.

Dosimetric verification of radiotherapy treatment planning systems in Serbia: national audit

Background

Independent external audits play an important role in quality assurance programme in radiation oncology. The audit supported by the IAEA in Serbia was designed to review the whole chain of activities in 3D conformal radiotherapy (3D-CRT) workflow, from patient data acquisition to treatment planning and dose delivery. The audit was based on the IAEA recommendations and focused on dosimetry part of the treatment planning and delivery processes.

Methods

The audit was conducted in three radiotherapy departments of Serbia. An anthropomorphic phantom was scanned with a computed tomography unit (CT) and treatment plans for eight different test cases involving various beam configurations suggested by the IAEA were prepared on local treatment planning systems (TPSs). The phantom was irradiated following the treatment plans for these test cases and doses in specific points were measured with an ionization chamber. The differences between the measured and calculated doses were reported.

Results

The measurements were conducted for different photon beam energies and TPS calculation algorithms. The deviation between the measured and calculated values for all test cases made with advanced algorithms were within the agreement criteria, while the larger deviations were observed for simpler algorithms. The number of measurements with results outside the agreement criteria increased with the increase of the beam energy and decreased with TPS calculation algorithm sophistication. Also, a few errors in the basic dosimetry data in TPS were detected and corrected.

Conclusions

The audit helped the users to better understand the operational features and limitations of their TPSs and resulted in increased confidence in dose calculation accuracy using TPSs. The audit results indicated the shortcomings of simpler algorithms for the test cases performed and, therefore the transition to more advanced algorithms is highly desirable.

Accurate heterogeneous dose calculation for lung cancer patients without high-resolution CT densities

The aim of this study was to investigate the relative accuracy of megavoltage photon-beam dose calculations employing either 5 bulk densities or independent voxel densities determined by calibration of the CT Houndsfield number. Full-resolution CT and bulk density treatment plans were generated for 70 lung or esophageal cancer tumors (66 cases) using a commercial treatment planning system with an adaptive convolution dose calculation algorithm (Pinnacle3, Philips Medicals Systems). Bulk densities were applied to segmented regions. Individual and population average densities were compared to the full-resolution plan for each case. Monitor units were kept constant and no normalizations were employed. Dose volume histograms (DVH) and dose difference distributions were examined for all cases. The average densities of the segmented air, lung, fat, soft tissue, and bone for the entire set were found to be 0.14, 0.26, 0.89, 1.02, and 1.12 g/cc, respectively. In all cases, the normal tissue DVH agreed to better than 2% in dose. In 62 of 70 target DVHs, agreement to better than 3% in dose was observed. Six cases demonstrated emphysema, one with bullous formations and one with a hiatus hernia having a large volume of gas. They required the additional assignment of density to the emphysemic lung and inflammatory changes to the lung, the regions of collapsed lung, the bullous formations, and the hernia gas. Bulk tissue density dose calculation provides an accurate method of heterogeneous dose calculation. However, patients with advanced emphysema may require high-resolution CT studies for accurate treatment planning.

Inhomogeneity correction and the analytic anisotropic algorithm

The ability of the analytic anisotropic algorithm (AAA), a superposition- convolution algorithm implemented in the Eclipse (Varian Medical Systems, Palo Alto, CA) treatment planning system (TPS), to accurately account for the presence of inhomogeneities in simple geometries is examined. The goal of 2% accuracy, as set out by the American Association of Physicists in Medicine Task Group 65, serves as a useful benchmark against which to evaluate the inhomogeneity correction capabilities of this treatment planning algorithm. A planar geometry phantom consisting of upper and lower layers of Solid Water (Gammex rmi, Middleton, WI) separated by a heterogeneity region of variable thickness, is modeled within the Eclipse TPS. Results obtained with the AAA are compared with experimental measurements. Seven different materials, spanning the range from air to aluminum, constitute the inhomogeneity layer. In general, the AAA overpredicts dose beyond low-density regions and underpredicts dose distal to volumes of high density. In many cases, the deviation between the AAA and experimental results exceeds the Task Group 65 target of 2%. The source of these deviations appears to arise from an inability of the AAA to correctly account for altered attenuation along primary ray paths.

Open low-field magnetic resonance imaging for target definition, dose calculations and set-up verification during three-dimensional CRT for glioblastoma multiforme

AIMS

To assess the effect on target delineation of using magnetic resonance simulation for planning of glioblastoma multiforme (GBM). Dose calculations derived from computed tomography- and magnetic resonance-derived plans were computed. The accuracy of set-up verification using magnetic resonance imaging (MRI)-based digital reconstructed radiographs (DRRs) was assessed.

MATERIALS AND METHODS

Ten patients with GBM were simulated using computed tomography and MRI. MRI was acquired with a low-field (0.23 T) MRI unit (SimMRI). Gross tumour volumes (GTVs) were delineated by two radiation oncologists on computed tomography and MRI. In total, 30 plans were generated using both the computed tomography, with (planbathoCT) and without (planCT) heterogeneity correction, and MRI data sets (planSimMRI). The minimum dose delivered (Dmin) to the GTV between computed tomography- and MRI-based plans was compared. The accuracy of set-up positioning using MRI DRRs was assessed by four radiation oncologists.

RESULTS

The mean GTVs delineated on computed tomography were significantly (P<0.001) larger than those contoured on MRI. The mean (+/-standard deviation) Dmin difference percentage was 0.3+/-0.8, 0.1+/-0.6 and -0.2+/-1.0% for the planCT/planbathoCT-, planCT/planSimMRI- and planbathoCT/planSimMRI-derived plans, respectively. The set-up differences observed with the computed tomography and MRI DRRs ranged from 1.0 to 4.0 mm (mean 1.5 mm; standard deviation+/-1.4).

CONCLUSIONS

GTVs defined on computed tomography were significantly larger than those delineated on MRI. Compared with computed tomography-derived plans, MRI-based dose calculations were accurate. The precision of set-up verifications based on computed tomography- and MRI-derived DRRs seemed similar. The use of MRI only for the planning of GBM should be further assessed.

An anatomically realistic lung model for Monte Carlo-based dose calculations

Treatment planning for disease sites with large variations of electron density in neighboring tissues requires an accurate description of the geometry. This self-evident statement is especially true for the lung, a highly complex organ having structures with a wide range of sizes that range from about 10(-4) to 1 cm. In treatment planning, the lung is commonly modeled by a voxelized geometry obtained using computed tomography (CT) data at various resolutions. The simplest such model, which is often used for QA and validation work, is the atomic mix or mean density model, in which the entire lung is homogenized and given a mean (volume-averaged) density. The purpose of this paper is (i) to describe a new heterogeneous random lung model, which is based on morphological data of the human lung, and (ii) use this model to assess the differences in dose calculations between an actual lung (as represented by our model) and a mean density (homogenized) lung. Eventually, we plan to use the random lung model to assess the accuracy of CT-based treatment plans of the lung. For this paper, we have used Monte Carlo methods to make accurate comparisons between dose calculations for the random lung model and the mean density model. For four realizations of the random lung model, we used a single photon beam, with two different energies (6 and 18 MV) and four field sizes (1 x 1, 5 x 5, 10 x 10, and 20 x 20 cm2). We found a maximum difference of 34% of D(max) with the 1 x 1, 18 MV beam along the central axis (CAX). A "shadow" region distal to the lung, with dose reduction up to 7% of D(max), exists for the same realization. The dose perturbations decrease for larger field sizes, but the magnitude of the differences in the shadow region is nearly independent of the field size. We also observe that, compared to the mean density model, the random structures inside the heterogeneous lung can alter the shape of the isodose lines, leading to a broadening or shrinking of the penumbra region. For small field sizes, the mean lung doses significantly depend on the structures' relative locations to the beam. In addition to these comparisons between the random lung and mean density models, we also provide a preliminary comparison between dose calculations for the random lung model and a voxelized version of this model at 0.4 x 0.4 x 0.4 cm3 resolution. Overall, this study is relevant to treatment planning for lung tumors, especially in situations where small field sizes are used. Our results show that for such situations, the mean density model of the lung is inadequate, and a more accurate CT model of the lung is required. Future work with our model will involve patient motion, setup errors, and recommendations for the resolution of CT models.

The impact of breathing motion versus heterogeneity effects in lung cancer treatment planning

The purpose of this study is to investigate the effects of tissue heterogeneity and breathing-induced motion/deformation on conformal treatment planning for pulmonary tumors and to compare the magnitude and the clinical importance of changes induced by these effects. Treatment planning scans were acquired at normal exhale/inhale breathing states for fifteen patients. The internal target volume (ITV) was defined as the union of exhale and inhale gross tumor volumes uniformly expanded by 5 mm. Anterior/posterior opposed beams (AP/PA) and three-dimensional (3D)-conformal plans were designed using the unit-density exhale ("static") dataset. These plans were further used to calculate (a) density-corrected ("heterogeneous") static dose and (b) heterogeneous cumulative dose, including breathing deformations. The DPM Monte Carlo code was used for dose computations. For larger than coin-sized tumors, relative to unit-density plans, tumor and lung doses increased in the heterogeneity-corrected plans. In comparing cumulative and static plans, larger normal tissue complication probability changes were observed for tumors with larger motion amplitudes and uncompensated breathing-induced hot/cold spots in lung. Accounting for tissue heterogeneity resulted in average increases of 9% and 7% in mean lung dose (MLD) for the 6 MV and 15 MV photon beams, respectively. Breathing-induced effects resulted in approximately 1% and 2% average decreases in MLD from the static value, for the 6 and 15 MV photon beams, respectively. The magnitude of these effects was not found to correlate with the treatment plan technique, i.e., AP/PA versus 3D-CRT. Given a properly designed ITV, tissue heterogeneity effects are likely to have a larger clinical significance on tumor and normal lung treatment evaluation metrics than four-dimensional respiratory-induced changes.

Predicting changes in dose distribution to tumor and normal tissue when correcting for heterogeneity in radiotherapy for lung cancer

PURPOSES

The purposes of this study were to examine dose alterations to gross tumor volume (GTV) and lung using heterogeneity corrections and to predict the magnitude of these changes.

METHODS

Three separate conformal plans were generated for 37 patients with lung cancer: plan 1 corrected for heterogeneity, plan 2 did not correct for heterogeneity, and plan 3 used identical beams and monitor units from plan 2 but with heterogeneous calculations. Plans 1 and 2 were normalized to the 95% isodose line. Mean dose (MeanDGTV), maximum dose (MaxDGTV), and minimum dose (MinDGTV) to GTV and V20 were compared between plans 1 and 3. For each patient, the amount of lung in all beam paths of plan 3 was quantified by a density correction factor and correlated with the percent change.

RESULTS

The median percent change in MeanDGTV, MaxDGTV, and MinDGTV between plan 3 and plan 1 was -4.7% (-0.1% to -19.1%, P < 0.0001), -5.59% (0.16% to -31.86%, P < 0.0001), and -4.88% (2.90% to -24.88%, P < 0.0001), respectively. The median V20 difference was -1% (1% to -8%). The density correction factor correlated with larger differences in MeanDGTV on univariate analysis.

CONCLUSIONS

Heterogeneity correction lowers the dose to GTV by 5%. This difference can be correlated with the density correction factor.

Retrospective monte carlo dose calculations with limited beam weight information

An important unresolved issue in outcomes analysis for lung complications is the effect of poor or completely lacking heterogeneity corrections in previously archived treatment plans. To estimate this effect, we developed a novel method based on Monte Carlo (MC) dose calculations which can be applied retrospectively to RTOG/AAPM-style archived treatment plans (ATP). We applied this method to 218 archived nonsmall cell lung cancer lung treatment plans that were originally calculated either without heterogeneity corrections or with primitive corrections. To retrospectively specify beam weights and wedges, beams were broken into Monte Carlo-generated beamlets, simulated using the VMC++ code, and mathematical optimization was used to match the archived water-based dose distributions. The derived beam weights (and any wedge effects) were then applied to Monte Carlo beamlets regenerated based on the patient computed tomography densities. Validation of the process was performed against five comparable lung treatment plans generated using a commercial convolution/superposition implementation. For the application here (normal lung, esophagus, and planning target volume dose distributions), the agreement was very good. Resulting MC and convolution/superposition values were similar when dose distributions without heterogeneity corrections or dose distributions with corrections were compared. When applied to the archived plans (218), the average absolute percent difference between water-based MC and water-based ATPs, for doses above 2.5% of the maximum dose was 1.8+/-0.6%. The average absolute percent difference between heterogeneity-corrected MC and water-based ATPs increased to 3.1+/-0.9%. The average absolute percent difference between the MC heterogeneity-corrected and the ATP heterogeneity-corrected dose distributions was 3.8+/-1.6% (available in 132/218 archives). The entire dose-volume-histograms for lung, tumor, and esophagus from the different calculation methods, as well as specific dose metrics, were compared. The average difference in maximum lung dose between water-based ATPs and heterogeneity-corrected MC dose distributions was -1.0+/-2.1 Gy. Potential errors in relying on primitive heterogeneity corrections are most evident from a comparison of maximum lung doses, for which the average MC heterogeneity-corrected values were 5.3+/-2.8 Gy less than the ATP heterogeneity-corrected values. We have demonstrated that recalculation of archived dose distributions, without explicit information about beam weights or wedges, is feasible using beamlet-based optimization methods. The method provides heterogeneity-corrected dose data consistent with convolution-superposition calculations and is one feasible approach for improving dosimetric data for outcomes analyses.

Current practice when treating lung cancer in Australasia

A multidisciplinary meeting was held by the radiation oncology department of South Western Sydney Area Cancer Services in March 2003. This meeting was advertised in all radiation oncology departments in Australia and New Zealand. As a precursor to this meeting, a survey was undertaken on the use of radiotherapy for treating lung cancer. All departments in Australia and New Zealand were asked to participate. The survey considered planning techniques, delivery set-up and prescription doses for non-small-cell and small-cell lung cancer and palliative and radical treatments. A wide range in the techniques used was seen across departments, particularly when prescription doses and fractionation were considered.

Monte Carlo dose calculations for a 6-MV photon beam in a thorax phantom

PURPOSE

In this study we evaluated the accuracy of the Monte Carlo (MC) and effective path length (EPL) methods for dose calculations in the inhomogeneous thorax phantom.

MATERIALS AND METHODS

The Philips SL 75/5 linear accelerator head was modeled using the MCNP4C Monte Carlo code. An anatomic inhomogeneous thorax phantom was irradiated with a 6-MV photon beam, and the doses along points of the central axis of the beam were measured by a small ionization chamber. The central axis relative dose was calculated by the MCNP4C code and the EPL method in a conventional treatment planning system. The results of calculations and measurements were compared.

RESULTS

For all measured points on the thorax phantom the results of the MC method were in agreement with the actual measurement (local difference was less than 2%). For the EPL method, the amount of error was dependent on the field size and the point location in the phantom. The maximum error was +19.5 and +26.8 for field sizes of 10 x 10 and 5 x 5 cm2 for lateral irradiation.

CONCLUSION

Our study showed large, unacceptable errors for EPL calculations in the lung for both field sizes. The accuracy of the MC method was better than the recommended value of 3%. Thus, application of this method is strongly recommended for lung dose calculations, especially for small field sizes.

Determination of CT-to-density conversion relationship for image-based treatment planning systems

The implementation of tissue inhomogeneity correction in image-based treatment planning will improve the accuracy of radiation dose calculations for patients undergoing external-beam radiotherapy. Before the tissue inhomogeneity correction can be applied, the relationship between the computed tomography (CT) value and density must be established. This tissue characterization relationship allows the conversion of CT value in each voxel of the CT images into density for use in the dose calculations. This paper describes the proper procedure of establishing the CT value to density conversion relationship. A tissue characterization phantom with 17 inserts made of different materials was scanned using a GE Lightspeed Plus CT scanner (120 kVp). These images were then downloaded into the Eclipse and Pinnacle treatment planning systems. At the treatment planning workstation, the axial images were retrieved to determine the CT value of the inserts. A region of interest was drawn on the central portion of the insert and the mean CT value and its standard deviation were determined. The mean CT value was plotted against the density of the tissue inserts and fitted with bilinear equations. A new set of CT values vs. densities was generated from the bilinear equations and then entered into the treatment planning systems. The need to obtain CT values through the treatment planning system is very clear. The 2 treatment planning systems use different CT value ranges, one from -1024 to 3071 and the other from 0 to 4096. If the range is correct, it would result in inappropriate use of the conversion curve. In addition to the difference in the range of CT values, one treatment planning system uses physical density, while the other uses relative electron density.

Phase III study of the Eastern Cooperative Oncology Group (ECOG 2597): induction chemotherapy followed by either standard thoracic radiotherapy or hyperfractionated accelerated radiotherapy for patients with unresectable stage IIIA and B non-small-cell lung cancer

PURPOSE

To compare once-daily radiation therapy (qdRT) with hyperfractionated accelerated radiation therapy (HART) after two cycles of induction chemotherapy.

PATIENTS AND METHODS

Eligible patients were treatment naive, and had stage IIIA and B unresectable non-small-cell lung cancer, Eastern Cooperative Oncology Group performance status 0/1, and normal organ function. Induction chemotherapy consisted of two cycles of carboplatin area under time-concentration curve 6 mg/mL . min plus paclitaxel 225 mg/m2 on day 1. RT consisted of arm 1 (qdRT), 64 Gy (2 Gy/d), versus arm 2 (HART), 57.6 Gy (1.5 Gy tid for 2.5 weeks). A total of 388 patients were needed to detect a 50% increase in median survival from 14 months of qdRT to 21 months of HART; accrual was not achieved and the study closed prematurely.

RESULTS

Of 141 patients enrolled, 83% were randomly assigned after chemotherapy to qdRT (n = 59) or HART (n = 60). Median survival was 20.3 and 14.9 months for HART and qdRT, respectively (P = .28). Overall response was 25% and 22% for HART and qdRT, respectively (P = .69). Two- and 3-year survival was 44% and 34% for HART, and 24% and 14% for qdRT, respectively. Grade > or = 3 toxicities included esophagitis in 14 v nine patients, and pneumonitis in 0 v 6 patients for HART and qdRT, respectively. Any subsequent trials of the HART regimen must address the issues that led to early closure, including slow accrual, logistics of HART, mucosal toxicity, and the fact that concurrent chemoradiotherapy now seems more effective than sequential treatment.

CONCLUSION

After two cycles of induction chemotherapy with carboplatin-paclitaxel, HART is feasible with an acceptable toxicity profile. Although statistical significance was not achieved and the study closed early, there was a positive statistical trend suggesting a survival advantage with the HART regimen.

Dose-volumetric parameters for predicting severe radiation pneumonitis after three-dimensional conformal radiation therapy for lung cancer

PURPOSE

To retrospectively evaluate dose-volumetric parameters for association with risk of severe (grade >/=3) radiation pneumonitis (RP) in patients after three-dimensional (3D) conformal radiation therapy for lung cancer.

MATERIALS AND METHODS

The study was approved by the institutional review board, which did not require informed consent. Data from 76 patients (66 men, 10 women; median age, 60 years; range, 35-79 years) with histologically proved lung cancer treated curatively with 3D conformal radiation therapy between August 2001 and October 2002 were retrospectively analyzed. Twenty patients underwent surgery before radiation therapy; 57 patients received chemotherapy. Median total radiation dose of 60 Gy (range, 54-66 Gy) was delivered in 30 (range, 27-33) fractions over 6 weeks. RP was scored by using Radiation Therapy Oncology Group criteria. Clinical parameters were analyzed. Dose-volumetric parameters analyzed were percentage of lung volume that received a dose of 20 Gy or more (V20), 30 Gy or more (V30), 40 Gy or more (V40), or 50 Gy or more (V50); mean lung dose (MLD); normal tissue complication probability (NTCP); and total dose. Fisher exact test was performed to compare clinical parameters between patients who developed severe RP and those who did not. Univariate and multivariate logistic regression analyses were performed to evaluate data for association between dose-volumetric parameters and severe RP. Pearson chi(2) test was used to assess data for correlations among dose-volumetric parameters. P < or = .05 was considered to indicate statistically significant difference.

RESULTS

Of 76 patients, 30 (39%) did not develop RP; 23 (30%) developed RP of grade 1; 11 (14%), grade 2; 11 (14%), grade 3; and 1 (1%), grade 4. None had grade 5 RP. Age (< 60 vs > or =60), sex, Karnofsky performance status (< 70 vs > or =70), forced expiratory volume in 1 second, presence of weight loss, preexisting lung disease, history of thoracic surgery, and history of chemotherapy did not significantly differ between patients who developed severe RP and those who did not. In univariate analyses, MLD, V20, V30, V40, V50, and NTCP were associated with severe RP (P < .05). In multivariate analysis, MLD was the only variable associated with severe RP.

CONCLUSION

MLD is a useful indicator of risk for development of severe RP after 3D conformal radiation therapy in patients with lung cancer.

Incorporating an improved dose-calculation algorithm in conformal radiotherapy of lung cancer: re-evaluation of dose in normal lung tissue

BACKGROUND AND PURPOSE

The low density of lung tissue causes a reduced attenuation of photons and an increased range of secondary electrons, which is inaccurately predicted by the algorithms incorporated in some commonly available treatment planning systems (TPSs). This study evaluates the differences in dose in normal lung tissue computed using a simple and a more correct algorithm. We also studied the consequences of these differences on the dose-effect relations for radiation-induced lung injury.

MATERIALS AND METHODS

The treatment plans of 68 lung cancer patients initially produced in a TPS using a calculation model that incorporates the equivalent-path length (EPL) inhomogeneity-correction algorithm, were recalculated in a TPS with the convolution-superposition (CS) algorithm. The higher accuracy of the CS algorithm is well-established. Dose distributions in lung were compared using isodoses, dose-volume histograms (DVHs), the mean lung dose (MLD) and the percentage of lung receiving >20 Gy (V20). Published dose-effect relations for local perfusion changes and radiation pneumonitis were re-evaluated.

RESULTS

Evaluation of isodoses showed a consistent overestimation of the dose at the lung/tumor boundary by the EPL algorithm of about 10%. This overprediction of dose was also reflected in a consistent shift of the EPL DVHs for the lungs towards higher doses. The MLD, as determined by the EPL and CS algorithm, differed on average by 17+/-4.5% (+/-1SD). For V20, the average difference was 12+/-5.7% (+/-1SD). For both parameters, a strong correlation was found between the EPL and CS algorithms yielding a straightforward conversion procedure. Re-evaluation of the dose-effect relations showed that lung complications occur at a 12-14% lower dose. The values of the TD(50) parameter for local perfusion reduction and radiation pneumonitis changed from 60.5 and 34.1 Gy to 51.1 and 29.2 Gy, respectively.

CONCLUSIONS

A simple tissue inhomogeneity-correction algorithm like the EPL overestimates the dose to normal lung tissue. Dosimetric parameters for lung injury (e.g. MLD, V20) computed using both algorithms are strongly correlated making an easy conversion feasible. Dose-effect relations should be refitted when more accurate dose data is available.

Treatment planning for lung cancer: traditional homogeneous point-dose prescription compared with heterogeneity-corrected dose-volume prescription

PURPOSE

To quantify the differences in doses to target volumes and critical thoracic structures calculated by traditional homogeneous point-dose prescription and heterogeneity-corrected volume-dose prescription.

METHODS AND MATERIALS

Between 1998 and 2001, 30 patients with inoperable Stage I/II non-small-cell lung cancer underwent radiation treatment planning at our institution. A commercially available convolution/superposition- based algorithm was used. Three treatment plans were calculated for each patient using identical beam geometries: one plan was generated by traditional homogeneous point-dose prescription, a second by the traditional method with heterogeneity correction, and a third by heterogeneity-corrected volume-dose prescription that would cover 95% of the planned target volume (PTV). Target volume coverage, isocenter dose, and dose uniformity in the second and third plans were compared.

RESULTS

The PTV, clinical target volume (CTV), and isocenter calculated by the heterogeneity-corrected volume-dose method were equivalent to those calculated by the traditional homogeneous point-dose method with heterogeneity correction. The fraction of the PTV covered by heterogeneity-corrected volume-dose prescription was significantly greater than the fraction covered by traditional homogeneous point-dose prescription with heterogeneity correction (p = 0.05). The dose prescribed using the traditional method would have been delivered to less than 90% of the PTV in 14 of 30 patients. There was no significant difference in the maximum and minimum doses to the PTV, the CTV, or the isocenter calculated by the traditional homogeneous method with heterogeneity correction and the heterogeneity-corrected volume-dose method. There was also no significant difference in the planned volume of lung receiving greater than 20 Gy as calculated by these two methods.

CONCLUSION

When compared with traditional homogeneous radiation treatment planning, heterogeneity-corrected methods produce equivalent PTV, CTV, and isocenter doses while providing superior PTV coverage.

Clinical implications of incorporating heterogeneity corrections in mantle field irradiation

PURPOSE

Patient dose calculations for mantle-field irradiation have traditionally been performed using homogeneous, water phantom data. The advent of computed tomography (CT)-based treatment planning now permits dose calculations to be corrected for actual patient density. Incorporation of full heterogeneity corrections is desirable, because calculations performed in this fashion more closely represent the actual dose delivered to the patient. In preparation for full clinical implementation of heterogeneity corrections in mantle irradiation, an evaluation of possible changes in dosimetry when transitioning from treatment plans generated without heterogeneity corrections to treatment plans that incorporated full heterogeneity corrections is presented.

MATERIALS AND METHODS

A retrospective analysis was performed of treatment plans with and without heterogeneity corrections for 15 consecutive patients who had undergone full mantle-field irradiation. Comparisons were made of the absolute delivered doses (in cGy per monitor unit) and the absolute volume (in cubic centimeters) enclosed by the isodose surface of the 30.6 Gy prescription line and the surface representing 90% of the prescribed dose. Dose-volume histograms (DVHs) were generated and studied to evaluate differences in the doses received by the lungs, heart, and spinal cord between corrected and uncorrected plans. Comparisons were made of the volumes of lung receiving at least 20 Gy, the volumes of heart receiving at least 25.2 Gy, and the maximum cord dose.

RESULTS

Dosimetric differences between heterogeneity-corrected and heterogeneity-uncorrected calculations were small. The mean total ratio of corrected-to-uncorrected dose per monitor unit was 1.01, with a standard deviation (SD) of 0.02. The mean corrected-to-uncorrected treated volume ratio (30.6 Gy) was 0.97, SD 0.14, and the mean corrected-to-uncorrected volume ratio of the 90% isodose surface was 0.99, SD 0.02. The ratio of the volume of lung receiving at least 20 Gy was 1.03, SD 0.02; the ratio of the volume of heart receiving at least 25.2 Gy was 1.01, SD 0.03; and the maximum spinal cord dose ratio was 1.02, SD 0.02.

CONCLUSIONS

In all patient treatment plans evaluated, no significant dosimetric differences were observed between heterogeneity-corrected and heterogeneity-uncorrected treatment plans. However, unpredictable differences in the prescription isodose (30.6 Gy) were observed. The differences in coverage at the 90% isodose volume were negligible. The dose administered to lung in heterogeneity-corrected plans demonstrates a higher dose overall, with the greatest increase occurring at volumes receiving at least 20 Gy. In light of these small dosimetric differences, we believe that heterogeneity corrections can be incorporated into full mantle-field treatment planning.

Evaluation of deep inspiration breath-hold lung treatment plans with Monte Carlo dose calculation

PURPOSE

To evaluate dosimetry of deep inspiration breath-hold (DIBH) relative to free breathing (FB) for three-dimensional conformal radiation therapy of lung cancer with 6-MV photons and Monte Carlo (MC) dose calculations.

METHODS AND MATERIALS

Static three-dimensional conformal radiation therapy, 6-MV plans, based on DIBH and FB CT images for five non-small-cell lung cancer patients, were generated on a clinical treatment planning system with equivalent path length tissue inhomogeneity correction. Margins of gross to planning target volume were not reduced for DIBH plans. Cord and lung toxicity determined the maximum treatment dose for each plan. Dose distributions were recalculated for the same beams with an MC dose calculation algorithm and electron density distributions derived from the CT images.

RESULTS

MC calculations showed decreased target coverage relative to treatment-planning system predictions. Lateral disequilibrium caused more degradation of target coverage for DIBH than for FB (approximately 4% worse than expected for FB vs. 8% for DIBH). However, with DIBH higher treatment doses could be delivered without violating normal tissue constraints, resulting in higher total doses to gross target volume and to >99% of planning target volume.

CONCLUSIONS

If DIBH enables prescription dose increases exceeding 10%, MC calculations indicate that, despite lateral disequilibrium, higher doses will be delivered to medium-to-large, partly mediastinal gross target volumes, providing that 6-MV photons are used and margins are not reduced.

CT based lung density correction verification with in vivo dosimetry using diodes

In vivo dose measurements with diodes are easy to perform. The first aim of our study was to show whether diode measurements of the patient exit doses are precise enough for verifying inhomogeneity corrections used for treatment planning. The second aim was to assess the precision of the modified Batho Law inhomogeneity correction of the CadPlan treatment planning system. For this purpose, entrance and lait doses were measured in the thoracic region of 115 patients. Diode measurements were sufficiently precise to verify the density corrections predicted by the treatment planning system (< 0.5% of ICRU dose). The measured doses were compared with calculations of the CadPlan treatment planning system. The mean deviation of the exit dose calculations within the measurements error was zero. The present results show that measurements of exit dose even in a small number of patients are sufficient to identify systematic errors in the dose calculation.

Impact of simple tissue inhomogeneity correction algorithms on conformal radiotherapy of lung tumours

BACKGROUND AND PURPOSE

Conformal radiotherapy requires accurate dose calculation at the dose specification point, at other points in the planning target volume (PTV) and in organs at risk. To assess the limitations of treatment planning of lung tumours, errors in dose values, calculated by some simple tissue inhomogeneity correction algorithms available in a number of currently applied treatment planning systems, have been quantified.

MATERIALS AND METHODS

Single multileaf collimator-shaped photon beams of 6, 8, 15 and 18 MV nominal energy were used to irradiate a 50 mm diameter spherical solid tumour, simulated by polystyrene, which was located centrally inside lung tissue, simulated by cork. The planned dose distribution was made conformal to the PTV, which was a 15 mm three-dimensional expansion of the tumour. Values of both the absolute dose at the International Commission on Radiation Units and Measurement (ICRU) reference point and relative dose distributions inside the PTV and in the lung were calculated using three inhomogeneity correction algorithms. The algorithms investigated in this study are the pencil beam algorithm with one-dimensional corrections, the modified Batho algorithm and the equivalent path length algorithm. The calculated data were compared with measurements for a simple beam set-up using radiographic film and ionization chambers.

RESULTS

For this specific configuration, deviations of up to 3.5% between calculated and measured values of the dose at the ICRU reference point were found. Discrepancies between measured and calculated beam fringe values (distance between the 50 and 90% isodose lines) of up to 14 mm have been observed. The differences in beam fringe and penumbra width (20-80%) increase with increasing beam energy. Our results demonstrate that an underdosage of the PTV up to 20% may occur if calculated dose values are used for treatment planning. The three algorithms predict a considerably higher dose in the lung, both along the central beam axis and in the lateral direction, compared with the actual delivered dose values.

CONCLUSIONS

The dose at the ICRU reference point of such a tumour in lung geometry is calculated with acceptable accuracy. Differences between calculated and measured dose distributions are primarily due to changes in electron transport in the lung, which are not adequately taken into account by the simple tissue inhomogeneity correction algorithms investigated in this study. Particularly for high photon beam energies, clinically unacceptable errors will be introduced in the choice of field sizes employed for conformal treatments, leading to underdosage of the PTV. In addition, the dose to the lung will be wrongly predicted which may influence the choice of the prescribed dose level in dose-escalation studies.

Analytical calculation of the portal scatter to primary dose ratio: an EGS4 Monte Carlo and experimental validation at large air gaps

An analytical approximation for the scatter to primary dose ratio (SPR) on the central axis was validated against Monte Carlo results and experimental measurements for homogeneous and inhomogeneous phantoms. The analytical approximation only included first-order Compton scatter. The contribution to the total SPR from first-order Compton scatter, multiply scattered photons and electron scatter was investigated using Monte Carlo simulation for homogeneous phantoms (up to 30 cm thick for 6 and 18 MV beams; source to detector distances from 150 to 230 cm) as well as for a neck, thorax and pelvis phantom. SPRs were measured on the central axis with an ionization chamber for water phantoms (up to 20 cm thick at 4 MV, 30 cm for 6 MV and 10 MV and 40 cm for 18 MV; source to detector distances of 185 and 200 cm) and for phantoms representing the neck, thorax and pelvis (for air gaps of 50 cm and larger). The mean difference between the experimental and analytical SPRs on the central axis for source to detector distances of 170 cm or greater was within: -0.003 (neck); -0.012 (thorax); -0.028 (pelvis, 10 MV) and 0.008 (pelvis, 18 MV) respectively.

A comparative description of three multipurpose phantoms (MPP) for external audits of photon beams in radiotherapy: the water MPP, the Umeå MPP and the EC MPP

AIM

To present a technical description and intercomparison of three multipurpose phantoms (MPP) developed for mailed dosimetry checks of therapeutic photon beams in reference and non-reference conditions.

MATERIALS

The W-MPP is a water MPP, whereas the Umeâ-MPP, made of perspex (PMMA, Plexiglas), and the EC-MPP, made of polystyrene, are solid MPPs. The W-MPP uses only thermoluminescent dosimeters (TLD) for dosimetry checks, the EC MPP uses film and TLD; the Umeâ phantom uses film and TLD, and offers in addition the possibility for ionization chamber measurements. Either using TLD or films, the MPPs have been designed to check on-axis and off-axis the following irradiation conditions: square and rectangular fields, asymmetric fields, wedged beams, oblique incidence and, for the solid MPPs, also the influence of inhomogeneities.

RESULTS AND DISCUSSION

The MPPs have been compared for different aspects going from their dosimetric performance (number of dosimetric parameters that can be checked) to some practical consideration in the use of the different MPPs (set-up time, stability, instruction sheets, etc.). From a comparison between the solid multi-purpose phantoms, it turns out that the EC-MPP is capable of checking the largest number of dosimetric parameters per beam, but has the longest set-up time ( approximately 2 h) per beam according to the users. The Umeå-MPP has a smaller number of set-ups (hence a smaller average time) and also includes some parameters not checked with the EC-MPP (e.g. SSD accuracy). The major drawback however of the Umeå-MPP is considered to be its high density (>1.1 g/cm(3)) which increases the difficulty of the analysis with the treatment planning system. The W-MPP checks the smallest number of parameters, but is the fastest in set-up time, the easiest for mailing, and is water equivalent, which is advantageous for the TPS checks. The major drawback of this MPP is however the inability to check complete dose distribution (film) or inhomogeneities.