Phantom validation of coregistration of PET and CT for image-guided radiotherapy

Radiotherapy treatment planning integrating positron emission tomography (PET) and computerized tomography (CT) is rapidly gaining acceptance in the clinical setting. Although hybrid systems are available, often the planning CT is acquired on a dedicated system separate from the PET scanner. A limiting factor to using PET data becomes the accuracy of the CT/PET registration. In this work, we use phantom and patient validation to demonstrate a general method for assessing the accuracy of CT/PET image registration and apply it to two multi-modality image registration programs. An IAEA (International Atomic Energy Association) brain phantom and an anthropomorphic head phantom were used. Internal volumes and externally mounted fiducial markers were filled with CT contrast and 18F-fluorodeoxyglucose (FDG). CT, PET emission, and PET transmission images were acquired and registered using two different image registration algorithms. CT/PET Fusion (GE Medical Systems, Milwaukee, WI) is commercially available and uses a semi-automated initial step followed by manual adjustment. Automatic Mutual Information-based Registration (AMIR), developed at our institution, is fully automated and exhibits no variation between repeated registrations. Registration was performed using distinct phantom structures; assessment of accuracy was determined from registration of the calculated centroids of a set of fiducial markers. By comparing structure-based registration with fiducial-based registration, target registration error (TRE) was computed at each point in a three-dimensional (3D) grid that spans the image volume. Identical methods were also applied to patient data to assess CT/PET registration accuracy. Accuracy was calculated as the mean with standard deviation of the TRE for every point in the 3D grid. Overall TRE values for the IAEA brain phantom are: CT/PET Fusion = 1.71 +/- 0.62 mm, AMIR = 1.13 +/- 0.53 mm; overall TRE values for the anthropomorphic head phantom are: CT/PET Fusion = 1.66 +/- 0.53 mm, AMIR = 1.15 +/- 0.48 mm. Precision (repeatability by a single user) measured for CT/PET Fusion: IAEA phantom = 1.59 +/- 0.67 mm and anthropomorphic head phantom = 1.63 +/- 0.52 mm. (AMIR has exact precision and so no measurements are necessary.) One sample patient demonstrated the following accuracy results: CT/PET Fusion = 3.89 +/- 1.61 mm, AMIR = 2.86 +/- 0.60 mm. Semi-automatic and automatic image registration methods may be used to facilitate incorporation of PET data into radiotherapy treatment planning in relatively rigid anatomic sites, such as head and neck. The overall accuracies in phantom and patient images are < 2 mm and < 4 mm, respectively, using either registration algorithm. Registration accuracy may decrease, however, as distance from the initial registration points (CT/PET fusion) or center of the image (AMIR) increases. Additional information provided by PET may improve dose coverage to active tumor subregions and hence tumor control. This study shows that the accuracy obtained by image registration with these two methods is well suited for image-guided radiotherapy.     

Evaluation of three types of reference image data for external beam radiotherapy target localization using digital tomosynthesis (DTS)

Digital tomosynthesis (DTS) is a fast, low-dose three-dimensional (3D) imaging approach which yields slice images with excellent in-plane resolution, though low plane-to-plane resolution. A stack of DTS slices can be reconstructed from a single limited-angle scan, with typical scan angles ranging from 10 degrees to 40 degrees and acquisition times of less than 10 s. The resulting DTS slices show soft tissue contrast approaching that of full cone-beam CT. External beam radiotherapy target localization using DTS requires the registration of on-board DTS images with corresponding reference image data. This study evaluates three types of reference volume: original reference CT, exact reference DTS (RDTS), and a more computationally efficient approximate reference DTS (RDTSapprox), as well as three different DTS scan angles (22 degrees, 44 degrees, and 65 degrees) for the DTS target localization task. Three-dimensional mutual information (MI) shared between reference and onboard DTS volumes was computed in a region surrounding the spine of a chest phantom, as translations spanning +/-5 mm and rotations spanning +/-5 degrees were simulated along each dimension in the reference volumes. The locations of the MI maxima were used as surrogates for registration accuracy, and the width of the MI peaks were used to characterize the registration robustness. The results show that conventional treatment planning CT volumes are inadequate reference volumes for direct registration with on-board DTS data. The efficient RDTSapprox method also appears insufficient for MI-based registration without further refinement of the technique, though it may be suitable for manual registration performed by a human observer. The exact RDTS volumes, on the other hand, delivered a 3D DTS localization accuracy of 0.5 mm and 0.50 along each axis, using only a single 44 degrees coronal on-board DTS scan of the chest phantom.     

