How does CT image noise affect 3D deformable image registration for image-guided radiotherapy planning?

An automated technique for global noise level measurement in CT image with a conjunction of image gradient

Automated assessment of noise level in clinical computed tomography (CT) images is a crucial technique for evaluating and ensuring the quality of these images. There are various factors that can impact CT image noise, such as statistical noise, electronic noise, structure noise, texture noise, artifact noise, etc. In this study, a method was developed to measure the global noise index (GNI) in clinical CT scans due to the fluctuation of x-ray quanta. Initially, a noise map is generated by sliding a 10 × 10 pixel for calculating Hounsfield unit (HU) standard deviation and the noise map is further combined with the gradient magnitude map. By employing Boolean operation, pixels with high gradients are excluded from the noise histogram generated with the noise map. By comparing the shape of the noise histogram from this method with Christianson's tissue-type global noise measurement algorithm, it was observed that the noise histogram computed in anthropomorphic phantoms had a similar shape with a close GNI value. In patient CT images, excluding the HU deviation due the structure change demonstrated to have consistent GNI values across the entire CT scan range with high heterogeneous tissue compared to the GNI values using Christianson's tissue-type method. The proposed GNI was evaluated in phantom scans and was found to be capable of comparing scan protocols between different scanners. The variation of GNI when using different reconstruction kernels in clinical CT images demonstrated a similar relationship between noise level and kernel sharpness as observed in uniform phantom: sharper kernel resulted in noisier images. This indicated that GNI was a suitable index for estimating the noise level in clinical CT images with either a smooth or grainy appearance. The study's results suggested that the algorithm can be effectively utilized to screen the noise level for a better CT image quality control.

Validation of deep learning-based CT image reconstruction for treatment planning

Deep learning-based CT image reconstruction (DLR) is a state-of-the-art method for obtaining CT images. This study aimed to evaluate the usefulness of DLR in radiotherapy. Data were acquired using a large-bore CT system and an electron density phantom for radiotherapy. We compared the CT values, image noise, and CT value-to-electron density conversion table of DLR and hybrid iterative reconstruction (H-IR) for various doses. Further, we evaluated three DLR reconstruction strength patterns (Mild, Standard, and Strong). The variations of CT values of DLR and H-IR were large at low doses, and the difference in average CT values was insignificant with less than 10 HU at doses of 100 mAs and above. DLR showed less change in CT values and smaller image noise relative to H-IR. The noise-reduction effect was particularly large in the low-dose region. The difference in image noise between DLR Mild and Standard/Strong was large, suggesting the usefulness of reconstruction intensities higher than Mild. DLR showed stable CT values and low image noise for various materials, even at low doses; particularly for Standard or Strong, the reduction in image noise was significant. These findings indicate the usefulness of DLR in treatment planning using large-bore CT systems.

Accurate and efficient internal deformation measurement of multiphase/porous materials via segmentation-aided digital volume correlation

In using the regular subvolume-based digital volume correlation (R-DVC) method, calculation points should be defined at the real material phase, and the local deformation within the interrogated subvolumes is assumed to be continuous. However, this basic assumption in R-DVC analysis is often violated when measuring the deformation near the interface when dealing with multiphase materials (including porous materials) or contact problems. This is because discontinuous deformation always presents in the calculation points located at the vicinity of interfaces of various phases. All these factors lead to increased measurement error and/or meaningless calculation burdens when using R-DVC. To address these issues, we propose a segmentation-aided DVC (S-DVC) for accuracy-enhanced internal deformation analysis near the interface. The presented S-DVC first divides the reference volume image into different portions according to the distinct gray scales within different material phases (or background) or objects. Based on the segmented reference volume image, we can ensure that subvolumes only contain the voxels from the same material phase/object and exclude other phases/objects. As such, the error due to undermatched shape function can be minimized and meaningless DVC calculation can be avoided. The accuracy, efficiency, and practicality of S-DVC over R-DVC are validated by a simulated compression test of nodular cast iron (multiphase material) and a real compression experiment of 3D printed polymer (porous material).

Clinical Assessment of a Novel Ring Gantry Linear Accelerator-Mounted Helical Fan-Beam kVCT System

Purpose

To assess clinically relevant image quality metrics (IQMs) of helical fan beam kilovoltage (kV) fan beam computed tomography (CT).

Methods and Materials

kVCT IQMs were evaluated on an Accuray Radixact unit equipped with helical fan beam kVCT to assess the capabilities of this newly available modality. kVCT IQMs were evaluated and compared to a kVCT simulator and linear accelerator-based cone beam CTs (CBCT) using a commercial CBCT image quality phantom. kVCTs were acquired on the Accuray Radixact for all combinations of kVp and mAs in fine mode using a 440-mm field of view (FOV). Evaluated IQMs were spatial resolution, overall uniformity, subject contrast, contrast-to-noise ratio (CNR), and effective slice thickness. Imaging dose was assessed for planar kV imaging.

Results

On this kVCT system spatial resolution and contrast were consistent across all settings with 0.28 ± 0.03 lp/mm and 9.8% ± 0.7% (both 95% confidence interval). CNR strongly depended on selected mode (views per rotation) and body size (mA per view) and ranged between 7.9 and 34.9. Overall uniformity was greater than 97% for all settings. Large FOV was not found to substantially affect the IQMs whereas small FOV affected IQMs due to its effect on pitch. Technique-matched CT simulator scans were comparable for uniformity and contrast, while spatial resolution was higher (0.43 ± 0.06 lp/mm), and CNR was between 4% (140 kVp) and 51% (100 kVp) lower. For kV-CBCT, spatial resolutions ranging from 0.37 to 0.44 lp/mm were achieved with comparable contrast, CNR, and uniformity to kVCT. All kVCT scans exhibit imaging artifacts due to helical acquisition. Clinical acquisitions of megavoltage (MV) CT, kV-CBCT, and kVCT on the same patient showed improved and comparable image quality of kVCT compared to MVCT and kV-CBCT, respectively.

Conclusions

Helical fan beam kVCT allows for daily image guidance for localization and setup verification with comparable performance to existing kV-CBCT systems. Scan parameters must be selected carefully to maximize image quality for the desired tasks. Due to the large effective slice thicknesses for all parameter combinations, kVCT scans should not be used for simulation or planning of stereotactic procedures. Finally, improved image quality over MVCT has the potential to greatly improve manual and automated adaptive monitoring and planning.

Past, present and future role of retinal imaging in neurodegenerative disease

