Autocalibration method for non-stationary CT bias correction

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.

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.

Automatic quantification of epicardial adipose tissue volume

PURPOSE

Epicardial fat is the adipose tissue between the serosal pericardial wall layer and the visceral layer. It is distributed mainly around the atrioventricular groove, atrial septum, ventricular septum and coronary arteries. Studies have shown that the density, thickness, volume and other characteristics of epicardial adipose tissue (EAT) are independently correlated with a variety of cardiovascular diseases. Given this association, the accurate determination of EAT volume is an essential aim of future research. Therefore, the purpose of this study was to establish a framework for fully automatic EAT segmentation and quantification in coronary computed tomography angiography (CCTA) scans.

METHODS

A set of 103 scans are randomly selected from our medical center. An automatic pipeline has been developed to segment and quantify the volume of EAT. First, a multi-slice deep neural network is used to simultaneously segment the pericardium in multiple adjacent slices. Then a deformable model is employed to reduce false positive and negative regions in the segmented binary pericardial images. Finally, the pericardium mask is used to define the region of interest (ROI) and the threshold method is utilized to extract the pixels ranging from -175 Hounsfield units (HU) to -15 HU for the segmentation of EAT.

RESULTS

The Dice indices of the pericardial segmentation using the proposed method with respect to the manual delineation results of two radiology experts were 97.1%  ±  0.7% and 96.9%  ±  0.6%, respectively. The inter-observer variability was also assessed, resulting in a Dice index of 97.0%  ±  0.7%. For the EAT segmentation results, the Dice indices between the proposed method and the two radiology experts were 93.4%  ±  1.5% and 93.3%  ±  1.3%, respectively, and the same measurement between the experts themselves was 93.6%  ±  1.9%. The Pearson's correlation coefficients between the EAT volumes computed from the results of the proposed method and the manual delineation by the two experts were 1.00 and 0.99 and the same coefficients between the experts was 0.99.

CONCLUSIONS

This work describes the development of a fully automatic EAT segmentation and quantification method from CCTA scans and the results compare favorably with the assessments of two independent experts. The proposed method is also packaged with a graphical user interface which can be found at https://github.com/MountainAndMorning/EATSeg.

Deep-learning-based direct inversion for material decomposition

PURPOSE

To develop a convolutional neural network (CNN) that can directly estimate material density distribution from multi-energy computed tomography (CT) images without performing conventional material decomposition.

METHODS

The proposed CNN (denoted as Incept-net) followed the general framework of encoder-decoder network, with an assumption that local image information was sufficient for modeling the nonlinear physical process of multi-energy CT. Incept-net was implemented with a customized loss function, including an in-house-designed image-gradient-correlation (IGC) regularizer to improve edge preservation. The network consisted of two types of customized multibranch modules exploiting multiscale feature representation to improve the robustness over local image noise and artifacts. Inserts with various densities of different materials [hydroxyapatite (HA), iodine, a blood-iodine mixture, and fat] were scanned using a research photon-counting detector (PCD) CT with two energy thresholds and multiple radiation dose levels. The network was trained using phantom image patches only, and tested with different-configurations of full field-of-view phantom and in vivo porcine images. Furthermore, the nominal mass densities of insert materials were used as the labels in CNN training, which potentially provided an implicit mass conservation constraint. The Incept-net performance was evaluated in terms of image noise, detail preservation, and quantitative accuracy. Its performance was also compared to common material decomposition algorithms including least-square-based material decomposition (LS-MD), total-variation regularized material decomposition (TV-MD), and U-net-based method.

RESULTS

Incept-net improved accuracy of the predicted mass density of basis materials compared with the U-net, TV-MD, and LS-MD: the mean absolute error (MAE) of iodine was 0.66, 1.0, 1.33, and 1.57 mgI/cc for Incept-net, U-net, TV-MD, and LS-MD, respectively, across all iodine-present inserts (2.0-24.0 mgI/cc). With the LS-MD as the baseline, Incept-net and U-net achieved comparable noise reduction (both around 95%), both higher than TV-MD (85%). The proposed IGC regularizer effectively helped both Incept-net and U-net to reduce image artifact. Incept-net closely conserved the total mass densities (i.e., mass conservation constraint) in porcine images, which heuristically validated the quantitative accuracy of its outputs in anatomical background. In general, Incept-net performance was less dependent on radiation dose levels than the two conventional methods; with approximately 40% less parameters, the Incept-net achieved relatively improved performance than the comparator U-net, indicating that performance gain by Incept-net was not achieved by simply increasing network learning capacity.

CONCLUSION

Incept-net demonstrated superior qualitative image appearance, quantitative accuracy, and lower noise than the conventional methods and less sensitive to dose change. Incept-net generalized and performed well with unseen image structures and different material mass densities. This study provided preliminary evidence that the proposed CNN may be used to improve the material decomposition quality in multi-energy CT.

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.

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.