Validating and improving CT ventilation imaging by correlating with ventilation 4D-PET/CT using 68Ga-labeled nanoparticles.

PURPOSE: CT ventilation imaging is a novel functional lung imaging modality based on deformable image registration. The authors present the first validation study of CT ventilation using positron emission tomography with (68)Ga-labeled nanoparticles (PET-Galligas). The authors quantify this agreement for different CT ventilation metrics and PET reconstruction parameters. METHODS: PET-Galligas ventilation scans were acquired for 12 lung cancer patients using a four-dimensional (4D) PET/CT scanner. CT ventilation images were then produced by applying B-spline deformable image registration between the respiratory correlated phases of the 4D-CT. The authors test four ventilation metrics, two existing and two modified. The two existing metrics model mechanical ventilation (alveolar air-flow) based on Hounsfield unit (HU) change (VHU) or Jacobian determinant of deformation (VJac). The two modified metrics incorporate a voxel-wise tissue-density scaling (ρVHU and ρVJac) and were hypothesized to better model the physiological ventilation. In order to assess the impact of PET image quality, comparisons were performed using both standard and respiratory-gated PET images with the former exhibiting better signal. Different median filtering kernels (σm = 0 or 3 mm) were also applied to all images. As in previous studies, similarity metrics included the Spearman correlation coefficient r within the segmented lung volumes, and Dice coefficient d20 for the (0 - 20)th functional percentile volumes. RESULTS: The best agreement between CT and PET ventilation was obtained comparing standard PET images to the density-scaled HU metric (ρVHU) with σm = 3 mm. This leads to correlation values in the ranges 0.22 ≤ r ≤ 0.76 and 0.38 ≤ d20 ≤ 0.68, with r = 0.42 ± 0.16 and d20 = 0.52 ± 0.09 averaged over the 12 patients. Compared to Jacobian-based metrics, HU-based metrics lead to statistically significant improvements in r and d20 (p < 0.05), with density scaled metrics also showing higher r than for unscaled versions (p < 0.02). r and d20 were also sensitive to image quality, with statistically significant improvements using standard (as opposed to gated) PET images and with application of median filtering. CONCLUSIONS: The use of modified CT ventilation metrics, in conjunction with PET-Galligas and careful application of image filtering has resulted in improved correlation compared to earlier studies using nuclear medicine ventilation. However, CT ventilation and PET-Galligas do not always provide the same functional information. The authors have demonstrated that the agreement can improve for CT ventilation metrics incorporating a tissue density scaling, and also with increasing PET image quality. CT ventilation imaging has clear potential for imaging regional air volume change in the lung, and further development is warranted

3-D lung deformation and function from respiratory-gated 4-D x-ray CT images : application to radiation treatment planning.

