Quantification of Emphysema Progression at CT Using Simultaneous Volume, Noise, and Bias Lung Density Correction.

Background CT attenuation is affected by lung volume, dosage, and scanner bias, leading to inaccurate emphysema progression measurements in multicenter studies. Purpose To develop and validate a method that simultaneously corrects volume, noise, and interscanner bias for lung density change estimation in emphysema progression at CT in a longitudinal multicenter study. Materials and Methods In this secondary analysis of the prospective Genetic Epidemiology of Chronic Obstructive Pulmonary Disease (COPDGene) study, lung function data were obtained from participants who completed baseline and 5-year follow-up visits from January 2008 to August 2017. CT emphysema progression was measured with volume-adjusted lung density (VALD) and compared with the joint volume-noise-bias-adjusted lung density (VNB-ALD). Reproducibility was studied under change of dosage protocol and scanner model with repeated acquisitions. Emphysema progression was visually scored in 102 randomly selected participants. A stratified analysis of clinical characteristics was performed that considered groups based on their combined lung density change measured by VALD and VNB-ALD. Results A total of 4954 COPDGene participants (mean age, 60 years ± 9 [SD]; 2511 male, 2443 female) were analyzed (1329 with repeated reduced-dose acquisition in the follow-up visit). Mean repeatability coefficients were 30 g/L ± 0.46 for VALD and 14 g/L ± 0.34 for VNB-ALD. VALD measurements showed no evidence of differences between nonprogressors and progressors (mean, -5.5 g/L ± 9.5 vs -8.6 g/L ± 9.6; P = .11), while VNB-ALD agreed with visual readings and showed a difference (mean, -0.67 g/L ± 4.8 vs -4.2 g/L ± 5.5; P < .001). Analysis of progression showed that VNB-ALD progressors had a greater decline in forced expiratory volume in 1 second (-42 mL per year vs -32 mL per year; Tukey-adjusted P = .002). Conclusion Simultaneously correcting volume, noise, and interscanner bias for lung density change estimation in emphysema progression at CT improved repeatability analyses and agreed with visual readings. It distinguished between progressors and nonprogressors and was associated with a greater decline in lung function metrics. Clinical trial registration no. NCT00608764 © RSNA, 2024 Supplemental material is available for this article. See also the editorial by Goo in this issue.

Quantitative CT Evaluation of Emphysema Progression over 10 Years in the COPDGene Study.

Background Long-term studies of chronic obstructive pulmonary disease (COPD) can evaluate emphysema progression. Adjustment for differences in equipment and scanning protocols of individual CT examinations have not been studied extensively. Purpose To evaluate emphysema progression in current and former smokers in the COPDGene cohort over three imaging points obtained at 5-year intervals accounting for individual CT parameters. Materials and Methods Current and former cigarette smokers enrolled between 2008 and 2011 from the COPDGene study were prospectively followed for 10 years between 2008 and 2020. Extent of emphysema as adjusted lung density (ALD) from quantitative CT was measured at baseline and at 5- and 10-year follow-up. Linear mixed models adjusted for CT technical characteristics were constructed to evaluate emphysema progression. Mean annual changes in ALD over consecutive 5-year study periods were estimated by smoking status and baseline emphysema. Results Of 8431 participants at baseline (mean age, 60 years ± 9 [SD]; 3905 female participants), 4913 were at 5-year follow-up and 1544 participants were at 10-year follow-up. There were 4134 (49%) participants who were current smokers, and 4449 (53%) participants had more than trace emphysema at baseline. Current smokers with more than trace emphysema showed the largest decline in ALD, with mean annual decreases of 1.4 g/L (95% CI: 1.2, 1.5) in the first 5 years and 0.9 g/L (95% CI: 0.7, 1.2) in the second 5 years. Accounting for CT noise, field of view, and scanner model improved model fit for estimation of emphysema progression (P < .001 by likelihood ratio test). Conclusion Evaluation at CT of emphysema progression in the COPDGene study showed that, during the span of 10 years, participants with pre-existing emphysema who continued smoking had the largest decline in ALD. Adjusting for CT equipment and protocol factors improved these longitudinal estimates. Clinical trial registration no. NCT00608764 © RSNA, 2023 Supplemental material is available for this article. See the editorial by Parraga and Kirby in this issue.

