A radiomics approach for lung nodule detection in thoracic CT images based on the dynamic patterns of morphological variation.

OBJECTIVES: To propose and evaluate a set of radiomic features, called morphological dynamics features, for pulmonary nodule detection, which were rooted in the dynamic patterns of morphological variation and needless precise lesion segmentation. MATERIALS AND METHODS: Two datasets were involved, namely, university hospital (UH) and LIDC datasets, comprising 72 CT scans (360 nodules) and 888 CT scans (2230 nodules), respectively. Each nodule was annotated by multiple radiologists. Denoted the category of nodules identified by at least k radiologists as ALk. A nodule detection algorithm, called CAD-MD algorithm, was proposed based on the morphological dynamics radiomic features, characterizing a lesion by ten sets of the same features with different values extracted from ten different thresholding results. Each nodule candidate was classified by a two-level classifier, including ten decision trees and a random forest, respectively. The CAD-MD algorithm was compared with a deep learning approach, the N-Net, using the UH dataset. RESULTS: On the AL1 and AL2 of the UH dataset, the AUC of the AFROC curves were 0.777 and 0.851 for the CAD-MD algorithm and 0.478 and 0.472 for the N-Net, respectively. The CAD-MD algorithm achieved the sensitivities of 84.4% and 91.4% with 2.98 and 3.69 FPs/scan and the N-Net 74.4% and 80.7% with 3.90 and 4.49 FPs/scan, respectively. On the LIDC dataset, the CAD-MD algorithm attained the sensitivities of 87.6%, 89.2%, 92.2%, and 95.0% with 4 FPs/scan for AL1-AL4, respectively. CONCLUSION: The morphological dynamics radiomic features might serve as an effective set of radiomic features for lung nodule detection. KEY POINTS: • Texture features varied with such CT system settings as reconstruction kernels of CT images, CT scanner models, and parameter settings, and so on. • Shape and first-order statistics were shown to be the most robust features against variation in CT imaging parameters. • The morphological dynamics radiomic features, which mainly characterized the dynamic patterns of morphological variation, were shown to be effective for lung nodule detection.

Longitudinal registration of thoracic CT images with radiation-induced lung diseases: A divide-and-conquer approach based on component structure wise registration using coherent point drift.

BACKGROUND AND OBJECTIVE: Registration of pulmonary computed tomography (CT) images with radiation-induced lung diseases (RILD) was essential to investigate the voxel-wise relationship between the formation of RILD and the radiation dose received by different tissues. Although various approaches had been developed for the registration of lung CTs, their performances remained clinically unsatisfactory for registration of lung CT images with RILD. The main difficulties arose from the longitudinal change in lung parenchyma, including RILD and volumetric change of lung cancers, after radiation therapy, leading to inaccurate registration and artifacts caused by erroneous matching of the RILD tissues. METHODS: To overcome the influence of the parenchymal changes, a divide-and-conquer approach rooted in the coherent point drift (CPD) paradigm was proposed. The proposed method was based on two kernel ideas. One was the idea of component structure wise registration. Specifically, the proposed method relaxed the intrinsic assumption of equal isotropic covariances in CPD by decomposing a lung and its surrounding tissues into component structures and independently registering the component structures pairwise by CPD. The other was the idea of defining a vascular subtree centered at a matched branch point as a component structure. This idea could not only provide a sufficient number of matched feature points within a parenchyma, but avoid being corrupted by the false feature points resided in the RILD tissues due to globally and indiscriminately sampling using mathematical operators. The overall deformation model was built by using the Thin Plate Spline based on all matched points. RESULTS: This study recruited 30 pairs of lung CT images with RILD, 15 of which were used for internal validation (leave-one-out cross-validation) and the other 15 for external validation. The experimental results showed that the proposed algorithm achieved a mean and a mean of maximum 1 % of average surface distances <2 and 8 mm, respectively, and a mean and a maximum target registration error <2 mm and 5 mm on both internal and external validation datasets. The paired two-sample t-tests corroborated that the proposed algorithm outperformed a recent method, the Stavropoulou's method, on the external validation dataset (p < 0.05). CONCLUSIONS: The proposed algorithm effectively reduced the influence of parenchymal changes, resulting in a reasonably accurate and artifact-free registration.

Automated marker-free longitudinal infrared breast image registration by GA-PSO.

The automated marker-free longitudinal Infrared (IR) breast image registration overcomes several challenges like no anatomic fiducial markers on the body surface, blurry boundaries, heat pattern variation by environmental and physiological factors, nonrigid deformation, etc, has the ability of quantitative pixel-wise analysis with the heat energy and patterns change in a time course study. To achieve the goal, scale-invariant feature transform, Harris corner, and Hessian matrix were employed to generate the feature points as anatomic fiducial markers, and hybrid genetic algorithm and particle swarm optimization minimizing the matching errors was used to find the appropriate corresponding pairs between the 1st IR image and thenth IR image. Moreover, the mechanism of the IR spectrogram hardware system has a high level of reproducibility. The performance of the proposed longitudinal image registration system was evaluated by the simulated experiments and the clinical trial. In the simulated experiments, the mean difference of our system is 1.64 mm, which increases 57.58% accuracy than manual determination and makes a 17.4% improvement than the previous study. In the clinical trial, 80 patients were captured several times of IR breast images during chemotherapy. Most of them were well aligned in the spatiotemporal domain. In the few cases with evident heat pattern dissipation and spatial deviation, it still provided a reliable comparison of vascular variation. Therefore, the proposed system is accurate and robust, which could be considered as a reliable tool for longitudinal approaches to breast cancer diagnosis.

Predictive Modelling of Lung Function using Emphysematous Density Distribution.

Target lung tissue selection remains a challenging task to perform for treating severe emphysema with lung volume reduction (LVR). In order to target the treatment candidate, the percentage of low attenuation volume (LAV%) representing the proportion of emphysema volume to whole lung volume is measured using computed tomography (CT) images. Although LAV% have shown to have a correlation with lung function in patients with chronic obstructive pulmonary disease (COPD), similar measurements of LAV% in whole lung or lobes may have large variations in lung function due to emphysema heterogeneity. The functional information of regional emphysema destruction is required for supporting the choice of optimal target. The purpose of this study is to develop an emphysema heterogeneity descriptor for the three-dimensional emphysematous bullae according to the size variations of emphysematous density (ED) and their spatial distribution. The second purpose is to derive a predictive model of airflow limitation based on the regional emphysema heterogeneity. Deriving the bullous representation and grouping them into four scales in the upper and lower lobes, a predictive model is computed using the linear model fitting to estimate the severity of lung function. A total of 99 subjects, 87 patients with mild to very severe COPD (Global Initiative for Chronic Obstructive Lung Disease (GOLD) stage I~IV) and 12 control participants with normal lung functions (forced expiratory volume in one second (FEV1)/forced vital capacity (FVC) > 0.7) were evaluated. The final model was trained with stratified cross-validation on randomly selected 75% of the dataset (n = 76) and tested on the remaining dataset (n = 23). The dispersed cases of LAV% inconsistent with their lung function outcome were evaluated, and the correlation study suggests that comparing to LAV of larger bullae, the widely spread smaller bullae with equivalent LAV has a larger impact on lung function. The testing dataset has the correlation of r = -0.76 (p < 0.01) between the whole lung LAV% and FEV1/FVC, whereas using two ED % of scales and location-dependent variables to predict the emphysema-associated FEV1/FVC, the results shows their correlation of 0.82 (p < 0.001) with clinical FEV1/FVC.