Many lung diseases or injuries can cause biomechanical or material property changes that can alter lung function. While the mechanical changes associated with the change of the material properties originate at a regional level, they remain largely asymptomatic and are invisible to global measures of lung function until they have advanced significantly and have aggregated. In the realm of external beam radiation therapy of patients suffering from lung cancer, determination of patterns of pre- and post-treatment motion, and measures of regional and global lung elasticity and function are clinically relevant. In this dissertation, we demonstrate that 4-D CT derived ventilation images, including mechanical strain, provide an accurate and physiologically relevant assessment of regional pulmonary function which may be incorporated into the treatment planning process. Our contributions are as follows: (i) A new volumetric deformable image registration technique based on 3-D optical flow (MOFID) has been designed and implemented which permits the possibility of enforcing physical constraints on the numerical solutions for computing motion field from respiratory-gated 4-D CT thoracic images. The proposed optical flow framework is an accurate motion model for the thoracic CT registration problem. (ii) A large displacement landmark-base elastic registration method has been devised for thoracic CT volumetric image sets containing large deformations or changes, as encountered for example in registration of pre-treatment and post-treatment images or multi-modality registration. (iii) Based on deformation maps from MOFIO, a novel framework for regional quantification of mechanical strain as an index of lung functionality has been formulated for measurement of regional pulmonary function. (iv) In a cohort consisting of seven patients with non-small cell lung cancer, validation of physiologic accuracy of the 4-0 CT derived quantitative images including Jacobian metric of ventilation, Vjac, and principal strains, (V?1, V?2, V?3, has been performed through correlation of the derived measures with SPECT ventilation and perfusion scans. The statistical correlations with SPECT have shown that the maximum principal strain pulmonary function map derived from MOFIO, outperforms all previously established ventilation metrics from 40-CT. It is hypothesized that use of CT -derived ventilation images in the treatment planning process will help predict and prevent pulmonary toxicity due to radiation treatment. It is also hypothesized that measures of regional and global lung elasticity and function obtained during the course of treatment may be used to adapt radiation treatment. Having objective methods with which to assess pre-treatment global and regional lung function and biomechanical properties, the radiation treatment dose can potentially be escalated to improve tumor response and local control

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

Comparison of hyperpolarised gas MRI and CT-based surrogates of ventilation

Background: Non-contrast CT-based surrogates of regional ventilation derived from pulmonary images acquired at multiple inflation levels have been proposed as alternatives to established modalities. However, their physiological accuracy has yet to be validated prior to clinical translation. Purpose: To address the hypothesis that these surrogates can provide information comparable to a direct measure of ventilation from hyperpolarised gas MRI ventilation via: i. development of a methodology for registering CT and gas MRI. ii. comparison of these surrogates with gas MRI at the lobar level. iii. evaluation of the impact of inflation levels when comparing gas MRI and ventilation CT. iv. development of an image acquisition and analysis framework to facilitate spatial correlations of both techniques. v. assessment of the effect of using different gases on the correlation. Methods: i. A method to indirectly register gas MRI to CT via same-breath 1H-structural MR images was developed and its accuracy was assessed. ii. A ventilation model based on expansion of lobar CT segmentations was compared with gas MRI lobar ventilation measurements. iii. The spatial overlap of ventilation CT was compared to gas MRI acquired at two different inflation levels. iv. An image acquisition protocol was designed to minimise differences in acquisition settings between scans such as posture and breathing manoeuvre and analysis methods were developed to enable direct regional and voxel level correlations. v. The effect of using two different noble gases, namely, 3He and 129Xe, on correlation with ventilation CT was assessed. Results: i. The indirect method of registration was more accurate than direct registration. ii. Despite subtle differences, lobar ventilation measurements derived from CT and hyperpolarised gas MRI were comparable. iii. Comparison of ventilation CT and gas MRI varied with inflation state. iv. The spatial correlation between ventilation CT and gas MRI increased at coarser levels. v. A marked improvement in correlation was observed for 3He and 129Xe MRI in contrast to when ventilation CT was compared with either 3He and 129Xe MRI. Conclusion: Although CT-based surrogates of ventilation show promise for replacing established ventilation modalities such as hyperpolarised gas MRI, particularly at coarser levels, they cannot be assumed to be equivalent to the techniques they purport to replace

Anatomical and Functional Lung Imaging with Volumetric Computed Tomography in Non-Small Cell Lung Cancer

Non-small cell lung cancer (NSCLC) is one of the most diagnosed cancers in Canada, and the leading cause of cancer deaths. A significant challenge in treating NSCLC is balancing aggressive treatment with the potentially severe side effects. In radiation therapy, the management of respiratory motion and the risks of radiation-induced lung injury (RILI) pose significant challenges. 4-dimensional computed tomography (4D-CT) is an important part of motion management, but images often suffer from motion-induced artifacts. Volumetric CT scanners with wide axial field-of-view (aFOV) may reduce these artifacts and present an opportunity to advance CT-based functional lung imaging. Chapter 2 presents a phantom imaging study to investigate the suitability of a 256-slice volumetric CT (vCT) scanner for radiotherapy treatment planning. The density of the highest density materials was under-estimated by the scanner, which can be addressed with the use of an appropriate relative electron density (RED) curve. An average RED curve for all aFOV settings may be used. Chapter 3 presents a study of phantom and NSCLC patient 4D-CT images acquired on a clinical scanner and a vCT scanner. The v4D-CT images were re-sampled to simulate a conventional acquisition using a narrow aFOV clinical scanner. The phantom images demonstrated that target contouring variability decreased in v4D-CT imaging as compared to clinical 4D-CT. In the patient images, mean Hausdorff distance between organs at risk (OAR) contours was significantly correlated to respiratory phase, indicating that motion artifacts contribute to this variability. Chapter 4 presents a novel acquisition and analysis pipeline to image lung ventilation (V), perfusion (Q) and V/Q ratio in a single volumetric CT scan. In a porcine study, these images of V and Q were significantly correlated to standard Xe-enhanced ventilation and PET perfusion images in voxel-wise analysis. In a NSCLC patient study, the images were sensitive to changes in V and Q between baseline imaging and follow-up 6 weeks after radiotherapy. In this thesis, I demonstrate that volumetric CT scanners are suitable for use in radiation therapy simulation and treatment planning, and detail two scanning protocols which may reduce the challenges posed by respiratory motion and RILI risk in NSCLC

Texture Analysis and Machine Learning to Predict Pulmonary Ventilation from Thoracic Computed Tomography

Chronic obstructive pulmonary disease (COPD) leads to persistent airflow limitation, causing a large burden to patients and the health care system. Thoracic CT provides an opportunity to observe the structural pathophysiology of COPD, whereas hyperpolarized gas MRI provides images of the consequential ventilation heterogeneity. However, hyperpolarized gas MRI is currently limited to research centres, due to the high cost of gas and polarization equipment. Therefore, I developed a pipeline using texture analysis and machine learning methods to create predicted ventilation maps based on non-contrast enhanced, single-volume thoracic CT. In a COPD cohort, predicted ventilation maps were qualitatively and quantitatively related to ground-truth MRI ventilation, and both maps were related to important patient lung function and quality-of-life measures. This study is the first to demonstrate the feasibility of predicting hyperpolarized MRI-based ventilation from single-volume, breath-hold thoracic CT, which has potential to translate pulmonary ventilation information to widely available thoracic CT imaging

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

Quantifying the reproducibility of lung ventilation images between 4-Dimensional Cone Beam CT and 4-Dimensional CT.

PURPOSE: Computed tomography ventilation imaging derived from four-dimensional cone beam CT (CTVI4DCBCT ) can complement existing 4DCT-based methods (CTVI4DCT ) to track lung function changes over a course of lung cancer radiation therapy. However, the accuracy of CTVI4DCBCT needs to be assessed since anatomic 4DCBCT has demonstrably poor image quality and small field of view (FOV) compared to treatment planning 4DCT. We perform a direct comparison between short interval CTVI4DCBCT and CTVI4DCT pairs to understand the patient specific image quality factors affecting the intermodality CTVI reproducibility in the clinic. METHODS AND MATERIALS: We analysed 51 pairs of 4DCBCT and 4DCT scans acquired within 1 day of each other for nine lung cancer patients. To assess the impact of image quality, CTVIs were derived from 4DCBCT scans reconstructed using both standard Feldkamp-Davis-Kress backprojection (CTVIFDK4DCBCT) and an iterative McKinnon-Bates Simultaneous Algebraic Reconstruction Technique (CTVIMKBSART4DCBCT). Also, the influence of FOV was assessed by deriving CTVIs from 4DCT scans that were cropped to a similar FOV as the 4DCBCT scans (CTVIcrop4DCT), or uncropped (CTVIuncrop4DCT). All CTVIs were derived by performing deformable image registration (DIR) between the exhale and inhale phases and evaluating the Jacobian determinant of deformation. Reproducibility between corresponding CTVI4DCBCT and CTVI4DCT pairs was quantified using the voxel-wise Spearman rank correlation and the Dice similarity coefficient (DSC) for ventilation defect regions (identified as the lower quartile of ventilation values). Mann-Whitney U-tests were applied to determine statistical significance of each reconstruction and cropping condition. RESULTS: The (mean ± SD) Spearman correlation between CTVIFDK4DCBCT and CTVIuncrop4DCT was 0.60 ± 0.23 (range -0.03-0.88) and the DSC was 0.64 ± 0.12 (0.34-0.83). By comparison, correlations between CTVIMKBSART4DCBCT and CTVIuncrop4DCT showed a small but statistically significant improvement with = 0.64 ± 0.20 (range 0.06-0.90, P = 0.03) and DSC = 0.66 ± 0.13 (0.31-0.87, P = 0.02). Intermodal correlations were noted to decrease with an increasing fraction of lung truncation in 4DCBCT relative to 4DCT, albeit not significantly (Pearson correlation R = 0.58, P = 0.002). CONCLUSIONS: This study demonstrates that DIR based CTVIs derived from 4DCBCT can exhibit reasonable to good voxel-level agreement with CTVIs derived from 4DCT. These correlations outperform previous cross-modality comparisons between 4DCT-based ventilation and nuclear medicine. The use of 4DCBCT scans with iterative reconstruction and minimal lung truncation is recommended to ensure better reproducibility between 4DCBCT- and 4DCT-based CTVIs

Pulmonary Image Segmentation and Registration Algorithms: Towards Regional Evaluation of Obstructive Lung Disease

Pulmonary imaging, including pulmonary magnetic resonance imaging (MRI) and computed tomography (CT), provides a way to sensitively and regionally measure spatially heterogeneous lung structural-functional abnormalities. These unique imaging biomarkers offer the potential for better understanding pulmonary disease mechanisms, monitoring disease progression and response to therapy, and developing novel treatments for improved patient care. To generate these regional lung structure-function measurements and enable broad clinical applications of quantitative pulmonary MRI and CT biomarkers, as a first step, accurate, reproducible and rapid lung segmentation and registration methods are required. In this regard, we first developed a 1H MRI lung segmentation algorithm that employs complementary hyperpolarized 3He MRI functional information for improved lung segmentation. The 1H-3He MRI joint segmentation algorithm was formulated as a coupled continuous min-cut model and solved through convex relaxation, for which a dual coupled continuous max-flow model was proposed and a max-flow-based efficient numerical solver was developed. Experimental results on a clinical dataset of 25 chronic obstructive pulmonary disease (COPD) patients ranging in disease severity demonstrated that the algorithm provided rapid lung segmentation with high accuracy, reproducibility and diminished user interaction. We then developed a general 1H MRI left-right lung segmentation approach by exploring the left-to-right lung volume proportion prior. The challenging volume proportion-constrained multi-region segmentation problem was approximated through convex relaxation and equivalently represented by a max-flow model with bounded flow conservation conditions. This gave rise to a multiplier-based high performance numerical implementation based on convex optimization theories. In 20 patients with mild- to-moderate and severe asthma, the approach demonstrated high agreement with manual segmentation, excellent reproducibility and computational efficiency. Finally, we developed a CT-3He MRI deformable registration approach that coupled the complementary CT-1H MRI registration. The joint registration problem was solved by exploring optical-flow techniques, primal-dual analyses and convex optimization theories. In a diverse group of patients with asthma and COPD, the registration approach demonstrated lower target registration error than single registration and provided fast regional lung structure-function measurements that were strongly correlated with a reference method. Collectively, these lung segmentation and registration algorithms demonstrated accuracy, reproducibility and workflow efficiency that all may be clinically-acceptable. All of this is consistent with the need for broad and large-scale clinical applications of pulmonary MRI and CT

Functional imaging for assessing regional lung ventilation in preclinical and clinical research

Dynamic heterogeneity in lung ventilation is an important measure of pulmonary function and may be characteristic of early pulmonary disease. While standard indices like spirometry, body plethysmography, and blood gases have been utilized to assess lung function, they do not provide adequate information on regional ventilatory distribution nor function assessments of ventilation during the respiratory cycle. Emerging technologies such as xenon CT, volumetric CT, functional MRI and X-ray velocimetry can assess regional ventilation using non-invasive radiographic methods that may complement current methods of assessing lung function. As a supplement to current modalities of pulmonary function assessment, functional lung imaging has the potential to identify respiratory disease phenotypes with distinct natural histories. Moreover, these novel technologies may offer an optimal strategy to evaluate the effectiveness of novel therapies and therapies targeting localized small airways disease in preclinical and clinical research. In this review, we aim to discuss the features of functional lung imaging, as well as its potential application and limitations to adoption in research.</p