Search

Quantitative imaging biomarkers (QIBs) provide medical image–derived intensity, texture, shape, and size features that may help characterize cancerous tumors and predict clinical outcomes. Successful clinical translation of QIBs depends on the robustness of their measurements. Biomarkers derived from positron emission tomography images are prone to measurement errors owing to differences in image processing factors such as the tumor segmentation method used to define volumes of interest over which to calculate QIBs. We illustrate a new Bayesian statistical approach to characterize the robustness of QIBs to different processing factors. Study data consist of 22 QIBs measured on 47 head and neck tumors in 10 positron emission tomography/computed tomography scans segmented manually and with semiautomated methods used by 7 institutional members of the NCI Quantitative Imaging Network. QIB performance is estimated and compared across institutions with respect to measurement errors and power to recover statistical associations with clinical outcomes. Analysis findings summarize the performance impact of different segmentation methods used by Quantitative Imaging Network members. Robustness of some advanced biomarkers was found to be similar to conventional markers, such as maximum standardized uptake value. Such similarities support current pursuits to better characterize disease and predict outcomes by developing QIBs that use more imaging information and are robust to different processing factors. Nevertheless, to ensure reproducibility of QIB measurements and measures of association with clinical outcomes, errors owing to segmentation methods need to be reduced.

Fractal biological structures are pervasive throughout the plant and animal kingdoms, with the mammalian lung being a quintessential example. The lung airway and vascular trees are generated during embryogenesis from a small set of building codes similar to Turing mechanisms that create robust trees ideally suited to their functions. Whereas the blood flow pattern generated by these fractal trees has been shown to be genetically determined, the geometry of the trees has not. We explored a newly established repository providing high-resolution bronchial trees from the four most commonly studied laboratory mice (B6C3F1, BALB/c, C57BL/6 and CD-1). The data fit a fractal model well for all animals with the fractal dimensions ranging from 1.54 to 1.67, indicating that the conducting airway of mice can be considered a self-similar and space-filling structure. We determined that the fractal dimensions of these airway trees differed by strain but not sex, reinforcing the concept that airway branching patterns are encoded within the DNA. The observations also highlight that future study design and interpretations may need to consider differences in airway geometry between mouse strains. NEW & NOTEWORTHY Similar to larger mammals such as humans, the geometries of the bronchial tree in mice are fractal structures that have repeating patterns from the trachea to the terminal branches. The airway geometries of the four most commonly studied mice are different and need to be considered when comparing results that employ different mouse strains. This variability in mouse airway geometries should be incorporated into computer models exploring toxicology and aerosol deposition in mouse models.

To facilitate computational toxicology, we developed an approach for generating high-resolution lung-anatomy and particle-deposition mouse models. Major processing steps of our method include mouse preparation, serial block-face cryomicrotome imaging, and highly automated image analysis for generating three-dimensional (3D) mesh-based models and volume-based models of lung anatomy (airways, lobes, sublobes, and near-acini structures) that are linked to local particle-deposition measurements. Analysis resulted in 34 mouse models covering 4 different mouse strains (B6C3F1: 8, BALB/C: 11, C57Bl/6: 8, and CD-1: 7) as well as both sexes (16 male and 18 female) and different particle sizes [2 μm ( n = 15), 1 μm ( n = 16), and 0.5 μm ( n = 3)]. On average, resulting mouse airway models had 1,616.9 ± 298.1 segments, a centerline length of 597.6 ± 59.8 mm, and 1,968.9 ± 296.3 outlet regions. In addition to 3D geometric lung models, matching detailed relative particle-deposition measurements are provided. All data sets are available online in the lapdMouse archive for download. The presented approach enables linking relative particle deposition to anatomical structures like airways. This will in turn improve the understanding of site-specific airflows and how they affect drug, environmental, or biological aerosol deposition. NEW & NOTEWORTHY Computer simulations of particle deposition in mouse lungs play an important role in computational toxicology. Until now, a limiting factor was the lack of high-resolution mouse lung models and measured local particle-deposition information, which are required for developing accurate modeling approaches (e.g., computational fluid dynamics). With the developed imaging and analysis approach, we address this issue and provide all of the raw and processed data in a publicly accessible repository.

