Reconstruction of a time-averaged midposition CT scan for radiotherapy planning of lung cancer patients using deformable registration.

PURPOSE: lower lobe lung tumors move with amplitudes of up to 2 cm due to respiration. To reduce respiration imaging artifacts in planning CT scans, 4D imaging techniques are used. Currently, we use a single (midventilation) frame of the 4D data set for clinical delineation of structures and radiotherapy planning. A single frame, however, often contains artifacts due to breathing irregularities, and is noisier than a conventional CT scan since the exposure per frame is lower. Moreover, the tumor may be displaced from the mean tumor position due to hysteresis. The aim of this work is to develop a framework for the acquisition of a good quality scan representing all scanned anatomy in the mean position by averaging transformed (deformed) CT frames, i.e., canceling out motion. A nonrigid registration method is necessary since motion varies over the lung. METHODS AND MATERIALS: 4D and inspiration breath-hold (BH) CT scans were acquired for 13 patients. An iterative multiscale motion estimation technique was applied to the 4D CT scan, similar to optical flow but using image phase (gray-value transitions from bright to dark and vice versa) instead. From the (4D) deformation vector field (DVF) derived, the local mean position in the respiratory cycle was computed and the 4D DVF was modified to deform all structures of the original 4D CT scan to this mean position. A 3D midposition (MidP) CT scan was then obtained by (arithmetic or median) averaging of the deformed 4D CT scan. Image registration accuracy, tumor shape deviation with respect to the BH CT scan, and noise were determined to evaluate the image fidelity of the MidP CT scan and the performance of the technique. RESULTS: Accuracy of the used deformable image registration method was comparable to established automated locally rigid registration and to manual landmark registration (average difference to both methods < 0.5 mm for all directions) for the tumor region. From visual assessment, the registration was good for the clearly visible features (e.g., tumor and diaphragm). The shape of the tumor, with respect to that of the BH CT scan, was better represented by the MidP reconstructions than any of the 4D CT frames (including MidV; reduction of "shape differences" was 66%). The MidP scans contained about one-third the noise of individual 4D CT scan frames. CONCLUSIONS: We implemented an accurate method to estimate the motion of structures in a 4D CT scan. Subsequently, a novel method to create a midposition CT scan (time-weighted average of the anatomy) for treatment planning with reduced noise and artifacts was introduced. Tumor shape and position in the MidP CT scan represents that of the BH CT scan better than MidV CT scan and, therefore, was found to be appropriate for treatment planning.

Fusion of respiration-correlated PET and CT scans: correlated lung tumour motion in anatomical and functional scans.

Lower lobe lung tumours in particular can move up to 2 cm in the cranio-caudal direction during the respiration cycle. This breathing motion causes image artefacts in conventional free-breathing computed tomography (CT) and positron emission tomography (PET) scanning, rendering delineation of structures for radiotherapy inaccurate. The purpose of this study was to develop a method for four-dimensional (4D) respiration-correlated (RC) acquisition of both CT and PET scans and to develop a framework to fuse these modalities. The breathing signal was acquired using a thermometer in the breathing airflow of the patient. Using this breathing signal, the acquired CT and PET data were grouped to the corresponding respiratory phases, thereby obtaining 4D CT and PET scans. Tumour motion curves were assessed in both image modalities. From these tumour motion curves, the deviation with respect to the mean tumour position was calculated for each phase. The absolute position of the centre of the tumour, relative to the bony anatomy, in the RCCT and gated PET scans was determined. This 4D acquisition and 4D fusion methodology was performed for five patients with lower lobe tumours. The peak-to-peak amplitude range in this sample group was 1-2 cm. The 3D tumour motion curve differed less than 1 mm between PET and CT for all phases. The mean difference in amplitude was less than 1 mm. The position of the centre of the tumour (relative to the bony anatomy) in the RCCT and gated PET scan was similar (difference <1 mm) when no atelectasis was present. Based on these results, we conclude that the method described in this study allows for accurate quantification of tumour motion in CT and PET scans and yields accurate respiration-correlated 4D anatomical and functional information on the tumour region.

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.

