Prediction of DVH parameter changes due to setup errors for breast cancer treatment based on 2D portal dosimetry.

Electronic portal imaging devices (EPIDs) are increasingly used for portal dosimetry applications. In our department, EPIDs are clinically used for two-dimensional (2D) transit dosimetry. Predicted and measured portal dose images are compared to detect dose delivery errors caused for instance by setup errors or organ motion. The aim of this work is to develop a model to predict dose-volume histogram (DVH) changes due to setup errors during breast cancer treatment using 2D transit dosimetry. First, correlations between DVH parameter changes and 2D gamma parameters are investigated for different simulated setup errors, which are described by a binomial logistic regression model. The model calculates the probability that a DVH parameter changes more than a specific tolerance level and uses several gamma evaluation parameters for the planning target volume (PTV) projection in the EPID plane as input. Second, the predictive model is applied to clinically measured portal images. Predicted DVH parameter changes are compared to calculated DVH parameter changes using the measured setup error resulting from a dosimetric registration procedure. Statistical accuracy is investigated by using receiver operating characteristic (ROC) curves and values for the area under the curve (AUC), sensitivity, specificity, positive and negative predictive values. Changes in the mean PTV dose larger than 5%, and changes in V90 and V95 larger than 10% are accurately predicted based on a set of 2D gamma parameters. Most pronounced changes in the three DVH parameters are found for setup errors in the lateral-medial direction. AUC, sensitivity, specificity, and negative predictive values were between 85% and 100% while the positive predictive values were lower but still higher than 54%. Clinical predictive value is decreased due to the occurrence of patient rotations or breast deformations during treatment, but the overall reliability of the predictive model remains high. Based on our predictive model, 2D transit dosimetry measurements can now directly be translated in clinically more relevant DVH parameter changes for the PTV during conventional breast treatment. In this way, the possibility to design decision protocols based on extracted DVH changes is created instead of undertaking elaborate actions such as repeated treatment planning or 3D dose reconstruction for a large group of patients.

A global calibration model for a-Si EPIDs used for transit dosimetry.

Electronic portal imaging devices (EPIDs) are not only applied for patient setup verification and detection of organ motion but are also increasingly used for dosimetric verification. The aim of our work is to obtain accurate dose distributions from a commercially available amorphous silicon (a-Si) EPID for transit dosimetry applications. For that purpose, a global calibration model was developed, which includes a correction procedure for ghosting effects, field size dependence and energy dependence of the a-Si EPID response. In addition, the long-term stability and additional buildup material for this type of EPID were determined. Differences in EPID response due to photon energy spectrum changes have been measured for different absorber thicknesses and field sizes, yielding off-axis spectrum correction factors based on transmission measurements. Dose measurements performed with an ionization chamber in a water tank were used as reference data, and the accuracy of the dosimetric calibration model was determined for a large range of treatment conditions. Gamma values using 3% as dose-difference criterion and 3 mm as distance-to-agreement criterion were used for evaluation. The field size dependence of the response could be corrected by a single kernel, fulfilling the gamma evaluation criteria in case of virtual wedges and intensity modulated radiation therapy fields. Differences in energy spectrum response amounted up to 30%-40%, but could be reduced to less than 3% using our correction model. For different treatment fields and (in)homogeneous phantoms, transit dose distributions satisfied in almost all situations the gamma criteria. We have shown that a-Si EPIDs can be accurately calibrated for transit dosimetry purposes.

Phased attenuation correction in respiration correlated computed tomography/positron emitted tomography.

The motion of lung tumors with respiration causes difficulties in the imaging with computed tomography (CT) and positronemitted tomography (PET). Since an accurate knowledge of the position of the tumor and the surrounding tissues is needed for radiation treatment planning, it is important to improve CT/PET image acquisition. The purpose of this study was to evaluate the potential to improve image acquisition using phased attenuation correction in respiration correlated CT/PET, where data of both modalities were binned retrospectively. Respiration correlated scans were made on a Siemens Biograph Sensation 16 CT/PET scanner which was modified to make a low pitch CT scan and list mode PET scan possible. A lollipop phantom was used in the experiments. The sphere with a diameter of 3.1 cm was filled with approximately 20 MBq 18F-FDG. Three longitudinal movement amplitudes were tested: 2.5, 3.9, and 4.8 cm. After collection of the raw CT data, list mode PET data, and the respiratory signal CT/PET images were binned to ten phases with the help of in-house-built software. Each PET phase was corrected for attenuation with CT data of the corresponding phase. For comparison, the attenuation correction was also performed with nonrespiration correlated (non-RC) CT data. The volume and the amplitude of the movement were calculated for every phaseof both the CT and PET data (with phased attenuation correction). Maximum and average activity concentrations were compared between the phased and nonphased attenuation corrected PET. With a standard non-RC CT/PET scan, the volume was underestimated by as much as 46% in CT and the PET volume was overestimated to 370%. The volumes found with RC-CT/PET scanning had average deviations of 1.9% (+/- 4.8%) and 1.5% (+/- 3.4%) from the actual volume, for the CT and PET volumes, respectively. Evaluation of the maximum activity concentration showed a clear displacement in the images with non-RC attenuation correction, and activity values were on average14% (+/- 12%) lower than with phased attenuation correction. The standard deviation of the maximum activity values found in the different phases was a factor of 10 smaller when phased attenuation correction was applied. In this phantom study, we have shown that a combination of respiration correlated CT/PET scanning with application of phased attenuation correction can improve the imaging of moving objects and can lead to improved volume estimation and a more precise localization and quantification of the activity.