Radiation therapy (RT), along with surgery and medical therapies, are the fundamental methods used to treat cancers, as well as a wide range of other diseases. RT is delivered in multiple forms, including external beam therapy, brachytherapy, and radiopharmaceutical therapy. RT is a completely image-guided treatment paradigm, and benefits from advances made in quantitative imaging (QI). The therapeutic effects of radiation are proven, and have improved with each advance, enabling more precisely delivered radiation dose to a tumor target and avoidance of normal tissues. Advances in QI enable improved target and normal tissue definitions, and advances in computer-based algorithmic tools enable enhanced consistency, efficiency, and depth in utilization of the rich information within QI. RT benefits from application of these enhanced tools to imaging to ultimately improve therapy.

Radiomics is an image analysis approach for extracting large amounts of quantitative information from medical images using a variety of computational methods. Our goal was to evaluate the utility of radiomic feature analysis from 18F-fluorothymidine positron emission tomography (FLT PET) obtained at baseline in prediction of treatment response in patients with head and neck cancer. Thirty patients with advanced-stage oropharyngeal or laryngeal cancer, treated with definitive chemoradiation therapy, underwent FLT PET imaging before treatment. In total, 377 radiomic features of FLT uptake and feature variants were extracted from volumes of interest, these features variants were defined by either the primary tumor or the total lesion burden, which consisted of the primary tumor and all FLT-avid nodes. Feature variants included normalized measurements of uptake, which were calculated by dividing lesion uptake values by the mean uptake value in the bone marrow. Feature reduction was performed using clustering to remove redundancy, leaving 172 representative features. Effects of these features on progression-free survival were modeled with Cox regression and P-values corrected for multiple comparisons. In total, 9 features were considered significant. Our results suggest that smaller, more homogenous lesions at baseline were associated with better prognosis. In addition, features extracted from total lesion burden had a higher concordance index than primary tumor features for 8 of the 9 significant features. Furthermore, total lesion burden features showed lower interobserver variability.

Medical imaging is utilized in a wide range of clinical applications. To enable a detailed quantitative analysis, medical images must often be segmented to label (delineate) structures of interest; for example, a tumor. Frequently, manual segmentation is utilized in clinical practice (e.g., in radiation oncology) to define such structures of interest. However, it can be quite time consuming and subject to substantial between-, and within-reader variability. A more reproducible, less variable, and more time efficient segmentation approach is likely to improve medical treatment. This potential has spurred the development of segmentation algorithms which harness computational power. Segmentation algorithms’ widespread use is limited due to difficulty in quantifying their performance relative to manual segmentation, which itself is subject to variation. This paper presents a statistical model which simultaneously estimates segmentation method accuracy, and between- and within-reader variability. The model is simultaneously fit for multiple segmentation methods within a unified Bayesian framework. The Bayesian model is compared to other methods used in literature via a simulation study, and application to head and neck cancer PET/CT data. The modeling framework is flexible and can be employed in numerous comparison applications. Several alternate applications are discussed in the paper.

Purpose Quality control plays an increasingly important role in quantitative PET imaging and is typically performed using phantoms. The purpose of this work was to develop and validate a fully-automated analysis method for two common PET/CT quality assurance phantoms: the NEMA NU-2 IQ and SNMMI/CTN oncology phantom. The algorithm was designed to only utilize the PET scan to enable the analysis of phantoms with thin-walled inserts. Methods We introduce a model-based method for automated analysis of phantoms with spherical inserts. Models are first constructed for each type of phantom to be analyzed. A robust insert detection algorithm uses the model to locate all inserts inside the phantom. First, candidates for inserts are detected using a scale-space detection approach. Second, candidates are given an initial label using a score-based optimization algorithm. Third, a robust model fitting step aligns the phantom model to the initial labeling and fixes incorrect labels. Finally, the detected insert locations are refined and measurements are taken for each insert and several background regions. In addition, an approach for automated selection of NEMA and CTN phantom models is presented. The method was evaluated on a diverse set of 15 NEMA and 20 CTN phantom PET/CT scans. NEMA phantoms were filled with radioactive tracer solution at 9.7:1 activity ratio over background, and CTN phantoms were filled with 4:1 and 2:1 activity ratio over background. For quantitative evaluation, an independent reference standard was generated by two experts using PET/CT scans of the phantoms. In addition, the automated approach was compared against manual analysis, which represents the current clinical standard approach, of the PET phantom scans by four experts. Results The automated analysis method successfully detected and measured all inserts in all test phantom scans. It is a deterministic algorithm (zero variability), and the insert detection RMS error (i.e., bias) was 0.97, 1.12, and 1.48 mm for phantom activity ratios 9.7:1, 4:1, and 2:1, respectively. For all phantoms and at all contrast ratios, the average RMS error was found to be significantly lower for the proposed automated method compared to the manual analysis of the phantom scans. The uptake measurements produced by the automated method showed high correlation with the independent reference standard (R2 ≥ 0.9987). In addition, the average computing time for the automated method was 30.6 seconds and was found to be significantly lower (p ≪ 0.001) compared to manual analysis (mean: 247.8 seconds). Conclusions The proposed automated approach was found to have less error when measured against the independent reference than the manual approach. It can be easily adapted to other phantoms with spherical inserts. In addition, it eliminates inter- and intra-operator variability in PET phantom analysis and is significantly more time efficient, and therefore, represents a promising approach to facilitate and simplify PET standardization and harmonization efforts. This article is protected by copyright. All rights reserved.