Retinal imaging technology is rapidly advancing and can provide ever-increasing amounts of information about the structure, function and molecular composition of retinal tissue in humans in vivo. Most importantly, this information can be obtained rapidly, non-invasively and in many cases using Food and Drug Administration-approved devices that are commercially available. Technologies such as optical coherence tomography have dramatically changed our understanding of retinal disease and in many cases have significantly improved their clinical management. Since the retina is an extension of the brain and shares a common embryological origin with the central nervous system, there has also been intense interest in leveraging the expanding armamentarium of retinal imaging technology to understand, diagnose and monitor neurological diseases. This is particularly appealing because of the high spatial resolution, relatively low-cost and wide availability of retinal imaging modalities such as fundus photography or OCT compared to brain imaging modalities such as magnetic resonance imaging or positron emission tomography. The purpose of this article is to review and synthesize current research about retinal imaging in neurodegenerative disease by providing examples from the literature and elaborating on limitations, challenges and future directions. We begin by providing a general background of the most relevant retinal imaging modalities to ensure that the reader has a foundation on which to understand the clinical studies that are subsequently discussed. We then review the application and results of retinal imaging methodologies to several prevalent neurodegenerative diseases where extensive work has been done including sporadic late onset Alzheimer's Disease, Parkinson's Disease and Huntington's Disease. We also discuss Autosomal Dominant Alzheimer's Disease and cerebrovascular small vessel disease, where the application of retinal imaging holds promise but data is currently scarce. Although cerebrovascular disease is not generally considered a neurodegenerative process, it is both a confounder and contributor to neurodegenerative disease processes that requires more attention. Finally, we discuss ongoing efforts to overcome the limitations in the field and unmet clinical and scientific needs.

Validation of accuracy deformable image registration contour propagation using a benchmark virtual HN phantom dataset

Virtual anatomic phantoms offer precise voxel mapping of the variation of anatomy with ground truth deformation vector fields (DVFs). Dice similarity coefficient (DSC) and mean distance to agreement (MDA) are the standard metrics for evaluating geometric contour congruence when testing deformable registration (DIR) algorithms. A HN virtual patient phantom data set was used for a kVCT‐kVCT automatic propagation contour validation study employing the Accuray DIR algorithm. Furthermore, since TomoTherapy uses MVCT images of the relevant anatomy for adaptive monitoring, the kVCT image data set quality was transformed to an MVCT image data set quality to study intermodal kVCT‐MVCT DIR accuracy. The results of the study indicate that the Accuray DIR algorithm can be expected to autopropagate HN contours adequately, on average, within tolerances recommended by TG‐132 (DSC 0.8‐0.9, MDA within voxel width). However, contours critical to dosimetric planning should always be visually proofed for accuracy. Using standard reconstruction MVCT image quality causes slightly less, but acceptable, agreement with ground truth contours.

An international survey on the clinical use of rigid and deformable image registration in radiotherapy

OBJECTIVES

Rigid image registration (RIR) and deformable image registration (DIR) are widely used in radiotherapy. This project aims to capture current international approaches to image registration.

METHODS

A survey was designed to identify variations in use, resources, implementation, and decision-making criteria for clinical image registration. This was distributed to radiotherapy centers internationally in 2018.

RESULTS

There were 57 responses internationally, from the Americas (46%), Australia/New Zealand (32%), Europe (12%), and Asia (10%). Rigid image registration and DIR were used clinically for computed tomography (CT)-CT registration (96% and 51%, respectively), followed by CT-PET (81% and 47%), CT-CBCT (84% and 19%), CT-MR (93% and 19%), MR-MR (49% and 5%), and CT-US (9% and 0%). Respondent centers performed DIR using dedicated software (75%) and treatment planning systems (29%), with 84% having some form of DIR software. Centers have clinically implemented DIR for atlas-based segmentation (47%), multi-modality treatment planning (65%), and dose deformation (63%). The clinical use of DIR for multi-modality treatment planning and accounting for retreatments was considered to have the highest benefit-to-risk ratio (69% and 67%, respectively).

CONCLUSIONS

This survey data provides useful insights on where, when, and how image registration has been implemented in radiotherapy centers around the world. DIR is mainly in clinical use for CT-CT (51%) and CT-PET (47%) for the head and neck (43-57% over all use cases) region. The highest benefit-risk ratio for clinical use of DIR was for multi-modality treatment planning and accounting for retreatments, which also had higher clinical use than for adaptive radiotherapy and atlas-based segmentation.

Voxel-based analysis in radiation oncology: A methodological cookbook

Recently, 2D or 3D methods for dose distribution analysis have been proposed as evolutions of the Dose Volume Histogram (DVH) approaches. Those methods, collectively referred to as pixel- or voxel-based (VB) methods, evaluate local dose response patterns and go beyond the organ-based philosophy of Normal Tissue Complication Probability (NTCP) modelling. VB methods have been introduced in the context of radiation oncology in the very last years following the virtuous example of neuroimaging experience. In radiation oncology setting, dose mapping is a suitable scheme to compare spatial patterns of local dose distributions between patients who develop toxicity and who do not. In this critical review, we present the methods that include spatial dose distribution information for evaluating different toxicity endpoints after radiation therapy. The review addresses two main topics. First, the critical aspects in dose map building, namely the spatial normalization of the dose distributions from different patients. Then, the issues related to the actual dose map comparison, i.e. the viable options for a robust VB statistical analysis and the potential pitfalls related to the adopted solutions. To elucidate the different theoretical and technical issues, the covered topics are illustrated in relation to practical applications found in the existing literature. We conclude the overview on the VB philosophy in radiation oncology by introducing new phenomenological approaches to NTCP modelling that accounts for inhomogeneous organ radiosensitivity.

Computed Tomography to Cone Beam Computed Tomography Deformable Image Registration for Contour Propagation Using Head and Neck, Patient-Based Computational Phantoms: A Multicenter Study

PURPOSE

To investigate the performance of various algorithms for deformable image registration (DIR) for propagating regions of interest (ROIs) using multiple commercial platforms, from computed tomography to cone beam computed tomography (CBCT) and megavoltage computed tomography.

METHODS AND MATERIALS

Fourteen institutions participated in the study using 5 commercial platforms: RayStation (RaySearch Laboratories, Stockholm, Sweden), MIM (Cleveland, OH), VelocityAI and SmartAdapt (Varian Medical Systems, Palo Alto, CA), and ABAS (Elekta AB, Stockholm, Sweden). Algorithms were tested on synthetic images generated with the ImSimQA (Oncology Systems Limited, Shrewsbury, UK) package by applying 2 specific deformation vector fields (DVF) to real head and neck patient datasets. On-board images from 3 systems were used: megavoltage computed tomography from Tomotherapy and 2 kinds of CBCT from a clinical linear accelerator. Image quality of the system was evaluated. The algorithms' accuracy was assessed by comparing the DIR-mapped ROIs returned by each center with those of the reference, using the Dice similarity coefficient and mean distance to conformity metrics. Statistical inference on the validation results was carried out to identify the prognostic factors of DIR performance.

RESULTS

Analyzing 840 DIR-mapped ROIs returned by the centers, it was demonstrated that DVF intensity and image quality were significant prognostic factors of DIR performance. The accuracy of the propagated contours was generally high, and acceptable DIR performance can be obtained with lower-dose CBCT image protocols.