Registration of electronic portal images for patient set-up verification

Images acquired from an electronic portal imaging device are aligned with digitally reconstructed radiographs (DRRs) or other portal images to verify patient positioning during radiation therapy. Most of the currently available computer aided registration methods are based on the manual placement of corresponding landmarks. The purpose of the paper is twofold: (a) the establishment of a methodology for patient set-up verification during radiotherapy based on the registration of electronic portal images, and (b) the evaluation of the proposed methodology in a clinical environment. The estimation of set-up errors, using the proposed methodology, can be accomplished by matching the portal image of the current fraction of the treatment with the portal image of the baseline treatment (reference portal image) using a nearly automated technique. The proposed registration method is tested on a number of phantom data as well as on data from four patients. The phantom data included portal images that corresponded to various positions of the phantom on the treatment couch. For each patient, a set of 30 portal images was used. For the phantom data (for both transverse and lateral portal images), the maximum absolute deviations of the translational shifts were within 1.5 mm, whereas the in-plane rotation angle error was less than 0.5 degrees. The two-way Anova revealed no statistical significant variability both within observer and between-observer measurements (P > 0.05). For the patient data, the mean values obtained with manual and the proposed registration methods were within 0.5 mm. In conclusion, the proposed registration method has been incorporated within a system, called ESTERR-PRO. Its image registration capability achieves high accuracy and both intra- and inter-user reproducibility. The system is fully operational within the Radiotherapy Department of 'HYGEIA' Hospital in Athens and it could be easily installed in any other clinical environment since it requires standardized hardware specifications and minimal human intervention.     

Automatic registration between reference and on-board digital tomosynthesis images for positioning verification

The authors developed a hybrid multiresolution rigid-body registration technique to automatically register reference digital tomosynthesis (DTS) images with on-board DTS images to guide patient positioning in radiation therapy. This hybrid registration technique uses a faster but less accurate static method to achieve an initial registration, followed by a slower but more accurate adaptive method to fine tune the registration. A multiresolution scheme is employed in the registration to further improve the registration accuracy, robustness, and efficiency. Normalized mutual information is selected as the criterion for the similarity measure and the downhill simplex method is used as the search engine. This technique was tested using image data both from an anthropomorphic chest phantom and from eight head-and-neck cancer patients. The effects of the scan angle and the region-of-interest (ROI) size on the registration accuracy and robustness were investigated. The necessity of using the adaptive registration method in the hybrid technique was validated by comparing the results of the static method and the hybrid method. With a 44 degrees scan angle and a large ROI covering the entire DTS volume, the average of the registration capture ranges in single-axis simulations was between -31 and +34 deg for rotations and between -89 and +78 mm for translations in the phantom study, and between -38 and +38 deg for rotations and between -58 and +65 mm for translations in the patient study. Decreasing the DTS scan angle from 44 degrees to 22 degrees mainly degraded the registration accuracy and robustness for the out-of-plane rotations. Decreasing the ROI size from the entire DTS volume to the volume surrounding the spinal cord reduced the capture ranges to between -23 and +18 deg for rotations and between -33 and +43 mm for translations in the phantom study, and between -18 and +25 deg for rotations and between -35 and +39 mm for translations in the patient study. Results also showed that the hybrid registration technique had much larger capture ranges than the static method alone in registering the out-of-plane rotations.     

Validation of an accelerated 'demons' algorithm for deformable image registration in radiation therapy

