No abstract

Radiomic features have a wide range of clinical applications, but variability due to image acquisition factors can affect their performance. The harmonization tool ComBat is a promising solution but is limited by inability to harmonize multimodal distributions, unknown imaging parameters, and multiple imaging parameters. In this study, we propose two methods for addressing these limitations. We propose a sequential method that allows for harmonization of radiomic features by multiple imaging parameters (Nested ComBat). We also employ a Gaussian Mixture Model (GMM)-based method (GMM ComBat) where scans are split into groupings based on the shape of the distribution used for harmonization as a batch effect and subsequent harmonization by a known imaging parameter. These two methods were evaluated on features extracted with CapTK and PyRadiomics from two public lung computed tomography datasets. We found that Nested ComBat exhibited similar performance to standard ComBat in reducing the percentage of features with statistically significant differences in distribution attributable to imaging parameters. GMM ComBat improved harmonization performance over standard ComBat (− 11%, − 10% for Lung3/CAPTK, Lung3/PyRadiomics harmonizing by kernel resolution). Features harmonized with a variant of the Nested method and the GMM split method demonstrated similar c-statistics and Kaplan–Meier curves when used in survival analyses.

IntroductionQuantitative CT (QCT) imaging-based metrics have quantified disease alterations in asthma and chronic obstructive pulmonary disease (COPD), respectively. We seek to characterise the similarity and disparity between these groups using QCT-derived airway and parenchymal metrics.MethodsAsthma and COPD subjects (former-smoker status) were selected with a criterion of post-bronchodilator FEV1 <80%. Healthy non-smokers were included as a control group. Inspiratory and expiratory QCT images of 75 asthmatic, 215 COPD and 94 healthy subjects were evaluated. We compared three segmental variables: airway circularity, normalised wall thickness and normalised hydraulic diameter, indicating heterogeneous airway shape, wall thickening and luminal narrowing, respectively. Using an image registration, we also computed six lobar variables including per cent functional small-airway disease, per cent emphysema, tissue fraction at inspiration, fractional-air-volume change, Jacobian and functional metric characterising anisotropic deformation.ResultsCompared with healthy subjects, both asthma and COPD subjects demonstrated a decreased airway circularity especially in large and upper lobar airways, and a decreased normalised hydraulic diameter in segmental airways. Besides, COPD subjects had more severe emphysema and small-airway disease, as well as smaller regional tissue fraction and lung deformation, compared with asthmatic subjects. The difference of emphysema, small-airway disease and tissue fraction between asthma and COPD was more prominent in upper and middle lobes.ConclusionsPatients with asthma and COPD, with a persistent FEV1 <80%, demonstrated similar alterations in airway geometry compared with controls, but different degrees of alterations in parenchymal regions. Density-based metrics measured at upper and middle lobes were found to be discriminant variables between patients with asthma and COPD.

Background Quantitative computed tomographic (QCT) imaging-based metrics enable to quantify smoking induced disease alterations and to identify imaging-based clusters for current smokers. We aimed to derive clinically meaningful sub-groups of former smokers using dimensional reduction and clustering methods to develop a new way of COPD phenotyping. Methods An imaging-based cluster analysis was performed for 406 former smokers with a comprehensive set of imaging metrics including 75 imaging-based metrics. They consisted of structural and functional variables at 10 segmental and 5 lobar locations. The structural variables included lung shape, branching angle, airway-circularity, airway-wall-thickness, airway diameter; the functional variables included regional ventilation, emphysema percentage, functional small airway disease percentage, Jacobian (volume change), anisotropic deformation index (directional preference in volume change), and tissue fractions at inspiration and expiration. Results We derived four distinct imaging-based clusters as possible phenotypes with the sizes of 100, 80, 141, and 85, respectively. Cluster 1 subjects were asymptomatic and showed relatively normal airway structure and lung function except airway wall thickening and moderate emphysema. Cluster 2 subjects populated with obese females showed an increase of tissue fraction at inspiration, minimal emphysema, and the lowest progression rate of emphysema. Cluster 3 subjects populated with older males showed small airway narrowing and a decreased tissue fraction at expiration, both indicating air-trapping. Cluster 4 subjects populated with lean males were likely to be severe COPD subjects showing the highest progression rate of emphysema. Conclusions QCT imaging-based metrics for former smokers allow for the derivation of statistically stable clusters associated with unique clinical characteristics. This approach helps better categorization of COPD sub-populations; suggesting possible quantitative structural and functional phenotypes. Electronic supplementary material The online version of this article (10.1186/s12931-019-1121-z) contains supplementary material, which is available to authorized users.

BackgroundClassification of COPD is usually based on the severity of airflow, which may not sensitively differentiate subpopulations. Using a multiscale imaging-based cluster analysis (MICA), we aim to identify subpopulations for current smokers with COPD.MethodsAmong the SPIROMICS subjects, we analyzed computed tomography images at total lung capacity (TLC) and residual volume (RV) of 284 current smokers. Functional variables were derived from registration of TLC and RV images, e.g. functional small airways disease (fSAD%). Structural variables were assessed at TLC images, e.g. emphysema and airway wall thickness and diameter. We employed an unsupervised method for clustering.ResultsFour clusters were identified. Cluster 1 had relatively normal airway structures; Cluster 2 had an increase of fSAD% and wall thickness; Cluster 3 exhibited a further increase of fSAD% but a decrease of wall thickness and airway diameter; Cluster 4 had a significant increase of fSAD% and emphysema. Clinically, Cluster 1 showed normal FEV1/FVC and low exacerbations. Cluster 4 showed relatively low FEV1/FVC and high exacerbations. While Cluster 2 and Cluster 3 showed similar exacerbations, Cluster 2 had the highest BMI among all clusters.ConclusionsAssociation of imaging-based clusters with existing clinical metrics suggests the sensitivity of MICA in differentiating subpopulations.Electronic supplementary materialThe online version of this article (10.1186/s12931-018-0888-7) contains supplementary material, which is available to authorized users.