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.