Purpose The purpose of this work was to assess the potential of deep convolutional neural networks in automated measurement of cerebellum tracer uptake in F-18 fluorodeoxyglucose (FDG) positron emission tomography (PET) scans. Methods Three different three-dimensional (3D) convolutional neural network architectures (U-Net, V-Net, and modified U-Net) were implemented and compared regarding their performance in 3D cerebellum segmentation in FDG PET scans. For network training and testing, 134 PET scans with corresponding manual volumetric segmentations were utilized. For segmentation performance assessment, a fivefold cross-validation was used, and the Dice coefficient as well as signed and unsigned distance errors were calculated. In addition, standardized uptake value (SUV) uptake measurement performance was assessed by means of a statistical comparison to an independent reference standard. Furthermore, a comparison to a previously reported active-shape-model-based approach was performed. Results Out of the three convolutional neural networks investigated, the modified U-Net showed significantly better segmentation performance. It achieved a Dice coefficient of 0.911 ± 0.026, a signed distance error of 0.220 ± 0.103 mm, and an unsigned distance error of 1.048 ± 0.340 mm. When compared to the independent reference standard, SUV uptake measurements produced with the modified U-Net showed no significant error in slope and intercept. The estimated reduction in total SUV measurement error was 95.1%. Conclusions The presented work demonstrates the potential of deep convolutional neural networks in automated SUV measurement of reference regions. While it focuses on the cerebellum, utilized methods can be generalized to other reference regions like the liver or aortic arch. Future work will focus on combining lesion and reference region analysis into one approach.

Background In the era of precision oncology and publicly available datasets, the amount of information available for each patient case has dramatically increased. From clinical variables and PET-CT radiomics measures to DNA-variant and RNA expression profiles, such a wide variety of data presents a multitude of challenges. Large clinical datasets are subject to sparsely and/or inconsistently populated fields. Corresponding sequencing profiles can suffer from the problem of high-dimensionality, where making useful inferences can be difficult without correspondingly large numbers of instances. In this paper we report a novel deployment of machine learning techniques to handle data sparsity and high dimensionality, while evaluating potential biomarkers in the form of unsupervised transformations of RNA data. We apply preprocessing, MICE imputation, and sparse principal component analysis (SPCA) to improve the usability of more than 500 patient cases from the TCGA-HNSC dataset for enhancing future oncological decision support for Head and Neck Squamous Cell Carcinoma (HNSCC). Results Imputation was shown to improve prognostic ability of sparse clinical treatment variables. SPCA transformation of RNA expression variables reduced runtime for RNA-based models, though changes to classifier performance were not significant. Gene ontology enrichment analysis of gene sets associated with individual sparse principal components (SPCs) are also reported, showing that both high- and low-importance SPCs were associated with cell death pathways, though the high-importance gene sets were found to be associated with a wider variety of cancer-related biological processes. Conclusions MICE imputation allowed us to impute missing values for clinically informative features, improving their overall importance for predicting two-year recurrence-free survival by incorporating variance from other clinical variables. Dimensionality reduction of RNA expression profiles via SPCA reduced both computation cost and model training/evaluation time without affecting classifier performance, allowing researchers to obtain experimental results much more quickly. SPCA simultaneously provided a convenient avenue for consideration of biological context via gene ontology enrichment analysis.

