Registration of electronic portal images for patient set-up verification

Computerized estimation of patient setup errors in portal images based on localized pelvic templates for prostate cancer radiotherapy

We have developed a computerized method for estimating patient setup errors in portal images based on localized pelvic templates for prostate cancer radiotherapy. The patient setup errors were estimated based on a template-matching technique that compared the portal image and a localized pelvic template image with a clinical target volume produced from a digitally reconstructed radiography (DRR) image of each patient. We evaluated the proposed method by calculating the residual error between the patient setup error obtained by the proposed method and the gold standard setup error determined by consensus between two radiation oncologists. Eleven training cases with prostate cancer were used for development of the proposed method, and then we applied the method to 10 test cases as a validation test. As a result, the residual errors in the anterior-posterior, superior-inferior and left-right directions were smaller than 2 mm for the validation test. The mean residual error was 2.65 ± 1.21 mm in the Euclidean distance for training cases, and 3.10 ± 1.49 mm for the validation test. There was no statistically significant difference in the residual error between the test for training cases and the validation test (P = 0.438). The proposed method appears to be robust for detecting patient setup error in the treatment of prostate cancer radiotherapy.

Validation for 2D/3D registration. I: A new gold standard data set

PURPOSE

In this article, the authors propose a new gold standard data set for the validation of two-dimensional/three-dimensional (2D/3D) and 3D/3D image registration algorithms.

METHODS

A gold standard data set was produced using a fresh cadaver pig head with attached fiducial markers. The authors used several imaging modalities common in diagnostic imaging or radiotherapy, which include 64-slice computed tomography (CT), magnetic resonance imaging using T1, T2, and proton density sequences, and cone beam CT imaging data. Radiographic data were acquired using kilovoltage and megavoltage imaging techniques. The image information reflects both anatomy and reliable fiducial marker information and improves over existing data sets by the level of anatomical detail, image data quality, and soft-tissue content. The markers on the 3D and 2D image data were segmented using ANALYZE 10.0 (AnalyzeDirect, Inc., Kansas City, KN) and an in-house software.

RESULTS

The projection distance errors and the expected target registration errors over all the image data sets were found to be less than 2.71 and 1.88 mm, respectively.

CONCLUSIONS

The gold standard data set, obtained with state-of-the-art imaging technology, has the potential to improve the validation of 2D/3D and 3D/3D registration algorithms for image guided therapy.

DRR and portal image registration for automatic patient positioning in radiotherapy treatment

Image processing turns out to be essential in the planning and verification of radiotherapy treatments. Before applying a radiotherapy treatment, a dosimetry planning must be performed. Usually, the planning is done by means of an X-ray volumetric analysis using computerized tomography, where the area to be radiated is marked out. During the treatment phase, it is necessary to place the patient under the particle accelerator exactly as considered in the dosimetry stage. Coarse alignment is achieved using fiduciary markers placed over the patient's skin as external references. Later, fine alignment is provided by comparing a digitally reconstructed radiography (DRR) from the planning stage and a portal image captured by the accelerator in the treatment stage. The preprocessing of DRR and portal images, as well as the minimization of the non-shared information between both kinds of images, is mandatory for the correct operation of the image registration algorithm. With this purpose, mathematical morphology and image processing techniques have been used. The present work describes a fully automatic method to calculate more accurately the necessary displacement of the couch to place the patient exactly at the planned position. The proposed method to achieve the correct positioning of the patient is based on advanced image registration techniques. Preliminary results show a perfect match with the displacement estimated by the physician.

A quality assurance phantom for electronic portal imaging devices