CONCLUSIONS

The performance of the systems proved to be image quality specific, depending on the DVF type and only partially on the platforms. All systems proved to be robust against image artifacts and noise, except the demon-based software.

Usefulness of hybrid deformable image registration algorithms in prostate radiation therapy

To evaluate the accuracy of commercially available hybrid deformable image registration (DIR) algorithms when using planning CT (pCT) and daily cone-beam computed tomography (CBCT) in radiation therapy for prostate cancer. The hybrid DIR algorithms in RayStation and MIM Maestro were evaluated. Contours of the prostate, bladder, rectum, and seminal vesicles (SVs) were used as region-of-interest (ROIs) to guide image deformation in the hybrid DIR and to compare the DIR accuracy. To evaluate robustness of the hybrid DIR for prostate cancer patients with organs with volume that vary on a daily basis, such as the bladder and rectum, the DIR algorithms were performed on ten pairs of CT volumes from ten patients who underwent prostate intensity-modulated radiation therapy or volumetric modulated arc therapy. In a visual evaluation, MIM caused unrealistic image deformation in soft tissues, organs, and pelvic bones. The mean dice similarity coefficient (DSC) ranged from 0.46 to 0.90 for the prostate, bladder, rectum, and SVs; the SVs had the lowest DSC. Target registration error (TRE) at the centroid of the ROIs was about 2 mm for the prostate and bladder, and about 6 mm for the rectum and SVs. RayStation did not cause unrealistic image deformation, and could maintain the shape of pelvic bones in most cases. The mean DSC and TRE at the centroid of the ROIs were about 0.9 and within 5 mm generally. In both software programs, the use of ROIs to guide image deformation had the possibility to reduce any unrealistic image deformation and might be effective to keep the DIR physically reasonable. The pCT/CBCT DIR for the prostate cancer did not reduce the DIR accuracy because of the use of ROIs to guide the image deformation.

Evaluation of automatic contour propagation in T2-weighted 4DMRI for normal-tissue motion assessment using internal organ-at-risk volume (IRV)

Purpose

The purpose of this study was to evaluate the quality of automatically propagated contours of organs at risk (OARs) based on respiratory‐correlated navigator‐triggered four‐dimensional magnetic resonance imaging (RC‐4DMRI) for calculation of internal organ‐at‐risk volume (IRV) to account for intra‐fractional OAR motion.

Methods and Materials

T2‐weighted RC‐4DMRI images were of 10 volunteers acquired and reconstructed using an internal navigator‐echo surrogate and concurrent external bellows under an IRB‐approved protocol. Four major OARs (lungs, heart, liver, and stomach) were delineated in the 10‐phase 4DMRI. Two manual‐contour sets were delineated by two clinical personnel and two automatic‐contour sets were propagated using free‐form deformable image registration. The OAR volume variation within the 10‐phase cycle was assessed and the IRV was calculated as the union of all OAR contours. The OAR contour similarity between the navigator‐triggered and bellows‐rebinned 4DMRI was compared. A total of 2400 contours were compared to the most probable ground truth with a 95% confidence level (S95) in similarity, sensitivity, and specificity using the simultaneous truth and performance level estimation (STAPLE) algorithm.

Results

Visual inspection of automatically propagated contours finds that approximately 5–10% require manual correction. The similarity, sensitivity, and specificity between manual and automatic contours are indistinguishable (P > 0.05). The Jaccard similarity indexes are 0.92 ± 0.02 (lungs), 0.89 ± 0.03 (heart), 0.92 ± 0.02 (liver), and 0.83 ± 0.04 (stomach). Volume variations within the breathing cycle are small for the heart (2.6 ± 1.5%), liver (1.2 ± 0.6%), and stomach (2.6 ± 0.8%), whereas the IRV is much larger than the OAR volume by: 20.3 ± 8.6% (heart), 24.0 ± 8.6% (liver), and 47.6 ± 20.2% (stomach). The Jaccard index is higher in navigator‐triggered than bellows‐rebinned 4DMRI by 4% (P < 0.05), due to the higher image quality of navigator‐based 4DMRI.

Conclusion

Automatic and manual OAR contours from Navigator‐triggered 4DMRI are not statistically distinguishable. The navigator‐triggered 4DMRI image provides higher contour quality than bellows‐rebinned 4DMRI. The IRVs are 20–50% larger than OAR volumes and should be considered in dose estimation.

Feasibility study of range-based registration using daily cone beam CT for intensity-modulated proton therapy

PURPOSE

Proton dose coverage is sensitive to proton beam range. The current practice of CT number-based registration for patient positioning focuses on matching the target and is not sufficient for proton therapy because the proton range depends on the medium traversed by the beam. Patient body deformations and anatomical changes result in range deviation in the target. We propose proton range-based registration to minimize the range deviation.

METHODS

The range was calculated from cone beam-computed tomography (CBCT) of the patient on couch, and the range deviation was the difference of the calculated range from that on the initial (day 1) CBCT. In the investigated prostate cases in which the main cause of range deviation was the rotation of femur bones, and in the investigated abdomen cases in which the main cause of range deviation was body growth and anatomic change, our range-based registration was used to obtain the optimal beam angle by minimizing the range deviation. The new angle was limited to be ±5° from that planned to prevent potentially increased dose to the organs at risk. To demonstrate the benefit of range-based registration, we investigated the range at the voxels on the surface of the target volume. The calculation error of range deviation due to CBCT scatter was investigated by using solid water phantoms with different thicknesses. Range-based registration using both CBCTs and CTs was performed in cases of two patients with pelvic rhabdomyosarcoma and one patient with upper abdominal tumor. The range was represented by the water-equivalent thickness to shorten the computation for online application purposes.

RESULTS

In the phantom study, the calculation error of range deviation due to CBCT scatter was within 2 mm for a 1-cm thickness change (the mean range deviation was 0.8 mm). In the CT study of the prostate cases, the range deviation (mean ± root-mean-square deviation) on the contour in each slice was efficiently reduced from 3.6 ± 2.8 mm to 2.1 ± 1.4 mm, with most slices being within 3 mm; in the CT study of the abdomen cases, the range deviation of the whole set was reduced from 4.4 ± 1.9 mm to 3.5 ± 2.1 mm. Both the mean and root-mean-square deviation of the range deviation on each treatment day were decreased. The dose coverage on the target was improved and the dose on the OARs was only slightly changed.

CONCLUSION

Range-based registration can efficiently mitigate range deviation due to patient positioning and anatomical changes. It can shorten patient positioning time and reduce the patient's dose from CBCT.

Optimization of an on-board imaging system for extremely rapid radiation therapy

PURPOSE

