Mass-preserving Image Registration Research Articles

CT-derived 3D-diaphragm motion in emphysema and IPF compared to normal subjects

Image registration-based local displacement analysis enables evaluation of respiratory motion between two computed tomography-captured lung volumes. The objective of this study was to compare diaphragm movement among emphysema, idiopathic pulmonary fibrosis (IPF) and normal subjects. 29 normal, 50 emphysema, and 51 IPF subjects were included. A mass preserving image registration technique was used to compute displacement vectors of local lung regions at an acinar scale. Movement of the diaphragm was assumed to be equivalent to movement of the basal lung within 5 mm from the diaphragm. Magnitudes and directions of displacement vectors were compared between the groups. Three-dimensional (3D) and apico-basal displacements were smaller in emphysema than normal subjects (P = 0.003, P = 0.002). Low lung attenuation area on expiration scan showed significant correlations with decreased 3D and apico-basal displacements (r = − 0.546, P < 0.0001; r = − 0.521, P < 0.0001) in emphysema patients. Dorsal–ventral displacement was smaller in IPF than normal subjects (P < 0.0001). The standard deviation of the displacement angle was greater in both emphysema and IPF patients than normal subjects (P < 0.0001). In conclusion, apico-basal movement of the diaphragm is reduced in emphysema while dorsal–ventral movement is reduced in IPF. Image registration technique to multi-volume CT scans provides insight into the pathophysiology of limited diaphragmatic motion in emphysema and IPF.

Relative Regional Air Volume Change Maps at the Acinar Scale Reflect Variable Ventilation in Low Lung Attenuation of COPD patients

The purpose of this study was to investigate regional air volume changes at the acinar scale of the lung in chronic obstructive pulmonary disease (COPD) patients using an image registration technique. Thirty-four emphysema patients and 24 subjects with normal chest CT and pulmonary function test (PFT) results were included in this retrospective study for which informed consent was waived by the institutional review board. After lung segmentation, a mass-preserving image registration technique was used to compute relative regional air volume changes (RRAVCs) between inspiration and expiration CT scans. After determining the appropriate thresholds of RRAVCs for low ventilation areas (LVAs), they were displayed and analyzed using color maps on the background inspiration CT image, and compared with the low attenuation area (LAA) map. Correlations between quantitative CT parameters and PFTs were assessed using Pearson's correlation test, and parameters were compared between emphysema and normal-CT patients using the Student's t-test. LVA percentage with an RRAVC threshold of 0.5 (%LVA0.5) showed the strongest correlations with FEV1/FVC (r = -0.566), FEV1 (r = -0.534), %LAA-950insp (r = 0.712), and %LAA-856exp (r = 0.775). %LVA0.5 was significantly higher (P < 0.001) in COPD patients than normal subjects. Despite the identical appearance of emphysematous lesions on the LAA-950insp map, the RRAVC map depicted a wide range of ventilation differences between these LAA clusters. RRAVC-based %LVA0.5 correlated well with FEV1/FVC, FEV1, %LAA-950insp and %LAA-856exp. RRAVC holds the potential for providing additional acinar scale functional information for emphysematous LAAs in inspiratory CT images, providing the basis for a novel set for emphysematous phenotypes.

Efficient methods for implementation of multi-level nonrigid mass-preserving image registration on GPUs and multi-threaded CPUs

Background and objectiveFaster and more accurate methods for registration of images are important for research involved in conducting population-based studies that utilize medical imaging, as well as improvements for use in clinical applications. We present a novel computation- and memory-efficient multi-level method on graphics processing units (GPU) for performing registration of two computed tomography (CT) volumetric lung images. MethodsWe developed a computation- and memory-efficient Diffeomorphic Multi-level B-Spline Transform Composite (DMTC) method to implement nonrigid mass-preserving registration of two CT lung images on GPU. The framework consists of a hierarchy of B-Spline control grids of increasing resolution. A similarity criterion known as the sum of squared tissue volume difference (SSTVD) was adopted to preserve lung tissue mass. The use of SSTVD consists of the calculation of the tissue volume, the Jacobian, and their derivatives, which makes its implementation on GPU challenging due to memory constraints. The use of the DMTC method enabled reduced computation and memory storage of variables with minimal communication between GPU and Central Processing Unit (CPU) due to ability to pre-compute values. The method was assessed on six healthy human subjects. ResultsResultant GPU-generated displacement fields were compared against the previously validated CPU counterpart fields, showing good agreement with an average normalized root mean square error (nRMS) of 0.044±0.015. Runtime and performance speedup are compared between single-threaded CPU, multi-threaded CPU, and GPU algorithms. Best performance speedup occurs at the highest resolution in the GPU implementation for the SSTVD cost and cost gradient computations, with a speedup of 112 times that of the single-threaded CPU version and 11 times over the twelve-threaded version when considering average time per iteration using a Nvidia Tesla K20X GPU. ConclusionsThe proposed GPU-based DMTC method outperforms its multi-threaded CPU version in terms of runtime. Total registration time reduced runtime to 2.9min on the GPU version, compared to 12.8min on twelve-threaded CPU version and 112.5min on a single-threaded CPU. Furthermore, the GPU implementation discussed in this work can be adapted for use of other cost functions that require calculation of the first derivatives.

