Automatic lung segmentation for accurate quantitation of volumetric X-ray CT images

Estimation of lung shunt fraction from simultaneous fluoroscopic and nuclear images

Radioembolisation with yttrium-90 (⁹⁰Y) is increasingly used as a treatment of unresectable liver malignancies. For safety, a scout dose of technetium-99m macroaggregated albumin (^99mTc-MAA) is used prior to the delivery of the therapeutic activity to mimic the deposition of ⁹⁰Y. One-day procedures are currently limited by the lack of nuclear images in the intervention room. To cope with this limitation, an interventional simultaneous fluoroscopic and nuclear imaging device is currently being developed. The purpose of this simulation study was to evaluate the accuracy of estimating the lung shunt fraction (LSF) of the scout dose in the intervention room with this device and compare it against current clinical methods.

METHODS

A male and female XCAT phantom, both with two respiratory profiles, were used to simulate various LSFs resulting from a scout dose of 150 MBq ^99mTc-MAA. Hybrid images were Monte Carlo simulated for breath-hold (5 s) and dynamic breathing (10 frames of 0.5 s) acquisitions. Nuclear images were corrected for attenuation with the fluoroscopic image and for organ overlap effects using a pre-treatment CT-scan. For comparison purposes, planar scintigraphy and mobile gamma camera images (both 300 s acquisition time) were simulated. Estimated LSFs were evaluated for all methods and compared to the phantom ground truth.

RESULTS

In the clinically relevant range of 10-20% LSF, hybrid imaging overestimated LSF with approximately 2 percentage points (pp) and 3 pp for the normal and irregular breathing phantoms, respectively. After organ overlap correction, LSF was estimated with a more constant error. Errors in planar scintigraphy and mobile gamma camera imaging were more dependent on LSF, body shape and breathing profile.

CONCLUSION

LSF can be estimated with a constant minor error with a hybrid imaging device. Estimated LSF is highly dependent on true LSF, body shape and breathing pattern when estimated with current clinical methods. The hybrid imaging device is capable of accurately estimating LSF within a few seconds in an interventional setting.

Multivariate modeling using quantitative CT metrics may improve accuracy of diagnosis of bronchiolitis obliterans syndrome after lung transplantation

BACKGROUND

To assess how quantitative CT (qCT) metrics compare to pulmonary function testing (PFT) and semi-quantitative image scores (SQS) to diagnose bronchiolitis obliterans syndrome (BOS), manifestation of chronic lung allograft dysfunction after lung transplantation (LTx), according to the type of LTx (unilateral or bilateral).

METHODS

Paired inspiratory-expiratory CT scans and PFTs of 176 LTx patients were analyzed retrospectively, and separated into BOS (78) and non-BOS (98) cohorts. SQS were assessed by 2 radiologists and graded (0-3) for features including mosaic attenuation and bronchiectasis. qCT metrics included lung volumes and air trapping volumes. Multivariate logistic regression (MVLR) and support vector machines (SVM) were used for the classification task.

RESULTS

MVLR and SVM models using PFT metrics demonstrated highest accuracy for bilateral LTx (max AUC 0.771), whereas models using qCT metrics-only outperformed models using SQS or PFTs in unilateral LTx (max AUC 0.817), to diagnose BOS. Adding PC (principal components) from qCT on top of PFT improved model diagnostic accuracy for all transplant types.

CONCLUSIONS

Combinations of qCT metrics augment the diagnostic performance of PFTs, are superior to SQS to predict BOS status, and outperform PFTs in the unilateral LTx group. This suggests that latent information on paired volumetric CT may allow early diagnosis of BOS in LTx patients, particularly in unilateral LTx.

Automated mediastinal lymph node detection from CT volumes based on intensity targeted radial structure tensor analysis

This paper presents a local intensity structure analysis based on an intensity targeted radial structure tensor (ITRST) and the blob-like structure enhancement filter based on it (ITRST filter) for the mediastinal lymph node detection algorithm from chest computed tomography (CT) volumes. Although the filter based on radial structure tensor analysis (RST filter) based on conventional RST analysis can be utilized to detect lymph nodes, some lymph nodes adjacent to regions with extremely high or low intensities cannot be detected. Therefore, we propose the ITRST filter, which integrates the prior knowledge on detection target intensity range into the RST filter. Our lymph node detection algorithm consists of two steps: (1) obtaining candidate regions using the ITRST filter and (2) removing false positives (FPs) using the support vector machine classifier. We evaluated lymph node detection performance of the ITRST filter on 47 contrast-enhanced chest CT volumes and compared it with the RST and Hessian filters. The detection rate of the ITRST filter was 84.2% with 9.1 FPs/volume for lymph nodes whose short axis was at least 10 mm, which outperformed the RST and Hessian filters.

Size variation and collapse of emphysema holes at inspiration and expiration CT scan: evaluation with modified length scale method and image co-registration

A novel approach of size-based emphysema clustering has been developed, and the size variation and collapse of holes in emphysema clusters are evaluated at inspiratory and expiratory computed tomography (CT). Thirty patients were visually evaluated for the size-based emphysema clustering technique and a total of 72 patients were evaluated for analyzing collapse of the emphysema hole in this study. A new approach for the size differentiation of emphysema holes was developed using the length scale, Gaussian low-pass filtering, and iteration approach. Then, the volumetric CT results of the emphysema patients were analyzed using the new method, and deformable registration was carried out between inspiratory and expiratory CT. Blind visual evaluations of EI by two readers had significant correlations with the classification using the size-based emphysema clustering method (r-values of reader 1: 0.186, 0.890, 0.915, and 0.941; reader 2: 0.540, 0.667, 0.919, and 0.942). The results of collapse of emphysema holes using deformable registration were compared with the pulmonary function test (PFT) parameters using the Pearson's correlation test. The mean extents of low-attenuation area (LAA), E1 (<1.5 mm), E2 (<7 mm), E3 (<15 mm), and E4 (≥15 mm) were 25.9%, 3.0%, 11.4%, 7.6%, and 3.9%, respectively, at the inspiratory CT, and 15.3%, 1.4%, 6.9%, 4.3%, and 2.6%, respectively at the expiratory CT. The extents of LAA, E2, E3, and E4 were found to be significantly correlated with the PFT parameters (r=-0.53, -0.43, -0.48, and -0.25), with forced expiratory volume in 1 second (FEV₁; -0.81, -0.62, -0.75, and -0.40), and with diffusing capacity of the lungs for carbon monoxide (cDLco), respectively. The fraction of emphysema that shifted to the smaller subgroup showed a significant correlation with FEV₁, cDLco, forced expiratory flow at 25%-75% of forced vital capacity, and residual volume (RV)/total lung capacity (r=0.56, 0.73, 0.40, and -0.58). A detailed assessment of the size variation and collapse of emphysema holes may be useful for understanding the dynamic collapse of emphysema and its functional relation.