Field size reduction enables iso-NTCP escalation of tumor control probability for irradiation of lung tumors.

PURPOSE: With the mean lung dose (MLD) as an estimator for the normal tissue complication probability (NTCP) of the lung, we assessed whether the probability of tumor control of lung tumors might be increased by dose escalation in combination with a reduction of field sizes, thus increasing target dose inhomogeneity while maintaining a constant MLD. METHODS AND MATERIALS: An 8-MV AP-PA irradiation of a lung tumor, located in a cylindrically symmetric lung-equivalent phantom, was modeled using numerical simulation. Movement of the clinical target volume (CTV) due to patient breathing and setup errors was simulated. The probability of tumor control, expressed as the equivalent uniform dose (EUD) of the CTV, was assessed as a function of field size, under the constraint of a constant MLD. The approach was tested for a treatment of a non-small cell lung cancer (NSCLC) patient using the beam directions of the clinically applied treatment plan. RESULTS: In the phantom simulation it was shown that by choosing field sizes that ensured a minimum dose of 95% in the CTV ("conventional" plan) taking into account setup errors and tumor motion, an EUD of the CTV of 43.8 Gy can be obtained for a prescribed dose of 44.2 Gy. By reducing the field size and thus shifting the 95% isodose surface inwards, the EUD increases to a maximum of 68.3 Gy with a minimum dose in the CTV of 55.2 Gy. This increase in EUD is caused by the fact that field size reduction enables escalation of the prescribed dose while maintaining a constant MLD. Further reduction of the field size results in decrease of the EUD because the minimum dose in the CTV becomes so low that it has a predominant effect on the EUD, despite further escalation of the prescribed dose. For the NSCLC patient, the EUD could be increased from an initial 62.2 Gy for the conventional plan, to 83.2 Gy at maximum. In this maximum, the prescribed dose is 88.1 Gy, and the minimum dose in the CTV is 67.4 Gy. In this case, the 95% isodose surface is conformed closely to the "static" CTV during treatment planning. CONCLUSIONS: Iso-NTCP escalation of the probability of tumor control is possible for lung tumors by reducing field sizes and allowing a larger dose inhomogeneity in the CTV. Optimum field sizes can be derived, having the highest EUD and highest minimum dose in the CTV under condition of a constant NTCP of the lungs. We conclude that the concept of homogeneous dose in the target volume is not the best approach to reach the highest probability of tumor control for lung tumors.

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.

The effect of breathing and set-up errors on the cumulative dose to a lung tumor.

BACKGROUND AND PURPOSE: To assess the impact of both set-up errors and respiration-induced tumor motion on the cumulative dose delivered to a clinical target volume (CTV) in lung, for an irradiation based on current clinically applied field sizes. MATERIALS AND METHODS: A cork phantom, having a 50 mm spherically shaped polystyrene insertion to simulate a gross tumor volume (GTV) located centrally in a lung was irradiated with two parallel opposed beams. The planned 95% isodose surface was conformed to the planning target volume (PTV) using a multi leaf collimator. The resulting margin between the CTV and the field edge was 16 mm in beam's eye view. A dose of 70 Gy was prescribed. Dose area histograms (DAHs) of the central plane of the CTV (GTV+5 mm) were determined using radiographic film for different combinations of set-up errors and respiration-induced tumor motion. The DAHs were evaluated using the population averaged tumor control probability (TCP(pop)) and the equivalent uniform dose (EUD) model. RESULTS: Compared with dose volume histograms of the entire CTV, DAHs overestimate the impact of tumor motion on tumor control. Due to the choice of field sizes a large part of the PTV will receive a too low dose resulting in an EUD of the central plane of the CTV of 68.9 Gy for the static case. The EUD drops to 68.2, 66.1 and 51.1 Gy for systematic set-up errors of 5, 10 and 15 mm, respectively. For random set-up errors of 5, 10 and 15 mm (1 SD), the EUD decreases to 68.7, 67.4 and 64.9 Gy, respectively. For similar amplitudes of respiration-induced motion, the EUD decreases to 68.8, 68.5 and 67.7 Gy, respectively. For a clinically relevant scenario of 7.5 mm systematic set-up error, 3 mm random set-up error and 5 mm amplitude of breathing motion, the EUD is 66.7 Gy. This corresponds with a tumor control probability TCP(pop) of 41.7%, compared with 50.0% for homogeneous irradiation of the CTV to 70 Gy. CONCLUSION: Systematic set-up errors have a dominant effect on the cumulative dose to the CTV. The effect of breathing motion and random set-up errors is smaller. Therefore the gain of controlling breathing motion during irradiation is expected to be small and efforts should rather focus on minimizing systematic errors. For the current clinically applied field sizes and a clinically relevant combination of set-up errors and breathing motion, the EUD of the central plane of the CTV is reduced by 3.3 Gy, at maximum, relative to homogeneous irradiation of the CTV to 70 Gy, for our worst case scenario.