Electronic portal imaging device (EPID) plays an important role in radiation therapy portal imaging, geometric and dosimetric verification. Consistent image quality and stable radiation response is necessary for proper utilization that requires routine quality assurance (QA). A commercial ‘EPID QC’ phantom weighing 3.8 kg with a dimension of 25×25×4.8 cm3 is used for EPID QA. This device has five essential tools to measure the geometric accuracy, signal‐to‐noise ratio (SNR), dose linearity, and the low‐ and the high‐contrast resolutions. It is aligned with beam divergence to measure the imaging and geometric parameters in both X and Y directions, and can be used as a baseline check for routine QA. The low‐contrast tool consists of a series of holes with various diameters and depths in an aluminum slab, very similar to the Las Vegas phantom. The high‐resolution contrast tool provides the modulation transfer function (MTF) in both the x‐ and y‐dimensions to measure the focal spot of linear accelerator that is important for imaging and small field dosimetry. The device is tested in different institutions with various amorphous silicon imagers including Elekta, Siemens and Varian units. Images of the QA phantom were acquired at 95.2 cm source‐skin‐distance (SSD) in the range 1–15 MU for a 26×26 cm2 field and phantom surface is set normal to the beam direction when gantry is at 0° and 90°. The epidSoft is a software program provided with the EPID QA phantom for analysis of the data. The preliminary results using the phantom on the tested EPID showed very good low‐contrast resolution and high resolution, and an MTF (0.5) in the range of 0.3–0.4 lp/mm. All imagers also exhibit satisfactory geometric accuracy, dose linearity and SNR, and are independent of MU and spatial orientations. The epidSoft maintains an image analysis record and provides a graph of the temporal variations in imaging parameters. In conclusion, this device is simple to use and provides testing on basic and advanced imaging parameters for daily QA on any imager used in clinical practices.

PACS number: 87.57 C‐, 87.57 N‐

Evaluation of Inter-fractional Setup Shifts for Site-specific Helical Tomotherapy Treatments

This paper proposes to summarize and analyze the daily patient setup shifts based on megavoltage computed tomography (MVCT) image registration results for Helical TomoTherapy® (HT) treatment. One hundred and fifty-five consecutive treatment plans for a total of 137 patients delivered by the HT unit through one year were collected in this study. The patient data included pelvis (26%), abdomen (23%), lung (21%), head and neck (10%), prostate (8%), and others (12%). All the translational and roll rotational shifts made via auto MVCT and kilovoltage computed tomography (kVCT) image registration were recorded at each fraction. Manual fine-tuning was followed if automatic registration result was not satisfactory. The mean shift ± one standard deviation (1 SD) was calculated for each patient based on the entire treatment course. For each treatment site, the average shift was analyzed as well as displacement in 3D vector. Statistical tests were performed to analyze the relationship of patient-specific, tumor site-specific, and fraction number association with the patient setup shifts. For all the treatment sites, the largest average shift was found in the anterior-posterior direction. The population standard deviations were between 1.2 and 5.6 mm for the X, Y, and Z directions and ranged from 0.2 to 0.6 degrees for the roll rotational correction. The largest standard deviations of the setup reproducibility in X, Y, and Z directions were found in lung patients (4.2 mm), abdomen, lung and spine patients (4.4 mm), and prostate patients (5.6 mm), respectively. The maximum 3D displacement was 10.9 mm for prostate patients' setup. ANOVA tests demonstrated the setup shifts were statistically different between patients even for those that were treated at the same tumor site in the translational directions. No strong correlation between the setup and the fraction number was found. In conclusion, the MVCT guided function in the HT treatment enables us to generate relatively accurate daily setup through registration with KVCT data sets. Our results indicate that lung, prostate, and abdominal patients are more prone to setup uncertainty and should be carefully evaluated.

ESTERR-PRO: a setup verification software system using electronic portal imaging

The purpose of the paper is to present and evaluate the performance of a new software-based registration system for patient setup verification, during radiotherapy, using electronic portal images. The estimation of setup errors, using the proposed system, can be accomplished by means of two alternate registration methods. (a) The portal image of the current fraction of the treatment is registered directly with the reference image (digitally reconstructed radiograph (DRR) or simulator image) using a modified manual technique. (b) The portal image of the current fraction of the treatment is registered with the portal image of the first fraction of the treatment (reference portal image) by applying a nearly automated technique based on self-organizing maps, whereas the reference portal has already been registered with a DRR or a simulator image. The proposed system was tested on phantom data and on data from six patients. The root mean square error (RMSE) of the setup estimates was 0.8 ± 0.3 (mean value ± standard deviation) for the phantom data and 0.3 ± 0.3 for the patient data, respectively, by applying the two methodologies. Furthermore, statistical analysis by means of the Wilcoxon nonparametric signed test showed that the results that were obtained by the two methods did not differ significantly (P value > 0.05).

