Patient-specific CT dosimetry calculation: a feasibility study.

Current estimation of radiation dose from computed tomography (CT) scans on patients has relied on the measurement of Computed Tomography Dose Index (CTDI) in standard cylindrical phantoms, and calculations based on mathematical representations of "standard man". Radiation dose to both adult and pediatric patients from a CT scan has been a concern, as noted in recent reports. The purpose of this study was to investigate the feasibility of adapting a radiation treatment planning system (RTPS) to provide patient-specific CT dosimetry. A radiation treatment planning system was modified to calculate patient-specific CT dose distributions, which can be represented by dose at specific points within an organ of interest, as well as organ dose-volumes (after image segmentation) for a GE Light Speed Ultra Plus CT scanner. The RTPS calculation algorithm is based on a semi-empirical, measured correction-based algorithm, which has been well established in the radiotherapy community. Digital representations of the physical phantoms (virtual phantom) were acquired with the GE CT scanner in axial mode. Thermoluminescent dosimeter (TLDs) measurements in pediatric anthropomorphic phantoms were utilized to validate the dose at specific points within organs of interest relative to RTPS calculations and Monte Carlo simulations of the same virtual phantoms (digital representation). Congruence of the calculated and measured point doses for the same physical anthropomorphic phantom geometry was used to verify the feasibility of the method. The RTPS algorithm can be extended to calculate the organ dose by calculating a dose distribution point-by-point for a designated volume. Electron Gamma Shower (EGSnrc) codes for radiation transport calculations developed by National Research Council of Canada (NRCC) were utilized to perform the Monte Carlo (MC) simulation. In general, the RTPS and MC dose calculations are within 10% of the TLD measurements for the infant and child chest scans. With respect to the dose comparisons for the head, the RTPS dose calculations are slightly higher (10%-20%) than the TLD measurements, while the MC results were within 10% of the TLD measurements. The advantage of the algebraic dose calculation engine of the RTPS is a substantially reduced computation time (minutes vs. days) relative to Monte Carlo calculations, as well as providing patient-specific dose estimation. It also provides the basis for a more elaborate reporting of dosimetric results, such as patient specific organ dose volumes after image segmentation.

Correction of motion-induced misalignment in co-registered PET/CT and MRI (T1/T2/FLAIR) head images for stereotactic radiosurgery.

The purpose was to evaluate and correct the co-registration of diagnostic PET/CT and MRI/MRI images for stereotactic radiosurgery (SRS) using 3D volumetric image registration (3DVIR). The 3DVIR utilizes the homogeneity of color distribution over a volumetric anatomical landmark as the registration criterion with submillimeter accuracy. Fifty-three PET/CT and MRI (T1, T2 and FLAIR) image sets of patients with brain lesions were acquired sequentially from a hybrid PET/CT or an MRI scanner with common diagnostic head holding devices. Twenty-five sets of head 18F-FDG-PET/CT images were scanned over a 10-minute interval and 14 whole-body sets were scanned over a 30-minute interval. Fourteen sets of MRI images were acquired, and each 3-modal image set (T1, T2 and FLAIR) was scanned in sequence at time 0, ~5 and ~20 minutes. The misalignments in these "co-registered" images were evaluated and corrected using the 3DVIR. Using the head immobilization devices commonly found in diagnostic PET/CT and MRI/MRI imaging, 80%-100% of these "co-registered" images were identified as misaligned. For PET/CT, the magnitude of misalignment was 0.4° ± 0.5° and 0.7 ± 0.4 mm for 10-minute scans, and 0.8° ± 1.2° and 2.7 ± 1.7 mm for 30-minute scans. For MRI/MRI, the magnitude was 0.2° ± 0.4° and 0.3 ± 0.2 mm for 5-minute scan intervals, and 1.1° ± 0.7° and 1.2 ± 1.4 mm for 20-minute intervals. Small, but significant, misalignment is present in the co-registered diagnostic PET/CT and MRI/MRI images and can be corrected in SRS treatment planning using the volumetric image registration for improved target localization within the clinical error tolerance.

Quantitative prediction of respiratory tidal volume based on the external torso volume change: a potential volumetric surrogate.