Robust lung segmentation is challenging, especially when tens of thousands of lung CT scans need to be processed, as required by large multi-center studies. The goal of this work was to develop and assess a method for the fusion of segmentation results from two different methods to generate lung segmentations that have a lower failure rate than individual input segmentations. As basis for the fusion approach, lung segmentations generated with a region growing and model-based approach were utilized. The fusion result was generated by comparing input segmentations and selectively combining them using a trained classification system. The method was evaluated on a diverse set of 204 CT scans of normal and diseased lungs. The fusion approach resulted in a Dice coefficient of 0.9855 ± 0.0106 and showed a statistically significant improvement compared to both input segmentation methods. In addition, the failure rate at different segmentation accuracy levels was assessed. For example, when requiring that lung segmentations must have a Dice coefficient of better than 0.97, the fusion approach had a failure rate of 6.13%. In contrast, the failure rate for region growing and model-based methods was 18.14% and 15.69%, respectively. Therefore, the proposed method improves the quality of the lung segmentations, which is important for subsequent quantitative analysis of lungs. Also, to enable a comparison with other methods, results on the LOLA11 challenge test set are reported.

Purpose: The purpose of this work was to develop, validate, and compare a highly computer-aided method for the segmentation of hot lesions in head and neck 18F-FDG PET scans. Methods: A semiautomated segmentation method was developed, which transforms the segmentation problem into a graph-based optimization problem. For this purpose, a graph structure around a user-provided approximate lesion centerpoint is constructed and a suitable cost function is derived based on local image statistics. To handle frequently occurring situations that are ambiguous (e.g., lesions adjacent to each other versus lesion with inhomogeneous uptake), several segmentation modes are introduced that adapt the behavior of the base algorithm accordingly. In addition, the authors present approaches for the efficient interactive local and global refinement of initial segmentations that are based on the “just-enough-interaction” principle. For method validation, 60 PET/CT scans from 59 different subjects with 230 head and neck lesions were utilized. All patients had squamous cell carcinoma of the head and neck. A detailed comparison with the current clinically relevant standard manual segmentation approach was performed based on 2760 segmentations produced by three experts. Results: Segmentation accuracy measured by the Dice coefficient of the proposed semiautomated and standard manual segmentation approach was 0.766 and 0.764, respectively. This difference was not statistically significant (p = 0.2145). However, the intra- and interoperator standard deviations were significantly lower for the semiautomated method. In addition, the proposed method was found to be significantly faster and resulted in significantly higher intra- and interoperator segmentation agreement when compared to the manual segmentation approach. Conclusions: Lack of consistency in tumor definition is a critical barrier for radiation treatment targeting as well as for response assessment in clinical trials and in clinical oncology decision-making. The properties of the authors approach make it well suited for applications in image-guided radiation oncology, response assessment, or treatment outcome prediction.

Chest wall strapping (CWS) induces breathing at low lung volumes but also increases parenchymal elastic recoil. In this study, we tested the hypothesis that CWS dilates airways via airway-parenchymal interdependence. In 11 subjects (6 healthy and 5 with mild to moderate COPD), pulmonary function tests and lung volumes were obtained in control (baseline) and the CWS state. Control and CWS-CT scans were obtained at 50% of control (baseline) total lung-capacity (TLC). CT lung volumes were analyzed by CT volumetry. If control and CWS-CT volumetry did not differ by more than 25%, airway dimensions were analyzed via automated airway segmentation. CWS-TLC was reduced on average to 71% of control-TLC in normal subjects and 79% of control-TLC in subjects with COPD. CWS increased expiratory airflow at 50% of control-TLC by 41% (3.50 ± 1.6 vs. 4.93 ± 1.9 l/s, P = 0.04) in normals and 316% in COPD(0.25 ± 0.05 vs 0.79 ± 0.39 l/s, P = 0.04). In 10 subjects (5 normals and 5 COPD), control and CWS-CT scans at 50% control-TLC did not differ more than 25% on CT volumetry and were included in the airway structure analysis. CWS increased the mean number of detectable airways with a diameter of ≤2 mm by 32.5% (65 ± 10 vs. 86 ± 124, P = 0.01) in normal subjects and by 79% (59 ± 19 vs. 104 ± 16, P = 0.01) in subjects with COPD. There was no difference in the number of detectable airways with diameters 2–4 mm and >4 mm in normal or in COPD subjects. In conclusion, CWS enhances the detection of small airways via automated CT airway segmentation and increases expiratory airflow in normal subjects as well as in subjects with mild to moderate COPD. NEW & NOTEWORTHY In normal and COPD subjects, chest wall strapping(CWS) increased the number of detectable small airways using automated CT airway segmentation. The concept of dysanapsis expresses the physiological variation in the geometry of the tracheobronchial tree and lung parenchyma based on development. We propose a dynamic concept to dysanapsis in which CWS leads to breathing at lower lung volumes with a corresponding increase in the size of small airways, a potentially novel, nonpharmacological treatment for COPD.