Next-generation extremely rapid radiation therapy systems could mitigate the need for motion management, improve patient comfort during the treatment, and increase patient throughput for cost effectiveness. Such systems require an on-board imaging system that is competitively priced, fast, and of sufficiently high quality to allow good registration between the image taken on the day of treatment and the image taken the day of treatment planning. In this study, three different detectors for a custom on-board CT system were investigated to select the best design for integration with an extremely rapid radiation therapy system.

METHODS

Three different CT detectors are proposed: low-resolution (all 4×4 mm pixels), medium-resolution (a combination of 4×4 mm pixels and 2×2 mm pixels), and high-resolution (all 1×1 mm pixels). An in-house program was used to generate projection images of a numerical anthropomorphic phantom and to reconstruct the projections into CT datasets, henceforth called "realistic" images. Scatter was calculated using a separate Monte Carlo simulation, and the model included an antiscatter grid and bowtie filter. Diagnostic-quality images of the phantom were generated to represent the patient scan at the time of treatment planning. Commercial deformable registration software was used to register the diagnostic-quality scan to images produced by the various on-board detector configurations. The deformation fields were compared against a "gold standard" deformation field generated by registering initial and deformed images of the numerical phantoms that were used to make the diagnostic and treatment-day images. Registrations of on-board imaging system data were judged by the amount their deformation fields differed from the corresponding gold standard deformation fields--the smaller the difference, the better the system. To evaluate the registrations, the pointwise distance between gold standard and realistic registration deformation fields was computed.

RESULTS

By most global metrics (e.g., mean, median, and maximum pointwise distance), the high-resolution detector had the best performance but the medium-resolution detector was comparable. For all medium- and high-resolution detector registrations, mean error between the realistic and gold standard deformation fields was less than 4 mm. By pointwise metrics (e.g., tracking a small lesion), the high- and medium-resolution detectors performed similarly. For these detectors, the smallest error between the realistic and gold standard registrations was 0.6 mm and the largest error was 3.6 mm.

CONCLUSIONS

The medium-resolution CT detector was selected as the best for an extremely rapid radiation therapy system. In essentially all test cases, data from this detector produced a significantly better registration than data from the low-resolution detector and a comparable registration to data from the high-resolution detector. The medium-resolution detector provides an appropriate compromise between registration accuracy and system cost.

Evaluation of the tool "Reg Refine" for user-guided deformable image registration

“Reg Refine” is a tool available in the MIM Maestro v6.4.5 platform (www.mimsoftware.com) that allows the user to actively participate in the deformable image registration process. The purpose of this work was to evaluate the efficacy of this tool and investigate strategies for how to apply it effectively. This was done by performing DIR on two publicly available ground‐truth models, the Pixel‐based Breathing Thorax Model (POPI) for lung, and the Deformable Image Registration Evaluation Project (DIREP) for head and neck. Image noise matched in both magnitude and texture to clinical CBCT scans was also added to each model to simulate the use case of CBCT–CT alignment. For lung, the results showed Reg Refine effective at improving registration accuracy when controlled by an expert user within the context of large lung deformation. CBCT noise was also shown to have no effect on DIR performance while using the MIM algorithm for this site. For head and neck, the results showed CBCT noise to have a large effect on the accuracy of registration, specifically for low‐contrast structures such as the brainstem and parotid glands. In these cases, the Reg Refine tool was able to improve the registration accuracy when controlled by an expert user. Several strategies for how to achieve these results have been outlined to assist other users and provide feedback for developers of similar tools.

PACS number(s): 87.44.Qr, 87.57.nj, 87.57.c

Cone Beam Computed Tomography: The Challenges and Strategies in Its Application for Dose Accumulation

Online image guidance using cone beam computed tomography (CBCT) has greatly improved the geometric precision of radiotherapy. Changes in anatomy are common during a course of fractionated treatment, resulting in dose deviation from the planned distribution. There is increased interest in performing dose accumulation to compute the actual delivered dose and to adapt the treatment when necessary. This can be achieved by delineating the volume of interest and by generating "dose of the day" through dose computation on the CBCT. However, the image quality and the accuracy of the CT number of CBCT are deemed to be inferior to fan beam CT, which increases the uncertainty associated in this process. A review of literature was conducted to assess the reliability of and to examine strategies for overcoming the challenges in using CBCT for volume delineation and dose computation. The review demonstrates that the uncertainty varies across body sites, and different strategies have been recommended to generate comparable results to images from CT simulators. This facilitates a better understanding of the potential and the limitation of using CBCT for dose accumulation.

The effects of computed tomography image characteristics and knot spacing on the spatial accuracy of B-spline deformable image registration in the head and neck geometry

OBJECTIVES

To explore the effects of computed tomography (CT) image characteristics and B-spline knot spacing (BKS) on the spatial accuracy of a B-spline deformable image registration (DIR) in the head-and-neck geometry.

METHODS

The effect of image feature content, image contrast, noise, and BKS on the spatial accuracy of a B-spline DIR was studied. Phantom images were created with varying feature content and varying contrast-to-noise ratio (CNR), and deformed using a known smooth B-spline deformation. Subsequently, the deformed images were repeatedly registered with the original images using different BKSs. The quality of the DIR was expressed as the mean residual displacement (MRD) between the known imposed deformation and the result of the B-spline DIR.Finally, for three patients, head-and-neck planning CT scans were deformed with a realistic deformation field derived from a rescan CT of the same patient, resulting in a simulated deformed image and an a-priori known deformation field. Hence, a B-spline DIR was performed between the simulated image and the planning CT at different BKSs. Similar to the phantom cases, the DIR accuracy was evaluated by means of MRD.

RESULTS

In total, 162 phantom registrations were performed with varying CNR and BKSs. MRD-values < 1.0 mm were observed with a BKS between 10-20 mm for image contrast ≥ ± 250 HU and noise < ± 200 HU. Decreasing the image feature content resulted in increased MRD-values at all BKSs. Using BKS = 15 mm for the three clinical cases resulted in an average MRD < 1.0 mm.

CONCLUSIONS

For synthetically generated phantoms and three real CT cases the highest DIR accuracy was obtained for a BKS between 10-20 mm. The accuracy decreased with decreasing image feature content, decreasing image contrast, and higher noise levels. Our results indicate that DIR accuracy in clinical CT images (typical noise levels < ± 100 HU) will not be effected by the amount of image noise.

Validation of Automatic Contour Propagation for 4D Treatment Planning Using Multiple Metrics

