Search

Accurate assessment of pulmonary emphysema is crucial to assess disease severity and subtype, to monitor disease progression, and to predict lung cancer risk. However, visual assessment is time-consuming and subject to substantial inter-rater variability while standard densitometry approaches to quantify emphysema remain inferior to visual scoring. We explore if machine learning methods that learn from a large dataset of visually assessed CT scans can provide accurate estimates of emphysema extent and if methods that learn from emphysema extent scoring can outperform algorithms that learn only from emphysema presence scoring. Four Multiple Instance Learning classifiers, trained on emphysema presence labels, and five Learning with Label Proportions classifiers, trained on emphysema extent labels, are compared. Performance is evaluated on 600 low-dose CT scans from the Danish Lung Cancer Screening Trial and we find that learning from emphysema presence labels, which are much easier to obtain, gives equally good performance to learning from emphysema extent labels. The best performing Multiple Instance Learning and Learning with Label Proportions classifiers, achieve intra-class correlation coefficients around 0.90 and average overall agreement with raters of 78% and 79% compared to an inter-rater agreement of 83.

Background: Accurate estimation of the malignancy risk of pulmonary nodules at chest CT is crucial for optimizing management in lung cancer screening. Purpose: To develop and validate a deep learning (DL) algorithm for malignancy risk estimation of pulmonary nodules detected at screening CT. Materials and Methods: In this retrospective study, the DL algorithm was developed with 16 077 nodules (1249 malignant) collected between 2002 and 2004 from the National Lung Screening Trial. External validation was performed in the following three cohorts collected between 2004 and 2010 from the Danish Lung Cancer Screening Trial: a full cohort containing all 883 nodules (65 malignant) and two cancer-enriched cohorts with size matching (175 nodules, 59 malignant) and without size matching (177 nodules, 59 malignant) of benign nodules selected at random. Algorithm performance was measured by using the area under the receiver operating characteristic curve (AUC) and compared with that of the Pan-Canadian Early Detection of Lung Cancer (PanCan) model in the full cohort and a group of 11 clinicians composed of four thoracic radiologists, five radiology residents, and two pulmonologists in the cancer-enriched cohorts. Results: The DL algorithm significantly outperformed the PanCan model in the full cohort (AUC, 0.93 [95% CI: 0.89, 0.96] vs 0.90 [95% CI: 0.86, 0.93]; P = .046). The algorithm performed comparably to thoracic radiologists in cancer-enriched cohorts with both random benign nodules (AUC, 0.96 [95% CI: 0.93, 0.99] vs 0.90 [95% CI: 0.81, 0.98]; P = .11) and size-matched benign nodules (AUC, 0.86 [95% CI: 0.80, 0.91] vs 0.82 [95% CI: 0.74, 0.89]; P = .26). Conclusion: The deep learning algorithm showed excellent performance, comparable to thoracic radiologists, for malignancy risk estimation of pulmonary nodules detected at screening CT. This algorithm has the potential to provide reliable and reproducible malignancy risk scores for clinicians, which may help optimize management in lung cancer screening.

An observer study with seven radiologists showed that using a dedicated CT lung screening viewer (with computer-aided detection, a segmentation algorithm, and automatic Lung Imaging Reporting and Data System [Lung-RADS] categorization) improved interobserver agreement for Lung-RADS categorization while substantially reducing the reading time. Purpose. To compare the inter- and intraobserver agreement and reading times achieved when assigning Lung Imaging Reporting and Data System (Lung-RADS) categories to baseline and follow-up lung cancer screening studies by using a dedicated CT lung screening viewer with integrated nodule detection and volumetric support with those achieved by using a standard picture archiving and communication system (PACS)-like viewer. Materials and Methods. Data were obtained from the National Lung Screening Trial (NLST). By using data recorded by NLST radiologists, scans were assigned to Lung-RADS categories. For each Lung-RADS category (1 or 2, 3, 4A, and 4B), 40 CT scans (20 baseline scans and 20 follow-up scans) were randomly selected for 160 participants (median age, 61 years; interquartile range, 58-66 years; 61 women) in total. Seven blinded observers independently read all CT scans twice in a randomized order with a 2-week washout period: once by using the standard PACS-like viewer and once by using the dedicated viewer. Observers were asked to assign a Lung-RADS category to each scan and indicate the risk-dominant nodule. Inter- and intraobserver agreement was analyzed by using Fleiss k values and Cohen weighted k values, respectively. Reading times were compared by using a Wilcoxon signed rank test. Results. The interobserver agreement was moderate for the standard viewer and substantial for the dedicated viewer, with Fleiss k values of 0.58 (95% CI: 0.55, 0.60) and 0.66 (95% CI: 0.64, 0.68), respectively. The intraobserver agreement was substantial, with a mean Cohen weighted k value of 0.67. The median reading time was significantly reduced from 160 seconds with the standard viewer to 86 seconds with the dedicated viewer (P < .001). Conclusion. Lung-RADS interobserver agreement increased from moderate to substantial when using the dedicated CT lung screening viewer. The median reading time was substantially reduced when scans were read by using the dedicated CT lung screening viewer.

