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 registration of megavoltage to kilovoltage CT images in helical tomotherapy: an evaluation of the setup verification process for the special case of a rigid head phantom

Precise daily target localization is necessary to achieve highly conformal radiation delivery. In helical tomotherapy, setup verification may be accomplished just prior to delivering each fraction by acquiring a megavoltage CT scan of the patient in the treatment position. This daily image set may be manually or automatically registered to the image set on which the treatment plan was calculated, in order to determine any needed adjustments. The system was tested by acquiring 104 MVCT scans of an anthropomorphic head phantom to which translational displacements had been introduced with respect to the planning image set. Registration results were compared against an independent, optically guided positioning system. The total experimental uncertainty was within approximately 1 mm. Although the registration of phantom images is not fully analogous to the registration of patient images, this study confirms that the system is capable of phantom localization with sub-voxel accuracy. In seven registration problems considered, expert human observers were able to perform manual registrations with comparable or inferior accuracy to automatic registration by mutual information. The time to compute an automatic registration is considerably shorter than the time required for manual registration. However, human evaluation of automatic results is necessary in order to identify occasional outliers, and to ensure that the registration is clinically acceptable, especially in the case of deformable patient anatomy.     

Intrathoracic tumour motion estimation from CT imaging using the 3D optical flow method

The purpose of this work was to develop and validate an automated method for intrathoracic tumour motion estimation from breath-hold computed tomography (BH CT) imaging using the three-dimensional optical flow method (3D OFM). A modified 3D OFM algorithm provided 3D displacement vectors for each voxel which were used to map tumour voxels on expiration BH CT onto inspiration BH CT images. A thoracic phantom and simulated expiration/inspiration BH CT pairs were used for validation. The 3D OFM was applied to the measured inspiration and expiration BH CT images from one lung cancer and one oesophageal cancer patient. The resulting displacements were plotted in histogram format and analysed to provide insight regarding the tumour motion. The phantom tumour displacement was measured as 1.20 and 2.40 cm with full-width at tenth maximum (FWTM) for the distribution of displacement estimates of 0.008 and 0.006 cm, respectively. The maximum error of any single voxel's motion estimate was 1.1 mm along the z-dimension or approximately one-third of the z-dimension voxel size. The simulated BH CT pairs revealed an rms error of less than 0.25 mm. The displacement of the oesophageal tumours was nonuniform and up to 1.4 cm, this was a new finding. A lung tumour maximum displacement of 2.4 cm was found in the case evaluated. In conclusion, 3D OFM provided an accurate estimation of intrathoracic tumour motion, with estimated errors less than the voxel dimension in a simulated motion phantom study. Surprisingly, oesophageal tumour motion was large and nonuniform, with greatest motion occurring at the gastro-oesophageal junction.     

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.     

Effects of x-ray and CT image enhancements on the robustness and accuracy of a rigid 3D/2D image registration

A rigid body three-dimensional/two-dimensional (3D/2D) registration method has been implemented using mutual information, gradient ascent, and 3D texturemap-based digitally reconstructed radiographs. Nine combinations of commonly used x-ray and computed tomography (CT) image enhancement methods, including window leveling, histogram equalization, and adaptive histogram equalization, were examined to assess their effects on accuracy and robustness of the registration method. From a set of experiments using an anthropomorphic chest phantom, we were able to draw several conclusions. First, the CT and x-ray preprocessing combination with the widest attraction range was the one that linearly stretched the histograms onto the entire display range on both CT and x-ray images. The average attraction ranges of this combination were 71.3 mm and 61.3 deg in the translation and rotation dimensions, respectively, and the average errors were 0.12 deg and 0.47 mm. Second, the combination of the CT image with tissue and bone information and the x-ray images with adaptive histogram equalization also showed subvoxel accuracy, especially the best in the translation dimensions. However, its attraction ranges were the smallest among the examined combinations (on average 36 mm and 19 deg). Last the bone-only information on the CT image did not show convergency property to the correct registration.     

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.     

Calibration of three-dimensional ultrasound images for image-guided radiation therapy

