Comparison of Image Registration Based Measures of Regional Lung Ventilation from Dynamic Spiral CT with Xe-CT

Purpose: Regional lung volume change as a function of lung inflation serves as an index of parenchymal and airway status as well as an index of regional ventilation and can be used to detect pathologic changes over time. In this article, we propose a new regional measure of lung mechanics --- the specific air volume change by corrected Jacobian. Methods: 4DCT and Xe-CT data sets from four adult sheep are used in this study. Nonlinear, 3D image registration is applied to register an image acquired near end inspiration to an image acquired near end expiration. Approximately 200 annotated anatomical points are used as landmarks to evaluate registration accuracy. Three different registration-based measures of regional lung mechanics are derived and compared: the specific air volume change calculated from the Jacobian (SAJ); the specific air volume change calculated by the corrected Jacobian (SACJ); and the specific air volume change by intensity change (SAI). Results: After registration, the mean registration error is on the order of 1 mm. For cubical ROIs in cubes with size 20 mm $\times$ 20 mm $\times$ 20 mm, the SAJ and SACJ measures show significantly higher correlation (linear regression, average $r^2=0.75$ and $r^2=0.82$) with the Xe-CT based measure of specific ventilation (sV) than the SAI measure. For ROIs in slabs along the ventral-dorsal vertical direction with size of 150 mm $\times$ 8 mm $\times$ 40 mm, the SAJ, SACJ, and SAI all show high correlation (linear regression, average $r^2=0.88$, $r^2=0.92$ and $r^2=0.87$) with the Xe-CT based sV without significant differences when comparing between the three methods. Conclusion: Given a deformation field by an image registration algorithm, significant differences between the SAJ, SACJ, and SAI measures were found at a regional level compared to the Xe-CT sV in four sheep that were studied

Inverse-Consistent Determination of Young\u27s Modulus of Human Lung

Human lung undergoes respiration-induced deformation due to sequential inhalation and exhalation. Accurate determination of lung deformation is crucial for tumor localization and targeted radiotherapy in patients with lung cancer. Numerical modeling of human lung dynamics based on underlying physics and physiology enables simulation and virtual visualization of lung deformation. Dynamical modeling is numerically complicated by the lack of information on lung elastic behavior, structural heterogeneity as well as boundary constrains. This study integrates physics-based modeling and image-based data acquisition to develop the patient-specific biomechanical model and consequently establish the first consistent Young\u27s modulus (YM) of human lung. This dissertation has four major components: (i) develop biomechanical model for computation of the flow and deformation characteristics that can utilize subject-specific, spatially-dependent lung material property; (ii) develop a fusion algorithm to integrate deformation results from a deformable image registration (DIR) and physics-based modeling using the theory of Tikhonov regularization; (iii) utilize fusion algorithm to establish unique and consistent patient specific Young\u27s modulus and; (iv) validate biomechanical model utilizing established patient-specific elastic property with imaging data. The simulation is performed on three dimensional lung geometry reconstructed from four-dimensional computed tomography (4DCT) dataset of human subjects. The heterogeneous Young\u27s modulus is estimated from a linear elastic deformation model with the same lung geometry and 4D lung DIR. The biomechanical model adequately predicts the spatio-temporal lung deformation, consistent with data obtained from imaging. The accuracy of the numerical solution is enhanced through fusion with the imaging data beyond the classical comparison of the two sets of data. Finally, the fused displacement results are used to establish unique and consistent patient-specific elastic property of the lung

What can computed tomography and magnetic resonance imaging tell us about ventilation?

This review provides a summary of pulmonary functional imaging approaches for determining pulmonary ventilation, with a specific focus on multi-detector x-ray computed tomography and magnetic resonance imaging (MRI). We provide the important functional definitions of pulmonary ventilation typically used in medicine and physiology and discuss the fact that some of the imaging literature describes gas distribution abnormalities in pulmonary disease that may or may not be related to the physiological definition or clinical interpretation of ventilation. We also review the current state-of-the-field in terms of the key physiological questions yet unanswered related to ventilation and gas distribution in lung disease. Current and emerging imaging research methods are described, including their strengths and the challenges that remain to translate these methods to more wide-spread research and clinical use. We also examine how computed tomography and MRI might be used in the future to gain more insight into gas distribution and ventilation abnormalities in pulmonary disease

Imaging Biomarkers of Pulmonary Structure and Function