Pulmonary Lobe Segmentation With Probabilistic Segmentation of the Fissures and a Groupwise Fissure Prior

A fully automated, unsupervised lobe segmentation algorithm is presented based on a probabilistic segmentation of the fissures and the simultaneous construction of a populationmodel of the fissures. A two-class probabilistic segmentation segments the lung into candidate fissure voxels and the surrounding parenchyma. This was combined with anatomical information and a groupwise fissure prior to drive non-parametric surface fitting to obtain the final segmentation. The performance of our fissure segmentation was validated on 30 patients from the chronic obstructive pulmonary disease COPDGene cohort, achieving a high median F₁ -score of 0.90 and showed general insensitivity to filter parameters. We evaluated our lobe segmentation algorithm on the Lobe and Lung Analysis 2011 dataset, which contains 55 cases at varying levels of pathology. We achieved the highest score of 0.884 of the automated algorithms. Our method was further tested quantitatively and qualitatively on 80 patients from the COPDgene study at varying levels of functional impairment. Accurate segmentation of the lobes is shown at various degrees of fissure incompleteness for 96% of all cases. We also show the utility of including a groupwise prior in segmenting the lobes in regions of grossly incomplete fissures.

Fuzzy-C-Means Clustering Based Segmentation and CNN-Classification for Accurate Segmentation of Lung Nodules

Objective: Accurate segmentation of abnormal and healthy lungs is very crucial for a steadfast computer-aided disease diagnostics. Methods: For this purpose a stack of chest CT scans are processed. In this paper, novel methods are proposed for segmentation of the multimodal grayscale lung CT scan. In the conventional methods using Markov–Gibbs Random Field (MGRF) model the required regions of interest (ROI) are identified. Result: The results of proposed FCM and CNN based process are compared with the results obtained from the conventional method using MGRF model. The results illustrate that the proposed method can able to segment the various kinds of complex multimodal medical images precisely. Conclusion: However, in this paper, to obtain an exact boundary of the regions, every empirical dispersion of the image is computed by Fuzzy C-Means Clustering segmentation. A classification process based on the Convolutional Neural Network (CNN) classifier is accomplished to distinguish the normal tissue and the abnormal tissue. The experimental evaluation is done using the Interstitial Lung Disease (ILD) database.

A novel lung cancer detection algorithm for CADs based on SSP and Level Set

The fuzzy degree of lung nodule boundary is the most important cue to judge the lung cancer in CT images. Based on this feature, the paper proposes a novel lung cancer detection method for CT images based on the super-pixels and the level set segmentation methods. In the proposed methods, the super-pixels method is used to segment the lung region and the suspected lung cancer lesion region in the CT image. The super-pixels method and a level set method are used to segment the suspected lung cancer lesion region simultaneously. Finally, the cancer is determined by the difference between results of the two segmentation methods. Experimental results show that the proposed algorithm has a high accuracy for lung cancer detection in CT images. For gross glass nodule, pleural nodule, the vascular nodules and solitary nodules, the sensitivity of the detection algorithm are respectively 91.3%, 96.3%, 80.9% and 82.3%.

Novel Assessment of Interstitial Lung Disease Using the "Computer-Aided Lung Informatics for Pathology Evaluation and Rating" (CALIPER) Software System in Idiopathic Inflammatory Myopathies

PURPOSE

To evaluate the correlation between measurements from quantitative thoracic high-resolution CT (HRCT) analysis with "Computer-Aided Lung Informatics for Pathology Evaluation and Rating" (CALIPER) software and measurements from pulmonary function tests (PFTs) in patients with idiopathic inflammatory myopathies (IIM)-associated interstitial lung disease (ILD).

METHODS

A cohort of patients with IIM-associated ILD seen at Mayo Clinic was identified from medical record review. Retrospective analysis of HRCT data and PFTs at baseline and 1 year was performed. The abnormalities in HRCT were quantified using CALIPER software.

RESULTS

A total of 110 patients were identified. At baseline, total interstitial abnormalities as measured by CALIPER, both by absolute volume and by percentage of total lung volume, had a significant negative correlation with diffusing capacity for carbon monoxide (DLCO), total lung capacity (TLC), and oxygen saturation. Analysis by subtype of interstitial abnormality revealed significant negative correlations between ground glass opacities (GGO) and reticular density (RD) with DLCO and TLC. At one year, changes of total interstitial abnormalities compared with baseline had a significant negative correlation with changes of TLC and oxygen saturation. A negative correlation between changes of total interstitial abnormalities and DLCO was also observed, but it was not statistically significant. Analysis by subtype of interstitial abnormality revealed negative correlations between changes of GGO and RD and changes of DLCO, TLC, and oxygen saturation, but most of the correlations did not achieve statistical significance.

CONCLUSION

CALIPER measurements correlate well with functional measurements in patients with IIM-associated ILD.

An integrated segmentation and shape-based classification scheme for distinguishing adenocarcinomas from granulomas on lung CT

PURPOSE

Distinguishing between benign granulmoas and adenocarcinomas is confounded by their similar visual appearance on routine CT scans. Unfortunately, owing to the inability to discriminate these lesions radigraphically, many patients with benign granulomas are subjected to unnecessary surgical wedge resections and biopsies for pathologic confirmation of cancer presence or absence. This suggests the need for improved computerized characterization of these nodules in order to distinguish between these two classes of lesions on CT scans. While there has been substantial interest in the use of textural analysis for radiomic characterization of lung nodules, relatively less work has been done in shape based characterization of lung nodules, particularly with respect to granulmoas and adenocarcinomas. The primary goal of this study is to evaluate the role of 3D shape features for discrimination of benign granulomas from malignant adenocarcinomas on lung CT images. Towards this end we present an integrated framework for segmentation, feature characterization and classification of these nodules on CT.

METHODS

The nodule segmentation method starts with separation of lung regions from the surrounding lung anatomy. Next, the lung CT scans are projected into and represented in a three dimensional spectral embedding (SE) space, allowing for better determination of the boundaries of the nodule. This then enables the application of a gradient vector flow active contour (SEGvAC) model for nodule boundary extraction. A set of 24 shape features from both 2D slices and 3D surface of the segmented nodules are extracted, including features pertaining to the angularity, spiculation, elongation and nodule compactness. A feature selection scheme, PCA-VIP, is employed to identify the most discriminating set of features to distinguish granulmoas from adenocarcinomas within a learning set of 82 patients. The features thus identified were then combined with a support vector machine classifier and independently validated on a distinct test set comprising 67 patients. The performance of the classifier for both of the training and validation cohorts was evaluated by the area under receiver characteristic curve (ROC).

RESULTS