The aim of this work is to provide insights into multiple metrics clinical validation of deformable image registration and contour propagation methods in 4D lung radiotherapy planning. The following indices were analyzed and compared: Volume Difference (VD), Dice Similarity Coefficient (DSC), Positive Predictive Value (PPV) and Surface Distances (SD). The analysis was performed on three patient datasets, using as reference a ground-truth volume generated by means of Simultaneous Truth And Performance Level Estimation (STAPLE) algorithm from the outlines of five experts. Significant discrepancies in the quality assessment provided by the different metrics in all the examined cases were found. Metrics sensitivity was more evident in presence of image artifacts and particularly for tubular anatomical structures, such as esophagus or spinal cord. Volume Differences did not account for position and DSC exhibited criticalities due to its intrinsic symmetry ( i.e. over- and under-estimation of the reference contours cannot be discriminated) and dependency on the total volume of the structure. PPV analysis showed more robust performance, as each voxel concurs to the classification of the propagation, but was not able to detect inclusion of propagated and ground-truth volumes. Mesh distances could interpret the actual shape of the structures, but might report higher mismatches in case of large local differences in the contour surfaces. According to our study, the combination of VD and SD for the validation of contour propagation algorithms in 4D could provide the necessary failure detection accuracy.

Using patient-specific phantoms to evaluate deformable image registration algorithms for adaptive radiation therapy

The quality of adaptive treatment planning depends on the accuracy of its underlying deformable image registration (DIR). The purpose of this study is to evaluate the performance of two DIR algorithms, B‐spline‐based deformable multipass (DMP) and deformable demons (Demons), implemented in a commercial software package. Evaluations were conducted using both computational and physical deformable phantoms. Based on a finite element method (FEM), a total of 11 computational models were developed from a set of CT images acquired from four lung and one prostate cancer patients. FEM generated displacement vector fields (DVF) were used to construct the lung and prostate image phantoms. Based on a fast‐Fourier transform technique, image noise power spectrum was incorporated into the prostate image phantoms to create simulated CBCT images. The FEM‐DVF served as a gold standard for verification of the two registration algorithms performed on these phantoms. The registration algorithms were also evaluated at the homologous points quantified in the CT images of a physical lung phantom. The results indicated that the mean errors of the DMP algorithm were in the range of 1.0~3.1mm for the computational phantoms and 1.9 mm for the physical lung phantom. For the computational prostate phantoms, the corresponding mean error was 1.0–1.9 mm in the prostate, 1.9–2.4 mm in the rectum, and 1.8–2.1 mm over the entire patient body. Sinusoidal errors induced by B‐spline interpolations were observed in all the displacement profiles of the DMP registrations. Regions of large displacements were observed to have more registration errors. Patient‐specific FEM models have been developed to evaluate the DIR algorithms implemented in the commercial software package. It has been found that the accuracy of these algorithms is patient‐dependent and related to various factors including tissue deformation magnitudes and image intensity gradients across the regions of interest. This may suggest that DIR algorithms need to be verified for each registration instance when implementing adaptive radiation therapy.

PACS numbers: 87.10.Kn, 87.55.km, 87.55.Qr, 87.57.nj

Effects of quantum noise in 4D-CT on deformable image registration and derived ventilation data

Quantum noise is common in CT images and is a persistent problem in accurate ventilation imaging using 4D-CT and deformable image registration (DIR). This study focuses on the effects of noise in 4D-CT on DIR and thereby derived ventilation data. A total of six sets of 4D-CT data with landmarks delineated in different phases, called point-validated pixel-based breathing thorax models (POPI), were used in this study. The DIR algorithms, including diffeomorphic morphons (DM), diffeomorphic demons (DD), optical flow and B-spline, were used to register the inspiration phase to the expiration phase. The DIR deformation matrices (DIRDM) were used to map the landmarks. Target registration errors (TRE) were calculated as the distance errors between the delineated and the mapped landmarks. Noise of Gaussian distribution with different standard deviations (SD), from 0 to 200 Hounsfield Units (HU) in amplitude, was added to the POPI models to simulate different levels of quantum noise. Ventilation data were calculated using the ΔV algorithm which calculates the volume change geometrically based on the DIRDM. The ventilation images with different added noise levels were compared using Dice similarity coefficient (DSC). The root mean square (RMS) values of the landmark TRE over the six POPI models for the four DIR algorithms were stable when the noise level was low (SD <150 HU) and increased with added noise when the level is higher. The most accurate DIR was DD with a mean RMS of 1.5 ± 0.5 mm with no added noise and 1.8 ± 0.5 mm with noise (SD = 200 HU). The DSC values between the ventilation images with and without added noise decreased with the noise level, even when the noise level was relatively low. The DIR algorithm most robust with respect to noise was DM, with mean DSC = 0.89 ± 0.01 and 0.66 ± 0.02 for the top 50% ventilation volumes, as compared between 0 added noise and SD = 30 and 200 HU, respectively. Although the landmark TRE were stable with low noise, the differences between ventilation images increased with noise level, even when the noise was low, indicating ventilation imaging from 4D-CT was sensitive to image noise. Therefore, high quality 4D-CT is essential for accurate ventilation images.

A study of the anatomic changes and dosimetric consequences in adaptive CRT of non-small-cell lung cancer using deformable CT and CBCT image registration

The aim of this study is to evaluate anatomic lung tumor changes and dosimetric consequences utilizing the deformable daily kilovolt (KV) cone-beam computer tomography (CBCT) image registration. Five patients diagnosed with NSCLC were treated with three-dimensional conformal radiotherapy (3D CRT) and 10 daily KV CBCT image sets were acquired for each patient. Each CBCT image and plan CT were imported into the deformable image registration (DIR) system. The plan CT image was deformed by the DIR system and a new contour on CBCT was obtained by using the auto-contouring function of the DIR. These contours were individually marked as CBCT f1, CBCT f2,..., and CBCT f10, and imported into a treatment planning system (TPS). The daily CBCT plan was individually generated with the same planning criteria based on new contours. These plans were individually marked as CBCTp1, CBCTp2,..., and CBCTp10, followed by generating a dose accumulation plan (DA plan) in original pCT image contour sets by adding all CBCT plans using Varian Eclipse TPS. The maximum, minimum and mean doses to the plan target volume (PTV) in the 5 DA plans were the same with the CT plans. However, the volume of radiation 5, 10, 20, 30, and 50 Gy of the total lungs in DA plans were less than those of the CT plans. The maximum dose of the spinal cord in the DA plans were average 27.96% less than the CT plans. The mean dose for the left, right, and total lungs in the DA plans were reduced by 13.80%, 23.65%, and 12.96%, respectively. The adaptive 3D CRT based on the deformable registration can reduce the dose to the lung and the spinal cord with the same PTV dose coverage. Moreover, it provides a method for further adaptive radiotherapy exploration.

A fully automated method for CT-on-rails-guided online adaptive planning for prostate cancer intensity modulated radiation therapy

PURPOSE

This study was designed to validate a fully automated adaptive planning (AAP) method which integrates automated recontouring and automated replanning to account for interfractional anatomical changes in prostate cancer patients receiving adaptive intensity modulated radiation therapy (IMRT) based on daily repeated computed tomography (CT)-on-rails images.

METHODS AND MATERIALS