Planning, computer optimization, and dosimetric verification of a segmented irradiation technique for prostate cancer.

PURPOSE: To develop and verify a multisegment technique for prostate irradiation that results in better sparing of the rectal wall compared to a conventional three-field technique, for patients with a concave-shaped planning target volume (PTV) overlapping the rectal wall. METHODS AND MATERIALS: Five patients have been selected with various degrees of overlap between PTV and rectal wall. The planned dose to the ICRU reference point is 78 Gy. The new technique consists of five beams, each having an open segment covering the entire PTV and several smaller segments in which the rectum is shielded. Segment weights are computer-optimized using an algorithm based on simulated annealing. The score function to be minimized consists of dose-volume constraints for PTV, rectal wall, and femoral heads. The resulting dose distribution is verified for each patient by using point measurements and line scans made with an ionization chamber in a water tank and by using film in a cylindrical polystyrene phantom. RESULTS: The final number of segments in the five-field technique ranges from 7 to 9 after optimization. Compared to the standard three-field technique, the maximum dose to the rectal wall decreases by approximately 3 Gy for patients with a large overlap and 1 Gy for patients with no overlap, resulting in a reduction of the normal tissue complication probability (NTCP) by a factor of 1.3 and 1.2, respectively. The mean dose to the PTV is the same for the two techniques, but the dose distribution is slightly less homogeneous with the five-field technique (Average standard deviation of five patients is 1.1 Gy and 1.7 Gy for the three-field and five-field technique, respectively). Ionization chamber measurements show that in the PTV, the calculated dose is in general within 1% of the measured dose. Outside the PTV, systematic dose deviations of up to 3% exist. Film measurements show that for the complete treatment, the position of the isodose lines in sagittal and coronal planes is calculated fairly accurately, the maximum distance between measured and calculated isodoses being 4 mm. CONCLUSIONS: We developed a relatively simple multisegment "step-and-shoot" technique that can be delivered within an acceptable time frame at the treatment machine (Extra time needed is approximately 3 minutes). The technique results in better sparing of the rectal wall compared to the conventional three-field technique. The technique can be planned and optimized relatively easily using automated procedures and a predefined score function. Dose calculation is accurate and can be verified for each patient individually.

Changes in local pulmonary injury up to 48 months after irradiation for lymphoma and breast cancer.