Quantitative biomarkers derived from medical images are being used increasingly to help diagnose disease, guide treatment, and predict clinical outcomes. Measurement of quantitative imaging biomarkers is subject to bias and variability from multiple sources, including the scanner technologies that produce images, the approaches for identifying regions of interest in images, and the algorithms that calculate biomarkers from regions. Moreover, these sources may differ within and between the quantification methods employed by institutions, thus making it difficult to develop and implement multi-institutional standards. We present a Bayesian framework for assessing bias and variability in imaging biomarkers derived from different quantification methods, comparing agreement to a reference standard, studying prognostic performance, and estimating sample size for future clinical studies. The statistical methods are illustrated with data obtained from a positron emission tomography challenge conducted by members of the NCI's Quantitative Imaging Network program, in which tumor volumes were measured manually and with seven different semi-automated segmentation algorithms. Estimates and comparisons of bias and variability in the resulting measurements are provided along with an R software package for the technical performance analysis and an online web application for sample size and power analysis.

We present a graph-based framework for airway tree reconstruction from computerized tomography (CT) scans and evaluate the performance of different feature categories and their combinations on five lung cohorts. The approach consists of two main processing steps. First, potential airway branch and connection candidates are identified and represented by a graph structure with weighted nodes and edges, respectively. Second, an optimization algorithm is utilized for generating an airway detection result by selecting a subset of airway branches and connections based on graph weights derived from image features. The performance of the algorithm with different feature categories and their combinations was assessed on a set of 50 lung CT scans from five different cohorts, including normal and diseased lungs. Results show trade-offs between feature categories/combinations in terms of correctly (true positive) and incorrectly (false positive) identified airways. Also, the performance of features in dependence of lung cohort was analyzed. Across all cohorts, a good trade-off with high true positive rate (TPR) and low false positive rate (FPR) was achieved by a combination of gray-value, local shape, and structural features. This combination enabled extracting 91.80% of reference airways (TPR) in combination with a low FPR of 1.00%. In addition, this variant was evaluated on the public EXACT'09 test set, and a comparison with other airway detection approaches is provided. One of the main advantages of the presented method is that it is robust against local disturbances/artifacts or other ambiguities that are frequently occurring in lung CT scans.

Segmentation of pulmonary lobes in inspiration and expiration chest CT scan pairs is an important prerequisite for lobe-based quantitative disease assessment. Conventional methods process each CT scan independently, resulting typically in lower segmentation performance at expiration compared to inspiration. To address this issue, we present an approach, which utilizes CT scans at both respiratory states. It consists of two main parts: a base method that processes a single CT scan and an extended method that utilizes the segmentation result obtained on the inspiration scan as a subject-specific prior for segmentation of the expiration scan. We evaluated the methods on a diverse set of 40 CT scan pairs. In addition, we compare the performance of our method to a registration-based approach. On inspiration scans, the base method achieved an average distance error of 0.59, 0.64, and 0.91 mm for the left oblique, right oblique, and right horizontal fissures, respectively, when compared with expert-based reference tracings. On expiration scans, the base method's errors were 1.54, 3.24, and 3.34 mm, respectively. In comparison, utilizing proposed subject-specific priors for segmentation of expiration scans allowed decreasing average distance errors to 0.82, 0.79, and 1.04 mm, which represents a significant improvement ([Formula: see text]) compared with all other methods investigated.