Nine prostate cancer patients treated at our institution were randomly selected. For the AAP method, contours on each repeat CT image were automatically generated by mapping the contours from the simulation CT image using deformable image registration. An in-house automated planning tool incorporated into the Pinnacle treatment planning system was used to generate the original and the adapted IMRT plans. The cumulative dose-volume histograms (DVHs) of the target and critical structures were calculated based on the manual contours for all plans and compared with those of plans generated by the conventional method, that is, shifting the isocenters by aligning the images based on the center of the volume (COV) of prostate (prostate COV-aligned).

RESULTS

The target coverage from our AAP method for every patient was acceptable, while 1 of the 9 patients showed target underdosing from prostate COV-aligned plans. The normalized volume receiving at least 70 Gy (V70), and the mean dose of the rectum and bladder were reduced by 8.9%, 6.4 Gy and 4.3%, 5.3 Gy, respectively, for the AAP method compared with the values obtained from prostate COV-aligned plans.

CONCLUSIONS

The AAP method, which is fully automated, is effective for online replanning to compensate for target dose deficits and critical organ overdosing caused by interfractional anatomical changes in prostate cancer.

A protocol to extend the longitudinal coverage of on-board cone-beam CT

The longitudinal coverage of a LINAC‐mounted CBCT scan is limited to the corresponding dimensional limits of its flat panel detector, which is often shorter than the length of the treatment field. These limits become apparent when fields are designed to encompass wide regions, as when providing nodal coverage. Therefore, we developed a novel protocol to acquire double orbit CBCT images using a commercial system, and combine the images to extend the longitudinal coverage for image‐guided adaptive radiotherapy (IGART). The protocol acquires two CBCT scans with a couch shift similar to the “step‐and‐shoot” cine CT acquisition, allowing a small longitudinal overlap of the two reconstructed volumes. An in‐house DICOM reading/writing software was developed to combine the two image sets into one. Three different approaches were explored to handle the possible misalignment between the two image subsets: simple stacking, averaging the overlapped volumes, and a 3D‐3D image registration with the three translational degrees of freedom. Using thermoluminescent dosimeters and custom‐designed holders for a CTDI phantom set, dose measurements were carried out to assess the resultant imaging dose of the technique and its geometric distribution. Deformable registration was tested on patient images generated with the double‐orbit protocol, using both the planning FBCT and the artificially deformed CBCT as source images. The protocol was validated on phantoms and has been employed clinically for IRB‐approved IGART studies for head and neck and prostate cancer patients.

PACS number: 87.57.nj

A method to estimate the effect of deformable image registration uncertainties on daily dose mapping

PURPOSE

To develop a statistical sampling procedure for spatially-correlated uncertainties in deformable image registration and then use it to demonstrate their effect on daily dose mapping.

METHODS

Sequential daily CT studies are acquired to map anatomical variations prior to fractionated external beam radiotherapy. The CTs are deformably registered to the planning CT to obtain displacement vector fields (DVFs). The DVFs are used to accumulate the dose delivered each day onto the planning CT. Each DVF has spatially-correlated uncertainties associated with it. Principal components analysis (PCA) is applied to measured DVF error maps to produce decorrelated principal component modes of the errors. The modes are sampled independently and reconstructed to produce synthetic registration error maps. The synthetic error maps are convolved with dose mapped via deformable registration to model the resulting uncertainty in the dose mapping. The results are compared to the dose mapping uncertainty that would result from uncorrelated DVF errors that vary randomly from voxel to voxel.

RESULTS

The error sampling method is shown to produce synthetic DVF error maps that are statistically indistinguishable from the observed error maps. Spatially-correlated DVF uncertainties modeled by our procedure produce patterns of dose mapping error that are different from that due to randomly distributed uncertainties.

CONCLUSIONS

Deformable image registration uncertainties have complex spatial distributions. The authors have developed and tested a method to decorrelate the spatial uncertainties and make statistical samples of highly correlated error maps. The sample error maps can be used to investigate the effect of DVF uncertainties on daily dose mapping via deformable image registration. An initial demonstration of this methodology shows that dose mapping uncertainties can be sensitive to spatial patterns in the DVF uncertainties.

Deformable registration using edge-preserving scale space for adaptive image-guided radiation therapy

Incorporating of daily cone-beam computer tomography (CBCT) image into online radiation therapy process can achieve adaptive image-guided radiation therapy (AIGRT). Registration of planning CT (PCT) and daily CBCT are the key issues in this process. In our work, a new multiscale deformable registration method is proposed by combining edge-preserving scale space with the multilevel free-form deformation (FFD) grids for CBCT-based AIGRT system. The edge-preserving scale space, which is able to select edges and contours of images according to their geometric size, is derived from the total variation model with the L1 norm (TV-L1). At each scale, despite the noise and contrast resolution differences between the PCT and CBCT, the selected edges and contours are sufficiently strong to drive the deformation using the FFD grid, and the edge-preserving property ensures more meaningful spatial information for mutual information (MI)-based registration. At last, the deformation fields are gained by a coarse to fine manner. Furthermore, in consideration of clinical application we designed an optimal estimation of the TV-L1 parameters by minimizing the defined offset function for automated registration. Six types of patients are studied in our work, including rectum, prostate, lung, H&N (head and neck), breast, and chest cancer patients. The experiment results demonstrate the significance of the proposed method both quantitatively with ground truth known and qualitatively with ground truth unknown. The applications for AIGRT, including adaptive deformable recontouring and redosing, and DVH (dose volume histogram) analysis in the course of radiation therapy are also studied.

Reconstruction of a cone-beam CT image via forward iterative projection matching

PURPOSE

To demonstrate the feasibility of reconstructing a cone-beam CT (CBCT) image by deformably altering a prior fan-beam CT (FBCT) image such that it matches the anatomy portrayed in the CBCT projection data set.

METHODS

A prior FBCT image of the patient is assumed to be available as a source image. A CBCT projection data set is obtained and used as a target image set. A parametrized deformation model is applied to the source FBCT image, digitally reconstructed radiographs (DRRs) that emulate the CBCT projection image geometry are calculated and compared to the target CBCT projection data, and the deformation model parameters are adjusted iteratively until the DRRs optimally match the CBCT projection data set. The resulting deformed FBCT image is hypothesized to be an accurate representation of the patient's anatomy imaged by the CBCT system. The process is demonstrated via numerical simulation. A known deformation is applied to a prior FBCT image and used to create a synthetic set of CBCT target projections. The iterative projection matching process is then applied to reconstruct the deformation represented in the synthetic target projections; the reconstructed deformation is then compared to the known deformation. The sensitivity of the process to the number of projections and the DRR/CBCT projection mismatch is explored by systematically adding noise to and perturbing the contrast of the target projections relative to the iterated source DRRs and by reducing the number of projections.

RESULTS