We used 82 and 67 studies from two different institutions respectively for training and independent validation of the model and the shape features. The Dice coefficient between automatically segmented nodules by SEGvAC and the manual delineations by expert radiologists (readers) was 0.84± 0.04 whereas inter-reader segmentation agreement was 0.79± 0.12. We also identified a set of consistent features (Roughness, Convexity and Spherecity) that were found to be strongly correlated across both manual and automated nodule segmentations (R > 0.80, p < 0.0001) and capture the marginal smoothness and 3D compactness of the nodules. On the independent validation set of 67 studies our classifier yielded a ROC AUC of 0.72 and 0.64 for manually- and automatically segmented nodules respectively. On a subset of 20 studies, the AUCs for the two expert radiologists and 1 pulmonologist were found to be 0.82, 0.68 and 0.58 respectively.

CONCLUSIONS

The major finding of this study was that certain shape features appear to differentially express between granulomas and adenocarcinomas and thus computer extracted shape cues could be used to distinguish these radiographically similar pathologies.

Chronic Obstructive Pulmonary Disease: Lobar Analysis with Hyperpolarized 129Xe MR Imaging

Purpose To compare lobar ventilation and apparent diffusion coefficient (ADC) values obtained with hyperpolarized xenon 129 (129Xe) magnetic resonance (MR) imaging to quantitative computed tomography (CT) metrics on a lobar basis and pulmonary function test (PFT) results on a whole-lung basis in patients with chronic obstructive pulmonary disease (COPD). Materials and Methods The study was approved by the National Research Ethics Service Committee; written informed consent was obtained from all patients. Twenty-two patients with COPD (Global Initiative for Chronic Obstructive Lung Disease stage II-IV) underwent hyperpolarized 129Xe MR imaging at 1.5 T, quantitative CT, and PFTs. Whole-lung and lobar 129Xe MR imaging parameters were obtained by using automated segmentation of multisection hyperpolarized 129Xe MR ventilation images and hyperpolarized 129Xe MR diffusion-weighted images after coregistration to CT scans. Whole-lung and lobar quantitative CT-derived metrics for emphysema and bronchial wall thickness were calculated. Pearson correlation coefficients were used to evaluate the relationship between imaging measures and PFT results. Results Percentage ventilated volume and average ADC at lobar 129Xe MR imaging showed correlation with percentage emphysema at lobar quantitative CT (r = -0.32, P < .001 and r = 0.75, P < .0001, respectively). The average ADC at whole-lung 129Xe MR imaging showed moderate correlation with PFT results (percentage predicted transfer factor of the lung for carbon monoxide [Tlco]: r = -0.61, P < .005) and percentage predicted functional residual capacity (r = 0.47, P < .05). Whole-lung quantitative CT percentage emphysema also showed statistically significant correlation with percentage predicted Tlco (r = -0.65, P < .005). Conclusion Lobar ventilation and ADC values obtained from hyperpolarized 129Xe MR imaging demonstrated correlation with quantitative CT percentage emphysema on a lobar basis and with PFT results on a whole-lung basis. © RSNA, 2016.

Electrical impedance tomography for lung ventilation monitoring of the dog

BACKGROUND

Electrical impedance tomography (EIT) is a radiation free technique which takes advantage of the different electrical conductivities of different tissues. Its main field of application is lung ventilation monitoring. The aim of this prospective study was to evaluate the feasibility of collecting EIT information on a sample of dogs with different thoracic shapes under clinical conditions by connecting an electrode belt without fur clipping.

MATERIAL AND METHODS

Fifteen pulmonary healthy dogs were anaesthetized, positioned in sternal recumbency and ventilated in a pressure-controlled mode at three different positive end-expiratory pressure levels (PEEP) of 0, 5 and 10 cmH₂O for five breaths each, with a peak inspiratory pressure of 15 cmH₂O. The impedance changes were recorded with a commercial EIT device applied around the thorax. Subsequently, the ventilation regime was repeated and a computed tomography scan (CT) of the same thoracic segment was performed for each PEEP level. The tidal volume (V_t) was recorded. For the collection of EIT data the sum of regional impedance changes was recorded. The impedance value of the entire lung (global) was recorded and the ventilated area was quartered into four regions of interest (ROI). In a CT image with the fewest adjacent organs, lung tissue was selected to obtain the mean value of lung radiodensitiy in Hounsfield-Units (HU) for the entire lung and for the four ROIs.

RESULTS

EIT recordings via the electrode belt were possible without clipping. There was a significant correlation for the parameters of aeration as measured by EIT and CT for both the entire ventilated lung and the corresponding ROIs. The increasing PEEP resulted in a proportional increase of the impedance, and there was a negative correlation between EIT and V_t. The better ventilated dorsal ROIs could be identified using both EIT and CT. An intra-assay coefficient of variation showed a good reproducibility for lung ventilation in anaesthetized dogs in the EIT.

DISCUSSION

The results show that EIT is a reliable method for evaluating the ventilation of dogs in a clinical setting. The accuracy of EIT might be improved by using a mesh corresponding to the different thoracic shapes of the dogs.

Isotropic Total Variation Regularization of Displacements in Parametric Image Registration

Spatial regularization is essential in image registration, which is an ill-posed problem. Regularization can help to avoid both physically implausible displacement fields and local minima during optimization. Tikhonov regularization (squared l₂ -norm) is unable to correctly represent non-smooth displacement fields, that can, for example, occur at sliding interfaces in the thorax and abdomen in image time-series during respiration. In this paper, isotropic Total Variation (TV) regularization is used to enable accurate registration near such interfaces. We further develop the TV-regularization for parametric displacement fields and provide an efficient numerical solution scheme using the Alternating Directions Method of Multipliers (ADMM). The proposed method was successfully applied to four clinical databases which capture breathing motion, including CT lung and MR liver images. It provided accurate registration results for the whole volume. A key strength of our proposed method is that it does not depend on organ masks that are conventionally required by many algorithms to avoid errors at sliding interfaces. Furthermore, our method is robust to parameter selection, allowing the use of the same parameters for all tested databases. The average target registration error (TRE) of our method is superior (10% to 40%) to other techniques in the literature. It provides precise motion quantification and sliding detection with sub-pixel accuracy on the publicly available breathing motion databases (mean TREs of 0.95 mm for DIR 4D CT, 0.96 mm for DIR COPDgene, 0.91 mm for POPI databases).

A Novel Mouse Segmentation Method Based on Dynamic Contrast Enhanced Micro-CT Images

With the development of hybrid imaging scanners, micro-CT is widely used in locating abnormalities, studying drug metabolism, and providing structural priors to aid image reconstruction in functional imaging. Due to the low contrast of soft tissues, segmentation of soft tissue organs from mouse micro-CT images is a challenging problem. In this paper, we propose a mouse segmentation scheme based on dynamic contrast enhanced micro-CT images. With a homemade fast scanning micro-CT scanner, dynamic contrast enhanced images were acquired before and after injection of non-ionic iodinated contrast agents (iohexol). Then the feature vector of each voxel was extracted from the signal intensities at different time points. Based on these features, the heart, liver, spleen, lung, and kidney could be classified into different categories and extracted from separate categories by morphological processing. The bone structure was segmented using a thresholding method. Our method was validated on seven BALB/c mice using two different classifiers: a support vector machine classifier with a radial basis function kernel and a random forest classifier. The results were compared to manual segmentation, and the performance was assessed using the Dice similarity coefficient, false positive ratio, and false negative ratio. The results showed high accuracy with the Dice similarity coefficient ranging from 0.709 ± 0.078 for the spleen to 0.929 ± 0.006 for the kidney.