Fast DRR splat rendering using common consumer graphics hardware

Digitally rendered radiographs (DRR) are a vital part of various medical image processing applications such as 2D/3D registration for patient pose determination in image-guided radiotherapy procedures. This paper presents a technique to accelerate DRR creation by using conventional graphics hardware for the rendering process. DRR computation itself is done by an efficient volume rendering method named wobbled splatting. For programming the graphics hardware, NVIDIAs C for Graphics (Cg) is used. The description of an algorithm used for rendering DRRs on the graphics hardware is presented, together with a benchmark comparing this technique to a CPU-based wobbled splatting program. Results show a reduction of rendering time by about 70%-90% depending on the amount of data. For instance, rendering a volume of 2 x 10(6) voxels is feasible at an update rate of 38 Hz compared to 6 Hz on a common Intel-based PC using the graphics processing unit (GPU) of a conventional graphics adapter. In addition, wobbled splatting using graphics hardware for DRR computation provides higher resolution DRRs with comparable image quality due to special processing characteristics of the GPU. We conclude that DRR generation on common graphics hardware using the freely available Cg environment is a major step toward 2D/3D registration in clinical routine.

Analysis of the rigid and deformable component of setup inaccuracies on portal images in head and neck radiotherapy

The issue of setup errors consisting of translation, rotation and deformation components in head and neck radiotherapy is addressed with a piecewise registration of small independent regions on a portal image to their reference position. These rectangular regions are termed featurelets as they contain relevant anatomical features. The resulting displacement vectors of each featurelet reflect both the center-of-mass (COM), i.e. the rigid, and the non-rigid component of the setup error. The displacement vectors of a series of daily portal images were subjected to a principal component analysis. In addition to the mean, systematic displacement of each featurelet, this analysis yields correlated patterns of anatomical deformations. Hence, the physiological movements of an individual patient can be obtained without a biomechanical model. It is shown that in the presence of setup errors that are due to rotations or deformations a correction by the COM displacement may deteriorate the error of parts of the anatomy further. The featurelet analysis can be used to refine setup correction protocols, tune spatially variable setup margins in treatment planning and optimize patient immobilization devices.

Evaluation of the precision of portal-image-guided head-and-neck localization: An intra- and interobserver study

There is increasing evidence that, for some patients, image-guided intensity-modulated radiation therapy (IMRT) for head-and-neck cancer patients may maintain target dose coverage and critical organ (e.g., parotids) dose closer to the planned doses than setup using lasers alone. We investigated inter- and intraobserver uncertainties in patient setup in head-and-neck cancer patients. Twenty-two sets of orthogonal digital portal images (from five patients) were selected from images used for daily localization of head-and-neck patients treated with IMRT. To evaluate interobserver variations, five radiation therapists compared the portal images with the plan digitally reconstructed radiographs and reported shifts for the isocenter (approximately C2) and for a supraclavicular reference point. One therapist repeated the procedure a month later to evaluate intraobserver variations. The procedure was then repeated with teams of two therapists. The frequencies for which agreement between the shift reported by the observer and the daily mean shift (average of all observers for a given image set) were less than 1.5 and 2.5 mm were calculated. Standard errors of measurement for the intra- and interobserver uncertainty (SEMintra and SEMinter) for the individual and teams were calculated. The data showed that there was very little difference between individual therapists and teams. At isocenter, 80%-90% of all reported shifts agreed with the daily average within 1.5 mm, showing consistency in the ways both individuals and teams interpret the images (SEMinter approximately 1 mm). This dropped to 65% for the supraclavicular point (SEMinter approximately 1.5 mm). Uncertainties increased for larger setup errors. In conclusion, image-guided patient positioning allows head-and-neck patients to be controlled within 3-4 mm. This is similar to the setup uncertainties found for most head-and-neck patients, but may provide some improvement for the subset of patients with larger setup uncertainties.