An external respiratory surrogate that not only highly correlates with but also quantitatively predicts internal tidal volume should be useful in guiding four-dimensional computed tomography (4DCT), as well as 4D radiation therapy (4DRT). A volumetric surrogate should have advantages over external fiducial point(s) for monitoring respiration-induced motion of the torso, which deforms in synchronization with a patient-specific breathing pattern. This study establishes a linear relationship between the external torso volume change (TVC) and lung air volume change (AVC) by validating a proposed volume conservation hypothesis (TVC = AVC) throughout the respiratory cycle using 4DCT and spirometry. Fourteen patients' torso 4DCT images and corresponding spirometric tidal volumes were acquired to examine this hypothesis. The 4DCT images were acquired using dual surrogates in ciné mode and amplitude-based binning in 12 respiratory stages, minimizing residual motion artifacts. Torso and lung volumes were calculated using threshold-based segmentation algorithms and volume changes were calculated relative to the full-exhalation stage. The TVC and AVC, as functions of respiratory stages, were compared, showing a high correlation (r = 0.992 +/- 0.005, p < 0.0001) as well as a linear relationship (slope = 1.027 +/- 0.061, R(2) = 0.980) without phase shift. The AVC was also compared to the spirometric tidal volumes, showing a similar linearity (slope = 1.030 +/- 0.092, R(2) = 0.947). In contrast, the thoracic and abdominal heights measured from 4DCT showed relatively low correlation (0.28 +/- 0.44 and 0.82 +/- 0.30, respectively) and location-dependent phase shifts. This novel approach establishes the foundation for developing an external volumetric respiratory surrogate.

A novel analytical approach to the prediction of respiratory diaphragm motion based on external torso volume change.

An analytical approach to predict respiratory diaphragm motion should have advantages over a correlation-based method, which cannot adapt to breathing pattern changes without re-calibration for a changing correlation and/or linear coefficient. To quantitatively calculate the diaphragm motion, a new expandable 'piston' respiratory (EPR) model was proposed and tested using 4DCT torso images of 14 patients. The EPR model allows two orthogonal lung motions (with a few volumetric constraints): (1) the lungs expand (DeltaV(EXP)) with the same anterior height variation as the thoracic surface, and (2) the lungs extend (DeltaV(EXT)) with the same inferior distance as the volumetrically equivalent 'piston' diaphragm. A volume conservation rule (VCR) established previously (Li et al 2009 Phys. Med. Biol. 54 1963-78) was applied to link the external torso volume change (TVC) to internal lung volume change (LVC) via lung air volume change (AVC). As the diaphragm moves inferiorly, the vacant space above the diaphragm inside the rib cage should be filled by lung tissue with a volume equal to DeltaV(EXT) (=LVC-DeltaV(EXP)), while the volume of non-lung tissues in the thoracic cavity should conserve. It was found that DeltaV(EXP) accounted for 3-24% of the LVC in these patients. The volumetric shape of the rib cage, characterized by the variation of cavity volume per slice over the piston motion range, deviated from a hollow cylinder by -1.1% to 6.0%, and correction was made iteratively if the variation is >3%. The predictions based on the LVC and TVC (with a conversion factor) were compared with measured diaphragm displacements (averaged from six pivot points), showing excellent agreements (0.2 +/- 0.7 mm and 0.2 +/- 1.2 mm, respectively), which are within clinically acceptable tolerance. Assuming motion synchronization between the piston and points of interest along the diaphragm, point motion was estimated but at higher uncertainty ( approximately 10% +/- 4%). This analytical approach provides a patient-independent technique to calculate the patient-specific diaphragm motion, using the anatomical and respiratory volumetric constraints.

Accuracy of 3D volumetric image registration based on CT, MR and PET/CT phantom experiments.