A greyscale-based fully automatic deformable image registration algorithm, originally known as the 'demons' algorithm, was implemented for CT image-guided radiotherapy. We accelerated the algorithm by introducing an 'active force' along with an adaptive force strength adjustment during the iterative process. These improvements led to a 40% speed improvement over the original algorithm and a high tolerance of large organ deformations. We used three methods to evaluate the accuracy of the algorithm. First, we created a set of mathematical transformations for a series of patient's CT images. This provides a 'ground truth' solution for quantitatively validating the deformable image registration algorithm. Second, we used a physically deformable pelvic phantom, which can measure deformed objects under different conditions. The results of these two tests allowed us to quantify the accuracy of the deformable registration. Validation results showed that more than 96% of the voxels were within 2 mm of their intended shifts for a prostate and a head-and-neck patient case. The mean errors and standard deviations were 0.5 mm+/-1.5 mm and 0.2 mm+/-0.6 mm, respectively. Using the deformable pelvis phantom, the result showed a tracking accuracy of better than 1.5 mm for 23 seeds implanted in a phantom prostate that was deformed by inflation of a rectal balloon. Third, physician-drawn contours outlining the tumour volumes and certain anatomical structures in the original CT images were deformed along with the CT images acquired during subsequent treatments or during a different respiratory phase for a lung cancer case. Visual inspection of the positions and shapes of these deformed contours agreed well with human judgment. Together, these results suggest that the accelerated demons algorithm has significant potential for delineating and tracking doses in targets and critical structures during CT-guided radiotherapy.     

Measurement of cone beam CT coincidence with megavoltage isocentre and image sharpness using the QUASAR Penta-Guide phantom

For image-guided radiotherapy (IGRT) systems based on cone beam CT (CBCT) integrated into a linear accelerator, the reproducible alignment of imager to x-ray source is critical to the registration of both the x-ray-volumetric image with the megavoltage (MV) beam isocentre and image sharpness. An enhanced method of determining the CBCT to MV isocentre alignment using the QUASARtrade mark Penta-Guide phantom was developed which improved both precision and accuracy. This was benchmarked against our existing method which used software and a ball-bearing (BB) phantom provided by Elekta. Additionally, a method of measuring an image sharpness metric (MTF(50)) from the edge response function of a spherical air cavity within the Penta-Guide phantom was developed and its sensitivity was tested by simulating misalignments of the kV imager. Reproducibility testing of the enhanced Penta-Guide method demonstrated a systematic error of <0.2 mm when compared to the BB method with near equivalent random error (s = 0.15 mm). The mean MTF(50) for five measurements was 0.278 +/- 0.004 lp mm(-1) with no applied misalignment. Simulated misalignments exhibited a clear peak in the MTF(50) enabling misalignments greater than 0.4 mm to be detected. The Penta-Guide phantom can be used to precisely measure CBCT-MV coincidence and image sharpness on CBCT-IGRT systems.     

Optimization of megavoltage CT scan registration settings for thoracic cases on helical tomotherapy

This study aims to investigate the settings that provide optimum registration accuracy when registering megavoltage CT (MVCT) studies acquired on tomotherapy with planning kilovoltage CT (kVCT) studies of patients with lung cancer. For each experiment, the systematic difference between the actual and planned positions of the thorax phantom was determined by setting the phantom up at the planning isocenter, generating and registering an MVCT study. The phantom was translated by 5 or 10 mm, MVCT scanned, and registration was performed again. A root-mean-square equation that calculated the residual error of the registration based on the known shift and systematic difference was used to assess the accuracy of the registration process. The phantom study results for 18 combinations of different MVCT/kVCT registration options are presented and compared to clinical registration data from 17 lung cancer patients. MVCT studies acquired with coarse (6 mm), normal (4 mm) and fine (2 mm) slice spacings could all be registered with similar residual errors. No specific combination of resolution and fusion selection technique resulted in a lower residual error. A scan length of 6 cm with any slice spacing registered with the full image fusion selection technique and fine resolution will result in a low residual error most of the time. On average, large corrections made manually by clinicians to the automatic registration values are infrequent. Small manual corrections within the residual error averages of the registration process occur, but their impact on the average patient position is small. Registrations using the full image fusion selection technique and fine resolution of 6 cm MVCT scans with coarse slices have a low residual error, and this strategy can be clinically used for lung cancer patients treated on tomotherapy. Automatic registration values are accurate on average, and a quick verification on a sagittal MVCT slice should be enough to detect registration outliers.     