A highly automated method for the segmentation of airways in the serial block-face cryomicrotome images of rat lungs is presented. First, a point inside of the trachea is manually specified. Then, a set of candidate airway centerline points is automatically identified. By utilizing a novel path extraction method, a centerline path between the root of the airway tree and each point in the set of candidate centerline points is obtained. Local disturbances are robustly handled by a novel path extraction approach, which avoids the shortcut problem of standard minimum cost path algorithms. The union of all centerline paths is utilized to generate an initial airway tree structure, and a pruning algorithm is applied to automatically remove erroneous subtrees or branches. Finally, a surface segmentation method is used to obtain the airway lumen. The method was validated on five image volumes of Sprague-Dawley rats. Based on an expert-generated independent standard, an assessment of airway identification and lumen segmentation performance was conducted. The average of airway detection sensitivity was 87.4% with a 95% confidence interval (CI) of (84.9, 88.6)%. A plot of sensitivity as a function of airway radius is provided. The combined estimate of airway detection specificity was 100% with a 95% CI of (99.4, 100)%. The average number and diameter of terminal airway branches was 1179 and 159 μm, respectively. Segmentation results include airways up to 31 generations. The regression intercept and slope of airway radius measurements derived from final segmentations were estimated to be 7.22 μm and 1.005, respectively. The developed approach enables the quantitative studies of physiology and lung diseases in rats, requiring detailed geometric airway models.

Introduction18 F-fluorodeoxyglucose (FDG) positron emission tomography (PET) is now a standard diagnostic imaging test performed in patients with head and neck cancer for staging, re-staging, radiotherapy planning, and outcome assessment. Currently, quantitative analysis of FDG PET scans is limited to simple metrics like maximum standardized uptake value, metabolic tumor volume, or total lesion glycolysis, which have limited predictive value. The goal of this work was to assess the predictive potential of new (i.e., nonstandard) quantitative imaging features on head and neck cancer outcome.MethodsThis retrospective study analyzed fifty-eight pre- and post-treatment FDG PET scans of patients with head and neck squamous cell cancer to calculate five standard and seventeen new features at baseline and post-treatment. Cox survival regression was used to assess the predictive potential of each quantitative imaging feature on disease-free survival.ResultsAnalysis showed that the post-treatment change of the average tracer uptake in the rim background region immediately adjacent to the tumor normalized by uptake in the liver represents a novel PET feature that is associated with disease-free survival (HR 1.95; 95% CI 1.27, 2.99) and has good discriminative performance (c index 0.791).ConclusionThe reported findings define a promising new direction for quantitative imaging biomarker research in head and neck squamous cell cancer and highlight the potential role of new radiomics features in oncology decision making as part of precision medicine.

Air trapping and airflow obstruction are being increasingly identified in infants with cystic fibrosis. These findings are commonly attributed to airway infection, inflammation, and mucus buildup.To learn if air trapping and airflow obstruction are present before the onset of airway infection and inflammation in cystic fibrosis.On the day they are born, piglets with cystic fibrosis lack airway infection and inflammation. Therefore, we used newborn wild-type piglets and piglets with cystic fibrosis to assess air trapping, airway size, and lung volume with inspiratory and expiratory X-ray computed tomography scans. Micro-computed tomography scanning was used to assess more distal airway sizes. Airway resistance was determined with a mechanical ventilator. Mean linear intercept and alveolar surface area were determined using stereologic methods.On the day they were born, piglets with cystic fibrosis exhibited air trapping more frequently than wild-type piglets (75% vs. 12.5%, respectively). Moreover, newborn piglets with cystic fibrosis had increased airway resistance that was accompanied by luminal size reduction in the trachea, mainstem bronchi, and proximal airways. In contrast, mean linear intercept length, alveolar surface area, and lung volume were similar between both genotypes.The presence of air trapping, airflow obstruction, and airway size reduction in newborn piglets with cystic fibrosis before the onset of airway infection, inflammation, and mucus accumulation indicates that cystic fibrosis impacts airway development. Our findings suggest that early airflow obstruction and air trapping in infants with cystic fibrosis might, in part, be caused by congenital airway abnormalities.

Recently, the optimal surface finding (OSF) and layered optimal graph image segmentation of multiple objects and surfaces (LOGISMOS) approaches have been reported with applications to medical image segmentation tasks. While providing high levels of performance, these approaches may locally fail in the presence of pathology or other local challenges. Due to the image data variability, finding a suitable cost function that would be applicable to all image locations may not be feasible. This paper presents a new interactive refinement approach for correcting local segmentation errors in the automated OSF-based segmentation. A hybrid desktop/virtual reality user interface was developed for efficient interaction with the segmentations utilizing state-of-the-art stereoscopic visualization technology and advanced interaction techniques. The user interface allows a natural and interactive manipulation of 3-D surfaces. The approach was evaluated on 30 test cases from 18 CT lung datasets, which showed local segmentation errors after employing an automated OSF-based lung segmentation. The performed experiments exhibited significant increase in performance in terms of mean absolute surface distance errors (2.54±0.75mm prior to refinement vs. 1.11±0.43mm post-refinement, p ≪0.001). Speed of the interactions is one of the most important aspects leading to the acceptance or rejection of the approach by users expecting real-time interaction experience. The average algorithm computing time per refinement iteration was 150ms, and the average total user interaction time required for reaching complete operator satisfaction was about 2min per case. This time was mostly spent on human-controlled manipulation of the object to identify whether additional refinement was necessary and to approve the final segmentation result. The reported principle is generally applicable to segmentation problems beyond lung segmentation in CT scans as long as the underlying segmentation utilizes the OSF framework. The two reported segmentation refinement tools were optimized for lung segmentation and might need some adaptation for other application domains.