Registration of DRRs and portal images for verification of stereotactic body radiotherapy: a feasibility study in lung cancer treatment

Image guidance has become a pre-requisite for hypofractionated radiotherapy where the applied dose per fraction is increased. Particularly in stereotactic body radiotherapy (SBRT) for lung tumours, one has to account for set-up errors and intrafraction tumour motion. In our feasibility study, we compared digitally reconstructed radiographs (DRRs) of lung lesions with MV portal images (PIs) to obtain the displacement of the tumour before irradiation. The verification of the tumour position was performed by rigid intensity based registration and three different merit functions such as the sum of squared pixel intensity differences, normalized cross correlation and normalized mutual information. The registration process then provided a translation vector that defines the displacement of the target in order to align the tumour with the isocentre. To evaluate the registration algorithms, 163 test images were created and subsequently, a lung phantom containing an 8 cm(3) tumour was built. In a further step, the registration process was applied on patient data, containing 38 tumours in 113 fractions. To potentially improve registration outcome, two filter types (histogram equalization and display equalization) were applied and their impact on the registration process was evaluated. Generated test images showed an increase in successful registrations when applying a histogram equalization filter whereas the lung phantom study proved the accuracy of the selected algorithms, i.e. deviations of the calculated translation vector for all test algorithms were below 1 mm. For clinical patient data, successful registrations occurred in about 59% of anterior-posterior (AP) and 46% of lateral projections, respectively. When patients with a clinical target volume smaller than 10 cm(3) were excluded, successful registrations go up to 90% in AP and 50% in lateral projection. In addition, a reliable identification of the tumour position was found to be difficult for clinical target volumes at the periphery of the lung, close to backbone or diaphragm. Moreover, tumour movement during shallow breathing strongly influences image acquisition for patient positioning. Recapitulating, 2D/3D image registration for lung tumours is an attractive alternative compared to conventional CT verification of the tumour position. Nevertheless, size and location of the tumour are limiting parameters for an accurate registration process.

3D interfractional patient position verification using 2D-3D registration of orthogonal images

Reproducible positioning of the patient during fractionated external beam radiation therapy is imperative to ensure that the delivered dose distribution matches the planned one. In this paper, we expand on a 2D-3D image registration method to verify a patient's setup in three dimensions (rotations and translations) using orthogonal portal images and megavoltage digitally reconstructed radiographs (MDRRs) derived from CT data. The accuracy of 2D-3D registration was improved by employing additional image preprocessing steps and a parabolic fit to interpolate the parameter space of the cost function utilized for registration. Using a humanoid phantom, precision for registration of three-dimensional translations was found to be better than 0.5 mm (1 s.d.) for any axis when no rotations were present. Three-dimensional rotations about any axis were registered with a precision of better than 0.2 degrees (1 s.d.) when no translations were present. Combined rotations and translations of up to 4 degrees and 15 mm were registered with 0.4 degrees and 0.7 mm accuracy for each axis. The influence of setup translations on registration of rotations and vice versa was also investigated and mostly agrees with a simple geometric model. Additionally, the dependence of registration accuracy on three cost functions, angular spacing between MDRRs, pixel size, and field-of-view, was examined. Best results were achieved by mutual information using 0.5 degrees angular spacing and a 10 x 10 cm2 field-of-view with 140 x 140 pixels. Approximating patient motion as rigid transformation, the registration method is applied to two treatment plans and the patients' setup errors are determined. Their magnitude was found to be < or = 6.1 mm and < or = 2.7 degrees for any axis in all of the six fractions measured for each treatment plan.

A novel method for automatic detection of patient out-of-plane rotation by comparing a single portal image to a reference image