Standardized platform for coregistration of nonconcurrent diffuse optical and magnetic resonance breast images obtained in different geometries

We present a novel methodology for combining breast image data obtained at different times, in different geometries, and by different techniques. We combine data based on diffuse optical tomography (DOT) and magnetic resonance imaging (MRI). The software platform integrates advanced multimodal registration and segmentation algorithms, requires minimal user experience, and employs computationally efficient techniques. The resulting superposed 3-D tomographs facilitate tissue analyses based on structural and functional data derived from both modalities, and readily permit enhancement of DOT data reconstruction using MRI-derived a-priori structural information. We demonstrate the multimodal registration method using a simulated phantom, and we present initial patient studies that confirm that tumorous regions in a patient breast found by both imaging modalities exhibit significantly higher total hemoglobin concentration (THC) than surrounding normal tissues. The average THC in the tumorous regions is one to three standard deviations larger than the overall breast average THC for all patients.     

Evaluation of deformable image registration and a motion model in CT images with limited features

Deformable image registration (DIR) is increasingly used in radiotherapy applications and provides the basis for a previously described model of patient-specific respiratory motion. We examine the accuracy of a DIR algorithm and a motion model with respiration-correlated CT (RCCT) images of software phantom with known displacement fields, physical deformable abdominal phantom with implanted fiducials in the liver and small liver structures in patient images. The motion model is derived from a principal component analysis that relates volumetric deformations with the motion of the diaphragm or fiducials in the RCCT. Patient data analysis compares DIR with rigid registration as ground truth: the mean ± standard deviation 3D discrepancy of liver structure centroid positions is 2.0 ± 2.2 mm. DIR discrepancy in the software phantom is 3.8 ± 2.0 mm in lung and 3.7 ± 1.8 mm in abdomen; discrepancies near the chest wall are larger than indicated by image feature matching. Marker's 3D discrepancy in the physical phantom is 3.6 ± 2.8 mm. The results indicate that visible features in the images are important for guiding the DIR algorithm. Motion model accuracy is comparable to DIR, indicating that two principal components are sufficient to describe DIR-derived deformation in these datasets.     

Evaluation of cost functions for gray value matching of two-dimensional images in radiotherapy

In external beam radiotherapy, portal imaging is applied for verification of the patient setup. Current automatic methods for portal image registration, which are often based on segmentation of anatomical structures, are especially successful for images of the pelvic region. For portal images of more complicated anatomical structures, e.g., lung, these techniques are less successful. It is desirable to have a method for image registration that is applicable for a wide range of treatment sites. In this study, a registration method for two-dimensional (2D) registration of portal and reference images based on intensity values was tested on portal images of various anatomical sites. Tests were performed with and without preprocessing (unsharp mask filtering followed by histogram equalization) for 96 image pairs and six cost functions. The images were obtained from treatments of the rectum, salivary gland, brain, prostate, and lung. To get insight into the behavior of the various cost functions, cost function values were computed for each portal image for 20,000 transformations of the corresponding reference image, translating the reference image in a range of +/- 1 cm and rotating +/- 10 degrees with respect to the clinical match. The automatic match was defined as the transformation associated with the global minimum (found by an exhaustive search). Without preprocessing, the registration reliability was low (less than 27%). With preprocessing, about 90% of the matches were successful, with a difference with our gold standard (manual registration) of about 1 mm and 1 degree SD. All tested cost functions performed similarly. However, the number of local minima using mutual information was larger than for the other tested cost functions. A cost function based on the mean product of the corresponding pixel values had the least number of local minima. In conclusion, gray value based registration of portal images is applicable for a wide range of treatment sites. However, pre-processing of the images is essential.     

2D/3D registration of endoscopic ultrasound to CT volume data