PURPOSE: To assess the recovery from early local pulmonary injury after irradiation and to determine whether regional differences exist. METHODS: For 110 patients treated for breast cancer or malignant lymphoma, single photon emission computed tomography (SPECT) perfusion and ventilation scans and CT scans were made before, 3, 18, and 48 months after radiotherapy. Dose-effect relations for changes in local perfusion, ventilation, and density were determined for each individual patient using spatially correlated SPECT and CT data sets, for each follow-up period. Average dose-effect relations for both subgroups were determined, as well as dose-effect relations for different regions. RESULTS: In general, partial improvement of local pulmonary injury was observed between 3 and 18 months for each of the three endpoints. After 18 months, no further improvement was seen. Patients with breast cancer and malignant lymphoma showed a similar improvement (except for the perfusion parameter), which was attributed to a recovery from the early radiation response and could not be explained by contraction effects of fibrosis of lung parenchyma. No regional differences in radiosensitivity 18 months after treatment were observed, except for the dorsal versus ventral region. This difference was attributed to a gravity-related effect in the measuring procedure. CONCLUSION: For all patients, a partial recovery from early local perfusion, ventilation, and density changes, was seen between 3 and 18 months after radiotherapy. After 18 months, local lung function did not further improve (lymphoma patients).

Dose-effect relations for early local pulmonary injury after irradiation for malignant lymphoma and breast cancer.

PURPOSE: To quantify the influence of treatment- and patient-related factors on the severity of early local pulmonary injury and to establish whether regional differences are present for local dose-effect relations for early radiation-induced pulmonary injury. METHODS: Forty-two patients with malignant lymphoma and 40 breast cancer patients were examined prior to and 3 months after radiotherapy. The lymphoma patients were irradiated with mantle fields to an average dose of 38 Gy and the breast cancer patients were irradiated with internal mammary node fields with or without tangential breast fields to an average dose of 50 Gy. Dose-effect relations for local perfusion, ventilation and density changes were determined using correlated single photon emission computed tomography (SPECT) and CT data. A multivariate analysis was performed to study the influence of irradiated volume, chemotherapy (CMF and MOPP/ABV), smoking, age and gender. In addition, dose-effect relations for different regions in the lung were determined. RESULTS: A similar and almost linear increase of early functional changes as a function of radiation dose was observed for perfusion and ventilation, whereas the shape of the dose-effect relation and the magnitude of early structural changes were different for density. For the three end-points studied, regional differences in radiosensitivity could not be demonstrated. For the posterior lung region compared to the anterior lung region, however, a difference was observed, which could be attributed to a gravity-related effect in the measuring procedure. Local structural changes (density) were significantly smaller for smokers (P = 0.002) and young patients (P = 0.007), whereas the CMF chemotherapy regimen given after radiotherapy (P = 0.017) significantly increased the amount of functional changes (perfusion). The magnitude of local pulmonary changes was independent of the irradiated volume, the MOPP/ABV chemotherapy regimen and gender. CONCLUSION: The dose-effect relations for early radiation-induced local pulmonary changes were independent of the irradiated volume, MOPP/ABV, gender and lung region. CMF, smoking and age influenced the magnitude of early pulmonary changes and should be taken into account in dose-escalation protocols.

Evaluation of two dose-volume histogram reduction models for the prediction of radiation pneumonitis.

PURPOSE: To evaluate the similarities between the mean lung dose and two dose-volume histogram (DVH) reduction techniques of 3D dose distributions of the lung. PATIENTS AND METHODS: DVHs of the lungs were calculated from 3D dose distributions of patients treated for malignant lymphoma (44), breast cancer (42) and lung cancer (20). With a DVH reduction technique, a DVH is summarized by the equivalent uniform dose (EUD), a quantity which is directly related to the normal tissue complication probability (NTCP). Two DVH reduction techniques were used. The first was based on an empirical model proposed by Kutcher et al. (Kutcher, G.J., Burman, C., Brewster, M.S., Goitein, M. and Mohan, R. Histogram reduction method for calculating complication probabilities for three-dimensional treatment planning evaluations. Int. J. Radiat. Oncol. Biol. Phys. 21: 137-146, 1991), which needs a volume exponent n. Several values for n were tested. The second technique was based on a radiobiological model, the parallel functional subunit model developed by Niemierko et al. (Niemierko, A. and Goitein, M. Modeling of normal tissue response to radiation: the critical volume model. Int. J. Radiat. Oncol. Biol. Phys. 25: 135-145, 1993) and Jackson et al. (Jackson, A., Kutcher, G.J. and Yorke, E.D. Probability of radiation-induced complications for normal tissues with parallel architecture subject to non-uniform irradiation. Med. Phys. 20: 613-625, 1993), for which a local dose-effect relation needed to be specified. This relation was obtained from an analysis of perfusion and ventilation SPECT data. RESULTS: It can be shown analytically that the two DVH reduction techniques are identical, if the local dose-effect relation obeys a power-law relationship in the clinical dose range. Local dose-effect relations based on perfusion and ventilation SPECT data can indeed be fitted with a power-law relationship in the range 0-80 Gy, from which values of n = 0.8-0.9 were deduced. These correspond to the commonly used value of n = 0.87 for lung tissue and yielded EUDn=0.87 values which were almost identical to the mean lung doses. For other n values, for which no experimental data are present, differences exist between EUD and mean dose values. Six patients with malignant lymphoma (6/44) and none of the breast cancer patients (0/42) developed radiation pneumonitis. These cases occurred only at high values for the mean lung dose. CONCLUSION: The two DVH reduction techniques are identical for lung and are very similar to mean dose calculations. The two techniques are also relatively similar for other model parameter values.