When there is no noise or contrast mismatch in the CBCT projection images, a set of 64 projections allows the known deformed CT image to be reconstructed to within a nRMS error of 1% and the known deformation to within a nRMS error of 7%. A CT image nRMS error of less than 4% is maintained at noise levels up to 3% of the mean projection intensity, at which the deformation error is 13%. At 1% noise level, the number of projections can be reduced to 8 while maintaining CT image and deformation errors of less than 4% and 13%, respectively. The method is sensitive to contrast mismatch between the simulated projections and the target projections when the soft-tissue contrast in the projections is low.

CONCLUSIONS

By using prior knowledge available in a FBCT image, the authors show that a CBCT image can be iteratively reconstructed from a comparatively small number of projection images, thus saving acquisition time and reducing imaging dose. This will enable more frequent daily imaging during radiation therapy. Because the process preserves the CT numbers of the FBCT image, the resulting 3D image intensities will be more accurate than a CBCT image reconstructed via conventional backprojection methods. Reconstruction errors are insensitive to noise at levels beyond what would typically be found in CBCT projection data, but are sensitive to contrast mismatch errors between the CBCT projection data and the DRRs.

Automatic image guidance for prostate IMRT using low dose CBCT

PURPOSE

Varian's On-Board Imager is a linac-integrated cone-beam CT (CBCT) system used at the authors' institution to acquire images prior to delivering each fraction of prostate intensity modulated radiotherapy. The images are used to determine a couch shift that realigns the tumor with the position obtained in the planning CT. However, this manual image-guided radiotherapy (IGRT) technique is operator dependent, time consuming, offers limited degrees of freedom, and requires significant imaging dose over the course of treatment. To overcome these problems, the authors propose two fully automatic IGRT techniques that require significantly less imaging dose.

METHODS

Dose is reduced by lowering the x-ray tube mA s during CBCT acquisition at the cost of increasing image noise. In "forward" IGRT, the CBCT image is automatically registered to the planning CT to obtain the necessary couch shift. The "reverse" technique offers additional degrees of freedom as it involves nonrigid registration of the planning CT to the CBCT. Both techniques were evaluated using images of an anthropomorphic phantom with simulated motion and by retrospectively analyzing data from ten prostate cancer patients.

RESULTS

IGRT error for the phantom data at 100% relative imaging dose was 8.2 +/- 3.7, 3.5 +/- 1.2,, and 2.1 +/- 0.6 mm for setup only, forward, and reverse techniques, respectively. For patient images acquired at 100% relative imaging dose, the errors were 5.4 +/- 1.7, 5.0 +/- 1.6, 5.0 +/- 2.0, and 4.0 +/- 1.6 mm for setup only, manual forward (performed clinically), automatic forward, and reverse IGRT, respectively. Furthermore, imaging dose could be reduced to 20% without a significant loss in image guidance accuracy.

CONCLUSIONS

The presented image guidance methods are accurate while requiring only 20% of the standard imaging dose. The combination of low dose, automation, and accuracy enables frequent corrections during treatment, possibly leading to reduced margins and improved treatment outcomes.

Assessing the intrinsic precision of 3D/3D rigid image registration results for patient setup in the absence of a ground truth

PURPOSE

To assess the precision and robustness of patient setup corrections computed from 3D/3D rigid registration methods using image intensity, when no ground truth validation is possible.

METHODS

Fifteen pairs of male pelvic CTs were rigidly registered using four different in-house registration methods. Registration results were compared for different resolutions and image content by varying the image down-sampling ratio and by thresholding out soft tissue to isolate bony landmarks. Intrinsic registration precision was investigated by comparing the different methods and by reversing the source and the target roles of the two images being registered.

RESULTS

The translational reversibility errors for successful registrations ranged from 0.0 to 1.69 mm. Rotations were less than 1 degrees. Mutual information failed in most registrations that used only bony landmarks. The magnitude of the reversibility error was strongly correlated with the success/ failure of each algorithm to find the global minimum.

CONCLUSIONS

Rigid image registrations have an intrinsic uncertainty and robustness that depends on the imaging modality, the registration algorithm, the image resolution, and the image content. In the absence of an absolute ground truth, the variation in the shifts calculated by several different methods provides a useful estimate of that uncertainty. The difference observed by reversing the source and target images can be used as an indication of robust convergence.

Optimal slice thickness for cone-beam CT with on-board imager

PURPOSE

To find the optimal slice thickness (Δτ) setting for patient registration with kilovoltage cone-beam CT (kVCBCT) on the Varian On Board Imager (OBI) system by investigating the relationship of slice thickness to automatic registration accuracy and contrast-to-noise ratio.

MATERIALS AND METHOD

Automatic registration was performed on kVCBCT studies of the head and pelvis of a RANDO anthropomorphic phantom. Images were reconstructed with 1.0 ≤ Δτ (mm) ≤ 5.0 at 1.0 mm increments. The phantoms were offset by a known amount, and the suggested shifts were compared to the known shifts by calculating the residual error. A uniform cylindrical phantom with cylindrical inserts of various known CT numbers was scanned with kVCBCT at 1.0 ≤ Δτ (mm) ≤ 5.0 at increments of 0.5 mm. The contrast-to-noise ratios for the inserts were measured at each Δτ.

RESULTS

For the planning CT slice thickness used in this study, there was no significant difference in residual error below a threshold equal to the planning CT slice thickness. For Δτ > 3.0 mm, residual error increased for both the head and pelvis phantom studies. The contrast-to-noise ratio is proportional to slice thickness until Δτ = 2.5 mm. Beyond this point, the contrast-to-noise ratio was not affected by Δτ.

CONCLUSION

Automatic registration accuracy is greatest when 1.0 ≤ Δτ (mm) ≤ 3.0 is used. Contrast-to-noise ratio is optimal for the 2.5 ≤ Δτ (mm) ≤ 5.0 range. Therefore 2.5 ≤ Δτ (mm) ≤ 3.0 is recommended for kVCBCT patient registration where the planning CT is 3.0 mm.

Image-based modeling of tumor shrinkage in head and neck radiation therapy

PURPOSE

Understanding the kinetics of tumor growth/shrinkage represents a critical step in quantitative assessment of therapeutics and realization of adaptive radiation therapy. This article presents a novel framework for image-based modeling of tumor change and demonstrates its performance with synthetic images and clinical cases.

METHODS

Due to significant tumor tissue content changes, similarity-based models are not suitable for describing the process of tumor volume changes. Under the hypothesis that tissue features in a tumor volume or at the boundary region are partially preserved, the kinetic change was modeled in two steps: (1) Autodetection of homologous tissue features shared by two input images using the scale invariance feature transformation (SIFT) method; and (2) establishment of a voxel-to-voxel correspondence between the images for the remaining spatial points by interpolation. The correctness of the tissue feature correspondence was assured by a bidirectional association procedure, where SIFT features were mapped from template to target images and reversely. A series of digital phantom experiments and five head and neck clinical cases were used to assess the performance of the proposed technique.