A new technique of patient positioning for radiotherapy/radiosurgery of extracranial tumours using three-dimensional (3D) ultrasound images has been developed. The ultrasound probe position is tracked within the treatment room via infrared light emitting diodes (IRLEDs) attached to the probe. In order to retrieve the corresponding room position of the ultrasound image, we developed an initial ultrasound probe calibration technique for both 2D and 3D ultrasound systems. This technique is based on knowledge of points in both room and image coordinates. We first tested the performance of three algorithms in retrieving geometrical transformations using synthetic data with different noise levels. Closed form solution algorithms (singular value decomposition and Horn's quaternion algorithms) were shown to outperform the Hooke and Jeeves iterative algorithm in both speed and accuracy. Furthermore, these simulations show that for a random noise level of 2.5, 5, 7.5 and 10 mm, the number of points required for a transformation accuracy better than 1 mm is 25, 100, 200 and 500 points respectively. Finally, we verified the tracking accuracy of this system using a specially designed ultrasound phantom. Since ultrasound images have a high noise level, we designed an ultrasound phantom that provides a large number of points for the calibration. This tissue equivalent phantom is made of nylon wires, and its room position is optically tracked using IRLEDs. By obtaining multiple images through the nylon wires, the calibration technique uses an average of 300 points for 3D ultrasound volumes and 200 for 2D ultrasound images, and its stability is very good for both rotation (standard deviation: 0.4 degrees) and translation (standard deviation: 0.3 mm) transformations. After this initial calibration procedure, the position of any voxel in the ultrasound image volume can be determined in world space, thereby allowing real-time image guidance of therapeutic procedures. Finally, the overall tracking accuracy of our 3D ultrasound image-guided positioning system was measured to be on average 0.2 mm, 0.9 mm and 0.6 mm for the AP, lateral and axial directions respectively.     

Automated 2D-3D registration of portal images and CT data using line-segment enhancement

In prostate radiotherapy, setup errors with respect to the patient's bony anatomy can be reduced by aligning 2D megavoltage (MV) portal images acquired during treatment to a reference 3D kilovoltage (kV) CT acquired for treatment planning purposes. The purpose of this study was to evaluate a fully automated 2D-3D registration algorithm to quantify setup errors in 3D through the alignment of line-enhanced portal images and digitally reconstructed radiographs computed from the CT. The line-enhanced images were obtained by correlating the images with a filter bank of short line segments, or "sticks" at different orientations. The proposed methods were validated on (1) accurately collected gold-standard data consisting of a 3D kV cone-beam CT scan of an anthropomorphic phantom of the pelvis and 2D MV portal images in the anterior-posterior (AP) view acquired at 15 different poses and (2) a conventional 3D kV CT scan and weekly 2D MV AP portal images of a patient over 8 weeks. The mean (and standard deviation) of the absolute registration error for rotations around the right-lateral (RL), inferior-superior (IS), and posterior-anterior (PA) axes were 0.212 degree (0.214 degree), 0.055 degree (0.033 degree) and 0.041 degree (0.039 degree), respectively. The corresponding registration errors for translations along the RL, IS, and PA axes were 0.161 (0.131) mm, 0.096 (0.033) mm, and 0.612 (0.485) mm. The mean (and standard deviation) of the total registration error was 0.778 (0.543) mm. Registration on the patient images was successful in all eight cases as determined visually. The results indicate that it is feasible to automatically enhance features in MV portal images of the pelvis for use within a completely automated 2D-3D registration framework for the accurate determination of patient setup errors. They also indicate that it is feasible to estimate all six transformation parameters from a 3D CT of the pelvis and a single portal image in the AP view.     

Cone beam CT with zonal filters for simultaneous dose reduction, improved target contrast and automated set-up in radiotherapy

Cone beam CT (CBCT) using a zonal filter is introduced. The aims are reduced concomitant imaging dose to the patient, simultaneous control of body scatter for improved image quality in the tumour target zone and preserved set-up detail for radiotherapy. Aluminium transmission diaphragms added to the CBCT x-ray tube of the Elekta Synergytrade mark linear accelerator produced an unattenuated beam for a central "target zone" and a partially attenuated beam for an outer "set-up zone". Imaging doses and contrast noise ratios (CNR) were measured in a test phantom for transmission diaphragms 12 and 24 mm thick, for 5 and 10 cm long target zones. The effect on automatic registration of zonal CBCT to conventional CT was assessed relative to full-field and lead-collimated images of an anthropomorphic phantom. Doses along the axis of rotation were reduced by up to 50% in both target and set-up zones, and weighted dose (two thirds surface dose plus one third central dose) was reduced by 10-20% for a 10 cm long target zone. CNR increased by up to 15% in zonally filtered CBCT images compared to full-field images. Automatic image registration remained as robust as that with full-field images and was superior to CBCT coned down using lead-collimation. Zonal CBCT significantly reduces imaging dose and is expected to benefit radiotherapy through improved target contrast, required to assess target coverage, and wide-field edge detail, needed for robust automatic measurement of patient set-up error.     

Automatic CT-SPECT registration of livers treated with radioactive microspheres

We have developed an automated technique to accurately register the CT and SPECT scans of the liver of patients treated with radioactive microspheres for tumour targeting assessment. An anthropomorphic phantom was used to validate the accuracy of the registration algorithm. The phantom liver and three fiducial markers placed on its surface were filled with 99mTc. The phantom was scanned with CT and SPECT scanners in different positions. The liver volume was contoured on the CT scans from which a three-dimensional liver mask was created. By constraining the liver volume to the volume obtained from the CT scans, another liver mask was automatically created from the SPECT images. An adaptive simulated annealing algorithm was used to minimize the difference between the two volumes formed by the two sets of masks. The algorithm involved rigid transformation of the SPECT mask to reach the optimization goal. The accuracy of the algorithm was evaluated by the superposition of the fiducial markers on the CT and SPECT. The registered SPECT overlaid on the CT scan of the phantom showed congruence of the fiducial markers on CT and SPECT images within 1 mm. The technique was applied to a patient image set who received the microsphere infusion procedure. The registered CT and SPECT images of the patient showed that the majority of the activity was concentrated in the tumour, indicating a successful tumour targeting.     

Acceptance Test of a Commercially Available Software for Automatic Image Registration of Computed Tomography (CT), Magnetic Resonance Imaging (MRI) And 99mTc-methoxyisobutylisonitrile (MIBI) Single-Photon Emission Computed Tomography (SPECT) Brain Images

This note describes a method to characterize the performances of image fusion software (Syntegra) with respect to accuracy and robustness. Computed tomography (CT), magnetic resonance imaging (MRI), and single-photon emission computed tomography (SPECT) studies were acquired from two phantoms and 10 patients. Image registration was performed independently by two couples composed of one radiotherapist and one physicist by means of superposition of anatomic landmarks. Each couple performed jointly and saved the registration. The two solutions were averaged to obtain the gold standard registration. A new set of estimators was defined to identify translation and rotation errors in the coordinate axes, independently from point position in image field of view (FOV). Algorithms evaluated were local correlation (LC) for CT-MRI, normalized mutual information (MI) for CT-MRI, and CT-SPECT registrations. To evaluate accuracy, estimator values were compared to limiting values for the algorithms employed, both in phantoms and in patients. To evaluate robustness, different alignments between images taken from a sample patient were produced and registration errors determined. LC algorithm resulted accurate in CT-MRI registrations in phantoms, but exceeded limiting values in 3 of 10 patients. MI algorithm resulted accurate in CT-MRI and CT-SPECT registrations in phantoms; limiting values were exceeded in one case in CT-MRI and never reached in CT-SPECT registrations. Thus, the evaluation of robustness was restricted to the algorithm of MI both for CT-MRI and CT-SPECT registrations. The algorithm of MI proved to be robust: limiting values were not exceeded with translation perturbations up to 2.5 cm, rotation perturbations up to 10° and roto-translational perturbation up to 3 cm and 5°.     

Automated 2D-3D registration of a radiograph and a cone beam CT using line-segment enhancement

The objective of this study was to develop a fully automated two-dimensional (2D)-three-dimensional (3D) registration framework to quantify setup deviations in prostate radiation therapy from cone beam CT (CBCT) data and a single AP radiograph. A kilovoltage CBCT image and kilovoltage AP radiograph of an anthropomorphic phantom of the pelvis were acquired at 14 accurately known positions. The shifts in the phantom position were subsequently estimated by registering digitally reconstructed radiographs (DRRs) from the 3D CBCT scan to the AP radiographs through the correlation of enhanced linear image features mainly representing bony ridges. Linear features were enhanced by filtering the images with "sticks," short line segments which are varied in orientation to achieve the maximum projection value at every pixel in the image. The mean (and standard deviations) of the absolute errors in estimating translations along the three orthogonal axes in millimeters were 0.134 (0.096) AP(out-of-plane), 0.021 (0.023) ML and 0.020 (0.020) SI. The corresponding errors for rotations in degrees were 0.011 (0.009) AP, 0.029 (0.016) ML (out-of-plane), and 0.030 (0.028) SI (out-of-plane). Preliminary results with megavoltage patient data have also been reported. The results suggest that it may be possible to enhance anatomic features that are common to DRRs from a CBCT image and a single AP radiography of the pelvis for use in a completely automated and accurate 2D-3D registration framework for setup verification in prostate radiotherapy. This technique is theoretically applicable to other rigid bony structures such as the cranial vault or skull base and piecewise rigid structures such as the spine.     

Patient position verification using CT images.

The use of ions in the radiotherapy of cancer patients requires an accurate patient positioning in order to exploit its potential benefits. Using CT images as the basis for the setup verification offers the advantage of a high in-plane resolution in combination with a geometrically accurate, volumetric information. Before each fraction a single CT slice is acquired at the isocenter level after the positioning procedure. This single slice is registered to the planning CT cube using automated image registration algorithms. Thus any erreonous translation or rotation can be detected and quantified. The registration process involves the interpolation of the volumetric data, the calculation of an energy function, and the minimization of this energy function. Several data interpolation functions as well as minimization algorithms were compared. CT studies with a head phantom were performed in which defined translations and rotations were simulated by moving a motor-driven treatment chair. Different slice thicknesses and anatomical sites were studied to investigate their potential influence on the registration accuracy. The accuracy of the registration was found to be a fraction of a voxel size for suitable combinations of algorithms (typically better than 0.16 mm/deg). A significant dependancy of the registration accuracy on the CT slice thickness and the anatomical site was found (the accuracy ranges from 0.05 mm/deg to 0.16 mm/deg depending on the site). The calculation time is dependant on the used algorithms and the magnitude of the setup error. For the standard combination of algorithms as proposed by the authors (Downhill Simplex minimization with Trilinear interpolation) the typical calculation time is about 20 s for a Sun UltraSPARC processor. Taking into account the mechanical accuracy of the setup device (motor-driven chair) the registration of CT images is thus a useful tool for detecting and quantifying any significant error in the patient position.     

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.     

In-bore setup and software for 3T MRI-guided transperineal prostate biopsy

MRI-guided prostate biopsy in conventional closed-bore scanners requires transferring the patient outside the bore during needle insertion due to the constrained in-bore space, causing a safety hazard and limiting image feedback. To address this issue, we present our custom-made in-bore setup and software to support MRI-guided transperineal prostate biopsy in a wide-bore 3 T MRI scanner. The setup consists of a specially designed tabletop and a needle-guiding template with a Z-frame that gives a physician access to the perineum of the patient at the imaging position and allows the physician to perform MRI-guided transperineal biopsy without moving the patient out of the scanner. The software and Z-frame allow registration of the template, target planning and biopsy guidance. Initially, we performed phantom experiments to assess the accuracy of template registration and needle placement in a controlled environment. Subsequently, we embarked on our clinical trial (N = 10). The phantom experiments showed that the translational errors of the template registration along the right-left (RP) and anterior-posterior (AP) axes were 1.1 ± 0.8 and 1.4 ± 1.1 mm, respectively, while the rotational errors around the RL, AP and superior-inferior axes were (0.8 ± 1.0)°, (1.7 ± 1.6)° and (0.0 ± 0.0)°, respectively. The 2D root-mean-square (RMS) needle-placement error was 3 mm. The clinical biopsy procedures were safely carried out in all ten clinical cases with a needle-placement error of 5.4 mm (2D RMS). In conclusion, transperineal prostate biopsy in a wide-bore 3T scanner is feasible using our custom-made tabletop setup and software, which supports manual needle placement without moving the patient out of the magnet.     

The influence of beam model differences in the comparison of dose calculation algorithms for lung cancer treatment planning

In this study, we show that beam model differences play an important role in the comparison of does calculated with various algorithms for lung cancer treatment planning. These differences may impact the accurate correlation of dose with clinical outcome. To accomplish this, we modified the beam model penumbral parameters in an equivalent path length (EPL) algorithm and subsequently compared the EPL doses with those generated with Monte Carlo (MC). A single AP beam was used for beam fitting. Two different beam models were generated for EPL calculations: (1) initial beam model (init_fit) and (2) optimized beam model (best_fit) , with parameters optimized to produce the best agreement with MC calculated profiles at several depths in a water phantom. For the 6 MV, AP beam, EPL(init_fit) calculations were on average within 2%/2 mm (1.4 mm max.) agreement with MC; the agreement for EPL(best_fit) was 2%/1.0 mm (1.3 mm max.) for EPL(best_fit). Treatment planning was performed using a realistic lung phantom using 6 and 15 MV photons. In all homogeneous phantom plans, EPL(best_fit) calculations were in better agreement with MC. In the heterogeneous 6 MV plan, differences between EPL(best_fit and init_fit) and MC were significant for the tumour. The EPL(init_fit), unlike the EPL(best_fit) calculation, showed large differences in the lung relative to MC. For the 15 MV heterogeneous plan, clinically important differences were found between EPL(best_fit or init_fit) and MC for tumour and lung, suggesting that the algorithmic difference in inhomogeneous cases, differences between EPL(best_fit) and MC for lung tissues were smaller compared to those between EPL(init_fit) and MC. Although the extent to which beam model differences impact the dose comparisons will be dependent upon beam parameters (orientation, field size and energy), and the size and location of the tumour, this study shows that failing to correctly account for beam model differences will lead to biased comparisons between dose algorithms. This may ultimately hinder our ability to accurately correlate dose with clinical outcome.     

Assessment of rigid multi-modality image registration consistency using the multiple sub-volume registration (MSR) method

Registration of different imaging modalities such as CT, MRI, functional MRI (fMRI), positron (PET) and single photon (SPECT) emission tomography is used in many clinical applications. Determining the quality of any automatic registration procedure has been a challenging part because no gold standard is available to evaluate the registration. In this note we present a method, called the 'multiple sub-volume registration' (MSR) method, for assessing the consistency of a rigid registration. This is done by registering sub-images of one data set on the other data set, performing a crude non-rigid registration. By analysing the deviations (local deformations) of the sub-volume registrations from the full registration we get a measure of the consistency of the rigid registration. Registration of 15 data sets which include CT, MR and PET images for brain, head and neck, cervix, prostate and lung was performed utilizing a rigid body registration with normalized mutual information as the similarity measure. The resulting registrations were classified as good or bad by visual inspection. The resulting registrations were also classified using our MSR method. The results of our MSR method agree with the classification obtained from visual inspection for all cases (p < 0.02 based on ANOVA of the good and bad groups). The proposed method is independent of the registration algorithm and similarity measure. It can be used for multi-modality image data sets and different anatomic sites of the patient.     

Image registration of a moving target phantom with helical tomotherapy: effect of the CT acquisition technique and action level proposal

This study aims to quantify the effects of target motion and resultant motion artifacts in planning and megavoltage CT (MVCT) studies on the automatic registration processes of helical tomotherapy. Clinical and experimental data were used to derive an action level for patient repositioning on helical tomotherapy. Planning CT studies of a respiratory motion phantom were acquired using conventional and four-dimensional CT (4D CT) techniques. MVCT studies were acquired on helical tomotherapy in the presence and absence of target motion and were registered with different planning CT studies. The residual errors of the registration process were calculated from the registration values to quantify the ability of the process to detect 5 or 10 mm translations of the phantom in two directions. Twenty-seven registration combinations of MVCT inter-slice spacing, technique and resolution were investigated. The residual errors were used as an estimate of the localization error of the registration process, and the accuracy of couch repositioning was determined from couch position measurements during 866 treatment fractions. These two parameters were used to calculate the action level for patient repositioning on helical tomotherapy. Automatic registration of an MVCT study with 0% breathing phase, average intensity and maximum intensity 4D CT projections did not differ from that of an MVCT study with a conventional planning CT. Motion artifacts in the MVCT or planning CT studies changed the accuracy of the automatic registration process by less than 2.0%. The action level for patient repositioning using MVCT studies of 6 mm inter-slice spacing was determined to be 0.7, 1.1 and 0.6 mm in the x-, y- and z-directions, respectively. These action levels have the greatest effect on treatments for disease sites in the brain.     

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.     

Fast free-form deformable registration via calculus of variations

In this paper, we present a fully automatic, fast and accurate deformable registration technique. This technique deals with free-form deformation. It minimizes an energy functional that combines both similarity and smoothness measures. By using calculus of variations, the minimization problem was represented as a set of nonlinear elliptic partial differential equations (PDEs). A Gauss-Seidel finite difference scheme is used to iteratively solve the PDE. The registration is refined by a multi-resolution approach. The whole process is fully automatic. It takes less than 3 min to register two three-dimensional (3D) image sets of size 256 x 256 x 61 using a single 933 MHz personal computer. Extensive experiments are presented. These experiments include simulations, phantom studies and clinical image studies. Experimental results show that our model and algorithm are suited for registration of temporal images of a deformable body. The registration of inspiration and expiration phases of the lung images shows that the method is able to deal with large deformations. When applied to the daily CT images of a prostate patient, the results show that registration based on iterative refinement of displacement field is appropriate to describe the local deformations in the prostate and the rectum. Similarity measures improved significantly after the registration. The target application of this paper is for radiotherapy treatment planning and evaluation that incorporates internal organ deformation throughout the course of radiation therapy. The registration method could also be equally applied in diagnostic radiology.