Harmonization of in-plane resolution in CT using multiple reconstructions from single acquisitions.

PURPOSE: To providea methodology that removes the spatial variability of in-plane resolution using different CT reconstructions. The methodology does not require any training, sinogram, or specific reconstruction method. METHODS: The methodology is formulated as a reconstruction problem. The desired sharp image is modeled as an unobservable variable to be estimated from an arbitrary number of observations with spatially variant resolution. The methodology comprises three steps: (1) density harmonization, which removes the density variability across reconstructions; (2) point spread function (PSF) estimation, which estimates a spatially variant PSF with arbitrary shape; (3) deconvolution, which is formulated as a regularized least squares problem. The assessment was performed with CT scans of phantoms acquired with three different Siemens scanners (Definition AS, Definition AS+, Drive). Four low-dose acquisitions reconstructed with backprojection and iterative methods were used for the resolution harmonization. A sharp, high-dose (HD) reconstruction was used as a validation reference. The different factors affecting the in-plane resolution (radial, angular, and longitudinal) were studied with regression analysis of the edge decay (between 10% and 90% of the edge spread function (ESF) amplitude). RESULTS: Results showed that the in-plane resolution improves remarkably and the spatial variability is substantially reduced without compromising the noise characteristics. The modulated transfer function (MTF) also confirmed a pronounced increase in resolution. The resolution improvement was also tested by measuring the wall thickness of tubes simulating airways. In all scanners, the resolution harmonization obtained better performance than the HD, sharp reconstruction used as a reference (up to 50 percentage points). The methodology was also evaluated in clinical scans achieving a noise reduction and a clear improvement in thin-layered structures. The estimated ESF and MTF confirmed the resolution improvement. CONCLUSION: We propose a versatile methodology to reduce the spatial variability of in-plane resolution in CT scans by leveraging different reconstructions available in clinical studies. The methodology does not require any sinogram, training, or specific reconstruction, and it is not limited to a fixed number of input images. Therefore, it can be easily adopted in multicenter studies and clinical practice. The results obtained with our resolution harmonization methodology evidence its suitability to reduce the spatially variant in-plane resolution in clinical CT scans without compromising the reconstruction's noise characteristics. We believe that the resolution increase achieved by our methodology may contribute in more accurate and reliable measurements of small structures such as vasculature, airways, and wall thickness.

SNR-enhanced diffusion MRI with structure-preserving low-rank denoising in reproducing kernel Hilbert spaces.

PURPOSE: To introduce, develop, and evaluate a novel denoising technique for diffusion MRI that leverages nonlinear redundancy in the data to boost the SNR while preserving signal information. METHODS: We exploit nonlinear redundancy of the dMRI data by means of kernel principal component analysis (KPCA), a nonlinear generalization of PCA to reproducing kernel Hilbert spaces. By mapping the signal to a high-dimensional space, a higher level of redundant information is exploited, thereby enabling better denoising than linear PCA. We implement KPCA with a Gaussian kernel, with parameters automatically selected from knowledge of the noise statistics, and validate it on realistic Monte Carlo simulations as well as with in vivo human brain submillimeter and low-resolution dMRI data. We also demonstrate KPCA denoising on multi-coil dMRI data. RESULTS: SNR improvements up to 2.7 × were obtained in real in vivo datasets denoised with KPCA, in comparison to SNR gains of up to 1.8 × using a linear PCA denoising technique called Marchenko-Pastur PCA (MPPCA). Compared to gold-standard dataset references created from averaged data, we showed that lower normalized root mean squared error was achieved with KPCA compared to MPPCA. Statistical analysis of residuals shows that anatomical information is preserved and only noise is removed. Improvements in the estimation of diffusion model parameters such as fractional anisotropy, mean diffusivity, and fiber orientation distribution functions were also demonstrated. CONCLUSION: Nonlinear redundancy of the dMRI signal can be exploited with KPCA, which allows superior noise reduction/SNR improvements than the MPPCA method, without loss of signal information.

Relationship between Emphysema Progression at CT and Mortality in Ever-Smokers: Results from the COPDGene and ECLIPSE Cohorts.