Registration-based assessment of regional lung function via volumetric CT images of normal subjects vs. severe asthmatics

The purpose of this work was to explore the use of image registration-derived variables associated with computed tomographic (CT) imaging of the lung acquired at multiple volumes. As an evaluation of the utility of such an imaging approach, we explored two groups at the extremes of population ranging from normal subjects to severe asthmatics. A mass-preserving image registration technique was employed to match CT images at total lung capacity (TLC) and functional residual capacity (FRC) for assessment of regional air volume change and lung deformation between the two states. Fourteen normal subjects and thirty severe asthmatics were analyzed via image registration-derived metrics together with their pulmonary function test (PFT) and CT-based air-trapping. Relative to the normal group, the severely asthmatic group demonstrated reduced air volume change (consistent with air trapping) and more isotropic deformation in the basal lung regions while demonstrating increased air volume change associated with increased anisotropic deformation in the apical lung regions. These differences were found despite the fact that both PFT-derived TLC and FRC in the two groups were nearly 100% of predicted values. Data suggest that reduced basal-lung air volume change in severe asthmatics was compensated by increased apical-lung air volume change and that relative increase in apical-lung air volume change in severe asthmatics was accompanied by enhanced anisotropic deformation. These data suggest that CT-based deformation, assessed via inspiration vs. expiration scans, provides a tool for distinguishing differences in lung mechanics when applied to the extreme ends of a population range.

A multiscale MDCT image-based breathing lung model with time-varying regional ventilation

A novel algorithm is presented that links local structural variables (regional ventilation and deforming central airways) to global function (total lung volume) in the lung over three imaged lung volumes, to derive a breathing lung model for computational fluid dynamics simulation. The algorithm constitutes the core of an integrative, image-based computational framework for subject-specific simulation of the breathing lung. For the first time, the algorithm is applied to three multi-detector row computed tomography (MDCT) volumetric lung images of the same individual. A key technique in linking global and local variables over multiple images is an in-house mass-preserving image registration method. Throughout breathing cycles, cubic interpolation is employed to ensure C1 continuity in constructing time-varying regional ventilation at the whole lung level, flow rate fractions exiting the terminal airways, and airway deformation. The imaged exit airway flow rate fractions are derived from regional ventilation with the aid of a three-dimensional (3D) and one-dimensional (1D) coupled airway tree that connects the airways to the alveolar tissue. An in-house parallel large-eddy simulation (LES) technique is adopted to capture turbulent-transitional-laminar flows in both normal and deep breathing conditions. The results obtained by the proposed algorithm when using three lung volume images are compared with those using only one or two volume images. The three-volume-based lung model produces physiologically-consistent time-varying pressure and ventilation distribution. The one-volume-based lung model under-predicts pressure drop and yields un-physiological lobar ventilation. The two-volume-based model can account for airway deformation and non-uniform regional ventilation to some extent, but does not capture the non-linear features of the lung.

Mass preserving image registration for lung CT

This paper presents a mass preserving image registration algorithm for lung CT images. To account for the local change in lung tissue intensity during the breathing cycle, a tissue appearance model based on the principle of preservation of total lung mass is proposed. This model is incorporated into a standard image registration framework with a composition of a global affine and several free-form B-Spline transformations with increasing grid resolution. The proposed mass preserving registration method is compared to registration using the sum of squared intensity differences as a similarity function on four groups of data: 44 pairs of longitudinal inspiratory chest CT scans with small difference in lung volume; 44 pairs of longitudinal inspiratory chest CT scans with large difference in lung volume; 16 pairs of expiratory and inspiratory CT scans; and 5 pairs of images extracted at end exhale and end inhale phases of 4D-CT images. Registration errors, measured as the average distance between vessel tree centerlines in the matched images, are significantly lower for the proposed mass preserving image registration method in the second, third and fourth group, while there is no statistically significant difference between the two methods in the first group. Target registration error, assessed via a set of manually annotated landmarks in the last group, was significantly smaller for the proposed registration method.

Motion Correction in Dual Gated Cardiac PET Using Mass-Preserving Image Registration

Respiratory and cardiac motion leads to image degradation in positron emission tomography (PET) studies of the human heart. In this paper we present a novel approach to motion correction based on dual gating and mass-preserving hyperelastic image registration. Thereby, we account for intensity modulations caused by the highly nonrigid cardiac motion. This leads to accurate and realistic motion estimates which are quantitatively validated on software phantom data and carried over to clinically relevant data using a hardware phantom. For patient data, the proposed method is first evaluated in a high statistic (20 min scans) dual gating study of 21 patients. It is shown that the proposed approach properly corrects PET images for dual-cardiac as well as respiratory-motion. In a second study the list mode data of the same patients is cropped to a scan time reasonable for clinical practice (3 min). This low statistic study not only shows the clinical applicability of our method but also demonstrates its robustness against noise obtained by hyperelastic regularization.