Experimental verification of a portal dose prediction model.

Electronic portal imaging devices (EPIDs) can be used to measure a two-dimensional (2D) dose distribution behind a patient, thus allowing dosimetric treatment verification. For this purpose we experimentally assessed the accuracy of a 2D portal dose prediction model based on pencil beam scatter kernels. A straightforward derivation of these pencil beam scatter kernels for portal dose prediction models is presented based on phantom measurements. The model is able to predict the 2D portal dose image (PDI) behind a patient, based on a PDI without the patient in the beam in combination with the radiological thickness of the patient, which requires in addition a PDI with the patient in the beam. To assess the accuracy of portal dose and radiological thickness values obtained with our model, various types of homogeneous as well as inhomogeneous phantoms were irradiated with a 6 MV photon beam. With our model we are able to predict a PDI with an accuracy better than 2% (mean difference) if the radiological thickness of the object in the beam is symmetrically situated around the isocenter. For other situations deviations up to 3% are observed for a homogeneous phantom with a radiological thickness of 17 cm and a 9 cm shift of the midplane-to-detector distance. The model can extract the radiological thickness within 7 mm (maximum difference) of the actual radiological thickness if the object is symmetrically distributed around the isocenter plane. This difference in radiological thickness is related to a primary portal dose difference of 3%. It can be concluded that our model can be used as an easy and accurate tool for the 2D verification of patient treatments by comparing predicted and measured PDIs. The model is also able to extract the primary portal dose with a high accuracy, which can be used as the input for a 3D dose reconstruction method based on back-projection.

Verification of treatment parameter transfer by means of electronic portal dosimetry.

Electronic portal imaging devices (EPIDs) are mainly used for patient setup verification during treatment but other geometric properties like block shape and leaf positions are also determined. Electronic portal dosimetry allows dosimetric treatment verification. By combining geometric and dosimetric information, the data transfer between treatment planning system (TPS) and linear accelerator can be verified which in particular is important when this transfer is not carried out electronically. We have developed a pretreatment verification procedure of geometric and dosimetric treatment parameters of a 10 MV photon beam using an EPID. Measurements were performed with a CCD camera-based iView EPID, calibrated to convert a greyscale EPID image into a two-dimensional absolute dose distribution. Central field dose calculations, independent of the TPS, are made to predict dose values at a focus-EPID distance of 157.5 cm. In the same EPID image, the presence of a wedge, its direction, and the field size defined by the collimating jaws were determined. The accuracy of the procedure was determined for open and wedged fields for various field sizes. Ionization chamber measurements were performed to determine the accuracy of the dose values measured with the EPID and calculated by the central field dose calculation. The mean difference between ionization chamber and EPID dose at the center of the fields was 0.8 +/- 1.2% (1 s.d.). Deviations larger than 2.5% were found for half fields and fields with a jaw in overtravel. The mean difference between ionization chamber results and the independent dose calculation was -0.21 +/- 0.6% (1 s.d.). For all wedged fields, the presence of the wedge was detected and the mean difference in actual and measured wedge direction was 0 +/- 3 degrees (1 s.d.). The mean field size differences in X and Y directions were 0.1 +/- 0.1 cm and 0.0 +/- 0.1 cm (1 s.d.), respectively. Pretreatment monitor unit verification is possible with high accuracy and also geometric parameters can be verified using the same EPID image.

Three-dimensional dose reconstruction of breast cancer treatment using portal imaging.

In this study, we present an algorithm for three-dimensional (3-D) dose reconstruction using portal images obtained with an electronic portal imaging device (EPID). For this purpose an algorithm for 2-D dose reconstruction, which was previously developed in our institution, was adapted. The external contour of the patient was used to correct for absorption of primary photons, but the presence of inhomogeneities was not taken into account. The accuracy of the algorithm was determined by irradiating two anthropomorphic breast phantoms with 6 MV photons. The dose values derived from portal images were compared with results from 3-D dose calculations, which, in turn, were verified with data obtained with an ionization chamber and film dosimetry. It was found that the application of contour information significantly improves the accuracy of 2-D dose reconstruction. If the total dose at the isocenter plane resulting from all treatment beams is reconstructed, the average deviation from the planned dose is 0.1%+/-1.7% (1 SD). If contour information is not available, the differences increase up to +/-20% for the individual beams. In that case, the dose can only be reconstructed with reasonable accuracy when (nearly) opposing beams are used. The average deviation of the 3-D reconstructed dose from the planned dose in the irradiated volume is 1.4%+/-5.4% (1 SD). If the irradiated volume is enclosed by planes less than 5 cm distant from the isocenter plane, then the average deviation is only 0.5%+/-3.4% (1 SD). It can be concluded that the proposed algorithm for a 3-D dose reconstruction allows a determination of the dose at the isocenter plane and the dose-volume histogram with an accuracy acceptable for an independent verification of the treatment.