Volumetric imaging is an essential diagnostic tool for medical practitioners. The use of popular techniques such as convolutional neural networks (CNN) for analysis of volumetric images is constrained by the availability of detailed (with local annotations) training data and GPU memory. In this paper, the volumetric image classification problem is posed as a multiinstance classification problem and a novel method is proposed to adaptively select positive instances from positive bags during the training phase. This method uses the extreme value theory to model the feature distribution of the images without a pathology and use it to identify positive instances of an imaged pathology. The experimental results, on three separate image classification tasks (i.e. classify retinal OCT images according to the presence or absence of fluid build-ups, emphysema detection in pulmonary 3D-CT images and detection of cancerous regions in 2D histopathology images) show that the proposed method produces classifiers that have similar performance to fully supervised methods and achieves the state of the art performance in all examined test cases.

Contains fulltext : 181941.pdf (Publisher’s version ) (Open Access) OBJECTIVES: To compare the PanCan model, Lung-RADS and the 1.2016 National Comprehensive Cancer Network (NCCN) guidelines for discriminating malignant from benign pulmonary nodules on baseline screening CT scans and the impact diameter measurement methods have on performances. METHODS: From the Danish Lung Cancer Screening Trial database, 64 CTs with malignant nodules and 549 baseline CTs with benign nodules were included. Performance of the systems was evaluated applying the system's original diameter definitions: D(longest-C) (PanCan), D(meanAxial) (NCCN), both obtained from axial sections, and D(mean3D) (Lung-RADS). Subsequently all diameter definitions were applied uniformly to all systems. Areas under the ROC curves (AUC) were used to evaluate risk discrimination. RESULTS: PanCan performed superiorly to Lung-RADS and NCCN (AUC 0.874 vs. 0.813, p = 0.003; 0.874 vs. 0.836, p = 0.010), using the original diameter specifications. When uniformly applying D(longest-C), D(mean3D) and D(meanAxial), PanCan remained superior to Lung-RADS (p < 0.001 - p = 0.001) and NCCN (p < 0.001 - p = 0.016). Diameter definition significantly influenced NCCN's performance with D(longest-C) being the worst (D(longest-C) vs. D(mean3D), p = 0.005; D(longest-C) vs. D(meanAxial), p = 0.016). CONCLUSIONS: Without follow-up information, the PanCan model performs significantly superiorly to Lung-RADS and the 1.2016 NCCN guidelines for discriminating benign from malignant nodules. The NCCN guidelines are most sensitive to nodule size definition. KEY POINTS: * PanCan model outperforms Lung-RADS and 1.2016 NCCN guidelines in identifying malignant pulmonary nodules. * Nodule size definition had no significant impact on Lung-RADS and PanCan model. * 1.2016 NCCN guidelines were significantly superior when using mean diameter to longest diameter. * Longest diameter achieved lowest performance for all models. * Mean diameter performed equivalently when derived from axial sections and from volumetry.