Asthma and chronic obstructive pulmonary disease (COPD) are characterized by airflow limitations resulting from airway obstruction and/or tissue destruction. The diagnosis and monitoring of these pulmonary diseases is primarily performed using spirometry, specifically the forced expiratory volume in one second (FEV1), which measures global airflow obstruction and provides no regional information of the different underlying disease pathologies. The limitations of spirometry and current therapies for lung disease patients have motivated the development of pulmonary imaging approaches, such as computed tomography (CT) and magnetic resonance imaging (MRI). Inhaled hyperpolarized noble gas MRI, specifically using helium-3 (3He) and xenon-129 (129Xe) gases, provides a way to quantify pulmonary ventilation by visualizing lung regions accessed by gas during a breath-hold, and alternatively, regions that are not accessed - coined “ventilation defects.” Despite the strong foundation and many advantages hyperpolarized 3He MRI has to offer research and patient care, clinical translation has been inhibited in part due to the cost and need for specialized equipment, including multinuclear-MR hardware and polarizers, and personnel. Accordingly, our objective was to develop and evaluate imaging biomarkers of pulmonary structure and function using MRI and CT without the use of exogenous contrast agents or specialized equipment. First, we developed and compared CT parametric response maps (PRM) with 3He MR ventilation images in measuring gas-trapping and emphysema in ex-smokers with and without COPD. We observed that in mild-moderate COPD, 3He MR ventilation abnormalities were related to PRM gas-trapping whereas in severe COPD, ventilation abnormalities correlated with both PRM gas-trapping and PRM emphysema. We then developed and compared pulmonary ventilation abnormalities derived from Fourier decomposition of free-breathing proton (1H) MRI (FDMRI) with 3He MRI in subjects with COPD and bronchiectasis. This work demonstrated that FDMRI and 3He MRI ventilation defects were strongly related in COPD, but not in bronchiectasis subjects. In COPD only, FDMRI ventilation defects were spatially related with 3He MRI ventilation defects and emphysema. Based on the FDMRI biomarkers developed in patients with COPD and bronchiectasis, we then evaluated ventilation heterogeneity in patients with severe asthma, both pre- and post-salbutamol as well as post-methacholine challenge, using FDMRI and 3He MRI. FDMRI free-breathing ventilation abnormalities were correlated with but under-estimated 3He MRI static ventilation defects. Finally, based on the previously developed free-breathing MRI approach, we developed a whole-lung free-breathing pulmonary 1H MRI technique to measure regional specific-ventilation and evaluated both asthmatics and healthy volunteers. These measurements not only provided similar information as specific-ventilation measured using plethysmography, but also information about regional ventilation defects that were correlated with 3He MRI ventilation abnormalities. These results demonstrated that whole-lung free-breathing 1H MRI biomarker of specific-ventilation may reflect ventilation heterogeneity and/or gas-trapping in asthma. These important findings indicate that imaging biomarkers of pulmonary structure and function using MRI and CT have the potential to regionally reveal the different pathologies in COPD and asthma without the use of exogenous contrast agents. The development and validation of these clinically meaningful imaging biomarkers are critically required to accelerate pulmonary imaging translation from the research workbench to being a part of the clinical workflow, with the overall goal to improve patient outcomes

Improving the Accuracy of CT-derived Attenuation Correction in Respiratory-Gated PET/CT Imaging

The effect of respiratory motion on attenuation correction in Fludeoxyglucose (18F) positron emission tomography (FDG-PET) was investigated. Improvements to the accuracy of computed tomography (CT) derived attenuation correction were obtained through the alignment of the attenuation map to each emission image in a respiratory gated PET scan. Attenuation misalignment leads to artefacts in the reconstructed PET image and several methods were devised for evaluating the attenuation inaccuracies caused by this. These methods of evaluation were extended to finding the frame in the respiratory gated PET which best matched the CT. This frame was then used as a reference frame in mono-modality compensation for misalignment. Attenuation correction was found to affect the quantification of tumour volumes; thus a regional analysis was used to evaluate the impact of mismatch and the benefits of compensating for misalignment. Deformable image registration was used to compensate for misalignment, however, there were inaccuracies caused by the poor signal-to-noise ratio (SNR) in PET images. Two models were developed that were robust to a poor SNR allowing for the estimation of deformation from very noisy images. Firstly, a cross population model was developed by statistically analysing the respiratory motion in 10 4DCT scans. Secondly, a 1D model of respiration was developed based on the physiological function of respiration. The 1D approach correctly modelled the expansion and contraction of the lungs and the differences in the compressibility of lungs and surrounding tissues. Several additional models were considered but were ruled out based on their poor goodness of fit to 4DCT scans. Approaches to evaluating the developed models were also used to assist with optimising for the most accurate attenuation correction. It was found that the multimodality registration of the CT image to the PET image was the most accurate approach to compensating for attenuation correction mismatch. Mono-modality image registration was found to be the least accurate approach, however, incorporating a motion model improved the accuracy of image registration. The significance of these findings is twofold. Firstly, it was found that motion models are required to improve the accuracy in compensating for attenuation correction mismatch and secondly, a validation method was found for comparing approaches to compensating for attenuation mismatch