Automatic gallbladder and gallstone regions segmentation in ultrasound image

PURPOSE

As gallbladder diseases including gallstone and cholecystitis are mainly diagnosed by using ultra-sonographic examinations, we propose a novel method to segment the gallbladder and gallstones in ultrasound images.

METHODS

The method is divided into five steps. Firstly, a modified Otsu algorithm is combined with the anisotropic diffusion to reduce speckle noise and enhance image contrast. The Otsu algorithm separates distinctly the weak edge regions from the central region of the gallbladder. Secondly, a global morphology filtering algorithm is adopted for acquiring the fine gallbladder region. Thirdly, a parameter-adaptive pulse-coupled neural network (PA-PCNN) is employed to obtain the high-intensity regions including gallstones. Fourthly, a modified region-growing algorithm is used to eliminate physicians' labeled regions and avoid over-segmentation of gallstones. It also has good self-adaptability within the growth cycle in light of the specified growing and terminating conditions. Fifthly, the smoothing contours of the detected gallbladder and gallstones are obtained by the locally weighted regression smoothing (LOESS).

RESULTS

We test the proposed method on the clinical data from Gansu Provincial Hospital of China and obtain encouraging results. For the gallbladder and gallstones, average similarity percent of contours (EVA) containing metrics dice's similarity , overlap fraction and overlap value is 86.01 and 79.81%, respectively; position error is 1.7675 and 0.5414 mm, respectively; runtime is 4.2211 and 0.6603 s, respectively. Our method then achieves competitive performance compared with the state-of-the-art methods.

CONCLUSIONS

The proposed method is potential to assist physicians for diagnosing the gallbladder disease rapidly and effectively.

Accurate Lungs Segmentation on CT Chest Images by Adaptive Appearance-Guided Shape Modeling

To accurately segment pathological and healthy lungs for reliable computer-aided disease diagnostics, a stack of chest CT scans is modeled as a sample of a spatially inhomogeneous joint 3D Markov-Gibbs random field (MGRF) of voxel-wise lung and chest CT image signals (intensities). The proposed learnable MGRF integrates two visual appearance sub-models with an adaptive lung shape submodel. The first-order appearance submodel accounts for both the original CT image and its Gaussian scale space (GSS) filtered version to specify local and global signal properties, respectively. Each empirical marginal probability distribution of signals is closely approximated with a linear combination of discrete Gaussians (LCDG), containing two positive dominant and multiple sign-alternate subordinate DGs. The approximation is separated into two LCDGs to describe individually the lungs and their background, i.e., all other chest tissues. The second-order appearance submodel quantifies conditional pairwise intensity dependencies in the nearest voxel 26-neighborhood in both the original and GSS-filtered images. The shape submodel is built for a set of training data and is adapted during segmentation using both the lung and chest appearances. The accuracy of the proposed segmentation framework is quantitatively assessed using two public databases (ISBI VESSEL12 challenge and MICCAI LOLA11 challenge) and our own database with, respectively, 20, 55, and 30 CT images of various lung pathologies acquired with different scanners and protocols. Quantitative assessment of our framework in terms of Dice similarity coefficients, 95-percentile bidirectional Hausdorff distances, and percentage volume differences confirms the high accuracy of our model on both our database (98.4±1.0%, 2.2±1.0 mm, 0.42±0.10%) and the VESSEL12 database (99.0±0.5%, 2.1±1.6 mm, 0.39±0.20%), respectively. Similarly, the accuracy of our approach is further verified via a blind evaluation by the organizers of the LOLA11 competition, where an average overlap of 98.0% with the expert's segmentation is yielded on all 55 subjects with our framework being ranked first among all the state-of-the-art techniques compared.

Monitoring change in volume of calcifications in juvenile idiopathic inflammatory myopathy: a pilot study using low dose computed tomography

BACKGROUND

Dystrophic calcifications may occur in patients with J uvenile Idiopathic Inflammatory Myopathy (JIIM) as well as other connective tissue and metabolic diseases, but a reliable method of measuring the volume of these calcifications has not been established. The purpose of this study is to determine the feasibility of low dose, limited slice, Computed Tomography (CT) to measure objectively in-situ calcification volumes in patients with JIIM over time.

METHODS

Ten JIIM patients (eight JDM, two Overlap) with calcifications were prospectively recruited over a 2-year period to undergo two limited, low dose, four-slice CT scans. Calculation of the volume of calcifications used a CT post processing workstation. Additional patient data included: Disease Activity Scores (DAS), Childhood Myositis Assessment Scale (CMAS), myositis specific antibodies (MSA), and the TNFα-308 promoter region A/G polymorphism. Statistical analysis utilized the Pearson correlation coefficient, the paired t-test and descriptive statistics.

RESULTS

Ten JIIM, mean age 14.54 ± 4.54 years, had a duration of untreated disease of 8.68 ± 5.65 months  MSA status: U1RNP (1), PM-Scl (1), Ro (1, 4 indeterminate), p155/140 (2), MJ (3), Mi-2 indeterminate (1), negative (3). 4/8 JDM (50%) were TNF-α-308 A+. Overall, the calcification volumes tended to decrease from the first to the second CT study by 0.5 cm³ (from 2.79 ± 1.98 cm³ to 2.29 ± 2.25 cm³). The average effective radiation dose was 0.007 ± 0.002, 0.010 ± 0.005, and 0.245 mSv for the upper extremity, lower extremity and chest, respectively (compared to a standard chest x-ray-- 0.02mSV effective dosage).

CONCLUSION

We conclude: 1) the limited low dose CT technique provides objective data about volume of the calcifications in JIIM; 2) measuring the volume of calcifications in an extremity is associated with minimal radiation exposure; 3) This method may be useful to evaluate the efficacy of therapies for JIIM dystrophic calcification.

Anatomy-based algorithm for automatic segmentation of human diaphragm in noncontrast computed tomography images

In-depth understanding of the diaphragm's anatomy and physiology has been of great interest to the medical community, as it is the most important muscle of the respiratory system. While noncontrast four-dimensional (4-D) computed tomography (CT) imaging provides an interesting opportunity for effective acquisition of anatomical and/or functional information from a single modality, segmenting the diaphragm in such images is very challenging not only because of the diaphragm's lack of image contrast with its surrounding organs but also because of respiration-induced motion artifacts in 4-D CT images. To account for such limitations, we present an automatic segmentation algorithm, which is based on a priori knowledge of diaphragm anatomy. The novelty of the algorithm lies in using the diaphragm's easy-to-segment contacting organs-including the lungs, heart, aorta, and ribcage-to guide the diaphragm's segmentation. Obtained results indicate that average mean distance to the closest point between diaphragms segmented using the proposed technique and corresponding manual segmentation is [Formula: see text], which is favorable. An important feature of the proposed technique is that it is the first algorithm to delineate the entire diaphragm. Such delineation facilitates applications, where the diaphragm boundary conditions are required such as biomechanical modeling for in-depth understanding of the diaphragm physiology.