Automatic three-dimensional matching of CT-SPECT and CT-CT to localize lung damage after radiotherapy.

UNLABELLED: The aim of this study was to develop a fast and clinically robust automatic method to register SPECT and CT scans of the lungs. METHODS: CT and SPECT scans were acquired in the supine position from 20 patients with healthy lungs. After partial irradiation of the lungs by radiotherapy, the scans were repeated. Two matching methods were compared: a conventional method with external skin markers and a new method using chamfer matching of the lung contours. In the latter method, a unique value for the SPECT threshold, needed for segmentation of the SPECT lungs, was determined by iteratively applying the chamfer matching algorithm. RESULTS: The new technique for CT-SPECT matching could be implemented in a fully automatic manner and required less than 2 min. No large systematic shifts or rotations were present between the matches obtained with the marker method and the lung contour method for healthy or partially irradiated lungs. For healthy lungs, the number of ventilation SPECT counts outside the CT-defined lung was taken as a measure for a good match. This number of outside counts was slightly lower for the new method than for the conventional method, which indicates that the accuracy of the new method is at least comparable to the conventional method. For ventilation, a systematic difference between the results of the matching methods, a small translation in the anterior --> posterior direction, could be attributed to an inconsistency of the marker positions (2 mm). For perfusion, a somewhat larger anterior --> posterior shift was found, which was attributed to the gravity force. CT-CT correlation on the lung contours using chamfer matching was tested with the same dataset. For accurate matching, the CT slices encompassing the diaphragm had to be deleted. CONCLUSION: The new method based on lung contour matching is a fast, automatic procedure and allows accurate clinical follow-up.

Prediction of overall pulmonary function loss in relation to the 3-D dose distribution for patients with breast cancer and malignant lymphoma.

PURPOSE: To predict the changes in pulmonary function tests (PFTs) 3-4 months after radiotherapy based on the three-dimensional (3-D) dose distribution and taking into account patient- and treatment-related factors. METHODS: For 81 patients with malignant lymphoma and breast cancer, PFTs (VA, VC, FEV1 and TL,COc) were performed prior to and 3-4 months after irradiation and dose-effect relations for early changes in local perfusion, ventilation and air-filled fraction were determined using correlated CT and SPECT data. The 3-D dose distribution of each patient was converted into four different dose-volume parameters, i.e. the mean dose in the lung and three overall response parameters (ORPs, which represent the average local injury over the complete lung). ORPs were determined using the dose-effect relations for early changes in local perfusion, ventilation and air-filled fraction. Correlation coefficients were calculated between these dose-volume parameters and the changes in PFTs. In addition, the impact of the variables chemotherapy (MOPP/ABV and CMF), tamoxifen, smoking, age and gender on the relation between the mean lung dose and the relative changes in PFTs following radiotherapy was studied using multiple regression analysis. RESULTS: The mean lung dose proved to be the easiest parameter to predict the reduction in PFTs 3-4 months following radiotherapy. For all patients the relation between the mean lung dose and the changes in PFTs could be described with one regression line through the origin and a slope of 1% reduction in PFT for each increase of 1 Gy in mean lung dose. Smoking and CMF chemotherapy influenced the reduction in PFTs significantly for VA and TL,COc, respectively. Patients treated with MOPP/ABV prior to radiotherapy had lower pre-radiotherapy PFTs than other patient groups, but did not show further deterioration after radiotherapy (at 3-4 months). CONCLUSIONS: The relative reduction in VA, VC, FEV1 and TL,COc 3-4 months after radiotherapy for breast cancer and malignant lymphoma can be estimated before radiotherapy based on the mean lung dose of each individual patient and taking into account the use of chemotherapy and smoking habits of the patient.