Improving the Accuracy of CT-derived Attenuation Correction in Respiratory-Gated PET/CT Imaging

The effect of respiratory motion on attenuation correction in Fludeoxyglucose (18F) positron emission tomography (FDG-PET) was investigated. Improvements to the accuracy of computed tomography (CT) derived attenuation correction were obtained through the alignment of the attenuation map to each emission image in a respiratory gated PET scan. Attenuation misalignment leads to artefacts in the reconstructed PET image and several methods were devised for evaluating the attenuation inaccuracies caused by this. These methods of evaluation were extended to finding the frame in the respiratory gated PET which best matched the CT. This frame was then used as a reference frame in mono-modality compensation for misalignment. Attenuation correction was found to affect the quantification of tumour volumes; thus a regional analysis was used to evaluate the impact of mismatch and the benefits of compensating for misalignment. Deformable image registration was used to compensate for misalignment, however, there were inaccuracies caused by the poor signal-to-noise ratio (SNR) in PET images. Two models were developed that were robust to a poor SNR allowing for the estimation of deformation from very noisy images. Firstly, a cross population model was developed by statistically analysing the respiratory motion in 10 4DCT scans. Secondly, a 1D model of respiration was developed based on the physiological function of respiration. The 1D approach correctly modelled the expansion and contraction of the lungs and the differences in the compressibility of lungs and surrounding tissues. Several additional models were considered but were ruled out based on their poor goodness of fit to 4DCT scans. Approaches to evaluating the developed models were also used to assist with optimising for the most accurate attenuation correction. It was found that the multimodality registration of the CT image to the PET image was the most accurate approach to compensating for attenuation correction mismatch. Mono-modality image registration was found to be the least accurate approach, however, incorporating a motion model improved the accuracy of image registration. The significance of these findings is twofold. Firstly, it was found that motion models are required to improve the accuracy in compensating for attenuation correction mismatch and secondly, a validation method was found for comparing approaches to compensating for attenuation mismatch

Evaluation of deformable image registration for improved 4D CT-derived ventilation for image guided radiotherapy

Recent treatment planning studies have demonstrated the use of physiologic images in radiation therapy treatment planning to identify regions for functional avoidance. This image-guided radiotherapy (IGRT) strategy may reduce the injury and/or functional loss following thoracic radiotherapy. 4D computed tomography (CT), developed for radiotherapy treatment planning, is a relatively new imaging technique that allows the acquisition of a time-varying sequence of 3D CT images of the patient\u27s lungs through the respiratory cycle. Guerrero et al. developed a method to calculate ventilation imaging from 4D CT, which is potentially better suited and more broadly available for IGRT than the current standard imaging methods. The key to extracting function information from 4D CT is the construction of a volumetric deformation field that accurately tracks the motion of the patient\u27s lungs during the respiratory cycle. The spatial accuracy of the displacement field directly impacts the ventilation images; higher spatial registration accuracy will result in less ventilation image artifacts and physiologic inaccuracies. Presently, a consistent methodology for spatial accuracy evaluation of the DIR transformation is lacking. Evaluation of the 4D CT-derived ventilation images will be performed to assess correlation with global measurements of lung ventilation, as well as regional correlation of the distribution of ventilation with the current clinical standard SPECT. This requires a novel framework for both the detailed assessment of an image registration algorithm\u27s performance characteristics as well as quality assurance for spatial accuracy assessment in routine application. Finally, we hypothesize that hypo-ventilated regions, identified on 4D CT ventilation images, will correlate with hypo-perfused regions in lung cancer patients who have obstructive lesions. A prospective imaging trial of patients with locally advanced non-small-cell lung cancer will allow this hypothesis to be tested. These advances are intended to contribute to the validation and clinical implementation of CT-based ventilation imaging in prospective clinical trials, in which the impact of this imaging method on patient outcomes may be tested