Automatic Approach for Lung Segmentation with Juxta-Pleural Nodules from Thoracic CT Based on Contour Tracing and Correction

This paper presents a fully automatic framework for lung segmentation, in which juxta-pleural nodule problem is brought into strong focus. The proposed scheme consists of three phases: skin boundary detection, rough segmentation of lung contour, and pulmonary parenchyma refinement. Firstly, chest skin boundary is extracted through image aligning, morphology operation, and connective region analysis. Secondly, diagonal-based border tracing is implemented for lung contour segmentation, with maximum cost path algorithm used for separating the left and right lungs. Finally, by arc-based border smoothing and concave-based border correction, the refined pulmonary parenchyma is obtained. The proposed scheme is evaluated on 45 volumes of chest scans, with volume difference (VD) 11.15 ± 69.63 cm³, volume overlap error (VOE) 3.5057 ± 1.3719%, average surface distance (ASD) 0.7917 ± 0.2741 mm, root mean square distance (RMSD) 1.6957 ± 0.6568 mm, maximum symmetric absolute surface distance (MSD) 21.3430 ± 8.1743 mm, and average time-cost 2 seconds per image. The preliminary results on accuracy and complexity prove that our scheme is a promising tool for lung segmentation with juxta-pleural nodules.

Mortality prediction in idiopathic pulmonary fibrosis: evaluation of computer-based CT analysis with conventional severity measures

Computer-based computed tomography (CT) analysis can provide objective quantitation of disease in idiopathic pulmonary fibrosis (IPF). A computer algorithm, CALIPER, was compared with conventional CT and pulmonary function measures of disease severity for mortality prediction.

CT and pulmonary function variables (forced expiratory volume in 1 s, forced vital capacity, diffusion capacity of the lung for carbon monoxide, transfer coefficient of the lung for carbon monoxide and composite physiologic index (CPI)) of 283 consecutive patients with a multidisciplinary diagnosis of IPF were evaluated against mortality. Visual and CALIPER CT features included total extent of interstitial lung disease, honeycombing, reticular pattern, ground glass opacities and emphysema. In addition, CALIPER scored pulmonary vessel volume (PVV) while traction bronchiectasis and consolidation were only scored visually. A combination of mortality predictors was compared with the Gender, Age, Physiology model.

On univariate analyses, all visual and CALIPER-derived interstitial features and functional indices were predictive of mortality to a 0.01 level of significance. On multivariate analysis, visual CT parameters were discarded. Independent predictors of mortality were CPI (hazard ratio (95% CI) 1.05 (1.02–1.07), p<0.001) and two CALIPER parameters: PVV (1.23 (1.08–1.40), p=0.001) and honeycombing (1.18 (1.06–1.32), p=0.002). A three-group staging system derived from this model was powerfully predictive of mortality (2.23 (1.85–2.69), p<0.0001).

CALIPER-derived parameters, in particular PVV, are more accurate prognostically than traditional visual CT scores. Quantitative tools such as CALIPER have the potential to improve staging systems in IPF.

Supine to upright lung mechanics: do changes in lung shape influence lung tissue deformation?

In this study we analyze lung shape change between the upright and supine postures and the effect of this shape change on the deformation of lung tissue under gravity. We use supine computed tomography images along with upright tomosynthesis images obtained on the same day to show that there is significant diaphragmatic movement between postures. Using a continuum model of lung tissue deformation under gravity we show that the shape changes due to this diaphragmatic movement could result in different lung tissue expansion patterns between supine and upright lungs. This is an essential consideration when interpreting imaging data acquired in different postures or translating data acquired in supine imaging to upright function.

Automatic multiorgan segmentation in CT images of the male pelvis using region-specific hierarchical appearance cluster models

PURPOSE

Accurate segmentation of pelvic organs in CT images is of great importance in external beam radiotherapy for prostate cancer. The aim of this studying is to develop a novel method for automatic, multiorgan segmentation of the male pelvis.

METHODS

The authors' segmentation method consists of several stages. First, a pretreatment includes parameterization, principal component analysis (PCA), and an established process of region-specific hierarchical appearance cluster (RSHAC) model which was executed on the training dataset. After the preprocessing, online automatic segmentation of new CT images is achieved by combining the RSHAC model with the PCA-based point distribution model. Fifty pelvic CT from eight prostate cancer patients were used as the training dataset. From another 20 prostate cancer patients, 210 CT images were used for independent validation of the segmentation method.

RESULTS

In the training dataset, 15 PCA modes were needed to represent 95% of shape variations of pelvic organs. When tested on the validation dataset, the authors' segmentation method had an average Dice similarity coefficient and mean absolute distance of 0.751 and 0.371 cm, 0.783 and 0.303 cm, 0.573 and 0.604 cm for prostate, bladder, and rectum, respectively. The automated segmentation process took on average 5 min on a personal computer equipped with Core 2 Duo CPU of 2.8 GHz and 8 GB RAM.

CONCLUSIONS

The authors have developed an efficient and reliable method for automatic segmentation of multiple organs in the male pelvis. This method should be useful for treatment planning and adaptive replanning for prostate cancer radiotherapy. With this method, the physicist can improve the work efficiency and stability.

Alterations of diaphragm and rib cage morphometry in severe COPD patients by CT analysis

Although it is known that in patients with COPD acute hyperinflation determines shortening of the inspiratory muscles, its effects on both diaphragm and rib cage morphology are still to be investigated. In this preliminary study the relationships between hyperinflation, emphysema, diaphragm and rib cage geometry were studied in 5 severe COPD patients and 5 healthy subjects. An automatic software was developed to obtain the 3-D reconstruction of diaphragm and rib cage from CT scans taken at total lung capacity (TLC) and residual volume (RV). Dome surface area (Ado), radius of curvature, length (Ld) and position (referred to xiphoid level) of the diaphragm and antero-posterior (A-P) and transverse (T) diameters of rib cage were calculated at both volumes. Ado and Ld were similar in COPD and controls when compared at similar absolute lung volumes. Radius of curvature was significantly higher in COPD than in controls only at TLC. In COPD, the range of diaphragm position was invariantly below the xiphoid level, while in controls the top of diaphragm dome was always above it. Rib cage diameters were not different at TLC. A-P diameter was greater in COPD than in controls at RV, while T diameters were similar. In conclusion, in severe COPD diaphragm and rib cage geometry is altered at RV. The lower position of diaphragm is associated to smaller A-P but not transversal rib cage diameters, such that rib cage adopts a more circular shape.