An analysis of the effect of ovoid shields in a selectron-LDR cervical applicator on dose distributions in rectum and bladder.

PURPOSE: A disadvantage of ovoid shields in a Fletcher-type applicator is that these shields cause artifacts on postimplant CT images. CT images, however, make it possible to calculate the dose distribution in the rectum and the bladder. To be able to estimate the possible advantage of having CT information over the use of ovoid shields without having CT information, we investigated the influence of shielding segments in a Fletcher-type Selectron-LDR applicator on the dose distribution in rectum and bladder. METHODS AND MATERIALS: Contours of rectum and bladder were delineated on transaxial CT slices of 15 unshielded applications. Of the volumes contained within these structures dose-volume histograms (DVHs) were calculated. In a similar way, DVHs of simulated shielded applications were calculated. The reduction, due to shielding, of the dose to the 2 cm3 (D2) and 5 cm3 (D5) volume of the cumulative DVHs of rectum and bladder, were determined. An isodose pattern in the sagittal plane through the center of each applicator was plotted to compare the location of the shielded area with the location of maximum dose in rectum and bladder in the unshielded situation. In two cases local dose reductions to the rectal wall were determined by calculating the dose in points at 10-mm intervals on the rectal contours. RESULTS: For the rectum, the reduction of D2 ranged from 0 to 11.1%, with an average of 5.0%; the reduction of D5 ranged from 2.3 to 12.1%, with an average of 6.4%. The reduction of D2 and D5 for the bladder ranged from 0 to 11.9% and from 0 to 11.6%, with average values of 2.2 and 2.6%, respectively. In 8 out of 15 cases the rectal maximum dose was located inferior to the shielded area. In all cases except one the bladder maximum dose was located superior to the shielded area. Local dose reductions on the rectal wall can be as high as 30% or more in an optimally shielded area. CONCLUSIONS: Reductions of D2 and D5 to rectum and bladder due to shielding are rather small, because the shielded area does usually not coincide with the high dose region and even if it does, the shielded area is too small to result in large reductions of these values. Because local dose reductions vary largely, one should proceed with caution when calculating the dose in just one rectal or bladder reference point. Because large overall dose reductions cannot be achieved with shielding, it is safe to use an unshielded applicator when post implant CT images are used to realize optimized dose distributions.

Recovery of overall and local lung function loss 18 months after irradiation for malignant lymphoma.

