Expert landmark correspondences are widely reported for evaluating deformable image registration (DIR) spatial accuracy. In this report, we present a framework for objective evaluation of DIR spatial accuracy using large sets of expert-determined landmark point pairs. Large samples (>1100) of pulmonary landmark point pairs were manually generated for five cases. Estimates of inter- and intra-observer variation were determined from repeated registration. Comparative evaluation of DIR spatial accuracy was performed for two algorithms, a gradient-based optical flow algorithm and a landmark-based moving least-squares algorithm. The uncertainty of spatial error estimates was found to be inversely proportional to the square root of the number of landmark point pairs and directly proportional to the standard deviation of the spatial errors. Using the statistical properties of this data, we performed sample size calculations to estimate the average spatial accuracy of each algorithm with 95% confidence intervals within a 0.5 mm range. For the optical flow and moving least-squares algorithms, the required sample sizes were 1050 and 36, respectively. Comparative evaluation based on fewer than the required validation landmarks results in misrepresentation of the relative spatial accuracy. This study demonstrates that landmark pairs can be used to assess DIR spatial accuracy within a narrow uncertainty range.

Two calculation methods to produce ventilation images from four-dimensional computed tomography (4DCT) acquired without added contrast have been reported. We reported a method to obtain ventilation images using deformable image registration (DIR) and the underlying CT density information. A second method performs the ventilation image calculation from the DIR result alone, using the Jacobian determinant of the deformation field to estimate the local volume changes resulting from ventilation. For each of these two approaches, there are variations on their implementation. In this study, two implementations of the Jacobian-based methodology are evaluated, as well as a single density change-based model for calculating the physiologic specific ventilation from 4DCT. In clinical practice, (99m)Tc-labeled aerosol single photon emission computed tomography (SPECT) is the standard method used to obtain ventilation images in patients. In this study, the distributions of ventilation obtained from the CT-based ventilation image calculation methods are compared with those obtained from the clinical standard SPECT ventilation imaging. Seven patients with 4DCT imaging and standard (99m)Tc-labeled aerosol SPECT/CT ventilation imaging obtained on the same day as part of a prospective validation study were selected. The results of this work demonstrate the equivalence of the Jacobian-based methodologies for quantifying the specific ventilation on a voxel scale. Additionally, we found that both Jacobian- and density-change-based methods correlate well with global measurements of the resting tidal volume. Finally, correlation with the clinical SPECT was assessed using the Dice similarity coefficient, which showed statistically higher (p-value < 10(-4)) correlation between density-change-based specific ventilation and the clinical reference than did either Jacobian-based implementation.

A novel method for dynamic ventilation imaging of the full respiratory cycle from four-dimensional computed tomography (4D CT) acquired without added contrast is presented. Three cases with 4D CT images obtained with respiratory gated acquisition for radiotherapy treatment planning were selected. Each of the 4D CT data sets was acquired during resting tidal breathing. A deformable image registration algorithm mapped each (voxel) corresponding tissue element across the 4D CT data set. From local average CT values, the change in fraction of air per voxel (i.e. local ventilation) was calculated. A 4D ventilation image set was calculated using pairs formed with the maximum expiration image volume, first the exhalation then the inhalation phases representing a complete breath cycle. A preliminary validation using manually determined lung volumes was performed. The calculated total ventilation was compared to the change in contoured lung volumes between the CT pairs (measured volume). A linear regression resulted in a slope of 1.01 and a correlation coefficient of 0.984 for the ventilation images. The spatial distribution of ventilation was found to be case specific and a 30% difference in mass-specific ventilation between the lower and upper lung halves was found. These images may be useful in radiotherapy planning.

Online adaptive radiation therapy (ART) promises the ability to deliver an 20 optimal treatment in response to daily patient anatomic variation. A major technical barrier for the clinical implementation of online ART is the requirement of rapid image segmentation. Deformable image registration (DIR) has been used as an automated segmentation method to transfer tumor/organ contours from the planning image to daily images. However, the current 25 computational time of DIR is insufficient for online ART. In this work, this issue is addressed by using computer graphics processing units (GPUs). A greyscale based DIR algorithm called demons and five of its variants were implemented on GPUs using the Compute Unified Device Architecture (CUDA) programming environment. The spatial accuracy of these algorithms was 30 evaluated over five sets of pulmonary 4DCT images with an average size of 256×256×100 and more than 1,100 expert-determined landmark point pairs each. For all the testing scenarios presented in this paper, the GPU-based DIR computation required around 7 to 11 seconds to yield an average 3D error ranging from 1.5 to 1.8 mm. It is interesting to find out that the original passive 35 force demons algorithms outperform subsequently proposed variants based on the combination of accuracy, efficiency, and ease of implementation. X. Gu et al.2