Background The relationship between emphysema progression and long-term outcomes is unclear. Purpose To determine the relationship between emphysema progression at CT and mortality among participants with emphysema. Materials and Methods In a secondary analysis of two prospective observational studies, COPDGene (clinicaltrials.gov, NCT00608764) and Evaluation of Chronic Obstructive Pulmonary Disease Longitudinally to Identify Predictive Surrogate End-points (ECLIPSE; clinicaltrials.gov, NCT00292552), emphysema was measured at CT at two points by using the volume-adjusted lung density at the 15th percentile of the lung density histogram (hereafter, lung density perc15) method. The association between emphysema progression rate and all-cause mortality was analyzed by using Cox regression adjusted for ethnicity, sex, baseline age, pack-years, and lung density, baseline and change in smoking status, forced expiratory volume in 1 second, and 6-minute walk distance. In COPDGene, respiratory mortality was analyzed by using the Fine and Gray method. Results A total of 5143 participants (2613 men [51%]; mean age, 60 years ± 9 [standard deviation]) in COPDGene and 1549 participants (973 men [63%]; mean age, 62 years ± 8) in ECLIPSE were evaluated, of which 2097 (40.8%) and 1179 (76.1%) had emphysema, respectively. Baseline imaging was performed between January 2008 and December 2010 for COPDGene and January 2006 and August 2007 for ECLIPSE. Follow-up imaging was performed after 5.5 years ± 0.6 in COPDGene and 3.0 years ± 0.2 in ECLIPSE, and mortality was assessed over the ensuing 5 years in both. For every 1 g/L per year faster rate of decline in lung density perc15, all-cause mortality increased by 8% in COPDGene (hazard ratio [HR], 1.08; 95% CI: 1.01, 1.16; P = .03) and 6% in ECLIPSE (HR, 1.06; 95% CI: 1.00, 1.13; P = .045). In COPDGene, respiratory mortality increased by 22% (HR, 1.22; 95% CI: 1.13, 1.31; P < .001) for the same increase in the rate of change in lung density perc15. Conclusion In ever-smokers with emphysema, emphysema progression at CT was associated with increased all-cause and respiratory mortality. © RSNA, 2021 Online supplemental material is available for this article. See also the editorial by Lee and Park in this issue.

Automatic estimation of aortic and mitral valve displacements in dynamic CTA with 4D graph-cuts.

The location of the mitral and aortic valves in dynamic cardiac imaging is useful for extracting functional derived parameters such as ejection fraction, valve excursions, and global longitudinal strain, and when performing anatomical structures tracking using slice following or valve intervention's planning. Completely automatic segmentation methods are still challenging tasks because of their fast movements and the different positions that prevent good visibility of the leaflets along the full cardiac cycle. In this article, we propose a processing pipeline to track the displacement of the aortic and mitral valve annuli from high-resolution cardiac four-dimensional computed tomographic angiography (4D-CTA). The proposed method is based on the dynamic separation of left ventricle, left atrium and aorta using statistical shape modeling and an energy minimization algorithm based on graph-cuts and has been evaluated on a set of 15 electrocardiography-gated 4D-CTAs. We report a mean agreement distance between manual annotations and our proposed method of 2.52±1.06 mm for the mitral annulus and 2.00±0.69 mm for the aortic valve annulus based on valve locations detected from manual anatomical landmarks. In addition, we show the effect of detecting the valvular planes on derived functional parameters (ejection fraction, global longitudinal strain, and excursions of the mitral and aortic valves).

Smaller Left Ventricle Size at Noncontrast CT Is Associated with Lower Mortality in COPDGene Participants.