Registration is critical for image-based treatment planning and image-guided treatment delivery. Although automatic registration is available, manual, visual-based image fusion using three orthogonal planar views (3P) is always employed clinically to verify and adjust an automatic registration result. However, the 3P fusion can be time consuming, observer dependent, as well as prone to errors, owing to the incomplete 3-dimensional (3D) volumetric image representations. It is also limited to single-pixel precision (the screen resolution). The 3D volumetric image registration (3DVIR) technique was developed to overcome these shortcomings. This technique introduces a 4th dimension in the registration criteria beyond the image volume, offering both visual and quantitative correlation of corresponding anatomic landmarks within the two registration images, facilitating a volumetric image alignment, and minimizing potential registration errors. The 3DVIR combines image classification in real-time to select and visualize a reliable anatomic landmark, rather than using all voxels for alignment. To determine the detection limit of the visual and quantitative 3DVIR criteria, slightly misaligned images were simulated and presented to eight clinical personnel for interpretation. Both of the criteria produce a detection limit of 0.1 mm and 0.1 degree. To determine the accuracy of the 3DVIR method, three imaging modalities (CT, MR and PET/CT) were used to acquire multiple phantom images with known spatial shifts. Lateral shifts were applied to these phantoms with displacement intervals of 5.0+/-0.1 mm. The accuracy of the 3DVIR technique was determined by comparing the image shifts determined through registration to the physical shifts made experimentally. The registration accuracy, together with precision, was found to be: 0.02+/-0.09 mm for CT/CT images, 0.03+/-0.07 mm for MR/MR images, and 0.03+/-0.35 mm for PET/CT images. This accuracy is consistent with the detection limit, suggesting an absence of detectable systematic error. This 3DVIR technique provides a superior alternative to the 3P fusion method for clinical applications.

A novel 3D volumetric voxel registration technique for volume-view-guided image registration of multiple imaging modalities.

PURPOSE: To provide more clinically useful image registration with improved accuracy and reduced time, a novel technique of three-dimensional (3D) volumetric voxel registration of multimodality images is developed. METHODS AND MATERIALS: This technique can register up to four concurrent images from multi-modalities with volume view guidance. Various visualization effects can be applied, facilitating global and internal voxel registration. Fourteen computed tomography/magnetic resonance (CT/MR) image sets and two computed tomography/positron emission tomography (CT/PET) image sets are used. For comparison, an automatic registration technique using maximization of mutual information (MMI) and a three-orthogonal-planar (3P) registration technique are used. RESULTS: Visually sensitive registration criteria for CT/MR and CT/PET have been established, including the homogeneity of color distribution. Based on the registration results of 14 CT/MR images, the 3D voxel technique is in excellent agreement with the automatic MMI technique and is indicatory of a global positioning error (defined as the means and standard deviations of the error distribution) using the 3P pixel technique: 1.8 degrees +/- 1.2 degrees in rotation and 2.0 +/- 1.3 (voxel unit) in translation. To the best of our knowledge, this is the first time that such positioning error has been addressed. CONCLUSION: This novel 3D voxel technique establishes volume-view-guided image registration of up to four modalities. It improves registration accuracy with reduced time, compared with the 3P pixel technique. This article suggests that any interactive and automatic registration should be safe-guarded using the 3D voxel technique.

A novel four-dimensional radiotherapy planning strategy from a tumor-tracking beam's eye view.