This paper describes a computer-aided navigation system using image fusion to support endoscopic interventions such as the accurate collection of biopsy specimens. An endoscope provides the physician with real-time ultrasound (US) and a video image. An image slice that corresponds to the corresponding image from the US scan head is derived from a preoperative computed tomography (CT) or magnetic resonance image volume data set using oblique reformatting and displayed side by side with the US image. The position of the image acquired by the US scan head is determined by a miniaturized electromagnetic tracking system (EMTS) after calibrating the endoscope's scan head. The transformation between the patient coordinate system and the preoperative data set is calculated using a 2D/3D registration. This is achieved by calibrating an intraoperative interventional CT slice with an optical tracking system (OTS) using the same algorithm as for the US calibration. The slice is then used for 2D/3D registration with the coordinate system of the preoperative volume. The fiducial registration error (FRE) for the US calibration was 2.0 mm +/- 0.4 mm; the interventional CT FRE was 0.36 +/- 0.12 mm; and the 2D/3D registration target registration error (TRE) was 1.8 +/- 0.3 mm. The point-to-point registration between the OTS and the EMTS had an FRE of 0.9 +/- 0.4 mm. Finally, we found an overall TRE for the complete system to be 3.9 +/- 0.6 mm.     

A fast, accurate, and automatic 2D-3D image registration for image-guided cranial radiosurgery

The authors developed a fast and accurate two-dimensional (2D)-three-dimensional (3D) image registration method to perform precise initial patient setup and frequent detection and correction for patient movement during image-guided cranial radiosurgery treatment. In this method, an approximate geometric relationship is first established to decompose a 3D rigid transformation in the 3D patient coordinate into in-plane transformations and out-of-plane rotations in two orthogonal 2D projections. Digitally reconstructed radiographs are generated offline from a preoperative computed tomography volume prior to treatment and used as the reference for patient position. A multiphase framework is designed to register the digitally reconstructed radiographs with the x-ray images periodically acquired during patient setup and treatment. The registration in each projection is performed independently; the results in the two projections are then combined and converted to a 3D rigid transformation by 2D-3D geometric backprojection. The in-plane transformation and the out-of-plane rotation are estimated using different search methods, including multiresolution matching, steepest descent minimization, and one-dimensional search. Two similarity measures, optimized pattern intensity and sum of squared difference, are applied at different registration phases to optimize accuracy and computation speed. Various experiments on an anthropomorphic head-and-neck phantom showed that, using fiducial registration as a gold standard, the registration errors were 0.33 +/- 0.16 mm (s.d.) in overall translation and 0.29 degrees +/- 0.11 degrees (s.d.) in overall rotation. The total targeting errors were 0.34 +/- 0.16 mm (s.d.), 0.40 +/- 0.2 mm (s.d.), and 0.51 +/- 0.26 mm (s.d.) for the targets at the distances of 2, 6, and 10 cm from the rotation center, respectively. The computation time was less than 3 s on a computer with an Intel Pentium 3.0 GHz dual processor.     

Three-dimensional Elastic Image Registration Based on Strain Energy Minimization: Application to Prostate Magnetic Resonance Imaging

The use of magnetic resonance (MR) imaging in conjunction with an endorectal coil is currently the clinical standard for the diagnosis of prostate cancer because of the increased sensitivity and specificity of this approach. However, imaging in this manner provides images and spectra of the prostate in the deformed state because of the insertion of the endorectal coil. Such deformation may lead to uncertainties in the localization of prostate cancer during therapy. We propose a novel 3-D elastic registration procedure that is based on the minimization of a physically motivated strain energy function that requires the identification of similar features (points, curves, or surfaces) in the source and target images. The Gauss–Seidel method was used in the numerical implementation of the registration algorithm. The registration procedure was validated on synthetic digital images, MR images from prostate phantom, and MR images obtained on patients. The registration error, assessed by averaging the displacement of a fiducial landmark in the target to its corresponding point in the registered image, was 0.2 ± 0.1 pixels on synthetic images. On the prostate phantom and patient data, the registration errors were 1.0 ± 0.6 pixels (0.6 ± 0.4 mm) and 1.8 ± 0.7 pixels (1.1 ± 0.4 mm), respectively. Registration also improved image similarity (normalized cross-correlation) from 0.72 ± 0.10 to 0.96 ± 0.03 on patient data. Registration results on digital images, phantom, and prostate data in vivo demonstrate that the registration procedure can be used to significantly improve both the accuracy of localized therapies such as brachytherapy or external beam therapy and can be valuable in the longitudinal follow-up of patients after therapy.     