Background Smokers with chronic obstructive pulmonary disease (COPD) have smaller left ventricles (LVs) due to reduced preload. Skeletal muscle wasting is also common in COPD, but less is known about its contribution to LV size. Purpose To explore the relationships between CT metrics of emphysema, venous vascular volume, and sarcopenia with the LV epicardial volume (LVEV) (myocardium and chamber) estimated from chest CT images in participants with COPD and then to determine the clinical relevance of the LVEV in multivariable models, including sex and anthropomorphic metrics. Materials and Methods The COPDGene study (ClinicalTrials.gov identifier: NCT00608764) is an ongoing prospective longitudinal observational investigation that began in 2006. LVEV, distal pulmonary venous blood volume for vessels smaller than 5 mm2 in cross section (BV5), CT emphysema, and pectoralis muscle area were retrospectively extracted from 3318 nongated, unenhanced COPDGene CT scans. Multivariable linear and Cox regression models were used to explore the association between emphysema, venous BV5, pectoralis muscle area, and LVEV as well as the association of LVEV with health status using the St George's Respiratory Questionnaire, 6-minute walk distance, and all-cause mortality. Results The median age of the cohort was 64 years (interquartile range, 57-70 years). Of the 2423 participants, 1806 were men and 617 were African American. The median LVEV between Global Initiative for Chronic Obstructive Lung Disease (GOLD) 1 and GOLD 4 COPD was reduced by 13.9% in women and 17.7% in men (P < .001 for both). In fully adjusted models, higher emphysema percentage (ß = -4.2; 95% confidence interval [CI]: -5.0, -3.4; P < .001), venous BV5 (ß = 7.0; 95% CI: 5.7, 8.2; P < .001), and pectoralis muscle area (ß = 2.7; 95% CI: 1.2, 4.1; P < .001) were independently associated with reduced LVEV. Reductions in LVEV were associated with improved health status (ß = 0.3; 95% CI: 0.1, 0.4) and 6-minute walk distance (ß = -12.2; 95% CI: -15.2, -9.3). These effects were greater in women than in men. The effect of reduced LVEV on mortality (hazard ratio: 1.07; 95% CI: 1.05, 1.09) did not vary by sex. Conclusion In women more than men with chronic obstructive pulmonary disease, a reduction in the estimated left ventricle epicardial volume correlated with a loss of pulmonary venous vasculature, greater pectoralis muscle sarcopenia, and lower all-cause mortality. © RSNA, 2020 Online supplemental material is available for this article.

Harmonization of chest CT scans for different doses and reconstruction methods.

PURPOSE: To develop and validate a computed tomography (CT) harmonization technique by combining noise-stabilization and autocalibration methodologies to provide reliable densitometry measurements in heterogeneous acquisition protocols. METHODS: We propose to reduce the effects of spatially variant noise such as nonuniform patterns of noise and biases. The method combines the statistical characterization of the signal-to-noise relationship in the CT image intensities, which allows us to estimate both the signal and spatially variant variance of noise, with an autocalibration technique that reduces the nonuniform biases caused by noise and reconstruction techniques. The method is firstly validated with anthropomorphic synthetic images that simulate CT acquisitions with variable scanning parameters: different dosage, nonhomogeneous variance of noise, and various reconstruction methods. We finally evaluate these effects and the ability of our method to provide consistent densitometric measurements in a cohort of clinical chest CT scans from two vendors (Siemens, n = 54 subjects; and GE, n = 50 subjects) acquired with several reconstruction algorithms (filtered back-projection and iterative reconstructions) with high-dose and low-dose protocols. RESULTS: The harmonization reduces the effect of nonhomogeneous noise without compromising the resolution of the images (25% RMSE reduction in both clinical datasets). An analysis through hierarchical linear models showed that the average biases induced by differences in dosage and reconstruction methods are also reduced up to 74.20%, enabling comparable results between high-dose and low-dose reconstructions. We also assessed the statistical similarity between acquisitions obtaining increases of up to 30% points and showing that the low-dose vs high-dose comparisons of harmonized data obtain similar and even higher similarity than the observed for high-dose vs high-dose comparisons of nonharmonized data. CONCLUSION: The proposed harmonization technique allows to compare measures of low-dose with high-dose acquisitions without using a specific reconstruction as a reference. Since the harmonization does not require a precalibration with a phantom, it can be applied to retrospective studies. This approach might be suitable for multicenter trials for which a reference reconstruction is not feasible or hard to define due to differences in vendors, models, and reconstruction techniques.

Arterial Vascular Pruning, Right Ventricular Size, and Clinical Outcomes in Chronic Obstructive Pulmonary Disease. A Longitudinal Observational Study.