A novel method for detecting out-of-plane patient rotation by comparing a single portal image to its reference image is presented. Out-of-plane rotation results in an apparent distortion of the anatomy in a portal image. This distortion can be mathematically predicted with the magnification varying at each point in the image. While scaling of points at equal depth is invariant under in-plane rotation or translation, and changes equally in both dimensions for an axial shift of the patient, a change of scaling in only one dimension can be ascribed to an out-of-plane rotation. For the two conditions that are used in this study, it is shown that out-of-plane rotation yields a different scaling of the image in two perpendicular directions and therefore it is feasible to calculate the scale factors as a function of out-of-plane rotation. Conversely the recovery of scale factors in two different directions at the same time would enable the magnitude of the out-of-plane rotation to be recovered. The properties of the Fourier transform of the image are used to align the portal image with the reference image (a simulator image or first approved portal image) prior to the recovery of the scale factors. Correlating the Fourier transform of the portal image on a log-scale with that of the reference image enables the scale factors to be automatically extracted from a single portal image. In the two approaches investigated, out-of-plane rotations of up to 41 degrees and 20 degrees (respectively) have been recovered with a maximum error of 2.4 degrees. This technique could be used to automatically detect patient roll or tilt prior to or during a treatment session.

Clinical use of electronic portal imaging

Accurate and routine target localization is necessary for successful outcome in radiation therapy treatments. Electronic portal imaging devices (EPIDs) provide an advanced tool with digital technology to improve target localization and maintain clinical efficiency. EPIDs are ubiquitous in the radiation therapy clinic, and they provide a powerful and flexible tool to collect and process data in a quantitative manner to improve treatment accuracy for virtually any treatment site. This manuscript provides an overview of the clinical implementation process for effective use of EPIDs. It continues with a review of correction strategies and finally highlights numerous examples of effective clinical application of EPID.

Intracranial stereotactic positioning systems: Report of the American Association of Physicists in Medicine Radiation Therapy Committee Task Group No. 68

Intracranial stereotactic positioning systems (ISPSs) are used to position patients prior to precise radiation treatment of localized lesions of the brain. Often, the lesion is located in close proximity to critical anatomic features whose functions should be maintained. Many types of ISPSs have been described in the literature and are commercially available. These are briefly reviewed. ISPS systems provide two critical functions. The first is to establish a coordinate system upon which a guided therapy can be applied. The second is to provide a method to reapply the coordinate system to the patient such that the coordinates assigned to the patient's anatomy are identical from application to application. Without limiting this study to any particular approach to ISPSs, this report introduces nomenclature and suggests performance tests to quantify both the stability of the ISPS to map diagnostic data to a coordinate system, as well as the ISPS's ability to be realigned to the patient's anatomy. For users who desire to develop a new ISPS system, it may be necessary for the clinical team to establish the accuracy and precision of each of these functions. For commercially available systems that have demonstrated an acceptable level of accuracy and precision, the clinical team may need to demonstrate local ability to apply the system in a manner consistent with that employed during the published testing. The level of accuracy and precision required of an individual ISPS system is dependent upon the clinical protocol (e.g., fractionation, margin, pathology, etc.). Each clinical team should provide routine quality assurance procedures that are sufficient to support the assumptions of accuracy and precision used during the planning process. The testing of ISPS systems can be grouped into two broad categories, type testing, which occurs prior to general commercialization, and site testing, performed when a commercial system is installed at a clinic. Guidelines to help select the appropriate tests as well as recommendations to help establish the required frequency of testing are provided. Because of the broad scope of different systems, it is important that both the manufacturer and user rigorously critique the system and set QA tests appropriate to the particular device and its possible weaknesses. Major recommendations of the Task Group include: introduction of a new nomenclature for reporting repositioning accuracy; comprehensive analysis of patient characteristics that might adversely affect positioning accuracy; performance of testing immediately before each treatment to establish that there are no gross positioning errors; a general request to the Medical Physics community for improved QA tools; implementation of weekly portal imaging (perhaps cone beam CT in the future) as a method of tracking fractionated patients (as per TG 40); and periodic routine reviews of positioning accuracy.