Three-dimensional velocity mapping of lung motion using vessel bifurcation pattern matching

A hybrid proton and hyperpolarized gas tagging MRI technique for lung respiratory motion imaging: a feasibility study

The aim of this work was to develop a novel hybrid 3D hyperpolarized (HP) gas tagging MRI (t-MRI) technique and to evaluate it for lung respiratory motion measurement with comparison to deformable image registrations (DIR) methods. Three healthy subjects underwent a hybrid MRI which combines 3D HP gas t-MRI with a low resolution (Low-R, 4.5 mm isotropic voxels) 3D proton MRI (p-MRI), plus a high resolution (High-R, 2.5 mm isotropic voxels) 3D p-MRI, during breath-holds at the end-of-inhalation (EOI) and the end-of-exhalation (EOE). Displacement vector field (DVF) of the lung motion was determined from the t-MRI images by tracking tagging grids and from the High-R p-MRI using three DIR methods (B-spline based method implemented by Velocity, Free Form Deformation by MIM, and B-spline by an open source software Elastix: denoted as A, B, and C, respectively), labeled as tDVF and dDVF, respectively. The tDVF from the HP gas t-MRI was used as ground-truth reference to evaluate performance of the three DIR methods. Differences in both magnitude and angle between the tDVF and dDVFs were analyzed. The mean lung motion of the three subjects was 37.3 mm, 8.9 mm and 12.9 mm, respectively. Relatively large discrepancies were observed between the tDVF and the dDVFs as compared to previously reported DIR errors. The mean  ±  standard deviation (SD) DVF magnitude difference was 8.3  ±  5.6 mm, 9.2  ±  4.5 mm, and 9.3  ±  6.1 mm, and the mean  ±  SD DVF angular difference was 29.1  ±  12.1°, 50.1  ±  28.6°, and 39.0  ±  6.3°, for the DIR Methods A, B, and C, respectively. These preliminary results showed that the hybrid HP gas t-MRI technique revealed different lung motion patterns as compared to the DIR methods. It may provide unique perspectives in developing and evaluating DIR of the lungs. Novelty and Significance We designed a MRI protocol that includes a novel hybrid MRI technique (3D HP gas t-MRI with a low resolution 3D p-MRI) plus a high resolution 3D p-MRI. We tested the novel hybrid MRI technique on three healthy subjects for measuring regional lung respiratory motion with comparison to deformable image registrations (DIR) methods, and observed relatively large discrepancies in lung motion between HP gas t-MRI and DIR methods.

Tracking the motion trajectories of junction structures in 4D CT images of the lung

Respiratory motion poses a major challenge in lung radiotherapy. Based on 4D CT images, a variety of intensity-based deformable registration techniques have been proposed to study the pulmonary motion. However, the accuracy achievable with these approaches can be sub-optimal because the deformation is defined globally in space. Therefore, the accuracy of the alignment of local structures may be compromised. In this work, we propose a novel method to detect a large collection of natural junction structures in the lung and use them as the reliable markers to track the lung motion. Specifically, detection of the junction centers and sizes is achieved by analysis of local shape profiles on one segmented image. To track the temporal trajectory of a junction, the image intensities within a small region of interest surrounding the center are selected as its signature. Under the assumption of the cyclic motion, we describe the trajectory by a closed B-spline curve and search for the control points by maximizing a metric of combined correlation coefficients. Local extrema are suppressed by improving the initial conditions using random walks from pair-wise optimizations. Several descriptors are introduced to analyze the motion trajectories. Our method was applied to 13 real 4D CT images. More than 700 junctions in each case are detected with an average positive predictive value of greater than 90%. The average tracking error between automated and manual tracking is sub-voxel and smaller than the published results using the same set of data.

Automated detection of junctions structures and tracking of their trajectories in 4D images

Junction structures, as the natural anatomical markers, are useful to study the organ or tumor motion. However, detection and tracking of the junctions in four-dimensional (4D) images are challenging. The paper presents a novel framework to automate this task. Detection of their centers and sizes is first achieved by an analysis of local shape profiles on one segmented reference image. Junctions are then separately tracked by simultaneously using neighboring intensity features from all images. Defined by a closed B-spline space curve, the individual trajectory is assumed to be cyclic and obtained by maximizing the metric of combined correlation coefficients. Local extrema are suppressed by improving the initial conditions using random walks from pair-wise optimizations. Our approach has been applied to analyze the vessel junctions in five real 4D respiration-gated computed tomography (CT) image datasets with promising results. More than 500 junctions in the lung are detected with an average accuracy of greater than 85% and the mean error between the automated and the manual tracking is sub-voxel.

Analysis of deformable image registration accuracy using computational modeling

Computer aided modeling of anatomic deformation, allowing various techniques and protocols in radiation therapy to be systematically verified and studied, has become increasingly attractive. In this study the potential issues in deformable image registration (DIR) were analyzed based on two numerical phantoms: One, a synthesized, low intensity gradient prostate image, and the other a lung patient's CT image data set. Each phantom was modeled with region-specific material parameters with its deformation solved using a finite element method. The resultant displacements were used to construct a benchmark to quantify the displacement errors of the Demons and B-Spline-based registrations. The results show that the accuracy of these registration algorithms depends on the chosen parameters, the selection of which is closely associated with the intensity gradients of the underlying images. For the Demons algorithm, both single resolution (SR) and multiresolution (MR) registrations required approximately 300 iterations to reach an accuracy of 1.4 mm mean error in the lung patient's CT image (and 0.7 mm mean error averaged in the lung only). For the low gradient prostate phantom, these algorithms (both SR and MR) required at least 1600 iterations to reduce their mean errors to 2 mm. For the B-Spline algorithms, best performance (mean errors of 1.9 mm for SR and 1.6 mm for MR, respectively) on the low gradient prostate was achieved using five grid nodes in each direction. Adding more grid nodes resulted in larger errors. For the lung patient's CT data set, the B-Spline registrations required ten grid nodes in each direction for highest accuracy (1.4 mm for SR and 1.5 mm for MR). The numbers of iterations or grid nodes required for optimal registrations depended on the intensity gradients of the underlying images. In summary, the performance of the Demons and B-Spline registrations have been quantitatively evaluated using numerical phantoms. The results show that parameter selection for optimal accuracy is closely related to the intensity gradients of the underlying images. Also, the result that the DIR algorithms produce much lower errors in heterogeneous lung regions relative to homogeneous (low intensity gradient) regions, suggests that feature-based evaluation of deformable image registration accuracy must be viewed cautiously.