PET/CT image registration: preliminary tests for its application to clinical dosimetry in radiotherapy

The quality of dosimetry in radiotherapy treatment requires the accurate delimitation of the gross tumor volume. This can be achieved by complementing the anatomical detail provided by CT images through fusion with other imaging modalities that provide additional metabolic and physiological information. Therefore, use of multiple imaging modalities for radiotherapy treatment planning requires an accurate image registration method. This work describes tests carried out on a Discovery LS positron emission/computed tomography (PET/CT) system by General Electric Medical Systems (GEMS), for its later use to obtain images to delimit the target in radiotherapy treatment. Several phantoms have been used to verify image correlation, in combination with fiducial markers, which were used as a system of external landmarks. We analyzed the geometrical accuracy of two different fusion methods with the images obtained with these phantoms. We first studied the fusion method used by the PET/CT system by GEMS (hardware fusion) on the basis that there is satisfactory coincidence between the reconstruction centers in CT and PET systems; and secondly the fiducial fusion, a registration method, by means of least-squares fitting algorithm of a landmark points system. The study concluded with the verification of the centroid position of some phantom components in both imaging modalities. Centroids were estimated through a calculation similar to center-of-mass, weighted by the value of the CT number and the uptake intensity in PET. The mean deviations found for the hardware fusion method were: deltax/ +/-sigma = 3.3 mm +/- 1.0 mm and /deltax/ +/-sigma = 3.6 mm +/- 1.0 mm. These values were substantially improved upon applying fiducial fusion based on external landmark points: /deltax/ +/-sigma = 0.7 mm +/- 0.8 mm and /deltax/ +/-sigma = 0.3 mm 1.7 mm. We also noted that differences found for each of the fusion methods were similar for both the axial and helical CT image acquisition protocols.     

A deformable phantom for 4D radiotherapy verification: design and image registration evaluation

Motion of thoracic tumors with respiration presents a challenge for three-dimensional (3D) conformal radiation therapy treatment. Validation of techniques aimed at measuring and minimizing the effects of respiratory motion requires a realistic deformable phantom for use as a gold standard. The purpose of this study was to develop and study the characteristics of a reproducible, tissue equivalent, deformable lung phantom. The phantom consists of a Lucite cylinder filled with water containing a latex balloon stuffed with dampened natural sponges. The balloon is attached to a piston that mimics the human diaphragm. Nylon wires and Lucite beads, emulating vascular and bronchial bifurcations, were uniformly glued at various locations throughout the sponges. The phantom is capable of simulating programmed irregular breathing patterns with varying periods and amplitudes. A tissue equivalent tumor, suitable for holding radiochromic film for dose measurements was embedded in the sponge. To assess phantom motion, eight 3D computed tomography data sets of the static phantom were acquired for eight equally spaced positions of the piston. The 3D trajectories of 12 manually chosen point landmarks and the tumor center-of-mass were studied. Motion reproducibility tests of the deformed phantom were established on seven repeat scans of three different states of compression. Deformable image registration (DIR) of the extreme breathing phases was performed. The accuracy of the DIR was evaluated by visual inspection of image overlays and quantified by the distance-to-agreement (DTA) of manually chosen point landmarks and triangulated surfaces obtained from 3D contoured structures. In initial tests of the phantom, a 20-mm excursion of the piston resulted in deformations of the balloon of 20 mm superior-inferior, 4 mm anterior-posterior, and 5 mm left-right. The change in the phantom mean lung density ranged from 0.24 (0.12 SD) g/cm3 at peak exhale to 0.19 (0.12 SD) g/cm3 at peak inhale. The SI displacement of the landmarks varied between 94% and 3% of the piston excursion for positions closer and farther away from the piston, respectively. The reproducibility of the phantom deformation was within the image resolution (0.7 x 0.7 x 1.25 mm3). Vector average registration accuracy based on point landmarks was found to be 0.5 (0.4 SD) mm. The tumor and lung mean 3D DTA obtained from triangulated surfaces were 0.4 (0.1 SD) mm and 1.0 (0.8 SD) mm, respectively. This phantom is capable of reproducibly emulating the physically realistic lung features and deformations and has a wide range of potential applications, including four-dimensional (4D) imaging, evaluation of deformable registration accuracy, 4D planning and dose delivery.     