RESULTS

The proposed technique can faithfully identify the known changes introduced when constructing the digital phantoms. The subsequent feature-guided thin plate spline calculation reproduced the "ground truth" with accuracy better than 1.5 mm. For the clinical cases, the new algorithm worked reliably for a volume change as large as 30%.

CONCLUSIONS

An image-based tumor kinetic algorithm was developed to model the tumor response to radiation therapy. The technique provides a practical framework for future application in adaptive radiation therapy.

Clinical evaluation of soft tissue organ boundary visualization on cone-beam computed tomographic imaging

PURPOSE

Cone-beam computed tomographic images (CBCTs) are increasingly used for setup correction, soft tissue targeting, and image-guided adaptive radiotherapy. However, CBCT image quality is limited by low contrast and imaging artifacts. This analysis investigates the detectability of soft tissue boundaries in CBCT by performing a multiple-observer segmentation study.

METHODS AND MATERIALS

In four prostate cancer patients prostate, bladder and rectum were repeatedly delineated by five observers on CBCTs and fan-beam CTs (FBCTs). A volumetric analysis of contouring variations was performed by calculating coefficients of variation (COV: standard deviation/average volume). The topographical distribution of contouring variations was analyzed using an average surface mesh-based method.

RESULTS

Observer- and patient-averaged COVs for FBCT/CBCT were 0.09/0.19 for prostate, 0.05/0.08 for bladder, and 0.09/0.08 for rectum. Contouring variations on FBCT were significantly smaller than on CBCT for prostate (p < 0.03) and bladder (p < 0.04), but not for rectum (p < 0.37; intermodality differences). Intraobserver variations from repeated contouring of the same image set were not significant for either FBCT or CBCT (p < 0.05). Average standard deviations of individual observers' contour differences from average surface meshes on FBCT vs. CBCT were 1.5 vs. 2.1 mm for prostate, 0.7 vs. 1.4 mm for bladder, and 1.3 vs. 1.5 mm for rectum. The topographical distribution of contouring variations was similar for FBCT and CBCT.

CONCLUSION

Contouring variations were larger on CBCT than FBCT, except for rectum. Given the well-documented uncertainty in soft tissue contouring in the pelvis, improvement of CBCT image quality and establishment of well-defined soft tissue identification rules are desirable for image-guided radiotherapy.

A pseudoinverse deformation vector field generator and its applications

PURPOSE

To present, implement, and test a self-consistent pseudoinverse displacement vector field (PIDVF) generator, which preserves the location of information mapped back-and-forth between image sets.

METHODS

The algorithm is an iterative scheme based on nearest neighbor interpolation and a subsequent iterative search. Performance of the algorithm is benchmarked using a lung 4DCT data set with six CT images from different breathing phases and eight CT images for a single prostrate patient acquired on different days. A diffeomorphic deformable image registration is used to validate our PIDVFs. Additionally, the PIDVF is used to measure the self-consistency of two nondiffeomorphic algorithms which do not use a self-consistency constraint: The ITK Demons algorithm for the lung patient images and an in-house B-Spline algorithm for the prostate patient images. Both Demons and B-Spline have been QAed through contour comparison. Self-consistency is determined by using a DIR to generate a displacement vector field (DVF) between reference image R and study image S (DVF(R-S)). The same DIR is used to generate DVF(S-R). Additionally, our PIDVF generator is used to create PIDVF(S-R). Back-and-forth mapping of a set of points (used as surrogates of contours) using DVF(R-S) and DVF(S-R) is compared to back-and-forth mapping performed with DVF(R-S) and PIDVF(S-R). The Euclidean distances between the original unmapped points and the mapped points are used as a self-consistency measure.

RESULTS

Test results demonstrate that the consistency error observed in back-and-forth mappings can be reduced two to nine times in point mapping and 1.5 to three times in dose mapping when the PIDVF is used in place of the B-Spline algorithm. These self-consistency improvements are not affected by the exchanging of R and S. It is also demonstrated that differences between DVF(S-R) and PIDVF(S-R) can be used as a criteria to check the quality of the DVF.

CONCLUSIONS

Use of DVF and its PIDVF will improve the self-consistency of points, contour, and dose mappings in image guided adaptive therapy.

The influence of CT image noise on proton range calculation in radiotherapy planning

The purpose of this note is to evaluate the relationship between the stochastic errors in CT numbers and the standard deviation of the computed proton beam range in radiotherapy planning. The stochastic voxel-to-voxel variation in CT numbers called 'noise,' may be due to signal registration, processing and numerical image reconstruction technique. Noise in CT images may cause a deviation in the computed proton range from the physical proton range, even assuming that the error due to CT number-stopping power calibration is removed. To obtain the probability density function (PDF) of the computed proton range, we have used the continuing slowing down approximation (CSDA) and the uncorrelated white Gaussian noise along the proton path. The model of white noise was accepted because for the slice-based fan-beam CT scanner; the power-spectrum properties apply only to the axial (x, y) domain and the noise is uncorrelated in the z domain. However, the possible influence of the noise power spectrum on the standard deviation of the range should be investigated in the future. A random number generator was utilized for noise simulation and this procedure was iteratively repeated to obtain convergence of range PDF, which approached a Gaussian distribution. We showed that the standard deviation of the range, sigma, increases linearly with the initial proton energy, computational grid size and standard deviation of the voxel values. The 95% confidence interval width of the range PDF, which is defined as 4sigma, may reach 0.6 cm for the initial proton energy of 200 MeV, computational grid 0.25 cm and 5% standard deviation of CT voxel values. Our results show that the range uncertainty due to random errors in CT numbers may be significant and comparable to the uncertainties due to calibration of CT numbers.

Feature-based rectal contour propagation from planning CT to cone beam CT

The purpose of this work is to develop a novel feature-based registration strategy to automatically map the rectal contours from planning computed tomography (CT) (pCT) to cone beam CT (CBCT). The rectal contours were manually outlined on the pCT. A narrow band with the outlined contour as its interior surface was then constructed, so that we can exclude the volume inside the rectum in the registration process. The corresponding contour in the CBCT was found by using a feature-based registration algorithm, which consists of two steps: (1) automatically searching for control points in the pCT and CBCT based on the features of the surrounding tissue and matching the homologous control points using the scale invariance feature transformation; and (2) using the control points for a thin plate spline transformation to warp the narrow band and mapping the corresponding contours from pCT to CBCT. The proposed contour propagation technique is applied to digital phantoms and clinical cases and, in all cases, the contour mapping results are found to be clinically acceptable. For clinical cases, the method yielded satisfactory results even when there were significant rectal content changes between the pCT and CBCT scans. As a consequence, the accordance between the rectal volumes after deformable registration and the manually segmented rectum was found to be more than 90%. The proposed technique provides a powerful tool for adaptive radiotherapy of prostate, rectal, and gynecological cancers in the future.