PURPOSE: To determine the local and overall pulmonary injury 3 to 18 months after irradiation and to investigate whether the changes in overall lung function can be predicted using the three-dimensional (3-D) dose distribution in combination with dose-effect relations for local injury; and to study the influence of chemotherapy on the injury. PATIENTS AND METHODS: Local perfusion (Q), ventilation (V), and tissue density were measured in 25 patients treated for malignant lymphoma, before, 3 to 4 months after, and 18 months after irradiation. Dose-effect relations for local injury, calculated using correlated single-photon emission computed tomographic (SPECT) and computed tomographic (CT) data, were combined with the 3-D dose distribution, to calculate the estimated mean local changes over the complete lung for each patient. The result was correlated with the actual changes in pulmonary function. RESULTS: A dose-dependent increase with injury was observed at 3 to 4 months after irradiation, which at 18 months had recovered by approximately 50% to 60%. The estimated mean relative reduction of local Q predicted the change in overall lung function within 10% of the actually observed values in 63% to 73% of patients. Chemotherapy given before radiotherapy enhanced radiation-induced reduction of local Q significantly, with dose-modifying factors of 1.22 and 1.37 at 3 to 4 months and 18 months, respectively. CONCLUSION: Partial recovery of radiation-induced reduction of local and overall lung function was observed at 18 months after irradiation. The overall functional outcome of most patients could be well predicted, based on the estimated mean local injury over the complete lung. Chemotherapy given before radiotherapy enhanced the radiation-induced reduction of local Q.

Estimation of overall pulmonary function after irradiation using dose-effect relations for local functional injury.

PURPOSE: To predict the pulmonary function 3-4 months after irradiation for malignant lymphoma from the three-dimensional (3-D) dose distribution. METHODS: Dose-effect relations for the relative reduction of local perfusion (Q) and local ventilation (V), were calculated in 25 patients, using correlated SPECT (Single Photon Emission Computed Tomography) and CT data. By combining the 3-D dose distribution of an individual patient with the dose-effect relations averaged over all patients, the average reduction of local Q and V (i.e., the overall response parameters) in the whole lung was estimated for each patient. Correlation coefficients were calculated between these overall response parameters and the change in standard lung function tests. In addition, the relation between the overall response parameters and the incidence of radiation pneumonitis was determined. RESULTS: The overall response parameter for perfusion was correlated with the change in standard lung function tests, with correlation coefficients varying between 0.53 (p = 0.007) and 0.71 (p < 0.001) for the change of Vital Capacity and Forced Expiratory Volume at 1 s, respectively. For the overall response parameter for ventilation similar correlations were observed. Four out of the 25 patients developed radiation pneumonitis; in these four patients the overall response parameter for perfusion was on average somewhat higher (13.2 +/- 1.4% (1 standard error of the mean)) than in patients without radiation pneumonitis (10.5 +/- 1.0%), but this difference was not significant. A higher incidence of radiation pneumonitis was observed for larger values of the overall response parameter for perfusion; in patient groups with an overall response parameter for perfusion of 0-5%, 5-10%, 10-15%, and 15-20%, the incidence of radiation pneumonitis was 0 (0/1), 10 (1/10), 13 (1/8) and 33% (2/6), respectively. CONCLUSION: By combining the 3-D dose distribution with the average dose-effect relations for local perfusion or ventilation, an overall response parameter can be calculated prior to irradiation, which is predictive for the radiation-induced change in the overall pulmonary function, and possibly for the incidence of radiation pneumonitis, in this group of patients.

Dose-effect relations for local functional and structural changes of the lung after irradiation for malignant lymphoma.

PURPOSE: To estimate the dose-effect relations for local functional (ventilation and perfusion) and structural (density) changes of the lung, 3-4 months after irradiation. METHODS: Twenty-five patients with malignant lymphoma were irradiated with a (modified) mantle field to an average dose of 38 Gy, given in 21 fractions. Single photon emission computed tomography (SPECT) ventilation (V) and perfusion (Q) scans, and CT scans were performed before and 3-4 months after radiation treatment. The three-dimensional dose distribution was calculated using the CT data. After correlation of SPECT and CT data sets, the average post-treatment value of V, Q and lung density per voxel was calculated relative to the pre-treatment value, per dose interval of 4 Gy. Subsequently, the dose-effect relations in each patient were normalized to the average value per voxel in the dose interval of 0-12 Gy. In addition, in each dose interval of 4 Gy the fraction of patients with changes larger than 20% was calculated for all three parameters. The dose-effect relations for perfusion and ventilation normalized to the low-dose regions, and the dose-incidence curves for the fraction of patients with changes larger than 20% were fitted for all three parameters, using a logistic model. RESULTS: Marked changes in the distribution of V and Q were found after irradiation. Prior to normalization to the low-dose regions, a change in V and Q was found in most patients in the dose interval of 0-12 Gy, varying from an increase of 37% to a decrease of 10%, which was followed by a decreasing trend at higher doses. The increase in the low-dose regions indicated a redistribution phenomenon, the magnitude of which was dependent of the irradiated volume. The logistic fit of the dose-effect relations for Q and V, normalized to the low-dose regions, resulted in values for D50 of 51 Gy and 54 Gy (given in 21 fractions on average), respectively, and for the steepness parameter k of 4.2 and 4.0, respectively. The logistic fit for the dose-incidence curves for Q, V and lung density resulted in values for D50 and k of 38 Gy, 37 Gy, 44 Gy and 10.3, 7.8 and 9.4, respectively. CONCLUSIONS: With the combined use of SPECT and CT scans, we have obtained dose-effect relations for local functional and structural damage in the lung, 3-4 months after irradiation.

