Standardization in Quantitative Imaging: A Multicenter Comparison of Radiomic Features from Different Software Packages on Digital Reference Objects and Patient Data Sets.

Radiomic features are being increasingly studied for clinical applications. We aimed to assess the agreement among radiomic features when computed by several groups by using different software packages under very tightly controlled conditions, which included standardized feature definitions and common image data sets. Ten sites (9 from the NCI's Quantitative Imaging Network] positron emission tomography-computed tomography working group plus one site from outside that group) participated in this project. Nine common quantitative imaging features were selected for comparison including features that describe morphology, intensity, shape, and texture. The common image data sets were: three 3D digital reference objects (DROs) and 10 patient image scans from the Lung Image Database Consortium data set using a specific lesion in each scan. Each object (DRO or lesion) was accompanied by an already-defined volume of interest, from which the features were calculated. Feature values for each object (DRO or lesion) were reported. The coefficient of variation (CV), expressed as a percentage, was calculated across software packages for each feature on each object. Thirteen sets of results were obtained for the DROs and patient data sets. Five of the 9 features showed excellent agreement with CV < 1%; 1 feature had moderate agreement (CV < 10%), and 3 features had larger variations (CV ≥ 10%) even after attempts at harmonization of feature calculations. This work highlights the value of feature definition standardization as well as the need to further clarify definitions for some features.

Drivers, challenges and opportunities of forage technology adoption by smallholder cattle households in Cambodia.

Forage technology has been successfully introduced into smallholder cattle systems in Cambodia as an alternative feed source to the traditional rice straw and native pastures, improving animal nutrition and reducing labour requirements of feeding cattle. Previous research has highlighted the positive impacts of forage technology including improved growth rates of cattle and household time savings. However, further research is required to understand the drivers, challenges and opportunities of forage technology for smallholder cattle households in Cambodia to facilitate widespread adoption and identify areas for further improvement. A survey of forage-growing households (n = 40) in July-September 2016 examined forage technology adoption experiences, including reasons for forage establishment, use of inputs and labour requirements of forage plot maintenance and use of forages (feeding, fattening, sale of grass or seedlings and silage). Time savings was reported as the main driver of forage adoption with household members spending approximately 1 h per day maintaining forages and feeding it to cattle. Water availability was reported as the main challenge to this activity. A small number of households also reported lack of labour, lack of fencing, competition from natural grasses, cost of irrigation and lack of experience as challenges to forage growing. Cattle fattening and sale of cut forage grass and seedlings was not found to be a widespread activity by interviewed households, with 25 and 10% of households reporting use of forages for these activities, respectively. Currently, opportunities exist for these households to better utilise forages through expansion of forage plots and cattle activities, although assistance is required to support these households in addressing current constraints, particularly availability of water, if the sustainability of this feed technology for smallholder cattle household is to be established in Cambodia.

Aromatase inhibitor-induced modulation of breast density: clinical and genetic effects.

BACKGROUND: Change in breast density may predict outcome of women receiving adjuvant hormone therapy for breast cancer. We performed a prospective clinical trial to evaluate the impact of inherited variants in genes involved in oestrogen metabolism and signalling on change in mammographic percent density (MPD) with aromatase inhibitor (AI) therapy. METHODS: Postmenopausal women with breast cancer who were initiating adjuvant AI therapy were enrolled onto a multicentre, randomised clinical trial of exemestane vs letrozole, designed to identify associations between AI-induced change in MPD and single-nucleotide polymorphisms in candidate genes. Subjects underwent unilateral craniocaudal mammography before and following 24 months of treatment. RESULTS: Of the 503 enrolled subjects, 259 had both paired mammograms at baseline and following 24 months of treatment and evaluable DNA. We observed a statistically significant decrease in mean MPD from 17.1 to 15.1% (P<0.001), more pronounced in women with baseline MPD ≥20%. No AI-specific difference in change in MPD was identified. No significant associations between change in MPD and inherited genetic variants were observed. CONCLUSION: Subjects with higher baseline MPD had a greater average decrease in MPD with AI therapy. There does not appear to be a substantial effect of inherited variants in biologically selected candidate genes.

Treatment response assessment of head and neck cancers on CT using computerized volume analysis.

BACKGROUND AND PURPOSE: Head and neck cancer can cause substantial morbidity and mortality. Our aim was to evaluate the potential usefulness of a computerized system for segmenting lesions in head and neck CT scans and for estimation of volume change of head and neck malignant tumors in response to treatment. MATERIALS AND METHODS: CT scans from a pretreatment examination and a post 1-cycle chemotherapy examination of 34 patients with 34 head and neck primary-site cancers were collected. The computerized system was developed in our laboratory. It performs 3D segmentation on the basis of a level-set model and uses as input an approximate bounding box for the lesion of interest. The 34 tumors included tongue, tonsil, vallecula, supraglottic, epiglottic, and hard palate carcinomas. As a reference standard, 1 radiologist outlined full 3D contours for each of the 34 primary tumors for both the pre- and posttreatment scans and a second radiologist verified the contours. RESULTS: The correlation between the automatic and manual estimates for both the pre- to post-treatment volume change and the percentage volume change for the 34 primary-site tumors was 0.95, with an average error of -2.4 ± 8.5% by automatic segmentation. There was no substantial difference and specific trend in the automatic segmentation accuracy for the different types of primary head and neck tumors, indicating that the computerized segmentation performs relatively robustly for this application. CONCLUSIONS: The tumor size change in response to treatment can be accurately estimated by the computerized segmentation system relative to radiologists' manual estimations for different types of head and neck tumors.

Polymer planar waveguide device using inverted channel structure with upper liquid crystal cladding.

We propose a composite waveguide configuration based on an inverted polymer channel structure with upper nematic liquid crystal cladding. This configuration can achieve a more homogenous liquid crystal molecular alignment between the core and the liquid crystal material by minimizing the rubbing damage during preparation of the alignment layer. We demonstrated our idea with a variable optical attenuator which exhibited a 24 dB of attenuation range over a tuning peak voltage of 10 V at 1550 nm.

Chondrodermatitis nodularis chronica helicis in monozygotic twins.

Chondrodermatitis nodularis chronica helicis (CNCH) is a benign inflammatory nodule of the helix. Patients report severe tenderness upon pressure. Commonly seen in middle-aged men, there are no reports of this disease in twins. We report middle-aged male monozygotic twins who simultaneously developed CNCH. This suggests, but does not prove, the possibility of a hereditary factor in the pathogenesis of CNCH.

Skin decontamination of glyphosate from human skin in vitro.

This study compared three model decontaminant solutions (tap water, isotonic saline, and hypertonic saline) for their ability to remove a model herbicide (glyphosate) from an in vitro human skin model. Human cadaver skin was dosed (approximately 375microg) of [14C]-glyphosate on 3cm2 per skin. After each exposure time (1, 3, and 30min post-dosing, respectively), the surface skin was washed three times (4ml per time) with each solution. After washing, the skin was stripped twice with tape discs. Lastly, the wash solutions, strippings, receptor fluid, and remainder of skin were liquid scintillation analyzer counted to determine the amount of glyphosate. There were no statistical differences among these groups at any time points. The total mass balance recovery at three time exposure points was between 94.8% and 102.4%. The wash off rates (glyphosate in wash solutions) at three different exposure times is 79-101.2%. Thus the three tested decontaminants possess similar effectiveness in removing glyphosate from skin. This in vitro model is not only economic and rapid, but also provides quantitative data that may aid screening for optimal decontaminants.

Effects of antioxidants on pyridinoline cross-link formation in culture supernatants of fibroblasts from normal skin and hypertrophic scars.

Free radicals are normally generated in many metabolic pathways. They are closely associated with inflammatory diseases. The aim of this study was to investigate the effects of free radicals and their antioxidants on the formation of pyridinoline using human fibroblasts from normal skin and hypertrophic scars. The significance of the increase in pyridinoline cross-links is that large quantities have been found in hypertrophic scars formed post-burn than in normal skin, and that catalase was effective in reducing the pyridinoline cross-link formation in hypertrophic scars. The pyridinoline cross-link concentration expressed in nM/ micro g hydroxyl proline was found to be higher in the culture supernatants of the fibroblasts from hypertrophic scars (9.04 +/- 2.74) than that of normal skin (7.55 +/- 2.1). When the human fibroblasts from normal skin and hypertrophic scar were subject to hydroxyl radicals generated by the Fenton reaction, there was no significant increase in pyridinoline cross-link concentration (nm/ micro g hydroxyl proline) in the supernatants compared with the control. When the controls plus various treatments with free radicals were subject to different antioxidants, including superoxide dismutase, catalase, glutathione peroxidase, and desferrioxamine, it was found that catalase alone was most effective in scavenging hydroxyl radicals as determined by the decrease in pyridinoline cross-links.

Selection of an optimal neural network architecture for computer-aided detection of microcalcifications--comparison of automated optimization techniques.

Many computer-aided diagnosis (CAD) systems use neural networks (NNs) for either detection or classification of abnormalities. Currently, most NNs are "optimized" by manual search in a very limited parameter space. In this work, we evaluated the use of automated optimization methods for selecting an optimal convolution neural network (CNN) architecture. Three automated methods, the steepest descent (SD), the simulated annealing (SA), and the genetic algorithm (GA), were compared. We used as an example the CNN that classifies true and false microcalcifications detected on digitized mammograms by a prescreening algorithm. Four parameters of the CNN architecture were considered for optimization, the numbers of node groups and the filter kernel sizes in the first and second hidden layers, resulting in a search space of 432 possible architectures. The area Az under the receiver operating characteristic (ROC) curve was used to design a cost function. The SA experiments were conducted with four different annealing schedules. Three different parent selection methods were compared for the GA experiments. An available data set was split into two groups with approximately equal number of samples. By using the two groups alternately for training and testing, two different cost surfaces were evaluated. For the first cost surface, the SD method was trapped in a local minimum 91% (392/432) of the time. The SA using the Boltzman schedule selected the best architecture after evaluating, on average, 167 architectures. The GA achieved its best performance with linearly scaled roulette-wheel parent selection; however, it evaluated 391 different architectures, on average, to find the best one. The second cost surface contained no local minimum. For this surface, a simple SD algorithm could quickly find the global minimum, but the SA with the very fast reannealing schedule was still the most efficient. The same SA scheme, however, was trapped in a local minimum on the first cost surface. Our CNN study demonstrated that, if optimization is to be performed on a cost surface whose characteristics are not known a priori, it is advisable that a moderately fast algorithm such as a SA using a Boltzman cooling schedule be used to conduct an efficient and thorough search, which may offer a better chance of reaching the global minimum.

Improvement of mammographic mass characterization using spiculation meausures and morphological features.

We are developing new computer vision techniques for characterization of breast masses on mammograms. We had previously developed a characterization method based on texture features. The goal of the present work was to improve our characterization method by making use of morphological features. Toward this goal, we have developed a fully automated, three-stage segmentation method that includes clustering, active contour, and spiculation detection stages. After segmentation, morphological features describing the shape of the mass were extracted. Texture features were also extracted from a band of pixels surrounding the mass. Stepwise feature selection and linear discriminant analysis were employed in the morphological, texture, and combined feature spaces for classifier design. The classification accuracy was evaluated using the area Az under the receiver operating characteristic curve. A data set containing 249 films from 102 patients was used. When the leave-one-case-out method was applied to partition the data set into trainers and testers, the average test Az for the task of classifying the mass on a single mammographic view was 0.83 +/- 0.02, 0.84 +/- 0.02, and 0.87 +/- 0.02 in the morphological, texture, and combined feature spaces, respectively. The improvement obtained by supplementing texture features with morphological features in classification was statistically significant (p = 0.04). For classifying a mass as malignant or benign, we combined the leave-one-case-out discriminant scores from different views of a mass to obtain a summary score. In this task, the test Az value using the combined feature space was 0.91 +/- 0.02. Our results indicate that combining texture features with morphological features extracted from automatically segmented mass boundaries will be an effective approach for computer-aided characterization of mammographic masses.

Computerized image analysis: estimation of breast density on mammograms.

An automated image analysis tool is being developed for the estimation of mammographic breast density. This tool may be useful for risk estimation or for monitoring breast density change in prevention or intervention programs. In this preliminary study, a data set of 4-view mammograms from 65 patients was used to evaluate our approach. Breast density analysis was performed on the digitized mammograms in three stages. First, the breast region was segmented from the surrounding background by an automated breast boundary-tracking algorithm. Second, an adaptive dynamic range compression technique was applied to the breast image to reduce the range of the gray level distribution in the low frequency background and to enhance the differences in the characteristic features of the gray level histogram for breasts of different densities. Third, rule-based classification was used to classify the breast images into four classes according to the characteristic features of their gray level histogram. For each image, a gray level threshold was automatically determined to segment the dense tissue from the breast region. The area of segmented dense tissue as a percentage of the breast area was then estimated. To evaluate the performance of the algorithm, the computer segmentation results were compared to manual segmentation with interactive thresholding by five radiologists. A "true" percent dense area for each mammogram was obtained by averaging the manually segmented areas of the radiologists. We found that the histograms of 6% (8 CC and 8 MLO views) of the breast regions were misclassified by the computer, resulting in poor segmentation of the dense region. For the images with correct classification, the correlation between the computer-estimated percent dense area and the "truth" was 0.94 and 0.91, respectively, for CC and MLO views, with a mean bias of less than 2%. The mean biases of the five radiologists' visual estimates for the same images ranged from 0.1% to 11%. The results demonstrate the feasibility of estimating mammographic breast density using computer vision techniques and its potential to improve the accuracy and reproducibility of breast density estimation in comparison with the subjective visual assessment by radiologists.

Automated registration of breast lesions in temporal pairs of mammograms for interval change analysis--local affine transformation for improved localization.

Analysis of interval change is important for mammographic interpretation. The aim of this study is to evaluate the use of an automated registration technique for computer-aided interval change analysis in mammography. Previously we developed a regional registration technique for identifying masses on temporal pairs of mammograms. In the current study, we improved lesion registration by including a local alignment step. Initially, the lesion position on the prior mammogram was estimated based on the breast geometry. An initial fan-shaped search region was then defined on the prior mammogram. In the second stage, the location of the fan-shaped region on the prior mammogram was refined by warping, based on an affine transformation and simplex optimization in a local region. In the third stage, a search for the best match between the lesion template from the current mammogram and a structure on the prior mammogram was carried out within the search region. This technique was evaluated on 124 temporal pairs of mammograms containing biopsyproven masses. Eighty-seven percent of the estimated lesion locations resulted in an area overlap of at least 50% with the true lesion locations and an average distance of 2.4 +/- 2.1 mm between their centroids. The average distance between the estimated and the true centroid of the lesions on the prior mammogram over all 124 temporal pairs was 4.2 +/- 5.7 mm. The registration accuracy was improved in comparison with our previous study that used a data set of 74 temporal pairs of mammograms. This improvement in accuracy resulted from the improved geometry estimation and the local affine transformation.

Analysis of uncertainties in estimates of components of variance in multivariate ROC analysis.

RATIONALE AND OBJECTIVES: Solutions have previously been presented to the problem of estimating the components of variance in the general linear model used for multivariate receiver operating characteristic (ROC) analysis. The case where the variance components do not change across the modalities under comparison was first treated, followed by the case where they are permitted to change. No analysis of uncertainties in these estimates has been presented previously. MATERIALS AND METHODS: For the case where the variance components do not change across modalities, the "jackknife-after-bootstrap" resampling procedure can be used together with conventional linear propagation of variance to solve for the uncertainties in estimates of the components. For the case where the components are permitted to change across modalities, a slight elaboration of this procedure is presented. RESULTS: The approach was validated by Monte Carlo simulations, where uncertainties in estimates of the variance components calculated by the jackknife-after-bootstrap procedure were found to converge in the mean to the Monte Carlo results over many independent trials. The method is exemplified with data from a study of readers-with and without the aid of a computer-assist modality-given the task of discriminating benign from malignant masses in mammography. CONCLUSION: The present approach is relevant to a broad class of problems where estimates of multiple contributions to the variance observed in ROC assessment of diagnostic modalities are desired, in particular, for the assessment of multiple-reader studies of computer-aided diagnosis in radiology where the variance components may change across reading modalities (eg, unaided vs computer-aided reading).

Digital mammography: observer performance study of the effects of pixel size on the characterization of malignant and benign microcalcifications.

RATIONALE AND OBJECTIVES: The authors performed this study to evaluate the effects of pixel size on the characterization of mammographic microcalcifications by radiologists. MATERIALS AND METHODS: Two-view mammograms of 112 microcalcification clusters were digitized with a laser scanner at a pixel size of 35 microm. Images with pixel sizes of 70, 105, and 140 microm were derived from the 35-microm-pixel size images by averaging neighboring pixels. The malignancy or benignity of the microcalcifications had been determined with findings at biopsy or 2-year follow-up. Region-of-interest images containing the microcalcifications were printed with a laser imager. Seven radiologists participated in a receiver operating characteristic (ROC) study to estimate the likelihood of malignancy. The classification accuracy was quantified with the area under the ROC curve (Az). The statistical significance of the differences in the Az values for different pixel sizes was estimated with the Dorfman-Berbaum-Metz method and the Student paired t test. The variance components were analyzed with a bootstrap method. RESULTS: The higher-resolution images did not result in better classification; the average Az with a pixel size of 35 microm was lower than that with pixel sizes of 70 and 105 microm. The differences in Az between different pixel sizes did not achieve statistical significance. CONCLUSION: Pixel sizes in the range studied do not have a strong effect on radiologists' accuracy in the characterization of microcalcifications. The low specificity of the image features of microcalcifications and the large interobserver and intraobserver variabilities may have prevented small advantages in image resolution from being observed.

Computer-aided characterization of mammographic masses: accuracy of mass segmentation and its effects on characterization.

Mass segmentation is used as the first step in many computer-aided diagnosis (CAD) systems for classification of breast masses as malignant or benign. The goal of this paper was to study the accuracy of an automated mass segmentation method developed in our laboratory, and to investigate the effect of the segmentation stage on the overall classification accuracy. The automated segmentation method was quantitatively compared with manual segmentation by two expert radiologists (R1 and R2) using three similarity or distance measures on a data set of 100 masses. The area overlap measures between R1 and R2, the computer and R1, and the computer and R2 were 0.76 +/- 0.13, 0.74 +/- 0.11, and 0.74 +/- 0.13, respectively. The interobserver difference in these measures between the two radiologists was compared with the corresponding differences between the computer and the radiologists. Using three similarity measures and data from two radiologists, a total of six statistical tests were performed. The difference between the computer and the radiologist segmentation was significantly larger than the interobserver variability in only one test. Two sets of texture, morphological, and spiculation features, one based on the computer segmentation, and the other based on radiologist segmentation, were extracted from a data set of 249 films from 102 patients. A classifier based on stepwise feature selection and linear discriminant analysis was trained and tested using the two feature sets. The leave-one-case-out method was used for data sampling. For case-based classification, the area Az under the receiver operating characteristic (ROC) curve was 0.89 and 0.88 for the feature sets based on the radiologist segmentation and computer segmentation, respectively. The difference between the two ROC curves was not statistically significant.

Analysis of temporal changes of mammographic features: computer-aided classification of malignant and benign breast masses.

A new classification scheme was developed to classify mammographic masses as malignant and benign by using interval change information. The masses on both the current and the prior mammograms were automatically segmented using an active contour method. From each mass, 20 run length statistics (RLS) texture features, 3 speculation features, and 12 morphological features were extracted. Additionally, 20 difference RLS features were obtained by subtracting the prior RLS features from the corresponding current RLS features. The feature space consisted of the current RLS features, the difference RLS features, the current and prior speculation features, and the current and prior mass sizes. Stepwise feature selection and linear discriminant analysis classification were used to select and merge the most useful features. A leave-one-case-out resampling scheme was used to train and test the classifier using 140 temporal image pairs (85 malignant, 55 benign) obtained from 57 biopsy-proven masses (33 malignant, 24 benign) in 56 patients. An average of 10 features were selected from the 56 training subsets: 4 difference RLS features, 4 RLS features, and 1 speculation feature from the current image, and 1 speculation feature from the prior, were most often chosen. The classifier achieved an average training Az of 0.92 and a test Az of 0.88. For comparison, a classifier was trained and tested using features extracted from the 120 current single images. This classifier achieved an average training Az of 0.90 and a test Az of 0.82. The information on the prior image significantly (p = 0.015) improved the accuracy for classification of the masses.

Feature selection and classifier performance in computer-aided diagnosis: the effect of finite sample size.

In computer-aided diagnosis (CAD), a frequently used approach for distinguishing normal and abnormal cases is first to extract potentially useful features for the classification task. Effective features are then selected from this entire pool of available features. Finally, a classifier is designed using the selected features. In this study, we investigated the effect of finite sample size on classification accuracy when classifier design involves stepwise feature selection in linear discriminant analysis, which is the most commonly used feature selection algorithm for linear classifiers. The feature selection and the classifier coefficient estimation steps were considered to be cascading stages in the classifier design process. We compared the performance of the classifier when feature selection was performed on the design samples alone and on the entire set of available samples, which consisted of design and test samples. The area Az under the receiver operating characteristic curve was used as our performance measure. After linear classifier coefficient estimation using the design samples, we studied the hold-out and resubstitution performance estimates. The two classes were assumed to have multidimensional Gaussian distributions, with a large number of features available for feature selection. We investigated the dependence of feature selection performance on the covariance matrices and means for the two classes, and examined the effects of sample size, number of available features, and parameters of stepwise feature selection on classifier bias. Our results indicated that the resubstitution estimate was always optimistically biased, except in cases where the parameters of stepwise feature selection were chosen such that too few features were selected by the stepwise procedure. When feature selection was performed using only the design samples, the hold-out estimate was always pessimistically biased. When feature selection was performed using the entire finite sample space, the hold-out estimates could be pessimistically or optimistically biased, depending on the number of features available for selection, the number of available samples, and their statistical distribution. For our simulation conditions, these estimates were always pessimistically (conservatively) biased if the ratio of the total number of available samples per class to the number of available features was greater than five.

Adverse effects of increased body weight on quantitative measures of mammographic image quality.

OBJECTIVE: The purpose of our study was to show that compressed breast thickness on mammograms in overweight and obese women exceeds the thickness in normal-weight women and that increased thickness results in image degradation. SUBJECTS AND METHODS: Three hundred consecutive routine mammograms were reviewed. Patients were categorized according to body mass index. Compression thickness, compressive force, kilovoltage, and milliampere-seconds were recorded. Geometric unsharpness and contrast degradation were calculated for each body mass index category. RESULTS: Body mass index categories were lean (3%), normal (36%), overweight (36%), and obese (25%). Body mass index was directly correlated with compressed thickness. In the mediolateral oblique view, the mean thickness of the obese category exceeded normal thickness by 18 mm (p < 0.01), corresponding to a 32% increase in geometric unsharpness. Mean obese thickness exceeded lean thickness by 33 mm (p < 0.01), corresponding to a 79% increase in unsharpness. Similar trends were observed for the craniocaudal view. In the mediolateral oblique projection, there was an increase of 1.0 kVp (p < 0.01) for obese compared with normal and 1.7 kVp (p < 0.01) between lean and obese, corresponding, respectively, to a 16% and a 25% decrease in image contrast because of scatter and kilovoltage changes. Milliampere-seconds increased by 47% on the mediolateral oblique images in the obese category compared with normal body mass index. CONCLUSION: An increased body mass index was associated with greater compressed breast thickness, resulting in increased geometric unsharpness, decreased image contrast, and greater potential for motion unsharpness.