Model-independent, multimodality deformable image registration by local matching of anatomical features and minimization of elastic energy

With respect to the demands of adaptive and 4D-radiotherapy applications, an algorithm is proposed for a fully automatic, multimodality deformable registration that follows the concept of translational relocation of regularly distributed image subvolumes governed by local anatomical features. Thereby, the problem of global deformable registration is broken down to multiple independent local registration steps which allows for straightforward parallelization of the algorithm. In a subsequent step, possible local misregistrations are corrected for by minimization of the elastic energy of the displacement field under consideration of image information. The final displacement field results from interpolation of the subvolume shift vectors. The algorithm can employ as a similarity measure both the correlation coefficient and mutual information. The latter allows the application to intermodality deformable registration problems. The typical calculation time on a modern multiprocessor PC is well below 1 min, which facilitates almost-interactive, "online" usage. CT-to-MRI and CT-to-cone-beam-CT registrations of head-and-neck data sets are presented, as well as inhale-to-exhale registrations of lung CT data sets. For quantitative evaluation of registration accuracy, a virtual thorax phantom was developed; additionally, a landmark-based evaluation on four lung respiratory-correlated CT data sets was performed. This consistently resulted in average registration residuals on the order of the voxel size or less (3D-residuals approximately 1-2 mm). Summarizing, the presented algorithm allows an accurate multimodality deformable registration with calculation times well below 1 min, and thus bears promise as a versatile basic tool in adaptive and 4D-radiotherapy applications.     

Deformable structure registration of bladder through surface mapping

Cumulative dose distributions in fractionated radiation therapy depict the dose to normal tissues and therefore may permit an estimation of the risk of normal tissue complications. However, calculation of these distributions is highly challenging because of interfractional changes in the geometry of patient anatomy. This work presents an algorithm for deformable structure registration of the bladder and the verification of the accuracy of the algorithm using phantom and patient data. In this algorithm, the registration process involves conformal mapping of genus zero surfaces using finite element analysis, and guided by three control landmarks. The registration produces a correspondence between fractions of the triangular meshes used to describe the bladder surface. For validation of the algorithm, two types of balloons were inflated gradually to three times their original size, and several computerized tomography (CT) scans were taken during the process. The registration algorithm yielded a local accuracy of 4 mm along the balloon surface. The algorithm was then applied to CT data of patients receiving fractionated high-dose-rate brachytherapy to the vaginal cuff, with the vaginal cylinder in situ. The patients' bladder filling status was intentionally different for each fraction. The three required control landmark points were identified for the bladder based on anatomy. Out of an Institutional Review Board (IRB) approved study of 20 patients, 3 had radiographically identifiable points near the bladder surface that were used for verification of the accuracy of the registration. The verification point as seen in each fraction was compared with its predicted location based on affine as well as deformable registration. Despite the variation in bladder shape and volume, the deformable registration was accurate to 5 mm, consistently outperforming the affine registration. We conclude that the structure registration algorithm presented works with reasonable accuracy and provides a means of calculating cumulative dose distributions.     

Automated CT marker segmentation for image registration in radionuclide therapy.

In this paper a novel, automated CT marker segmentation technique for image registration is described. The technique, which is based on analysing each CT slice contour individually, treats the cross sections of the external markers as protrusions of the slice contour. Knowledge-based criteria, using the shape and dimensions of the markers, are defined to enable marker identification and segmentation. Following segmentation, the three-dimensional (3D) markers' centroids are localized using an intensity-weighted algorithm. Finally, image registration is performed using a least-squares fit algorithm. The technique was applied to both simulated and patient studies. The patients were undergoing 131I-mIBG radionuclide therapy with each study comprising several 99mTc single photon emission computed tomography (SPECT) scans and one CT marker scan. The mean residual 3D registration errors (+/- 1 SD) computed for the simulated and patient studies were 1.8 +/- 0.3 mm and 4.3 +/- 0.5 mm respectively.     