Free-breathing Pulmonary MR Imaging to Quantify Regional Ventilation

Purpose: To measure regional specific ventilation with free-breathing hydrogen 1 (1H) magnetic resonance (MR) imaging without exogenous contrast material and to investigate correlations with hyperpolarized helium 3 (3He) MR imaging and pulmonary function test measurements in healthy volunteers and patients with asthma. Materials and Methods: Subjects underwent free-breathing 1H and static breath-hold hyperpolarized 3He MR imaging as well as spirometry and plethysmography; participants were consecutively recruited between January and June 2017. Free-breathing 1H MR imaging was performed with an optimized balanced steady-state free-precession sequence; images were retrospectively grouped into tidal inspiration or tidal expiration volumes with exponentially weighted phase interpolation. MR imaging volumes were coregistered by using optical flow deformable registration to generate 1H MR imaging-derived specific ventilation maps. Hyperpolarized 3He MR imaging- and 1H MR imaging-derived specific ventilation maps were coregistered to quantify regional specific ventilation within hyperpolarized 3He MR imaging ventilation masks. Differences between groups were determined with the Mann-Whitney test and relationships were determined with Spearman (ρ) correlation coefficients. Statistical analyses were performed with software. Results: Thirty subjects (median age: 50 years; interquartile range [IQR]: 30 years), including 23 with asthma and seven healthy volunteers, were evaluated. Both 1H MR imaging-derived specific ventilation and hyperpolarized 3He MR imaging-derived ventilation percentage were significantly greater in healthy volunteers than in patients with asthma (specific ventilation: 0.14 [IQR: 0.05] vs 0.08 [IQR: 0.06], respectively, P \u3c .0001; ventilation percentage: 99% [IQR: 1%] vs 94% [IQR: 5%], P \u3c .0001). For all subjects, 1H MR imaging-derived specific ventilation correlated with plethysmography-derived specific ventilation (ρ = 0.54, P = .002) and hyperpolarized 3He MR imaging-derived ventilation percentage (ρ = 0.67, P \u3c .0001) as well as with forced expiratory volume in 1 second (FEV1) (ρ = 0.65, P = .0001), ratio of FEV1 to forced vital capacity (ρ = 0.75, P \u3c .0001), ratio of residual volume to total lung capacity (ρ = -0.68, P \u3c .0001), and airway resistance (ρ = -0.51, P = .004). 1H MR imaging-derived specific ventilation was significantly greater in the gravitational-dependent versus nondependent lung in healthy subjects (P = .02) but not in patients with asthma (P = .1). In patients with asthma, coregistered 1H MR imaging specific ventilation and hyperpolarized 3He MR imaging maps showed that specific ventilation was diminished in corresponding 3He MR imaging ventilation defects (0.05 ± 0.04) compared with well-ventilated regions (0.09 ± 0.05) (P \u3c .0001). Conclusion: 1H MR imaging-derived specific ventilation correlated with plethysmography-derived specific ventilation and ventilation defects seen by using hyperpolarized 3He MR imaging. © RSNA, 2018 Online supplemental material is available for this article

Evaluating Small Airways Disease in Asthma and COPD using the Forced Oscillation Technique and Magnetic Resonance Imaging

Obstructive lung disease, including asthma and chronic obstructive pulmonary disease (COPD), is characterized by heterogeneous ventilation. Unfortunately, the underlying structure-function relationships and the relationships between measurements of heterogeneity and patient quality-of-life in obstructive lung disease are not well understood. Hyperpolarized noble gas MRI is used to visualize and quantify ventilation distribution and the forced oscillation technique (FOT) applies a multi-frequency pressure oscillation at the mouth to measure respiratory impedance to airflow (including resistance and reactance). My objective was to use FOT, ventilation MRI and computational airway tree modeling to better understand ventilation heterogeneity in asthma and COPD. FOT-measured respiratory system impedance was correlated with MRI ventilation heterogeneity and both were related to quality-of-life in asthma and COPD. FOT-measurements and model-predictions of reactance and small-airways resistance were correlated in asthma and COPD respectively. This study is the first to demonstrate the relationships between FOT-measured impedance, MRI ventilation heterogeneity, and patient quality-of-life