To investigate the feasibility of four-dimensional radiotherapy (4DRT) planning from a tumor-tracking beam's eye view (ttBEV) with reliable gross tumor volume (GTV) delineation, realistic normal tissue representation, high planning accuracy and low clinical workload, we propose and validate a novel 4D conformal planning strategy based on a synthesized 3.5D computed tomographic (3.5DCT) image with a motion-compensated tumor. To recreate patient anatomy from a ttBEV in the moving tumor coordinate system for 4DRT planning (or 4D planning), the centers of delineated GTVs in all phase CT images of 4DCT were aligned, and then the aligned CTs were averaged to produce a new 3.5DCT image. This GTV-motion-compensated CT contains a motionless target (with motion artifacts minimized) and motion-blurred normal tissues (with a realistic temporal density average). Semi-automatic threshold-based segmentation of the tumor, lung and body was applied, while manual delineation was used for other organs at risk (OARs). To validate this 3.5DCT-based 4D planning strategy, five patients with peripheral lung lesions of small size (<5 cm(3)) and large motion range (1.2-3.5 cm) were retrospectively studied for stereotactic body radiotherapy (SBRT) using 3D conformal radiotherapy planning tools. The 3.5DCT-based 4D plan (3.5DCT plan) with 9-10 conformal beams was compared with the 4DCT-based 4D plan (4DCT plan). The 4DCT plan was derived from multiple 3D plans based on all phase CT images, each of which used the same conformal beam configuration but with an isocenter shift to aim at the moving tumor and a minor beam aperture and weighting adjustment to maintain plan conformality. The dose-volume histogram (DVH) of the 4DCT plan was created with two methods: one is an integrated DVH (iDVH(4D)), which is defined as the temporal average of all 3D-phase-plan DVHs, and the other (DVH(4D)) is based on the dose distribution in a reference phase CT image by dose warping from all phase plans using the displacement vector field (DVF) from a free-form deformable image registration (DIR). The DVH(3.5D) (for the 3.5DCT plan) was compared with both iDVH(4D) and DVH(4D). To quantify the DVH difference between the 3.5DCT plan and the 4DCT plan, two methods were used: relative difference (%) of the areas underneath the DVH curves and the volumes receiving more than 20% (V20) and 50% (V50) of prescribed dose of these 4D plans. The volume of the delineated GTV from different phase CTs varied dramatically from 24% to 112% among the five patients, whereas the GTV from 3.5DCT deviated from the averaged GTV in 4DCT by only -6%±6%. For planning tumor volume (PTV) coverage, the difference between the DVH(3.5D) and iDVH(4D) was negligible (<1% area), whereas the DVH(3.5D) and DVH(4D) were quite different, due to DIR uncertainty (∼2 mm), which propagates to PTV dose coverage with a pronounced uncertainty for small tumors (0.3-4.0 cm(3)) in stereotactic plans with sharp dose falloff around PTV. For OARs, such as the lung, heart, cord and esophagus, the three DVH curves (DVH(3.5D), DVH(4D) and iDVH(4D)) were found to be almost identical for the same patients, especially in high-dose regions. For the tumor-containing lung, the relative difference of the areas underneath the DVH curves was found to be small (5.3% area on average), of which 65% resulted from the low-dose region (D < 20%). The averaged V20 difference between the two 4D plans was 1.2% ± 0.8%. For the mean lung dose (MLD), the 3.5DCT plan differed from the 4DCT plan by -1.1%±1.3%. GTV-motion-compensated CT (3.5DCT) produces an accurate and reliable GTV delineation, which is close to the mean GTV from 4DCT. The 3.5DCT plan is equivalent to the 4DCT plan with <1% dose difference to the PTV and negligible dose difference in OARs. The 3.5DCT approach simplifies 4D planning and provides accurate dose calculation without a substantial increase of clinical workload for motion-tracking delivery to treat small peripheral lung tumors with large motion.

Intra- and inter-radiation therapist reproducibility of daily isocenter verification using prostatic fiducial markers.

BACKGROUND: We sought to determine the intra- and inter-radiation therapist reproducibility of a previously established matching technique for daily verification and correction of isocenter position relative to intraprostatic fiducial markers (FM). MATERIALS AND METHODS: With the patient in the treatment position, anterior-posterior and left lateral electronic images are acquired on an amorphous silicon flat panel electronic portal imaging device. After each portal image is acquired, the therapist manually translates and aligns the fiducial markers in the image to the marker contours on the digitally reconstructed radiograph. The distances between the planned and actual isocenter location is displayed. In order to determine the reproducibility of this technique, four therapists repeated and recorded this operation two separate times on 20 previously acquired portal image datasets from two patients. The data were analyzed to obtain the mean variability in the distances measured between and within observers. RESULTS: The mean and median intra-observer variability ranged from 0.4 to 0.7 mm and 0.3 to 0.6 mm respectively with a standard deviation of 0.4 to 1.0 mm. Inter-observer results were similar with a mean variability of 0.9 mm, a median of 0.6 mm, and a standard deviation of 0.7 mm. When using a 5 mm threshold, only 0.5% of treatments will undergo a table shift due to intra or inter-observer error, increasing to an error rate of 2.4% if this threshold were reduced to 3 mm. CONCLUSION: We have found high reproducibility with a previously established method for daily verification and correction of isocenter position relative to prostatic fiducial markers using electronic portal imaging.