We present a highly automated approach to obtain detailed structural models of airway trees from ex vivo porcine lung tissue imaged with a high resolution micro-CT scanner. Such information is an important prerequisite to systematically study models of lung disease that affect airway morphology. The method initially identifies all tubular airway-like structures in the lung. In a second processing step, these structures are grouped into a connected airway tree by utilizing prior knowledge about the airway trees branching pattern. The method was evaluated on 12 micro-CT scans from four tracheal lobes of piglets imaged at three different inflation levels. For this study, two control piglets and two cystic fibrosis piglets were used. For systematic validation of our approach, an airway nomenclature was developed for the pig airway tree. Out of more than 3500 airway tree segments assessed during evaluation, 88.45% were correctly identified by the method. No false positive airway branches were found. A detailed performance analysis for different airway tree hierarchy levels, lung inflation levels and piglets with/without cystic fibrosis is presented in the paper.

Purpose Radiomics utilizes a large number of image-derived features for quantifying tumor characteristics that can in turn be correlated with response and prognosis. Unfortunately, extraction and analysis of such image-based features is subject to measurement variability and bias. The challenge for radiomics is particularly acute in Positron Emission Tomography (PET) where limited resolution, a high noise component related to the limited stochastic nature of the raw data, and the wide variety of reconstruction options confound quantitative feature metrics. Extracted feature quality is also affected by tumor segmentation methods used to define regions over which to calculate features, making it challenging to produce consistent radiomics analysis results across multiple institutions that use different segmentation algorithms in their PET image analysis. Understanding each element contributing to these inconsistencies in quantitative image feature and metric generation is paramount for ultimate utilization of these methods in multi-institutional trials and clinical oncology decision making. Methods To assess segmentation quality and consistency at the multi-institutional level, we conducted a study of seven institutional members of the National Cancer Institute Quantitative Imaging Network. For the study, members were asked to segment a common set of phantom PET scans acquired over a range of imaging conditions as well as a second set of head and neck cancer (HNC) PET scans. Segmentations were generated at each institution using their preferred approach. In addition, participants were asked to repeat segmentations with a time interval between initial and repeat segmentation. This procedure resulted in overall 806 phantom insert and 641 lesion segmentations. Subsequently, the volume was computed from the segmentations and compared to the corresponding reference volume by means of statistical analysis. Results On the two test sets (phantom and HNC PET scans), the performance of the seven segmentation approaches was as follows. On the phantom test set, the mean relative volume errors ranged from 29.9 to 87.8% of the ground truth reference volumes, and the repeat difference for each institution ranged between -36.4 to 39.9%. On the HNC test set, the mean relative volume error ranged between -50.5 to 701.5%, and the repeat difference for each institution ranged between -37.7 to 31.5%. In addition, performance measures per phantom insert/lesion size categories are given in the paper. On phantom data, regression analysis resulted in coefficient of variation (CV) components of 42.5% for scanners, 26.8% for institutional approaches, 21.1% for repeated segmentations, 14.3% for relative contrasts, 5.3% for count statistics (acquisition times), and 0.0% for repeated scans. Analysis showed that the CV components for approaches and repeated segmentations were significantly larger on the HNC test set with increases by 112.7% and 102.4%, respectively. Conclusion Analysis results underline the importance of PET scanner reconstruction harmonization and imaging protocol standardization for quantification of lesion volumes. In addition, to enable a distributed multi-site analysis of FDG PET images, harmonization of analysis approaches and operator training in combination with highly automated segmentation methods seems to be advisable. Future work will focus on quantifying the impact of segmentation variation on radiomics system performance.