Background: Chest x-ray is central in screening and diagnosis of tuberculosis. However, sputum culture remains gold standard for diagnosis. Aim: To establish the rate of normal chest x-rays in tuberculosis patients found by spot sputum culture screening, and compare them to a group identified through passive case finding. Method: Chest x-rays from 39 culture-positive patients, identified by spot sputum culture screening in Copenhagen from 2012 to 2014, were included in the study (spot sputum culture group(SSC)). 39 normal chest x-rays from persons screened by mobile x-ray, and 39 chest x-rays from tuberculosis-patients identified through passive case finding(PCF) were anonymised and randomised. Two respiratory physicians and two radiologists assessed the chest x-rays. Results: The normal chest x-ray rate was higher in the non-tuberculosis control group (median = 32 (82.1%), range = 74.4% – 100%), compared to the SSC group (median = 7 (17.9%), range = 10.3% – 33.3%), and the PCF controls (median = 3(7.7%), range = 2.6% – 15.4%). In the SSC group 14 (35.9%) were categorized as normal by at least one study participant. Conclusion: A substantial minority of patients diagnosed with tuberculosis by spot sputum culture screening, and through passive case finding would not have been identified with chest x-ray alone, highlighting that a normal chest x-ray does not exclude pulmonary tuberculosis. Keywords: Tuberculosis, Chest x-ray changes, Passive case finding, Active case finding, Chest x-ray assessment, Normal chest x-ray

Lung-RADS represents a categorical system published by the American College of Radiology to standardise management in lung cancer screening. The purpose of the study was to quantify how well readers agree in assigning Lung-RADS categories to screening CTs; secondary goals were to assess causes of disagreement and evaluate its impact on patient management. For the observer study, 80 baseline and 80 follow-up scans were randomly selected from the NLST trial covering all Lung-RADS categories in an equal distribution. Agreement of seven observers was analysed using Cohen’s kappa statistics. Discrepancies were correlated with patient management, test performance and diagnosis of malignancy within the scan year. Pairwise interobserver agreement was substantial (mean kappa 0.67, 95% CI 0.58–0.77). Lung-RADS category disagreement was seen in approximately one-third (29%, 971) of 3360 reading pairs, resulting in different patient management in 8% (278/3360). Out of the 91 reading pairs that referred to scans with a tumour diagnosis within 1 year, discrepancies in only two would have resulted in a substantial management change. Assignment of lung cancer screening CT scans to Lung-RADS categories achieves substantial interobserver agreement. Impact of disagreement on categorisation of malignant nodules was low. • Lung-RADS categorisation of low-dose lung screening CTs achieved substantial interobserver agreement. • Major cause for disagreement was assigning a different nodule as risk-dominant. • Disagreement led to a different follow-up time in 8% of reading pairs.

As of April 2015, participants in the Danish Lung Cancer Screening Trial had been followed for at least 5 years since their last screening.Mortality, causes of death, and lung cancer findings are reported to explore the effect of computed tomography (CT) screening.A total of 4,104 participants aged 50-70 years at the time of inclusion and with a minimum 20 pack-years of smoking were randomized to have five annual low-dose CT scans (study group) or no screening (control group).Follow-up information regarding date and cause of death, lung cancer diagnosis, cancer stage, and histology was obtained from national registries. No differences between the two groups in lung cancer mortality (hazard ratio, 1.03; 95% confidence interval, 0.66-1.6; P = 0.888) or all-cause mortality (hazard ratio, 1.02; 95% confidence interval, 0.82-1.27; P = 0.867) were observed. More cancers were found in the screening group than in the no-screening group (100 vs. 53, respectively; P0.001), particularly adenocarcinomas (58 vs. 18, respectively; P0.001). More early-stage cancers (stages I and II, 54 vs. 10, respectively; P0.001) and stage IIIa cancers (15 vs. 3, respectively; P = 0.009) were found in the screening group than in the control group. Stage IV cancers were nonsignificantly more frequent in the control group than in the screening group (32 vs. 23, respectively; P = 0.278). For the highest-stage cancers (T4N3M1, 21 vs. 8, respectively; P = 0.025), this difference was statistically significant, indicating an absolute stage shift. Older participants, those with chronic obstructive pulmonary disease, and those with more than 35 pack-years of smoking had a significantly increased risk of death due to lung cancer, with nonsignificantly fewer deaths in the screening group.No statistically significant effects of CT screening on lung cancer mortality were found, but the results of post hoc high-risk subgroup analyses showed nonsignificant trends that seem to be in good agreement with the results of the National Lung Screening Trial. Clinical trial registered with www.clinicaltrials.gov (NCT00496977).