Quantifying local lung perfusion and ventilation using correlated SPECT and CT data.

UNLABELLED: A clinically applicable method for quantifying lung perfusion and ventilation on a subregional (local) scale from SPECT scans in order to estimate local lung function in patients with pre-existing pulmonary disease and to monitor local treatment effects was developed and evaluated. METHODS: SPECT 99mTc perfusion and 81mKr ventilation images were corrected for photon attenuation and scatter effect with a postreconstruction correction method incorporating a variable-effective linear-attenuation coefficient calculated from spatially-correlated CT data. A new algorithm was developed to quantify local ventilation from the SPECT data, which, in contrast with other algorithms, makes no assumptions on ventilation homogeneity over the lung. The quantification procedure was applied to clinical data from patients with a normal lung function and from patients suffering from radiation-induced pulmonary dysfunction. RESULTS: The calculated attenuation correction factors on the observed number of counts in the lung range from 2.0 to 3.0 and 2.3 to 3.5 for 81mKr and 99mTc, respectively, showing a systematic increase from the diaphragm to the lung apex. As a result of this correction, the values of local perfusion and ventilation differ 10%-15% from values calculated without attenuation correction. The calculated values of the local ventilation are 10%-50% lower than those found by quantification algorithms which assume homogeneous ventilation. CONCLUSIONS: The methods presented here are robust with respect to uncertainties in the input parameters and yield realistic values for perfusion and ventilation distribution in the lung with an intrinsic accuracy (largely determined by count statistics) of about 10%.

A new method to determine dose-effect relations for local lung-function changes using correlated SPECT and CT data.

PURPOSE: To determine dose-effect relations for regional lung-function changes after radiotherapy. METHODS: Single Photon Emission Computed Tomography (SPECT) was performed to quantify regional ventilation and perfusion. CT scans were used to calculate the three-dimensional (3-D) dose distribution. Both SPECT and CT scans were performed prior to radiotherapy and 5 months after the start of the treatment. To obtain combined 3-D information on ventilation, perfusion and dose, the SPECT data were correlated with the corresponding CT data. The relative changes in ventilation and perfusion were calculated in each SPECT voxel (voxel size about 6 x 6 x 6 mm) and related to the dose in that voxel. The average relative changes were determined per dose interval of 4 Gy. This procedure was evaluated using the data from five patients treated for Hodgkin's disease with mantle field irradiation with a prescribed total dose of 40-42 Gy. RESULTS: Dose-effect relations for perfusion were observed in all patients, while in four of the five patients, a dose-effect relation was found for ventilation. The maximal uncertainty of the calculated radiation dose was 11%: a difference between the position of the patient during treatment and during CT scanning caused a maximal dose uncertainty of 6%, while the accuracy of the dose calculation algorithm itself was estimated to be within 5%. CONCLUSION: The results indicate that the combined use of SPECT and CT information is an effective method for determining dose-effect relations for regional lung function parameters in each individual patient.