Rationale: Cor pulmonale (right ventricular [RV] dilation) and cor pulmonale parvus (RV shrinkage) are both described in chronic obstructive pulmonary disease (COPD). The identification of emphysema as a shared risk factor suggests that additional disease characterization is needed to understand these widely divergent cardiac processes.Objectives: To explore the relationship between computed tomography measures of emphysema and distal pulmonary arterial morphology with RV volume, and their association with exercise capacity and mortality in ever-smokers with COPD enrolled in the COPDGene Study.Methods: Epicardial (myocardium and chamber) RV volume (RVEV), distal pulmonary arterial blood vessel volume (arterial BV5: vessels <5 mm2 in cross-section), and objective measures of emphysema were extracted from 3,506 COPDGene computed tomography scans. Multivariable linear and Cox regression models and the log-rank test were used to explore the association between emphysema, arterial BV5, and RVEV with exercise capacity (6-min-walk distance) and all-cause mortality.Measurements and Main Results: The RVEV was approximately 10% smaller in Global Initiative for Chronic Obstructive Lung Disease stage 4 versus stage 1 COPD (P < 0.0001). In multivariable modeling, a 10-ml decrease in arterial BV5 (pruning) was associated with a 1-ml increase in RVEV. For a given amount of emphysema, relative preservation of the arterial BV5 was associated with a smaller RVEV. An increased RVEV was associated with reduced 6-minute-walk distance and in those with arterial pruning an increased mortality.Conclusions: Pulmonary arterial pruning is associated with clinically significant increases in RV volume in smokers with COPD and is related to exercise capacity and mortality in COPD.Clinical trial registered with www.clinicaltrials.gov (NCT00608764).

Scalar diffusion-MRI measures invariant to acquisition parameters: A first step towards imaging biomarkers.

An imaging biomarker is a biologic feature in an image that is relevant to a patient's diagnosis or prognosis. In order to qualify as a biomarker, a measure must be robust and reproducible. However, the usual scalar measures derived from diffusion tensor imaging are known to be highly dependent on the variation of the acquisition parameters, which prevents their possible use as biomarkers. In this work, we propose a new set of quantitative measures based on diffusion magnetic resonance imaging from single-shell acquisitions that are designed to be robust to the variations of several acquisition parameters (number of gradient directions, b-value and SNR) while keeping a high discrimination power on differences in the diffusion characteristics of the tissue. These new scalar measures are analytically obtained from a generic diffusion function that does not require the calculation of a diffusion tensor. This way, on one hand, we avoid the use of a specific diffusion model and, on the other hand, we make easier the statistical characterization of the measures. Accordingly, the analysis of the measures bias is carried out and it is used to minimize their dependency with respect to the acquisition noise for different SNRs. The robustness and discrimination power of the measures are tested for different number of gradients, b-values and SNRs using a realistic phantom and three real datasets: (1) 13 control subjects and different acquisition parameters; (2) a public data set from a single subject acquired using multiple shells and (3) 32 schizophrenia patients and 32 age and sex-matched healthy controls with a varying number of gradient directions. The proposed quantitative measures exhibit low variability to the changes of the acquisition parameters, while at the same time they preserve a discrimination power that is able to detect significant changes in the anisotropy of the diffusion.

NOVIFAST: A Fast Algorithm for Accurate and Precise VFA MRI Mapping.

In quantitative magnetic resonance mapping, the variable flip angle (VFA) steady state spoiled gradient recalled echo (SPGR) imaging technique is popular as it provides a series of high resolution weighted images in a clinically feasible time. Fast, linear methods that estimate maps from these weighted images have been proposed, such as DESPOT1 and iterative re-weighted linear least squares. More accurate, non-linear least squares (NLLS) estimators are in play, but these are generally much slower and require careful initialization. In this paper, we present NOVIFAST, a novel NLLS-based algorithm specifically tailored to VFA SPGR mapping. By exploiting the particular structure of the SPGR model, a computationally efficient, yet accurate and precise map estimator is derived. Simulation and in vivo human brain experiments demonstrate a twenty-fold speed gain of NOVIFAST compared with conventional gradient-based NLLS estimators while maintaining a high precision and accuracy. Moreover, NOVIFAST is eight times faster than the efficient implementations of the variable projection (VARPRO) method. Furthermore, NOVIFAST is shown to be robust against initialization.

Interstitial Features at Chest CT Enhance the Deleterious Effects of Emphysema in the COPDGene Cohort.

