Three-dimensional heart dose reconstruction to estimate normal tissue complication probability after breast irradiation using portal dosimetry

Irradiation of the heart is one of the major concerns during radiotherapy of breast cancer. Three-dimensional (3D) treatment planning would therefore be useful but cannot always be performed for left-sided breast treatments, because CT data may not be available. However, even if 3D dose calculations are available and an estimate of the normal tissue damage can be made, uncertainties in patient positioning may significantly influence the heart dose during treatment. Therefore, 3D reconstruction of the actual heart dose during breast cancer treatment using electronic imaging portal device (EPID) dosimetry has been investigated. A previously described method to reconstruct the dose in the patient from treatment portal images at the radiological midsurface was used in combination with a simple geometrical model of the irradiated heart volume to enable calculation of dose-volume histograms (DVHs), to independently verify this aspect of the treatment without using 3D data from a planning CT scan. To investigate the accuracy of our method, the DVHs obtained with full 3D treatment planning system (TPS) calculations and those obtained after resampling the TPS dose in the radiological midsurface were compared for fifteen breast cancer patients for whom CT data were available. In addition, EPID dosimetry as well as 3D dose calculations using our TPS, film dosimetry, and ionization chamber measurements were performed in an anthropomorphic phantom. It was found that the dose reconstructed using EPID dosimetry and the dose calculated with the TPS agreed within 1.5% in the lung/heart region. The dose-volume histograms obtained with EPID dosimetry were used to estimate the normal tissue complication probability (NTCP) for late excess cardiac mortality. Although the accuracy of these NTCP calculations might be limited due to the uncertainty in the NTCP model, in combination with our portal dosimetry approach it allows incorporation of the actual heart dose. For the anthropomorphic phantom, and for fifteen patients for whom CT data were available to test our method, the average difference between the NTCP values obtained with our method and those resulting from the dose distributions calculated with the TPS was 0.1% +/- 0.3% (1 SD). Most NTCP values were 1%-2% lower than those obtained using the method described by Hurkmans et al. [Radiother. Oncol. 62, 163-171 (2002)], using the maximum heart distance determined from a simulator image as a single pre-treatment parameter. A similar difference between the two methods was found for twelve patients using in vivo EPID dosimetry; the average NTCP value obtained with EPID dosimetry was 0.9%, whereas an average NTCP value of 2.2% was derived using the method of Hurkmans et al. The results obtained in this study show that EPID dosimetry is well suited for in vivo verification of the heart dose during breast cancer treatment, and can be used to estimate the NTCP for late excess cardiac mortality. To the best of our knowledge, this is the first study using portal dosimetry to calculate a DVH and NTCP of an organ at risk.