Emphysema is part of chronic obstructive pulmonary disease, a leading cause of mortality worldwide. Visual assessment of emphysema presence is useful for identifying subjects at risk and for research into disease development. We train a machine learning method to predict emphysema from visually assessed expert labels. We use a multiple instance learning approach to predict both scan-level and region-level emphysema presence. We evaluate performance on 600 low-dose CT scans from the Danish Lung Cancer Screening Study and achieve an AUC of 0.82 for scan-level prediction and AUCs between 0.76 and 0.88 for region-level prediction.

Screening for lung cancer should be limited to a high-risk-population, and abnormalities in low-dose computed tomography (CT) screening images may be relevant for predicting the risk of lung cancer. Our aims were to compare the occurrence of visually detected emphysema and interstitial abnormalities in subjects with and without lung cancer in a screening population of smokers. Low-dose chest CT examinations (baseline and latest possible) of 1990 participants from The Danish Lung Cancer Screening Trial were independently evaluated by two observers who scored emphysema and interstitial abnormalities. Emphysema (lung density) was also measured quantitatively. Emphysema was seen more frequently and its extent was greater among participants with lung cancer on baseline (odds ratio (OR), 1.8, p = 0.017 and p = 0.002) and late examinations (OR 2.6, p

We present a fast and robust atlas-based algorithm for labeling airway trees, using geodesic distances in a geometric tree-space. Possible branch label configurations for an unlabeled airway tree are evaluated using distances to a training set of labeled airway trees. In tree-space, airway tree topology and geometry change continuously, giving a natural automatic handling of anatomical differences and noise. A hierarchical approach makes the algorithm efficient, assigning labels from the trachea and downwards. Only the airway centerline tree is used, which is relatively unaffected by pathology. The algorithm is evaluated on 80 segmented airway trees from 40 subjects at two time points, labeled by three medical experts each, testing accuracy, reproducibility and robustness in patients with chronic obstructive pulmonary disease (COPD). The accuracy of the algorithm is statistically similar to that of the experts and not significantly correlated with COPD severity. The reproducibility of the algorithm is significantly better than that of the experts, and negatively correlated with COPD severity. Evaluation of the algorithm on a longitudinal set of 8724 trees from a lung cancer screening trial shows that the algorithm can be used in large scale studies with high reproducibility, and that the negative correlation of reproducibility with COPD severity can be explained by missing branches, for instance due to segmentation problems in COPD patients. We conclude that the algorithm is robust to COPD severity given equally complete airway trees, and comparable in performance to that of experts in pulmonary medicine, emphasizing the suitability of the labeling algorithm for clinical use.

PURPOSE: To compare human observers to a mathematically derived computer model for differentiation between malignant and benign pulmonary nodules detected on baseline screening computed tomography (CT) scans.METHODS: A case-cohort study design was chosen. The study group consisted of 300 chest CT scans from the Danish Lung Cancer Screening Trial (DLCST). It included all scans with proven malignancies (n = 62) and two subsets of randomly selected baseline scans with benign nodules of all sizes (n = 120) and matched in size to the cancers, respectively (n = 118). Eleven observers and the computer model (PanCan) assigned a malignancy probability score to each nodule. Performances were expressed by area under the ROC curve (AUC). Performance differences were tested using the Dorfman, Berbaum and Metz method. Seven observers assessed morphological nodule characteristics using a predefined list. Differences in morphological features between malignant and size-matched benign nodules were analyzed using chi-square analysis with Bonferroni correction. A significant difference was defined at p < 0.004.RESULTS: Performances of the model and observers were equivalent (AUC 0.932 versus 0.910, p = 0.184) for risk-assessment of malignant and benign nodules of all sizes. However, human readers performed superior to the computer model for differentiating malignant nodules from size-matched benign nodules (AUC 0.819 versus 0.706, p < 0.001). Large variations between observers were seen for ROC areas and ranges of risk scores. Morphological findings indicative of malignancy referred to border characteristics (spiculation, p < 0.001) and perinodular architectural deformation (distortion of surrounding lung parenchyma architecture, p < 0.001; pleural retraction, p = 0.002).CONCLUSIONS: Computer model and human observers perform equivalent for differentiating malignant from randomly selected benign nodules, confirming the high potential of computer models for nodule risk estimation in population based screening studies. However, computer models highly rely on size as discriminator. Incorporation of other morphological criteria used by human observers to superiorly discriminate size-matched malignant from benign nodules, will further improve computer performance.

The introduction of lung cancer screening programs will produce an unprecedented amount of chest CT scans in the near future, which radiologists will have to read in order to decide on a patient follow-up strategy. According to the current guidelines, the workup of screen-detected nodules strongly relies on nodule size and nodule type. In this paper, we present a deep learning system based on multi-stream multi-scale convolutional networks, which automatically classifies all nodule types relevant for nodule workup. The system processes raw CT data containing a nodule without the need for any additional information such as nodule segmentation or nodule size and learns a representation of 3D data by analyzing an arbitrary number of 2D views of a given nodule. The deep learning system was trained with data from the Italian MILD screening trial and validated on an independent set of data from the Danish DLCST screening trial. We analyze the advantage of processing nodules at multiple scales with a multi-stream convolutional network architecture, and we show that the proposed deep learning system achieves performance at classifying nodule type that surpasses the one of classical machine learning approaches and is within the inter-observer variability among four experienced human observers.

The introduction of lung cancer screening programs will produce an unprecedented amount of chest CT scans in the near future, which radiologists will have to read in order to decide on a patient follow-up strategy. According to the current guidelines, the workup of screen-detected nodules strongly relies on nodule size and nodule type. In this paper, we present a deep learning system based on multi-stream multi-scale convolutional networks, which automatically classifies all nodule types relevant for nodule workup. The system processes raw CT data containing a nodule without the need for any additional information such as nodule segmentation or nodule size and learns a representation of 3D data by analyzing an arbitrary number of 2D views of a given nodule. The deep learning system was trained with data from the Italian MILD screening trial and validated on an independent set of data from the Danish DLCST screening trial. We analyze the advantage of processing nodules at multiple scales with a multi-stream convolutional network architecture, and we show that the proposed deep learning system achieves performance at classifying nodule type that surpasses the one of classical machine learning approaches and is within the inter-observer variability among four experienced human observers., Published on Scientific Reports

To evaluate interobserver agreement and time-trend in chest CT assessment of emphysema, airways, and interstitial abnormalities in a lung cancer screening cohort.Visual assessment of baseline and fifth-year examination of 1990 participants was performed independently by two observers. Results were standardised by means of an electronic score sheet; kappa and time-trend analyses were performed.Interobserver agreement was substantial in early emphysema diagnosis; highly significant (p 0.001) time-trends in both emphysema presence and grading were found (higher prevalence and grade of emphysema in late CT examinations). Significant progression in emphysema was seen in continuous smokers, but not in former smokers. Agreement on centrilobular emphysema subtype was substantial; agreement on paraseptal subtype, moderate. Agreement on panlobular and mixed subtypes was only fair. Agreement was fair regarding airway analysis. Interstitial abnormalities were infrequent in the cohort, and agreement on these was fair to moderate. A highly significant time-trend was found regarding interstitial abnormalities, which were more frequent in late examinations.Visual scoring of chest CT is able to characterise the presence, pattern, and progression of early emphysema. Continuous smokers progress; former smokers do not.• Substantial interobserver consistency in determining early-stage emphysema in low-dose CT. • Longitudinal analyses show clear time-trends for emphysema presence and grading. • For continuous smokers, progression of emphysema was seen in all lung zones. • For former smokers, progression of emphysema was undetectable by visual assessment. • Onset and progression of interstitial abnormalities are visually detectable.

This paper introduces a graph construction method for multi-dimensional and multi-surface segmentation problems. Such problems can be solved by searching for the optimal separating surfaces given the space of graph columns defined by an initial coarse surface. Conventional straight graph columns are not well suited for surfaces with high curvature, we therefore propose to derive columns from properly generated, non-intersecting flow lines. This guarantees solutions that do not self-intersect. The method is applied to segment human airway walls in computed tomography images in three-dimensions. Phantom measurements show that the inner and outer radii are estimated with sub-voxel accuracy. Two-dimensional manually annotated cross-sectional images were used to compare the results with those of another recently published graph based method. The proposed approach had an average overlap of 89.3 ± 5.8 %, and was on average within 0.096 ± 0.097 mm of the manually annotated surfaces, which is significantly better than what the previously published approach achieved. A medical expert visually evaluated 499 randomly extracted cross-sectional images from 499 scans and preferred the proposed approach in 68.5%, the alternative approach in 11.2%, and in 20.3% no method was favoured. Airway abnormality measurements obtained with the method on 490 scan pairs from a lung cancer screening trial correlate significantly with lung function and are reproducible; repeat scan R 2 of measures of the airway lumen diameter and wall area percentage in the airways from generation 0 (trachea) to 5 range from 0.96 to 0.73.

Computed tomography (CT) is an obvious modality for subclassification of COPD. Traditionally, the pulmonary involvement of chronic obstructive pulmonary disease (COPD) in smokers is understood as a combination of deleterious effects of smoking on small airways (chronic bronchitis and small airways disease) and distal to the airways with destruction and loss of lung parenchyma (emphysema). However, segmentation of airways is still experimental; with contemporary high-resolution CT (HRCT) we can just see the "entrance" of small airways, and until now changes in airway morphology that have been observed in COPD are subtle. Furthermore, recent results indicate that emphysema may also be the essential pathophysiologic mechanism behind the airflow limitation of COPD. The definition of COPD excludes bronchiectasis as a symptomatic subtype of COPD, and CT findings in chronic bronchitis and exacerbations of COPD are rather unspecific. This leaves emphysema as the most obvious candidate for subclassification of COPD. Both chest radiologists and pulmonary physicians are quite familiar with the appearance of various patterns of emphysema on HRCT, such as centrilobular, panlobular, and paraseptal emphysema. However, it has not yet been possible to develop operational definitions of these patterns that can be used by computer software to automatically classify CT scans into distinct patterns. In conclusion, even though various emphysema patterns can be recognized visually, CT has not yet demonstrated a great potential for automated subclassification of COPD, and it is an open question whether it will ever be possible to achieve success equivalent to that obtained by HRCT in the area of interstitial lung diseases.

Contains fulltext : 164462.pdf (Publisher’s version ) (Open Access) We propose a novel Computer-Aided Detection (CAD) system for pulmonary nodules using multi-view convolutional networks (ConvNets), for which discriminative features are automatically learnt from the training data. The network is fed with nodule candidates obtained by combining three candidate detectors specifically designed for solid, subsolid, and large nodules. For each candidate, a set of 2-D patches from differently oriented planes is extracted. The proposed architecture comprises multiple streams of 2-D ConvNets, for which the outputs are combined using a dedicated fusion method to get the final classification. Data augmentation and dropout are applied to avoid overfitting. On 888 scans of the publicly available LIDCIDRI dataset, our method reaches high detection sensitivities of 85.4% and 90.1% at 1 and 4 false positives per scan, respectively. An additional evaluation on independent datasets from the ANODE09 challenge and DLCST is performed. We showed that the proposed multi-view ConvNets is highly suited to be used for false positive reduction of a CAD system.

It remains unknown whether non-electrocardiogram-gated coronary artery calcium (CAC) score in lung cancer screening provides incremental prognostic value. The aim of this study was to evaluate the prognostic value of CAC in the Danish Lung Cancer Screening Trial (DLCST), in addition to conducting a systematic review and meta-analysis including previously published studies regarding CAC in lung cancer screening.In DLCST, we measured Agatston CAC scores in 1,945 current and former smokers. Causes of death were extracted from the Danish National Death Registry. We used Cox proportional hazards model to determine hazard ratios (HRs) of CAC scores. A weighted fixed-effects model was used for the meta-analysis.Median follow-up in DLCST was 7.1 years, and 55% were men. Overall survival rates associated with CAC scores of 0, 1-400, and400 were 98%, 96%, and 92% (p0.001), respectively. Adjusted HR of cardiovascular death associated with CAC400 was 3.8 (1.0-15) (p0.05). The meta-analysis included 28,045 asymptomatic participants. A high non-gated CAC score was associated with fatal or non-fatal cardiovascular events (p0.0001).Assessment of non-electrocardiogram-gated CAC in lung cancer screening programs is a robust prognostic measure of fatal or non-fatal cardiovascular events in current and former smokers independent of traditional cardiovascular risk factors.

Computer-Aided Detection (CAD) has been shown to be a promising tool for automatic detection of pulmonary nodules from computed tomography (CT) images. However, the vast majority of detected nodules are benign and do not require any treatment. For effective implementation of lung cancer screening programs, accurate identification of malignant nodules is the key. We investigate strategies to improve the performance of a CAD system in detecting nodules with a high probability of being cancers. Two strategies were proposed: (1) combining CAD detections with a recently published lung cancer risk prediction model and (2) the combination of multiple CAD systems. First, CAD systems were used to detect the nodules. Each CAD system produces markers with a certain degree of suspicion. Next, the malignancy probability was automatically computed for each marker, given nodule characteristics measured by the CAD system. Last, CAD degree of suspicion and malignancy probability were combined using the product rule. We evaluated the method using 62 nodules which were proven to be malignant cancers, from 180 scans of the Danish Lung Cancer Screening Trial. The malignant nodules were considered as positive samples, while all other findings were considered negative. Using a product rule, the best proposed system achieved an improvement in sensitivity, compared to the best individual CAD system, from 41.9% to 72.6% at 2 false positives (FPs)/scan and from 56.5% to 88.7% at 8 FPs/scan. Our experiment shows that combining a nodule malignancy probability with multiple CAD systems can increase the performance of computerized detection of lung cancer.

We present a fully automatic method for the assessment of spiculation of pulmonary nodules in low-dose Computed Tomography (CT) images. Spiculation is considered as one of the indicators of nodule malignancy and an important feature to assess in order to decide on a patient-tailored follow-up procedure. For this reason, lung cancer screening scenario would benefit from the presence of a fully automatic system for the assessment of spiculation. The presented framework relies on the fact that spiculated nodules mainly differ from non-spiculated ones in their morphology. In order to discriminate the two categories, information on morphology is captured by sampling intensity profiles along circular patterns on spherical surfaces centered on the nodule, in a multi-scale fashion. Each intensity profile is interpreted as a periodic signal, where the Fourier transform is applied, obtaining a spectrum. A library of spectra is created by clustering data via unsupervised learning. The centroids of the clusters are used to label back each spectrum in the sampling pattern. A compact descriptor encoding the nodule morphology is obtained as the histogram of labels along all the spherical surfaces and used to classify spiculated nodules via supervised learning. We tested our approach on a set of nodules from the Danish Lung Cancer Screening Trial (DLCST) dataset. Our results show that the proposed method outperforms other 3-D descriptors of morphology in the automatic assessment of spiculation.

Contains fulltext : 154627.pdf (Publisher’s version ) (Open Access) We present a novel descriptor for the characterization of pulmonary nodules in computed tomography (CT) images. The descriptor encodes information on nodule morphology and has scale-invariant and rotation-invariant properties. Information on nodule morphology is captured by sampling intensity profiles along circular patterns on spherical surfaces centered on the nodule, in a multi-scale fashion. Each intensity profile is interpreted as a periodic signal, where the Fourier transform is applied, obtaining a spectrum. A library of spectra is created and labeled via unsupervised clustering, obtaining a Bag-of- Frequencies, which is used to assign each spectra a label. The descriptor is obtained as the histogram of labels along all the spheres. Additional contributions are a technique to estimate the nodule size, based on the sampling strategy, as well as a technique to choose the most informative plane to cut a 2-D view of the nodule in the 3-D image. We evaluate the descriptor on several nodule morphology classification problems, namely discrimination of nodules versus vascular structures and characterization of spiculation. We validate the descriptor on data from European screening trials NELSON and DLCST and we compare it with state-of-the-art approaches for 3-D shape description in medical imaging and computer vision, namely SPHARM and 3-D SIFT, outperforming them in all the considered experiments.

Lung cancer risk models should be externally validated to test generalizability and clinical usefulness. The Danish Lung Cancer Screening Trial (DLCST) is a population-based prospective cohort study, used to assess the discriminative performances of the PanCan models. From the DLCST database, 1,152 nodules from 718 participants were included. Parsimonious and full PanCan risk prediction models were applied to DLCST data, and also coefficients of the model were recalculated using DLCST data. Receiver operating characteristics (ROC) curves and area under the curve (AUC) were used to evaluate risk discrimination. AUCs of 0.826–0.870 were found for DLCST data based on PanCan risk prediction models. In the DLCST, age and family history were significant predictors (p = 0.001 and p = 0.013). Female sex was not confirmed to be associated with higher risk of lung cancer; in fact opposing effects of sex were observed in the two cohorts. Thus, female sex appeared to lower the risk (p = 0.047 and p = 0.040) in the DLCST. High risk discrimination was validated in the DLCST cohort, mainly determined by nodule size. Age and family history of lung cancer were significant predictors and could be included in the parsimonious model. Sex appears to be a less useful predictor. • High accuracy in logistic modelling for lung cancer risk stratification of nodules. • Lung cancer risk prediction is primarily based on size of pulmonary nodules. • Nodule spiculation, age and family history of lung cancer are significant predictors. • Sex does not appear to be a useful risk predictor.

Statistical analysis of anatomical trees is hard to perform due to differences in the topological structure of the trees. In this paper we define statistical properties of leaf-labeled anatomical trees with geometric edge attributes by considering the anatomical trees as points in the geometric space of leaf-labeled trees. This tree-space is a geodesic metric space where any two trees are connected by a unique shortest path, which corresponds to a tree deformation. However, tree-space is not a manifold, and the usual strategy of performing statistical analysis in a tangent space and projecting onto tree-space is not available. Using tree-space and its shortest paths, a variety of statistical properties, such as mean, principal component, hypothesis testing and linear discriminant analysis can be defined. For some of these properties it is still an open problem how to compute them; others (like the mean) can be computed, but efficient alternatives are helpful in speeding up algorithms that use means iteratively, like hypothesis testing. In this paper, we take advantage of a very large dataset (N = 8016) to obtain computable approximations, under the assumption that the data trees parametrize the relevant parts of tree-space well. Using the developed approximate statistics, we illustrate how the structure and geometry of airway trees vary across a population and show that airway trees with Chronic Obstructive Pulmonary Disease come from a different distribution in tree-space than healthy ones. Software is available from http://image.diku.dk/aasa/software.php.

To study the effect of inspiration on airway dimensions measured in voluntary inspiration breath-hold examinations. 961 subjects with normal spirometry were selected from the Danish Lung Cancer Screening Trial. Subjects were examined annually for five years with low-dose CT. Automated software was utilized to segment lungs and airways, identify segmental bronchi, and match airway branches in all images of the same subject. Inspiration level was defined as segmented total lung volume (TLV) divided by predicted total lung capacity (pTLC). Mixed-effects models were used to predict relative change in lumen diameter (ALD) and wall thickness (AWT) in airways of generation 0 (trachea) to 7 and segmental bronchi (R1-R10 and L1-L10) from relative changes in inspiration level. Relative changes in ALD were related to relative changes in TLV/pTLC, and this distensibility increased with generation (p

We present a fast and robust supervised algorithm for labeling anatomical airway trees, based on geodesic distances in a geometric tree-space. Possible branch label configurations for a given tree are evaluated based on distances to a training set of labeled trees. In tree-space, the tree topology and geometry change continuously, giving a natural way to automatically handle anatomical differences and noise. The algorithm is made efficient using a hierarchical approach, in which labels are assigned from the top down. We only use features of the airway centerline tree, which are relatively unaffected by pathology. A thorough leave-one-patient-out evaluation of the algorithm is made on 40 segmented airway trees from 20 subjects labeled by 2 medical experts. We evaluate accuracy, reproducibility and robustness in patients with chronic obstructive pulmonary disease (COPD). Performance is statistically similar to the inter- and intra-expert agreement, and we found no significant correlation between COPD stage and labeling accuracy.

Statistical analysis of anatomical trees is hard to perform due to differences in the topological structure of the trees. In this paper we define statistical properties of leaf-labeled anatomical trees with geometric edge attributes by considering the anatomical trees as points in the geometric space of leaf-labeled trees. This tree-space is a geodesic metric space where any two trees are connected by a unique shortest path, which corresponds to a tree deformation. However, tree-space is not a manifold, and the usual strategy of performing statistical analysis in a tangent space and projecting onto tree-space is not available. Using tree-space and its shortest paths, a variety of statistical properties, such as mean, principal component, hypothesis testing and linear discriminant analysis can be defined. For some of these properties it is still an open problem how to compute them; others (like the mean) can be computed, but efficient alternatives are helpful in speeding up algorithms that use means iteratively, like hypothesis testing. In this paper, we take advantage of a very large dataset (N = 8016) to obtain computable approximations, under the assumption that the data trees parametrize the relevant parts of tree-space well. Using the developed approximate statistics, we illustrate how the structure and geometry of airway trees vary across a population and show that airway trees with Chronic Obstructive Pulmonary Disease come from a different distribution in tree-space than healthy ones. Software is available from http://image.diku.dk/aasa/software.php.