Novel image registration quality evaluator (RQE) with an implementation for automated patient positioning in cranial radiation therapy

In external beam radiation therapy, digitally reconstructed radiographs (DRRs) and portal images are used to verify patient setup based either on a visual comparison or, less frequently, with automated registration algorithms. A registration algorithm can be trapped in local optima due to irregularity of patient anatomy, image noise and artifacts, and/or out-of-plane shifts, resulting in an incorrect solution. Thus, human observation, which is subjective, is still required to check the registration result. We propose to use a novel image registration quality evaluator (RQE) to automatically identify misregistrations as part of an algorithm-based decision-making process for verification of patient positioning. A RQE, based on an adaptive pattern classifier, is generated from a pair of reference and target images to determine the acceptability of a registration solution given an optimization process. Here we applied our RQE to patient positioning for cranial radiation therapy. We constructed two RQEs-one for the evaluation of intramodal registrations (i.e., portal-portal); the other for intermodal registrations (i.e., portal-DRR). Mutual information, because of its high discriminatory ability compared with other measures (i.e., correlation coefficient and partitioned intensity uniformity), was chosen as the test function for both RQEs. We adopted 1 mm translation and 1 degree rotation as the maximal acceptable registration errors, reflecting desirable clinical setup tolerances for cranial radiation therapy. Receiver operating characteristic analysis was used to evaluate the performance of the RQE, including computations of sensitivity and specificity. The RQEs showed very good performance for both intramodal and intermodal registrations using simulated and phantom data. The sensitivity and the specificity were 0.973 and 0.936, respectively, for the intramodal RQE using phantom data. Whereas the sensitivity and the specificity were 0.961 and 0.758, respectively, for the intermodal RQE using phantom data. Phantom experiments also indicated our RQEs detected out-of-plane deviations exceeding 2.5 mm and 2.50. A preliminary retrospective clinical study of the RQE on cranial portal imaging also yielded good sensitivity > or = 0.857) and specificity (> or = 0.987). Clinical implementation of a RQE could potentially reduce the involvement of the human observer for routine patient positioning verification, while increasing setup accuracy and reducing setup verification time.     

Automatic 3D-to-2D registration for CT and dual-energy digital radiography for calcification detection

We are investigating three-dimensional (3D) to two-dimensional (2D) registration methods for computed tomography (CT) and dual-energy digital radiography (DEDR). CT is an established tool for the detection of cardiac calcification. DEDR could be a cost-effective alternative screening tool. In order to utilize CT as the "gold standard" to evaluate the capability of DEDR images for the detection and localization of calcium, we developed an automatic, intensity-based 3D-to-2D registration method for 3D CT volumes and 2D DEDR images. To generate digitally reconstructed radiography (DRR) from the CT volumes, we developed several projection algorithms using the fast shear-warp method. In particular, we created a Gaussian-weighted projection for this application. We used normalized mutual information (NMI) as the similarity measurement. Simulated projection images from CT values were fused with the corresponding DEDR images to evaluate the localization of cardiac calcification. The registration method was evaluated by digital phantoms, physical phantoms, and clinical data sets. The results from the digital phantoms show that the success rate is 100% with a translation difference of less than 0.8 mm and a rotation difference of less than 0.2 degrees. For physical phantom images, the registration accuracy is 0.43 +/- 0.24 mm. Color overlay and 3D visualization of clinical images show that the two images registered well. The NMI values between the DRR and DEDR images improved from 0.21 +/- 0.03 before registration to 0.25 +/- 0.03 after registration. Registration errors measured from anatomic markers decreased from 27.6 +/- 13.6 mm before registration to 2.5 +/- 0.5 mm after registration. Our results show that the automatic 3D-to-2D registration is accurate and robust. This technique can provide a useful tool for correlating DEDR with CT images for screening coronary artery calcification.