Rheumatoid arthritis-associated autoantibodies and subclinical interstitial lung disease: the Multi-Ethnic Study of Atherosclerosis

BACKGROUND

Adults with interstitial lung disease (ILD) often have serologic evidence of autoimmunity of uncertain significance without overt autoimmune disease. We examined associations of rheumatoid arthritis (RA)-associated antibodies with subclinical ILD in community-dwelling adults.

METHODS

We measured serum rheumatoid factor (RF) and anticyclic citrullinated peptide antibody (anti-CCP) and high attenuation areas (HAAs; CT attenuation values between -600 and -250 Hounsfield units) on cardiac CT in 6736 community-dwelling US adults enrolled in the Multi-Ethnic Study of Atherosclerosis. We measured interstitial lung abnormalities (ILAs) in 2907 full-lung CTs at 9.5-year median follow-up. We used generalised linear and additive models to examine associations between autoantibodies and both HAA and ILA, and tested for effect modification by smoking.

RESULTS

In adjusted models, HAA increased by 0.49% (95% CI 0.11% to 0.86%) per doubling of RF IgM and by 0.95% (95% CI 0.50% to 1.40%) per RF IgA doubling. ILA prevalence increased by 11% (95% CI 3% to 20%) per RF IgA doubling. Smoking modified the associations of both RF IgM and anti-CCP with both HAA and ILA (interaction p values varied from 0.01 to 0.09). Among ever smokers, HAA increased by 0.81% (95% CI 0.33% to 1.30%) and ILA prevalence increased by 14% (95% CI 5% to 24%,) per RF IgM doubling; and HAA increased by 1.31% (95% CI 0.45% to 2.18%) and ILA prevalence increased by 13% (95% CI 2% to 24%) per anti-CCP doubling. Among never smokers, no meaningful associations were detected.

CONCLUSIONS

RA-related autoimmunity is associated with both quantitative and qualitative subclinical ILD phenotypes on CT, particularly among ever smokers.

Tissue segmentation of computed tomography images using a Random Forest algorithm: a feasibility study

There is a need for robust, fully automated whole body organ segmentation for diagnostic CT. This study investigates and optimizes a Random Forest algorithm for automated organ segmentation; explores the limitations of a Random Forest algorithm applied to the CT environment; and demonstrates segmentation accuracy in a feasibility study of pediatric and adult patients. To the best of our knowledge, this is the first study to investigate a trainable Weka segmentation (TWS) implementation using Random Forest machine-learning as a means to develop a fully automated tissue segmentation tool developed specifically for pediatric and adult examinations in a diagnostic CT environment. Current innovation in computed tomography (CT) is focused on radiomics, patient-specific radiation dose calculation, and image quality improvement using iterative reconstruction, all of which require specific knowledge of tissue and organ systems within a CT image. The purpose of this study was to develop a fully automated Random Forest classifier algorithm for segmentation of neck-chest-abdomen-pelvis CT examinations based on pediatric and adult CT protocols. Seven materials were classified: background, lung/internal air or gas, fat, muscle, solid organ parenchyma, blood/contrast enhanced fluid, and bone tissue using Matlab and the TWS plugin of FIJI. The following classifier feature filters of TWS were investigated: minimum, maximum, mean, and variance evaluated over a voxel radius of 2 (n) , (n from 0 to 4), along with noise reduction and edge preserving filters: Gaussian, bilateral, Kuwahara, and anisotropic diffusion. The Random Forest algorithm used 200 trees with 2 features randomly selected per node. The optimized auto-segmentation algorithm resulted in 16 image features including features derived from maximum, mean, variance Gaussian and Kuwahara filters. Dice similarity coefficient (DSC) calculations between manually segmented and Random Forest algorithm segmented images from 21 patient image sections, were analyzed. The automated algorithm produced segmentation of seven material classes with a median DSC of 0.86  ±  0.03 for pediatric patient protocols, and 0.85  ±  0.04 for adult patient protocols. Additionally, 100 randomly selected patient examinations were segmented and analyzed, and a mean sensitivity of 0.91 (range: 0.82-0.98), specificity of 0.89 (range: 0.70-0.98), and accuracy of 0.90 (range: 0.76-0.98) were demonstrated. In this study, we demonstrate that this fully automated segmentation tool was able to produce fast and accurate segmentation of the neck and trunk of the body over a wide range of patient habitus and scan parameters.

Semi-automatic classification of textures in thoracic CT scans

The textural patterns in the lung parenchyma, as visible on computed tomography (CT) scans, are essential to make a correct diagnosis in interstitial lung disease. We developed one automatic and two interactive protocols for classification of normal and seven types of abnormal lung textures. Lungs were segmented and subdivided into volumes of interest (VOIs) with homogeneous texture using a clustering approach. In the automatic protocol, VOIs were classified automatically by an extra-trees classifier that was trained using annotations of VOIs from other CT scans. In the interactive protocols, an observer iteratively trained an extra-trees classifier to distinguish the different textures, by correcting mistakes the classifier makes in a slice-by-slice manner. The difference between the two interactive methods was whether or not training data from previously annotated scans was used in classification of the first slice. The protocols were compared in terms of the percentages of VOIs that observers needed to relabel. Validation experiments were carried out using software that simulated observer behavior. In the automatic classification protocol, observers needed to relabel on average 58% of the VOIs. During interactive annotation without the use of previous training data, the average percentage of relabeled VOIs decreased from 64% for the first slice to 13% for the second half of the scan. Overall, 21% of the VOIs were relabeled. When previous training data was available, the average overall percentage of VOIs requiring relabeling was 20%, decreasing from 56% in the first slice to 13% in the second half of the scan.

Inter-observer Variability Analysis of Automatic Lung Delineation in Normal and Disease Patients