Purpose To determine if interstitial features at chest CT enhance the effect of emphysema on clinical disease severity in smokers without clinical pulmonary fibrosis. Materials and Methods In this retrospective cohort study, an objective CT analysis tool was used to measure interstitial features (reticular changes, honeycombing, centrilobular nodules, linear scar, nodular changes, subpleural lines, and ground-glass opacities) and emphysema in 8266 participants in a study of chronic obstructive pulmonary disease (COPD) called COPDGene (recruited between October 2006 and January 2011). Additive differences in patients with emphysema with interstitial features and in those without interstitial features were analyzed by using t tests, multivariable linear regression, and Kaplan-Meier analysis. Multivariable linear and Cox regression were used to determine if interstitial features modified the effect of continuously measured emphysema on clinical measures of disease severity and mortality. Results Compared with individuals with emphysema alone, those with emphysema and interstitial features had a higher percentage predicted forced expiratory volume in 1 second (absolute difference, 6.4%; P < .001), a lower percentage predicted diffusing capacity of lung for carbon monoxide (DLCO) (absolute difference, 7.4%; P = .034), a 0.019 higher right ventricular-to-left ventricular (RVLV) volume ratio (P = .029), a 43.2-m shorter 6-minute walk distance (6MWD) (P < .001), a 5.9-point higher St George's Respiratory Questionnaire (SGRQ) score (P < .001), and 82% higher mortality (P < .001). In addition, interstitial features modified the effect of emphysema on percentage predicted DLCO, RVLV volume ratio, 6WMD, SGRQ score, and mortality (P for interaction < .05 for all). Conclusion In smokers, the combined presence of interstitial features and emphysema was associated with worse clinical disease severity and higher mortality than was emphysema alone. In addition, interstitial features enhanced the deleterious effects of emphysema on clinical disease severity and mortality.

Disease Staging and Prognosis in Smokers Using Deep Learning in Chest Computed Tomography.

RATIONALE: Deep learning is a powerful tool that may allow for improved outcome prediction. OBJECTIVES: To determine if deep learning, specifically convolutional neural network (CNN) analysis, could detect and stage chronic obstructive pulmonary disease (COPD) and predict acute respiratory disease (ARD) events and mortality in smokers. METHODS: A CNN was trained using computed tomography scans from 7,983 COPDGene participants and evaluated using 1,000 nonoverlapping COPDGene participants and 1,672 ECLIPSE participants. Logistic regression (C statistic and the Hosmer-Lemeshow test) was used to assess COPD diagnosis and ARD prediction. Cox regression (C index and the Greenwood-Nam-D'Agnostino test) was used to assess mortality. MEASUREMENTS AND MAIN RESULTS: In COPDGene, the C statistic for the detection of COPD was 0.856. A total of 51.1% of participants in COPDGene were accurately staged and 74.95% were within one stage. In ECLIPSE, 29.4% were accurately staged and 74.6% were within one stage. In COPDGene and ECLIPSE, the C statistics for ARD events were 0.64 and 0.55, respectively, and the Hosmer-Lemeshow P values were 0.502 and 0.380, respectively, suggesting no evidence of poor calibration. In COPDGene and ECLIPSE, CNN predicted mortality with fair discrimination (C indices, 0.72 and 0.60, respectively), and without evidence of poor calibration (Greenwood-Nam-D'Agnostino P values, 0.307 and 0.331, respectively). CONCLUSIONS: A deep-learning approach that uses only computed tomography imaging data can identify those smokers who have COPD and predict who are most likely to have ARD events and those with the highest mortality. At a population level CNN analysis may be a powerful tool for risk assessment.

Statistical Framework for the Definition of Emphysema in CT Scans: Beyond Density Mask.

Lung parenchyma destruction (emphysema) is a major factor in the description of Chronic Obstructive Pulmonary Disease (COPD) and its prognosis. It is defined as an abnormal enlargement of air spaces distal to the terminal bronchioles and the destruction of alveolar walls. In CT imaging, the presence of emphysema is observed by a local decrease of the lung density and the diagnose is usually set as more than 5% of the lung below -950 HU, the so-called emphysema density mask. There is still debate, however, about the definition of this percentage and many researchers set it depending on the population under study. Additionally, the -950 HU threshold may vary depending on factors as the slice thickness or the respiratory phase of the acquisition. In this paper we propose (1) a statistical framework that provides an automatic definition of the density threshold based on the statistical characterization of air and lung parenchyma; (2) the definition of a statistical test for emphysema detection that accounts for the CT noise characteristics. Results show that this novel statistical framework improves the quantification of emphysema against a visual reference and improves the association of emphysema with the pulmonary function tests.