Human interaction has become almost mandatory for an automated medical system wishing to be accepted by clinical regulatory agencies such as Food and Drug Administration. Since this interaction causes variability in the gathered data, the inter-observer and intra-observer variability must be analyzed in order to validate the accuracy of the system. This study focuses on the variability from different observers that interact with an automated lung delineation system that relies on human interaction in the form of delineation of the lung borders. The database consists of High Resolution Computed Tomography (HRCT): 15 normal and 81 diseased patients' images taken retrospectively at five levels per patient. Three observers manually delineated the lungs borders independently and using software called ImgTracer™ (AtheroPoint™, Roseville, CA, USA) to delineate the lung boundaries in all five levels of 3-D lung volume. The three observers consisted of Observer-1: lesser experienced novice tracer who is a resident in radiology under the guidance of radiologist, whereas Observer-2 and Observer-3 are lung image scientists trained by lung radiologist and biomedical imaging scientist and experts. The inter-observer variability can be shown by comparing each observer's tracings to the automated delineation and also by comparing each manual tracing of the observers with one another. The normality of the tracings was tested using D'Agostino-Pearson test and all observers tracings showed a normal P-value higher than 0.05. The analysis of variance (ANOVA) test between three observers and automated showed a P-value higher than 0.89 and 0.81 for the right lung (RL) and left lung (LL), respectively. The performance of the automated system was evaluated using Dice Similarity Coefficient (DSC), Jaccard Index (JI) and Hausdorff (HD) Distance measures. Although, Observer-1 has lesser experience compared to Obsever-2 and Obsever-3, the Observer Deterioration Factor (ODF) shows that Observer-1 has less than 10% difference compared to the other two, which is under acceptable range as per our analysis. To compare between observers, this study used regression plots, Bland-Altman plots, two tailed T-test, Mann-Whiney, Chi-Squared tests which showed the following P-values for RL and LL: (i) Observer-1 and Observer-3 were: 0.55, 0.48, 0.29 for RL and 0.55, 0.59, 0.29 for LL; (ii) Observer-1 and Observer-2 were: 0.57, 0.50, 0.29 for RL and 0.54, 0.59, 0.29 for LL; (iii) Observer-2 and Observer-3 were: 0.98, 0.99, 0.29 for RL and 0.99, 0.99, 0.29 for LL. Further, CC and R-squared coefficients were computed between observers which came out to be 0.9 for RL and LL. All three observers however manage to show the feature that diseased lungs are smaller than normal lungs in terms of area.

Segmentation and Image Analysis of Abnormal Lungs at CT: Current Approaches, Challenges, and Future Trends

The computer-based process of identifying the boundaries of lung from surrounding thoracic tissue on computed tomographic (CT) images, which is called segmentation, is a vital first step in radiologic pulmonary image analysis. Many algorithms and software platforms provide image segmentation routines for quantification of lung abnormalities; however, nearly all of the current image segmentation approaches apply well only if the lungs exhibit minimal or no pathologic conditions. When moderate to high amounts of disease or abnormalities with a challenging shape or appearance exist in the lungs, computer-aided detection systems may be highly likely to fail to depict those abnormal regions because of inaccurate segmentation methods. In particular, abnormalities such as pleural effusions, consolidations, and masses often cause inaccurate lung segmentation, which greatly limits the use of image processing methods in clinical and research contexts. In this review, a critical summary of the current methods for lung segmentation on CT images is provided, with special emphasis on the accuracy and performance of the methods in cases with abnormalities and cases with exemplary pathologic findings. The currently available segmentation methods can be divided into five major classes: (a) thresholding-based, (b) region-based, (c) shape-based, (d) neighboring anatomy-guided, and (e) machine learning-based methods. The feasibility of each class and its shortcomings are explained and illustrated with the most common lung abnormalities observed on CT images. In an overview, practical applications and evolving technologies combining the presented approaches for the practicing radiologist are detailed.

Novel Logistic Regression Model of Chest CT Attenuation Coefficient Distributions for the Automated Detection of Abnormal (Emphysema or ILD) Versus Normal Lung

RATIONALE AND OBJECTIVES

We evaluated the role of automated quantitative computed tomography (CT) scan interpretation algorithm in detecting interstitial lung disease (ILD) and/or emphysema in a sample of elderly subjects with mild lung disease. We hypothesized that the quantification and distributions of CT attenuation values on lung CT, over a subset of Hounsfield units (HUs) range (-1000 HU, 0 HU), can differentiate early or mild disease from normal lung.

MATERIALS AND METHODS

We compared the results of quantitative spiral rapid end-exhalation (functional residual capacity, FRC) and end-inhalation (total lung capacity, TLC) CT scan analyses of 52 subjects with radiographic evidence of mild fibrotic lung disease to the results of 17 normal subjects. Several CT value distributions were explored, including (1) that from the peripheral lung taken at TLC (with peels at 15 or 65 mm), (2) the ratio of (1) to that from the core of lung, and (3) the ratio of (2) to its FRC counterpart. We developed a fused-lasso logistic regression model that can automatically identify sub-intervals of -1000 HU and 0 HU over which a CT value distribution provides optimal discrimination between abnormal and normal scans.

RESULTS

The fused-lasso logistic regression model based on (2) with 15-mm peel identified the relative frequency of CT values of over -1000 HU and -900 and those over -450 HU and -200 HU as a means of discriminating abnormal versus normal lung, resulting in a zero out-sample false-positive rate, and 15% false-negative rate of that was lowered to 12% by pooling information.

CONCLUSIONS

We demonstrated the potential usefulness of this novel quantitative imaging analysis method in discriminating ILD and/or emphysema from normal lungs.

Estimation of vocal fold plane in 3D CT images for diagnosis of vocal fold abnormalities

Vocal folds are the key body structures that are responsible for phonation and regulating air movement into and out of lungs. Various vocal fold disorders may seriously impact the quality of life. When diagnosing vocal fold disorders, CT of the neck is the commonly used imaging method. However, vocal folds do not align with the normal axial plane of a neck and the plane containing vocal cords and arytenoids does vary during phonation. It is therefore important to generate an algorithm for detecting the actual plane containing vocal folds. In this paper, we propose a method to automatically estimate the vocal fold plane using vertebral column and anterior commissure localization. Gray-level thresholding, connected component analysis, rule based segmentation and unsupervised k-means clustering were used in the proposed algorithm. The anterior commissure segmentation method achieved an accuracy of 85%, a good estimate of the expert assessment.

Neurotransmitter vesicle movement dynamics in living neurons

The communication between two neurons is established by endogenous chemical particles aggregated in vesicles that move along the axons. It is known that an abnormal transport of these vesicles is correlated with neurodegenerative diseases. The quantification of the dynamics of vesicles movement can therefore be a window to study early detection of such diseases. Nevertheless, most of the studies in the literature rely on manual tracking techniques. In this paper we present a novel methodology for quantifying neurotransmitter vesicle dynamics by using a combination of image tracking and classification algorithms. We use confocal microscopy videos of living neurons to detect and classify vesicles using support vector machine (SVM), while motion is extracted via global nearest neighbor (GNN) tracking approach. Results of the classification algorithm are presented and compared to a ground truth dataset defined by experts. Sensitivity of 90% and specificity of 97% were obtained at a much lower computational cost than an established method from the literature (0.24s/frame vs. 125s/frame). These preliminary results suggest the great potential of the method and tool we have been developing for single particle movement dynamics measure in living cells.

Automatic segmentation of the rima glottidis in 4D laryngeal CT scans in Parkinson's disease

Parkinson's disease (PD) is a progressive, incurable neuro-degenerative disease. Symptoms appear when approximately 70% of mid-brain dopaminergic neurons have died. Temporal analysis of the calculated area of the rima glottidis may give an indication of vocal impairment. In this paper, we present an automatic segmentation algorithm to segment the rima glottidis from 4D CT images using texture features and support vector machines (SVM). Automatic two dimensional region growing is then applied as a post processing step to segment the area accurately. The proposed segmentation algorithm resulted in accurate segmentation and we demonstrate a high correlation between the manually segmented area and automatic segmentation.

Lung Lesion Extraction Using a Toboggan Based Growing Automatic Segmentation Approach

The accurate segmentation of lung lesions from computed tomography (CT) scans is important for lung cancer research and can offer valuable information for clinical diagnosis and treatment. However, it is challenging to achieve a fully automatic lesion detection and segmentation with acceptable accuracy due to the heterogeneity of lung lesions. Here, we propose a novel toboggan based growing automatic segmentation approach (TBGA) with a three-step framework, which are automatic initial seed point selection, multi-constraints 3D lesion extraction and the final lesion refinement. The new approach does not require any human interaction or training dataset for lesion detection, yet it can provide a high lesion detection sensitivity (96.35%) and a comparable segmentation accuracy with manual segmentation (P > 0.05), which was proved by a series assessments using the LIDC-IDRI dataset (850 lesions) and in-house clinical dataset (121 lesions). We also compared TBGA with commonly used level set and skeleton graph cut methods, respectively. The results indicated a significant improvement of segmentation accuracy . Furthermore, the average time consumption for one lesion segmentation was under 8 s using our new method. In conclusion, we believe that the novel TBGA can achieve robust, efficient and accurate lung lesion segmentation in CT images automatically.

Airway Segmentation and Centerline Extraction from Thoracic CT - Comparison of a New Method to State of the Art Commercialized Methods

INTRODUCTION

Our motivation is increased bronchoscopic diagnostic yield and optimized preparation, for navigated bronchoscopy. In navigated bronchoscopy, virtual 3D airway visualization is often used to guide a bronchoscopic tool to peripheral lesions, synchronized with the real time video bronchoscopy. Visualization during navigated bronchoscopy, the segmentation time and methods, differs. Time consumption and logistics are two essential aspects that need to be optimized when integrating such technologies in the interventional room. We compared three different approaches to obtain airway centerlines and surface.

METHOD

CT lung dataset of 17 patients were processed in Mimics (Materialize, Leuven, Belgium), which provides a Basic module and a Pulmonology module (beta version) (MPM), OsiriX (Pixmeo, Geneva, Switzerland) and our Tube Segmentation Framework (TSF) method. Both MPM and TSF were evaluated with reference segmentation. Automatic and manual settings allowed us to segment the airways and obtain 3D models as well as the centrelines in all datasets. We compared the different procedures by user interactions such as number of clicks needed to process the data and quantitative measures concerning the quality of the segmentation and centrelines such as total length of the branches, number of branches, number of generations, and volume of the 3D model.

RESULTS

The TSF method was the most automatic, while the Mimics Pulmonology Module (MPM) and the Mimics Basic Module (MBM) resulted in the highest number of branches. MPM is the software which demands the least number of clicks to process the data. We found that the freely available OsiriX was less accurate compared to the other methods regarding segmentation results. However, the TSF method provided results fastest regarding number of clicks. The MPM was able to find the highest number of branches and generations. On the other hand, the TSF is fully automatic and it provides the user with both segmentation of the airways and the centerlines. Reference segmentation comparison averages and standard deviations for MPM and TSF correspond to literature.

CONCLUSION

The TSF is able to segment the airways and extract the centerlines in one single step. The number of branches found is lower for the TSF method than in Mimics. OsiriX demands the highest number of clicks to process the data, the segmentation is often sparse and extracting the centerline requires the use of another software system. Two of the software systems performed satisfactory with respect to be used in preprocessing CT images for navigated bronchoscopy, i.e. the TSF method and the MPM. According to reference segmentation both TSF and MPM are comparable with other segmentation methods. The level of automaticity and the resulting high number of branches plus the fact that both centerline and the surface of the airways were extracted, are requirements we considered particularly important. The in house method has the advantage of being an integrated part of a navigation platform for bronchoscopy, whilst the other methods can be considered preprocessing tools to a navigation system.

Supervised recursive segmentation of volumetric CT images for 3D reconstruction of lung and vessel tree

Three dimensional reconstruction of lung and vessel tree has great significance to 3D observation and quantitative analysis for lung diseases. This paper presents non-sheltered 3D models of lung and vessel tree based on a supervised semi-3D lung tissues segmentation method. A recursive strategy based on geometric active contour is proposed instead of the "coarse-to-fine" framework in existing literature to extract lung tissues from the volumetric CT slices. In this model, the segmentation of the current slice is supervised by the result of the previous one slice due to the slight changes between adjacent slice of lung tissues. Through this mechanism, lung tissues in all the slices are segmented fast and accurately. The serious problems of left and right lungs fusion, caused by partial volume effects, and segmentation of pleural nodules can be settled meanwhile during the semi-3D process. The proposed scheme is evaluated by fifteen scans, from eight healthy participants and seven participants suffering from early-stage lung tumors. The results validate the good performance of the proposed method compared with the "coarse-to-fine" framework. The segmented datasets are utilized to reconstruct the non-sheltered 3D models of lung and vessel tree.

Accurate and efficient separation of left and right lungs from 3D CT scans: A generic hysteresis approach

Separation of left and right lungs from binary segmentation is often necessary for quantitative image-based pulmonary disease evaluation. In this article, we present a new fully automated approach for accurate, robust, and efficient lung separation using 3-D CT scans. Our method follows a hysteresis setting that utilizes information from both lung regions and background gaps. First, original segmentation is separated by subtracting the gaps between left and right lungs, which are enhanced with Hessian filtering. Second, the 2-D separation manifold in 3-D image space is estimated based on the distance information from the two subsets. Finally, the separation manifold is projected back to the original segmentation in order to produce the separated lungs through optimization for addressing minor local variations. An evaluation on over 400 human and 100 small animal 3-D CT images with various abnormalities is performed. The proposed scheme successfully separated all connections on the candidate CT images. Using hysteresis mechanism, each phase is performed robustly and 3-D information is utilized to achieve a generic, efficient, and accurate solution.

Automatic Classification of Normal and Cancer Lung CT Images Using Multiscale AM-FM Features

Computer-aided diagnostic (CAD) systems provide fast and reliable diagnosis for medical images. In this paper, CAD system is proposed to analyze and automatically segment the lungs and classify each lung into normal or cancer. Using 70 different patients' lung CT dataset, Wiener filtering on the original CT images is applied firstly as a preprocessing step. Secondly, we combine histogram analysis with thresholding and morphological operations to segment the lung regions and extract each lung separately. Amplitude-Modulation Frequency-Modulation (AM-FM) method thirdly, has been used to extract features for ROIs. Then, the significant AM-FM features have been selected using Partial Least Squares Regression (PLSR) for classification step. Finally, K-nearest neighbour (KNN), support vector machine (SVM), naïve Bayes, and linear classifiers have been used with the selected AM-FM features. The performance of each classifier in terms of accuracy, sensitivity, and specificity is evaluated. The results indicate that our proposed CAD system succeeded to differentiate between normal and cancer lungs and achieved 95% accuracy in case of the linear classifier.