A CT Scan Harmonization Technique to Detect Emphysema and Small Airway Diseases.

Recent studies have suggested the central role of small airway destruction in the pathogenesis of COPD leading to further parenchymal destruction. This evidence has sparked the interest in in-vivo assessment of small airway disease overall at the early onset of the disease. The parametric response mapping (PRM) technique has been proposed to distinguish gas trapping due to small airway disease from low attenuation areas due to emphysema. Despite its success, the PRM technique shows some limitations that are precluding the interpretation of its results. The density value used to assess gas trapping highly depends on acquisition parameters, such as dose and reconstruction kernel, and changes in body size, that introduce inhomogeneous photon absorption patterns. In particular, many studies using PRM employ inspiratory and expiratory images that are obtained at different dose levels. Emphysema impact in early disease may be confounded with the gas trapping due to the noise introduced by differences in the acquisition during the PRM. In this work, we propose a CT harmonization technique to remove the nuisance factors to distinguish between small airway disease and emphysema. Our results show that the measurements based on CT harmonization provide an increase in the detection of both emphysema and airway disease, resulting in a statistically significant impact of both components and a better association with lung function measures.

Autocalibration method for non-stationary CT bias correction.

Computed tomography (CT) is a widely used imaging modality for screening and diagnosis. However, the deleterious effects of radiation exposure inherent in CT imaging require the development of image reconstruction methods which can reduce exposure levels. The development of iterative reconstruction techniques is now enabling the acquisition of low-dose CT images whose quality is comparable to that of CT images acquired with much higher radiation dosages. However, the characterization and calibration of the CT signal due to changes in dosage and reconstruction approaches is crucial to provide clinically relevant data. Although CT scanners are calibrated as part of the imaging workflow, the calibration is limited to select global reference values and does not consider other inherent factors of the acquisition that depend on the subject scanned (e.g. photon starvation, partial volume effect, beam hardening) and result in a non-stationary noise response. In this work, we analyze the effect of reconstruction biases caused by non-stationary noise and propose an autocalibration methodology to compensate it. Our contributions are: 1) the derivation of a functional relationship between observed bias and non-stationary noise, 2) a robust and accurate method to estimate the local variance, 3) an autocalibration methodology that does not necessarily rely on a calibration phantom, attenuates the bias caused by noise and removes the systematic bias observed in devices from different vendors. The validation of the proposed methodology was performed with a physical phantom and clinical CT scans acquired with different configurations (kernels, doses, algorithms including iterative reconstruction). The results confirmed the suitability of the proposed methods for removing the intra-device and inter-device reconstruction biases.

Statistical characterization of noise for spatial standardization of CT scans: Enabling comparison with multiple kernels and doses.

Computerized tomography (CT) is a widely adopted modality for analyzing directly or indirectly functional, biological and morphological processes by means of the image characteristics. However, the potential utilization of the information obtained from CT images is often limited when considering the analysis of quantitative information involving different devices, acquisition protocols or reconstruction algorithms. Although CT scanners are calibrated as a part of the imaging workflow, the calibration is circumscribed to global reference values and does not circumvent problems that are inherent to the imaging modality. One of them is the lack of noise stationarity, which makes quantitative biomarkers extracted from the images less robust and stable. Some methodologies have been proposed for the assessment of non-stationary noise in reconstructed CT scans. However, those methods focused on the non-stationarity only due to the reconstruction geometry and are mainly based on the propagation of the variance of noise throughout the whole reconstruction process. Additionally, the philosophy followed in the state-of-the-art methods is based on the reduction of noise, but not in the standardization of it. This means that, even if the noise is reduced, the statistics of the signal remain non-stationary, which is insufficient to enable comparisons between different acquisitions with different statistical characteristics. In this work, we propose a statistical characterization of noise in reconstructed CT scans that leads to a versatile statistical model that effectively characterizes different doses, reconstruction kernels, and devices. The statistical model is generalized to deal with the partial volume effect via a localized mixture model that also describes the non-stationarity of noise. Finally, we propose a stabilization scheme to achieve stationary variance. The validation of the proposed methodology was performed with a physical phantom and clinical CT scans acquired with different configurations (kernels, doses, algorithms including iterative reconstruction). The results confirmed its suitability to enable comparisons with different doses, and acquisition protocols.