Three dimensional texture analysis in MRI: a preliminary evaluation in gliomas

Magnetic resonance imaging (MRI)-based intratumoral and peritumoral radiomics for prognosis prediction in glioma patients

AIM

The purpose of this study was to identify robust radiological features from intratumoral and peritumoral regions, evaluate MRI protocols, and machine learning methods for overall survival stratification of glioma patients, and explore the relationship between radiological features and the tumour microenvironment.

MATERIAL AND METHODS

A retrospective analysis was conducted on 163 glioma patients, divided into a training set (n=113) and a testing set (n=50). For each patient, 2135 features were extracted from clinical MRI. Feature selection was performed using the Minimum Redundancy Maximum Relevance method and the Random Forest (RF) algorithm. Prognostic factors were assessed using the Cox proportional hazards model. Four machine learning models (RF, Logistic Regression, Support Vector Machine, and XGBoost) were trained on clinical and radiological features from tumour and peritumoral regions. Model evaluations on the testing set used receiver operating characteristic curves.

RESULTS

Among the 163 patients, 96 had an overall survival (OS) of less than three years postsurgery, while 67 had an OS of more than three years. Univariate Cox regression in the validation set indicated that age (p=0.003) and tumour grade (p<0.001) were positively associated with the risk of death within three years postsurgery. The final predictive model incorporated 13 radiological and 7 clinical features. The RF model, combining intratumor and peritumor radiomics, achieved the best predictive performance (AUC = 0.91; ACC = 0.86), outperforming single-region models.

CONCLUSION

Combined intratumoral and peritumoral radiomics can improve survival prediction and have potential as a practical imaging biomarker to guide clinical decision-making.

Persistent fatigue in post-acute COVID syndrome is associated with altered T1 MRI texture in subcortical structures: a preliminary investigation

Post-acute COVID syndrome (PACS) is a global health concern and is often associated with debilitating symptoms. Post-COVID fatigue is a particularly frequent and troubling issue, and its underlying mechanisms remain incompletely understood. One potential contributor is micropathological injury of subcortical and brainstem structures, as has been identified in other patient populations. Texture-based analysis (TA) may be used to measure such changes in anatomical MRI data. The present study develops a methodology of voxel-wise TA mapping in subcortical and brainstem regions, which is then applied to T1-weighted MRI data from a cohort of 48 individuals who had PACS (32 with and 16 without ongoing fatigue symptoms) and 15 controls who had cold and flu-like symptoms but tested negative for COVID-19. Both groups were assessed an average of 4-5 months post-infection. There were no significant differences between PACS and control groups, but significant differences were observed within the PACS groups, between those with and without fatigue symptoms. This included reduced texture energy and increased entropy, along with reduced texture correlation, cluster shade and profile in the putamen, pallidum, thalamus and brainstem. These findings provide new insights into the neurophysiological mechanisms that underlie PACS, with altered tissue texture as a potential biomarker of this debilitating condition.

Texture Estimation for Abnormal Tissue Segmentation in Brain MRI

This chapter discusses multifractal texture estimation and characterization of brain lesions (necrosis, edema, enhanced tumor, nonenhanced tumor, etc.) in magnetic resonance (MR) images. This work formulates the complex texture of tumor in MR images using a stochastic model known as multifractional Brownian motion (mBm). Mathematical derivations of the mBm model and corresponding algorithm to extract the spatially varying multifractal texture feature are discussed. Extracted multifractal texture feature is fused with other effective features to enhance the tissue characteristics. Segmentation of the tissues is performed using a feature-based classification method. The efficacy of the mBm texture feature in segmenting different abnormal tissues is demonstrated using a large-scale publicly available clinical dataset. Experimental results and performance of the methods confirm the efficacy of the proposed technique in an automatic segmentation of abnormal tissues in multimodal (T1, T2, Flair, and T1contrast) brain MRIs.

Preoperative prediction of lymphovascular space invasion in endometrioid adenocarcinoma: an MRI-based radiomics nomogram with consideration of the peritumoral region

BACKGROUND

Lymphovascular space invasion (LVSI) of endometrial cancer (EC) is a postoperative histological index, which is associated with lymph node metastases. A preoperative acknowledgement of LVSI status might aid in treatment decision-making.

PURPOSE

To explore the utility of multiparameter magnetic resonance imaging (MRI) and radiomic features obtained from intratumoral and peritumoral regions for predicting LVSI in endometrioid adenocarcinoma (EEA).

MATERIAL AND METHODS

A total of 334 EEA tumors were retrospectively analyzed. Axial T2-weighted (T2W) imaging and apparent diffusion coefficient (ADC) mapping were conducted. Intratumoral and peritumoral regions were manually annotated as the volumes of interest (VOIs). A support vector machine was applied to train the prediction models. Multivariate logistic regression analysis was used to develop a nomogram based on clinical and tumor morphological parameters and the radiomics score (RadScore). The predictive performance of the nomogram was assessed by the area under the receiver operator characteristic curve (AUC) in the training and validation cohorts.

RESULTS

Among the features obtained from different imaging modalities (T2W imaging and ADC mapping) and VOIs, the RadScore had the best performance in predicting LVSI classification (AUCtrain = 0.919, and AUCvalidation = 0.902). The nomogram based on age, CA125, maximum anteroposterior tumor diameter on sagittal T2W images, tumor area ratio, and RadScore was established to predict LVSI had AUC values in the training and validation cohorts of 0.962 (sensitivity 94.0%, specificity 86.0%) and 0.965 (sensitivity 90.0%, specificity 85.3%), respectively.

CONCLUSION

The intratumoral and peritumoral imaging features were complementary, and the MRI-based radiomics nomogram might serve as a non-invasive biomarker to preoperatively predict LVSI in patients with EEA.

Multi-modality radiomics nomogram based on DCE-MRI and ultrasound images for benign and malignant breast lesion classification

OBJECTIVE

To develop a multi-modality radiomics nomogram based on DCE-MRI, B-mode ultrasound (BMUS) and strain elastography (SE) images for classifying benign and malignant breast lesions.

MATERIAL AND METHODS

In this retrospective study, 345 breast lesions from 305 patients who underwent DCE-MRI, BMUS and SE examinations were randomly divided into training (n = 241) and testing (n = 104) datasets. Radiomics features were extracted from manually contoured images. The inter-class correlation coefficient (ICC), Mann-Whitney U test and the least absolute shrinkage and selection operator (LASSO) regression were applied for feature selection and radiomics signature building. Multivariable logistic regression was used to develop a radiomics nomogram incorporating radiomics signature and clinical factors. The performance of the radiomics nomogram was evaluated by its discrimination, calibration, and clinical usefulness and was compared with BI-RADS classification evaluated by a senior breast radiologist.

RESULTS

The All-Combination radiomics signature derived from the combination of DCE-MRI, BMUS and SE images showed better diagnostic performance than signatures derived from single modality alone, with area under the curves (AUCs) of 0.953 and 0.941 in training and testing datasets, respectively. The multi-modality radiomics nomogram incorporating the All-Combination radiomics signature and age showed excellent discrimination with the highest AUCs of 0.964 and 0.951 in two datasets, respectively, which outperformed all single modality radiomics signatures and BI-RADS classification. Furthermore, the specificity of radiomics nomogram was significantly higher than BI-RADS classification (both p < 0.04) with the same sensitivity in both datasets.

CONCLUSION

The proposed multi-modality radiomics nomogram based on DCE-MRI and ultrasound images has the potential to serve as a non-invasive tool for classifying benign and malignant breast lesions and reduce unnecessary biopsy.

Radiomics in Head and Neck Cancer Outcome Predictions

Head and neck cancer has great regional anatomical complexity, as it can develop in different structures, exhibiting diverse tumour manifestations and high intratumoural heterogeneity, which is highly related to resistance to treatment, progression, the appearance of metastases, and tumour recurrences. Radiomics has the potential to address these obstacles by extracting quantitative, measurable, and extractable features from the region of interest in medical images. Medical imaging is a common source of information in clinical practice, presenting a potential alternative to biopsy, as it allows the extraction of a large number of features that, although not visible to the naked eye, may be relevant for tumour characterisation. Taking advantage of machine learning techniques, the set of features extracted when associated with biological parameters can be used for diagnosis, prognosis, and predictive accuracy valuable for clinical decision-making. Therefore, the main goal of this contribution was to determine to what extent the features extracted from Computed Tomography (CT) are related to cancer prognosis, namely Locoregional Recurrences (LRs), the development of Distant Metastases (DMs), and Overall Survival (OS). Through the set of tumour characteristics, predictive models were developed using machine learning techniques. The tumour was described by radiomic features, extracted from images, and by the clinical data of the patient. The performance of the models demonstrated that the most successful algorithm was XGBoost, and the inclusion of the patients' clinical data was an asset for cancer prognosis. Under these conditions, models were created that can reliably predict the LR, DM, and OS status, with the area under the ROC curve (AUC) values equal to 0.74, 0.84, and 0.91, respectively. In summary, the promising results obtained show the potential of radiomics, once the considered cancer prognosis can, in fact, be expressed through CT scans.

Comparison of three image-analysis-based visual texture calculation methods: energy, entropy, and texture change index

BACKGROUND

Three image analysis methods to measure visual texture were applied to an image with much texture (scaled carp), and one with little texture (mirror carp). For each method, the effect of image rotation at 0°, 10°, 20°, 30°, 45°, 60°, 75° and 90° was evaluated.

RESULTS

Rotation did not change energy (E) and entropy (H) calculations using image histograms. Using co-occurrence matrices with different step size d (1-19 in increases of 2) and step angle θ (0°, 45°, 90° and 135°) showed that, as d increased, E decreased and H increased, and the number of legitimate pixel pairs decreased linearly. Averaging E and H at different θ values rendered the results rotation invariant. Theoretically, the 'texture primitives' method is not rotation independent. However, the variation in texture change index (TCI) with image rotation was negligible. Also, the increase in TCI between the less textured image and the more textured image was 5.3-11. In comparison, the E values from histograms for the images above were 0.0069-0.0081. For co-occurrence matrix-based calculations, at d = 1 and for all θ, E range was from 220 to 389 for scaled carp and from 232 to 412 for mirror carp.

CONCLUSION

The more than doubling of TCI for these images implies that it is a more precise method than energy and entropy to discern between visual texture levels. © 2022 Society of Chemical Industry.

An update on radiomics techniques in primary liver cancers

BACKGROUND

Radiomics is a progressing field of research that deals with the extraction of quantitative metrics from medical images. Radiomic features detention indirectly tissue features such as heterogeneity and shape and can, alone or in combination with demographic, histological, genomic, or proteomic data, be used for decision support system in clinical setting.

METHODS

This article is a narrative review on Radiomics in Primary Liver Cancers. Particularly, limitations and future perspectives are discussed.

RESULTS

In oncology, assessment of tissue heterogeneity is of particular interest: genomic analysis have demonstrated that the degree of tumour heterogeneity is a prognostic determinant of survival and an obstacle to cancer control. Therefore, that Radiomics could support cancer detection, diagnosis, evaluation of prognosis and response to treatment, so as could supervise disease status in hepatocellular carcinoma (HCC) and Intrahepatic Cholangiocarcinoma (ICC) patients. Radiomic analysis is a convenient radiological image analysis technique used to support clinical decisions as it is able to provide prognostic and / or predictive biomarkers that allow a fast, objective and repeatable tool for disease monitoring.

CONCLUSIONS

Although several studies have shown that this analysis is very promising, there is little standardization and generalization of the results, which limits the translation of this method into the clinical context. The limitations are mainly related to the evaluation of data quality, repeatability, reproducibility, overfitting of the model.

TRIAL REGISTRATION

Not applicable.

Differentiation of Cartilage Repair Techniques Using Texture Analysis from T2 Maps

OBJECTIVE

The aim of this study was to investigate texture features from T2 maps as a marker for distinguishing the maturation of repair tissue after 2 different cartilage repair procedures.

DESIGN

Seventy-nine patients, after either microfracture (MFX) or matrix-associated chondrocyte transplantation (MACT), were examined on a 3-T magnetic resonance (MR) scanner with morphological and quantitative (T2 mapping) MR sequences 2 years after surgery. Twenty-one texture features from a gray-level co-occurrence matrix (GLCM) were extracted. The texture feature difference between 2 repair types was assessed individually for the femoral condyle and trochlea/anterior condyle using linear regression models. The stability and reproducibility of texture features for focal cartilage were calculated using intra-observer variability and area under curve from receiver operating characteristics.

RESULTS

There was no statistical significance found between MFX and MACT for T2 values (P = 0.96). There was, however, found a statistical significance between MFX and MACT in femoral condyle in GLCM features autocorrelation (P < 0.001), sum of squares (P = 0.023), sum average (P = 0.005), sum variance (P = 0.0048), and sum entropy (P = 0.05); and in anterior condyle/trochlea homogeneity (P = 0.02) and dissimilarity (P < 0.001).

CONCLUSION

Texture analysis using GLCM provides a useful extension to T2 mapping for the characterization of cartilage repair tissue by increasing its sensitivity to tissue structure. Some texture features were able to distinguish between repair tissue after different cartilage repair procedures, as repair tissue texture (and hence, probably collagen organization) 24 months after MACT more closely resembled healthy cartilage than did MFX repair tissue.

Classification of neoplastic and inflammatory brain disease using MRI texture analysis in 119 dogs

Magnetic resonance imaging is the primary method used to diagnose canine glial cell neoplasia and noninfectious inflammatory meningoencephalitis. Subjective differentiation of these diseases can be difficult due to overlapping imaging characteristics. This study utilizes texture analysis (TA) of intra-axial lesions both as a means to quantitatively differentiate these broad categories of disease and to help identify glial tumor grade/cell type and specific meningoencephalitis subtype in a group of 119 dogs with histologically confirmed diagnoses. Fifty-nine dogs with gliomas and 60 dogs with noninfectious inflammatory meningoencephalitis were retrospectively recruited and randomly split into training (n = 80) and test (n = 39) cohorts. Forty-five of 120 texture metrics differed significantly between cohorts after correcting for multiple testing (false discovery rate < 0.05). After training the random forest algorithm, the classification accuracy for the test set was 85% (sensitivity 89%, specificity 81%). TA was only partially able to differentiate the inflammatory subtypes (granulomatous meningoencephalitis [GME], necrotizing meningoencephalitis [NME], and necrotizing leukoencephalitis [NLE]) (out-of-bag error rate of 35.0%) and was unable to identify metrics that could correctly classify glioma grade or cell type (out-of-bag error rate of 59.6% and 47.5%, respectively). Multiple demographic differences, such as patient age, sex, weight, and breed were identified between disease cohorts and subtypes which may be useful in prioritizing differential diagnoses. TA of MR images with a random forest algorithm provided classification accuracy of inflammatory and neoplastic brain disease approaching the accuracy of previously reported subjective radiologist evaluation.

Radiomics in hepatic metastasis by colorectal cancer

Background

Radiomics is an emerging field and has a keen interest, especially in the oncology field. The process of a radiomics study consists of lesion segmentation, feature extraction, consistency analysis of features, feature selection, and model building. Manual segmentation is one of the most critical parts of radiomics. It can be time-consuming and suffers from variability in tumor delineation, which leads to the reproducibility problem of calculating parameters and assessing spatial tumor heterogeneity, particularly in large or multiple tumors. Radiomic features provides data on tumor phenotype as well as cancer microenvironment. Radiomics derived parameters, when associated with other pertinent data and correlated with outcomes data, can produce accurate robust evidence based clinical decision support systems. The principal challenge is the optimal collection and integration of diverse multimodal data sources in a quantitative manner that delivers unambiguous clinical predictions that accurately and robustly enable outcome prediction as a function of the impending decisions.

Methods

The search covered the years from January 2010 to January 2021. The inclusion criterion was: clinical study evaluating radiomics of liver colorectal metastases. Exclusion criteria were studies with no sufficient reported data, case report, review or editorial letter.

Results

We recognized 38 studies that assessed radiomics in mCRC from January 2010 to January 2021. Twenty were on different tpics, 5 corresponded to most criteria; 3 are review, or letter to editors; so 10 articles were included.

Conclusions

In colorectal liver metastases radiomics should be a valid tool for the characterization of lesions, in the stratification of patients based on the risk of relapse after surgical treatment and in the prediction of response to chemotherapy treatment.

Quantitative MRI texture analysis in chronic active multiple sclerosis lesions

OBJECTIVE

Recently, there has been an increasing interest in "chronic enlarging" or "chronic active" multiple sclerosis (MS) lesions that are associated with clinical disability. However, investigation of dynamic lesion volume changes requires longitudinal MRI data from two or more time points. The aim of this study was to investigate the application of texture analysis (TA) on baseline T1-weighted 3D magnetization-prepared rapid acquisition gradient-echo (MPRAGE) images to differentiate chronic active from chronic stable MS lesions.

MATERIAL AND METHODS

To identify chronic active lesions as compared to non-enhancing stable lesions, two MPRAGE datasets acquired on a 3 T MRI at baseline and after 12 months follow-up were applied to the Voxel-Guided Morphometry (VGM) algorithm. TA was performed on the baseline MPRAGE images, 36 texture features were extracted for each lesion.

RESULTS

Overall, 374 chronic MS lesions (155 chronic active and 219 chronic stable lesions) from 60 MS patients were included in the final analysis. Multiple texture features including "DISCRETIZED_HISTO_Energy", "GLCM_Energy", "GLCM_Contrast" and "GLCM_Dissimilarity" were significantly higher in chronic active as compared to chronic stable lesions. Partial least squares regression yielded an area under the curve of 0.7 to differentiate both lesion types.

CONCLUSION

Our results suggest that multiple texture features extracted from MPRAGE images indicate higher intralesional heterogeneity, however they demonstrate only a fair accuracy to differentiate chronic active from chronic stable MS lesions.

Texture analysis for chemotherapy response evaluation in osteosarcoma using MR imaging

The efficacy of MRI-based statistical texture analysis (TA) in predicting chemotherapy response among patients with osteosarcoma was assessed. Forty patients (male: female = 31:9; age = 17.2 ± 5.7 years) with biopsy-proven osteosarcoma were analyzed in this prospective study. Patients were scheduled for three cycles of neoadjuvant chemotherapy (NACT) and diffusion-weighted MRI acquisition at three time points: at baseline (t0), after the first NACT (t1) and after the third NACT (t2) using a 1.5 T scanner. Eight patients (nonsurvivors) died during NACT while 34 patients (survivors) completed the NACT regimen followed by surgery. Histopathological evaluation was performed in the resected tumor to assess NACT response (responder [≤50% viable tumor] and nonresponder [>50% viable tumor]) and revealed nonresponder: responder = 20:12. Apparent diffusion coefficient (ADC) and intravoxel incoherent motion (IVIM) parameters, diffusion coefficient (D), perfusion coefficient (D*) and perfusion fraction (f) were evaluated. A total of 25 textural features were evaluated on ADC, D, D* and f parametric maps and structural T1-weighted (T1W) and T2-weighted (T2W) images in the entire tumor volume using 3D TA methods gray-level cooccurrence matrix (GLCM), neighborhood gray-tone-difference matrix (NGTDM) and run-length matrix (RLM). Receiver-operating-characteristic curve analysis was performed on the selected textural feature set to assess the role of TA features (a) as marker(s) of tumor aggressiveness leading to mortality at baseline and (b) in predicting the NACT response among survivors in the course of treatment. Findings showed that the NGTDM features coarseness, busyness and strength quantifying tumor heterogeneity in D, D* and f maps and T1W and T2W images were useful markers of tumor aggressiveness in identifying the nonsurvivor group (area-under-the-curve [AUC] = 0.82-0.88) at baseline. The GLCM features contrast and correlation, NGTDM features contrast and complexity and RLM feature short-run-low-gray-level-emphasis quantifying homogeneity/terogeneity in tumor were effective markers for predicting chemotherapeutic response using D (AUC = 0.80), D* (AUC = 0.80) and T2W (AUC = 0.70) at t0, and D* (AUC = 0.80) and f (AUC = 0.70) at t1. 3D statistical TA features might be useful as imaging-based markers for characterizing tumor aggressiveness and predicting chemotherapeutic response in patients with osteosarcoma.

Advanced intra-tumoural structural characterisation of hepatocellular carcinoma utilising FDG-PET/CT: a comparative study of radiomics and metabolic features in 3D and 2D

PURPOSE

The aim of our work is to evaluate the correlation of two-dimensional (2D) and three-dimensional (3D) radiomics and metabolic features of hepatocellular carcinoma (HCC) with tumour diameter, staging, and metabolic tumour volume (MTV).

MATERIAL AND METHODS

Thirty-three patients with HCC were studied using 18F-fluorodeoxyglucose positron-emission tomography with computed tomography (18F [FDG] PET/CT). The tumours were segmented from the PET images after CT correction. Metabolic parameters and 35 radiomics features were compared using 2D and 3D modes. The metabolic parameters and tumour morphology were compared using 2 different types of software. Tumour heterogeneity was studied in both metabolic parameters and radiomics features. Finally, the correlation between the metabolic and radiomics features in 3D mode, as well as tumour morphology and staging according to the American Joint Committee on Cancer (AJCC) staging were studied.

RESULTS

Most of the metabolic parameters and radiomics features are statically stable through the 2D and 3D modes. Most of the 3D mode features show a correlation with metabolic parameters; the total lesion glycolysis (TLG) shows the highest correlation, with a Spearman correlation coefficient (rs) of 0.9776. Also, the grey level run length matrix/run length non-uniformity (GLRLM_RLNU) from radiomics features exhibits a correlation with a Spearman correlation coefficient of 0.9733. Maximum tumour diameter is correlated with TLG and GLRLM_RLNU, with rs equal to 0.7461 and 0.7143, respectively. Regarding AJCC staging, some features show a medium but prognostic correlation. In the case of 2D-mode features, all metabolic and radiomics features show no significant correlation with MTV, AJCC staging, and tumour maximum diameter.

CONCLUSIONS

Most of the normal metabolic parameters and radiomics features are statistically stable through the 3D and 2D modes. 3D radiomics features are significantly correlated with tumour volume, maximum diameter, and staging. Conversely, 2D features have negligible correlation with the same parameters. Therefore, 3D mode features are preferable and can accurately evaluate tumour heterogeneity.

Toward automated classification of monolayer versus few-layer nanomaterials using texture analysis and neural networks

The need for a fast and robust method to characterize nanostructure thickness is growing due to the tremendous number of experiments and their associated applications. By automatically analyzing the microscopic image texture of MoS2 and WS2, it was possible to distinguish monolayer from few-layer nanostructures with high accuracy for both materials. Three methods of texture analysis (TA) were used: grey level histogram (GLH), grey levels co-occurrence matrix (GLCOM), and run-length matrix (RLM), which correspond to first, second, and higher-order statistical methods, respectively. The best discriminating features were automatically selected using the Fisher coefficient, for each method, and used as a base for classification. Two classifiers were used: artificial neural networks (ANN), and linear discriminant analysis (LDA). RLM with ANN was found to give high classification accuracy, which was 89% and 95% for MoS2 and WS2, respectively. The result of this work suggests that RLM, as a higher-order TA method, associated with an ANN classifier has a better ability to quantify and characterize the microscopic structure of nanolayers, and, therefore, categorize thickness to the proper class.

Texture analysis in susceptibility-weighted imaging may be useful to differentiate acute from chronic multiple sclerosis lesions

OBJECTIVES

To evaluate the diagnostic performance of texture analysis (TA) applied on non-contrast-enhanced susceptibility-weighted imaging (SWI) to differentiate acute (enhancing) from chronic (non-enhancing) multiple sclerosis (MS) lesions.

METHODS

We analyzed 175 lesions from 58 patients with relapsing-remitting MS imaged on a 3.0 T MRI scanner and applied TA on T2-w and SWI images to extract texture features. We evaluated the presence or absence of lesion enhancement on T1-w post-contrast images and performed a computational statistical analysis to assess if there was any significant correlation between the texture features and the presence of lesion activity. ROC curves and leave-one-out cross-validation were used to evaluate the performance of individual features and multiparametric models in the identification of active lesions.

RESULTS

Multiple TA features obtained from SWI images showed a significantly different distribution in acute and chronic lesions (AUC, 0.617-0.720). Multiparametric predictive models based on logistic ridge regression and partial least squares regression yielded an AUC of 0.778 and 0.808, respectively. Results from T2-w images did not show any significant predictive ability of neither individual features nor multiparametric models.

CONCLUSIONS

Texture analysis on SWI sequences may be useful to differentiate acute from chronic MS lesions. The good diagnostic performance could help to reduce the need of intravenous contrast agent administration in follow-up MRI studies.

KEY POINTS

• Texture analysis applied on SWI sequences may be useful to differentiate acute from chronic multiple sclerosis lesions • The good diagnostic performance could help to minimize the need of intravenous contrast agent administration in follow-up MRI studies.

Equilibrium CT Texture Analysis for the Evaluation of Hepatic Fibrosis: Preliminary Evaluation against Histopathology and Extracellular Volume Fraction

Background: Evaluate equilibrium contrast-enhanced CT (EQ-CT) texture analysis (EQ-CTTA) against histologically-quantified fibrosis, serum-based enhanced liver fibrosis panel (ELF) and imaging-based extracellular volume fraction (ECV) in chronic hepatitis. Methods: This study was a re-analysis of image data from a previous prospective study. Pre- and equilibrium-phase post-IV contrast CT datasets were collected from patients with chronic hepatitis with contemporaneous liver biopsy and serum ELF measurement between April 2011 and July 2013. Biopsy samples were analysed to derive collagen proportionate area (CPA). EQ-CTTA was performed with a filtration histogram technique using texture analysis software, with texture quantification using statistical and histogram-based metrics (mean, skewness, standard deviation, entropy, etc.). Association between pre-contrast and EQ-CTTA against CPA, ECV and ELF was evaluated using Spearman’s rank correlation coefficient (r_s). Results: Complete datasets collected in 29 patients (16 male; 13 female), mean age (range): 49 (22–66 years). Liver ECV, CPA and ELF had a median (interquartile range) of 0.26 (0.24–0.29); 5.0 (3.0–13.7) and 9.71 (8.39–10.92). Difference in segment VII hepatic CTTA (medium texture scale) between EQ-CT and pre-contrast images was significantly and positively associated with ELF score (mean: r_s = 0.69, p < 0.001; skewness: r_s = 0.57, p = 0.007). Significant negative associations were observed between pre-contrast and EQ-CT whole hepatic CTTA (coarse texture scale) with CPA (pre-contrast, SD: r_s = −0.66, p < 0.001) and ECV (EQ-CT, entropy: r_s = −0.58, p = 0.006). Conclusions: Hepatic EQ-CTTA demonstrates significant association with validated markers of liver fibrosis, suggesting a role in non-invasive quantification of severity in diffuse fibrosis.

Classification of radiographic lung pattern based on texture analysis and machine learning

This study evaluated the feasibility of using texture analysis and machine learning to distinguish radiographic lung patterns. A total of 1200 regions of interest (ROIs) including four specific lung patterns (normal, alveolar, bronchial, and unstructured interstitial) were obtained from 512 thoracic radiographs of 252 dogs and 65 cats. Forty-four texture parameters based on eight methods of texture analysis (first-order statistics, spatial gray-level-dependence matrices, gray-level-difference statistics, gray-level run length image statistics, neighborhood gray-tone difference matrices, fractal dimension texture analysis, Fourier power spectrum, and Law's texture energy measures) were used to extract textural features from the ROIs. The texture parameters of each lung pattern were compared and used for training and testing of artificial neural networks. Classification performance was evaluated by calculating accuracy and the area under the receiver operating characteristic curve (AUC). Forty texture parameters showed significant differences between the lung patterns. The accuracy of lung pattern classification was 99.1% in the training dataset and 91.9% in the testing dataset. The AUCs were above 0.98 in the training set and above 0.92 in the testing dataset. Texture analysis and machine learning algorithms may potentially facilitate the evaluation of medical images.

Texture Analysis to Detect Cerebral Degeneration in Amyotrophic Lateral Sclerosis

BACKGROUND

Evidence of cerebral degeneration is not apparent on routine brain MRI in amyotrophic lateral sclerosis (ALS). Texture analysis can detect change in images based on the statistical properties of voxel intensities. Our objective was to test the utility of texture analysis in detecting cerebral degeneration in ALS. A secondary objective was to determine whether the performance of texture analysis is dependent on image resolution.

METHODS

High-resolution (0.5×0.5 mm2 in-plane) coronal T2-weighted MRI of the brain were acquired from 12 patients with ALS and 19 healthy controls on a 4.7 Tesla MRI system. Image data sets at lower resolutions were created by down-sampling to 1×1, 2×2, 3×3, and 4×4 mm2. Texture features were extracted from a slice encompassing the corticospinal tract at the different resolutions and tested for their discriminatory power and correlations with clinical measures. Subjects were also classified by visual assessment by expert reviewers.

RESULTS

Texture features were different between ALS patients and healthy controls at 1×1, 2×2, and 3×3 mm2 resolutions. Texture features correlated with measures of upper motor neuron function and disability. Optimal classification performance was achieved when best-performing texture features were combined with visual assessment at 2×2 mm2 resolution (0.851 area under the curve, 83% sensitivity, 79% specificity).

CONCLUSIONS

Texture analysis can detect subtle abnormalities in MRI of ALS patients. The clinical yield of the method is dependent on image resolution. Texture analysis holds promise as a potential source of neuroimaging biomarkers in ALS.

A radiomics evaluation of 2D and 3D MRI texture features to classify brain metastases from lung cancer and melanoma

Brain metastases are occasionally detected before diagnosing their primary site of origin. In these cases, simple visual examination of medical images of the metastases is not enough to identify the primary cancer, so an extensive evaluation is needed. To avoid this procedure, a radiomics approach on magnetic resonance (MR) images of the metastatic lesions is proposed to classify two of the most frequent origins (lung cancer and melanoma). In this study, 50 T1-weighted MR images of brain metastases from 30 patients were analyzed: 27 of lung cancer and 23 of melanoma origin. A total of 43 statistical texture features were extracted from the segmented lesions in 2D and 3D. Five predictive models were evaluated using a nested cross-validation scheme. The best classification results were achieved using 3D texture features for all the models, obtaining an average AUC > 0.9 in all cases and an AUC = 0.947 ± 0.067 when using the best model (naïve Bayes).

Molecular classification of patients with grade II/III glioma using quantitative MRI characteristics

BACKGROUND

Molecular markers of WHO grade II/III glioma are known to have important prognostic and predictive implications and may be associated with unique imaging phenotypes. The purpose of this study is to determine whether three clinically relevant molecular markers identified in gliomas-IDH, 1p/19q, and MGMT status-show distinct quantitative MRI characteristics on FLAIR imaging.

METHODS

Sixty-one patients with grade II/III gliomas who had molecular data and MRI available prior to radiation were included. Quantitative MRI features were extracted that measured tissue heterogeneity (homogeneity and pixel correlation) and FLAIR border distinctiveness (edge contrast; EC). T-tests were conducted to determine whether patients with different genotypes differ across the features. Logistic regression with LASSO regularization was used to determine the optimal combination of MRI and clinical features for predicting molecular subtypes.

RESULTS

Patients with IDH wildtype tumors showed greater signal heterogeneity (p = 0.001) and lower EC (p = 0.008) within the FLAIR region compared to IDH mutant tumors. Among patients with IDH mutant tumors, 1p/19q co-deleted tumors had greater signal heterogeneity (p = 0.002) and lower EC (p = 0.005) compared to 1p/19q intact tumors. MGMT methylated tumors showed lower EC (p = 0.03) compared to the unmethylated group. The combination of FLAIR border distinctness, heterogeneity, and pixel correlation optimally classified tumors by IDH status.

CONCLUSION

Quantitative imaging characteristics of FLAIR heterogeneity and border pattern in grade II/III gliomas may provide unique information for determining molecular status at time of initial diagnostic imaging, which may then guide subsequent surgical and medical management.

Classifying brain metastases by their primary site of origin using a radiomics approach based on texture analysis: a feasibility study

OBJECTIVE

To examine the capability of MRI texture analysis to differentiate the primary site of origin of brain metastases following a radiomics approach.

METHODS

Sixty-seven untreated brain metastases (BM) were found in 3D T1-weighted MRI of 38 patients with cancer: 27 from lung cancer, 23 from melanoma and 17 from breast cancer. These lesions were segmented in 2D and 3D to compare the discriminative power of 2D and 3D texture features. The images were quantized using different number of gray-levels to test the influence of quantization. Forty-three rotation-invariant texture features were examined. Feature selection and random forest classification were implemented within a nested cross-validation structure. Classification was evaluated with the area under receiver operating characteristic curve (AUC) considering two strategies: multiclass and one-versus-one.

RESULTS

In the multiclass approach, 3D texture features were more discriminative than 2D features. The best results were achieved for images quantized with 32 gray-levels (AUC = 0.873 ± 0.064) using the top four features provided by the feature selection method based on the p-value. In the one-versus-one approach, high accuracy was obtained when differentiating lung cancer BM from breast cancer BM (four features, AUC = 0.963 ± 0.054) and melanoma BM (eight features, AUC = 0.936 ± 0.070) using the optimal dataset (3D features, 32 gray-levels). Classification of breast cancer and melanoma BM was unsatisfactory (AUC = 0.607 ± 0.180).

CONCLUSION

Volumetric MRI texture features can be useful to differentiate brain metastases from different primary cancers after quantizing the images with the proper number of gray-levels.

KEY POINTS

• Texture analysis is a promising source of biomarkers for classifying brain neoplasms. • MRI texture features of brain metastases could help identifying the primary cancer. • Volumetric texture features are more discriminative than traditional 2D texture features.

Cloud-Based NoSQL Open Database of Pulmonary Nodules for Computer-Aided Lung Cancer Diagnosis and Reproducible Research

Lung cancer is the leading cause of cancer-related deaths in the world, and its main manifestation is pulmonary nodules. Detection and classification of pulmonary nodules are challenging tasks that must be done by qualified specialists, but image interpretation errors make those tasks difficult. In order to aid radiologists on those hard tasks, it is important to integrate the computer-based tools with the lesion detection, pathology diagnosis, and image interpretation processes. However, computer-aided diagnosis research faces the problem of not having enough shared medical reference data for the development, testing, and evaluation of computational methods for diagnosis. In order to minimize this problem, this paper presents a public nonrelational document-oriented cloud-based database of pulmonary nodules characterized by 3D texture attributes, identified by experienced radiologists and classified in nine different subjective characteristics by the same specialists. Our goal with the development of this database is to improve computer-aided lung cancer diagnosis and pulmonary nodule detection and classification research through the deployment of this database in a cloud Database as a Service framework. Pulmonary nodule data was provided by the Lung Image Database Consortium and Image Database Resource Initiative (LIDC-IDRI), image descriptors were acquired by a volumetric texture analysis, and database schema was developed using a document-oriented Not only Structured Query Language (NoSQL) approach. The proposed database is now with 379 exams, 838 nodules, and 8237 images, 4029 of them are CT scans and 4208 manually segmented nodules, and it is allocated in a MongoDB instance on a cloud infrastructure.

Imaging prediction of nonalcoholic steatohepatitis using computed tomography texture analysis

OBJECTIVES

To determine if texture analysis of non-contrast-enhanced CT (NECT) images is able to predict nonalcoholic steatohepatitis (NASH).

METHODS

NECT images from 88 patients who underwent a liver biopsy for the diagnosis of suspected NASH were assessed and texture feature parameters were obtained without and with filtration. The patient population was divided into a predictive learning dataset and a validation dataset, and further divided into groups according to the prediction of liver fibrosis as assessed by hyaluronic acid levels. The reference standard was the histological result of a liver biopsy. A predictive model for NASH was developed using parameters derived from the learning dataset that demonstrated areas under the receiver operating characteristic curve (AUC) of >0.65. The resulting model was then applied to the validation dataset.

RESULTS

In patients without suspected fibrosis, the texture parameter mean without filter and skewness with a 2-mm filter were selected for the NASH prediction model. The AUC of the predictive model for the validation dataset was 0.94 and the accuracy was 94%. In patients with suspicion of fibrosis, the mean without filtration and kurtosis with a 4-mm filter were selected for the NASH prediction model. The AUC for the validation dataset was 0.60 and the accuracy was 42%.

CONCLUSIONS

In patients without suspicion of fibrosis, NECT texture analysis effectively predicted NASH.

KEY POINTS

• In patients without suspicion of fibrosis, NECT texture analysis effectively predicted NASH. • The mean without filtration and skewness with a 2-mm filter were modest predictors of NASH in patients without suspicion of liver fibrosis. • Hepatic fibrosis masks the characteristic texture features of NASH.

Magnetic Resonance Image Texture Analysis of the Periaqueductal Gray Matter in Episodic Migraine Patients without T2-Visible Lesions

Objective

The periaqueductal gray matter (PAG), a small midbrain structure, presents dysfunction in migraine. However, the precise neurological mechanism is still not well understood. Herein, the aim of this study was to investigate the texture characteristics of altered PAG in episodic migraine (EM) patients based on high resolution brain structural magnetic resonance (MR) images.

Materials and Methods

The brain structural MR images were obtained from 18 normal controls (NC), 18 EM patients and 16 chronic migraine (CM) patients using a 3T MR system. A PAG template was created using the International Consortium Brain Mapping 152 gray matter model, and the individual PAG segment was developed by applying the deformation field from the structural image segment to the PAG template. A grey level co-occurrence matrix was used to calculate the texture parameters including the angular second moment (ASM), contrast, correlation, inverse difference moment (IDM) and entropy.

Results

There was a significant difference for ASM, IDM and entropy in the EM group (998.629 ± 0.162 × 10⁻³, 999.311 ± 0.073 × 10⁻³, 916.354 ± 0.947 × 10⁻⁵) compared to that found in the NC group (998.760 ± 0.110 × 10⁻³, 999.358 ± 0.037 × 10⁻³ and 841.198 ± 0.575 × 10⁻⁵) (p < 0.05). The entropy was significantly lower among the patients with CM (864.116 ± 0.571 × 10⁻⁵) than that found among patients with EM (p < 0.05). The area under the receiver operating characteristic curve was 0.776 and 0.750 for ASM and entropy in the distinction of the EM from NC groups, respectively. ASM was negatively related to disease duration (DD) and the Migraine Disability Assessment Scale (MIDAS) scores in the EM group, and entropy was positively related to DD and MIDAS in the EM group (p < 0.05).

Conclusion

The present study identified altered MR image texture characteristics of the PAG in EM. The identified texture characteristics could be considered as imaging biomarkers for EM.

Lack of robustness of textural measures obtained from 3D brain tumor MRIs impose a need for standardization

Purpose

Textural measures have been widely explored as imaging biomarkers in cancer. However, their robustness under dynamic range and spatial resolution changes in brain 3D magnetic resonance images (MRI) has not been assessed. The aim of this work was to study potential variations of textural measures due to changes in MRI protocols.

Materials and methods

Twenty patients harboring glioblastoma with pretreatment 3D T1-weighted MRIs were included in the study. Four different spatial resolution combinations and three dynamic ranges were studied for each patient. Sixteen three-dimensional textural heterogeneity measures were computed for each patient and configuration including co-occurrence matrices (CM) features and run-length matrices (RLM) features. The coefficient of variation was used to assess the robustness of the measures in two series of experiments corresponding to (i) changing the dynamic range and (ii) changing the matrix size.

Results

No textural measures were robust under dynamic range changes. Entropy was the only textural feature robust under spatial resolution changes (coefficient of variation under 10% in all cases).

Conclusion

Textural measures of three-dimensional brain tumor images are not robust neither under dynamic range nor under matrix size changes. Standards should be harmonized to use textural features as imaging biomarkers in radiomic-based studies. The implications of this work go beyond the specific tumor type studied here and pose the need for standardization in textural feature calculation of oncological images.

Texture features of periaqueductal gray in the patients with medication-overuse headache

BACKGROUND

Periaqueductal gray (PAG) is the descending pain modulatory center, and PAG dysfunction had been recognized in migraine. Here we propose to investigate altered PAG texture features (quantitative approach for extracting texture descriptors for images) in the patients with medication-overuse headache (MOH) based on high resolution brain structural image to understand the MOH pathogenesis.

METHODS

The brain structural images were obtained from 32 normal controls (NC) and 44 MOH patients on 3.0 T MR system. PAG template was created based on the ICBM152 gray matter template, and the individual PAG segment was performed by applying the deformation field to the PAG template after structural image segment. Grey-level co-occurrence matrix (GLCM) was performed to measure the texture parameters including angular second moment (ASM), Contrast, Correlation, inverse difference moment (IDM) and Entropy.

RESULTS

Contrast was increased in MOH patients (9.28 ± 3.11) compared with that in NC (7.94 ± 0.65) (P < 0.05), and other texture features showed no significant difference between MOH and NC (P > 0.05). The area under the ROC curve was 0.697 for Contrast in the distinction of MOH from NC, and the cut-off value of Contrast was 8.11 with sensitivity 70.5% and specificity 62.5%. The contrast was negatively with the sleep scores (r = -0.434, P = 0.003).

CONCLUSION

Texture Contrast could be used to identify the altered MR imaging characteristics in MOH in understanding the MOH pathogenesis, and it could also be considered as imaging biomarker in for MOH diagnosis.

T2-weighted MRI-derived textural features reflect prostate cancer aggressiveness: preliminary results

PURPOSE

To evaluate the diagnostic relevance of T2-weighted (T2W) MRI-derived textural features relative to quantitative physiological parameters derived from diffusion-weighted (DW) and dynamic contrast-enhanced (DCE) MRI in Gleason score (GS) 3+4 and 4+3 prostate cancers.

MATERIALS AND METHODS

3T multiparametric-MRI was performed on 23 prostate cancer patients prior to prostatectomy. Textural features [angular second moment (ASM), contrast, correlation, entropy], apparent diffusion coefficient (ADC), and DCE pharmacokinetic parameters (K^trans and V_e) were calculated from index tumours delineated on the T2W, DW, and DCE images, respectively. The association between the textural features and prostatectomy GS and the MRI-derived parameters, and the utility of the parameters in differentiating between GS 3+4 and 4+3 prostate cancers were assessed statistically.

RESULTS

ASM and entropy correlated significantly (p < 0.05) with both GS and median ADC. Contrast correlated moderately with median ADC. The textural features correlated insignificantly with K^trans and V_e. GS 4+3 cancers had significantly lower ASM and higher entropy than 3+4 cancers, but insignificant differences in median ADC, K^trans, and V_e. The combined texture-MRI parameters yielded higher classification accuracy (91%) than the individual parameter sets.

CONCLUSION

T2W MRI-derived textural features could serve as potential diagnostic markers, sensitive to the pathological differences in prostate cancers.

KEY POINTS

• T2W MRI-derived textural features correlate significantly with Gleason score and ADC. • T2W MRI-derived textural features differentiate Gleason score 3+4 from 4+3 cancers. • T2W image textural features could augment tumour characterization.

Automatic weighing attribute to retrieve similar lung cancer nodules

BACKGROUND

Cancer is a disease characterized as an uncontrolled growth of abnormal cells that invades neighboring tissues and destroys them. Lung cancer is the primary cause of cancer-related deaths in the world, and it diagnosis is a complex task for specialists and it presents some big challenges as medical image interpretation process, pulmonary nodule detection and classification. In order to aid specialists in the early diagnosis of lung cancer, computer assistance must be integrated in the imaging interpretation and pulmonary nodule classification processes. Methods of Content-Based Image Retrieval (CBIR) have been described as one promising technique to computer-aided diagnosis and is expected to aid radiologists on image interpretation with a second opinion. However, CBIR presents some limitations: image feature extraction process and appropriate similarity measure. The efficiency of CBIR systems depends on calculating image features that may be relevant to the case similarity analysis. When specialists classify a nodule, they are supported by information from exams, images, etc. But each information has more or less weight over decision making about nodule malignancy. Thus, finding a way to measure the weight allows improvement of image retrieval process through the assignment of higher weights to that attributes that best characterize the nodules.

METHODS

In this context, the aim of this work is to present a method to automatically calculate attribute weights based on local learning to reflect the interpretation on image retrieval process. The process consists of two stages that are performed sequentially and cyclically: Evaluation Stage and Training Stage. At each iteration the weights are adjusted according to retrieved nodules. After some iterations, it is possible reach a set of attribute weights that optimize the recovery of similar nodes.

RESULTS

The results achieved by updated weights were promising because was possible increase precision by 10% to 6% on average to retrieve of benign and malignant nodules, respectively, with recall of 25% compared with tests without weights associated to attributes in similarity metric. The best result, we reaching values over 100% of precision average until thirtieth lung cancer nodule retrieved.

CONCLUSIONS

Based on the results, WED applied to the three vectors used attributes (3D TA, 3D MSA and InV), with weights adjusted by the process, always achieved better results than those found with ED. With the weights, the Precision was increased on average by 17.3% compared with using ED.

Applications and limitations of radiomics

Radiomics is an emerging field in quantitative imaging that uses advanced imaging features to objectively and quantitatively describe tumour phenotypes. Radiomic features have recently drawn considerable interest due to its potential predictive power for treatment outcomes and cancer genetics, which may have important applications in personalized medicine. In this technical review, we describe applications and challenges of the radiomic field. We will review radiomic application areas and technical issues, as well as proper practices for the designs of radiomic studies.

Radiomics: a new application from established techniques

The increasing use of biomarkers in cancer have led to the concept of personalized medicine for patients. Personalized medicine provides better diagnosis and treatment options available to clinicians. Radiological imaging techniques provide an opportunity to deliver unique data on different types of tissue. However, obtaining useful information from all radiological data is challenging in the era of "big data". Recent advances in computational power and the use of genomics have generated a new area of research termed Radiomics. Radiomics is defined as the high throughput extraction of quantitative imaging features or texture (radiomics) from imaging to decode tissue pathology and creating a high dimensional data set for feature extraction. Radiomic features provide information about the gray-scale patterns, inter-pixel relationships. In addition, shape and spectral properties can be extracted within the same regions of interest on radiological images. Moreover, these features can be further used to develop computational models using advanced machine learning algorithms that may serve as a tool for personalized diagnosis and treatment guidance.

MRI Texture Analysis Reveals Bulbar Abnormalities in Friedreich Ataxia

BACKGROUND AND PURPOSE

Texture analysis is an image processing technique that can be used to extract parameters able to describe meaningful features of an image or ROI. Texture analysis based on the gray level co-occurrence matrix gives a second-order statistical description of the image or ROI. In this work, the co-occurrence matrix texture approach was used to extract information from brain MR images of patients with Friedreich ataxia and a control group, to see whether texture parameters were different between these groups. A longitudinal analysis was also performed.

MATERIALS AND METHODS

Twenty patients and 21 healthy controls participated in the study. Both groups had 2 sets of T1-weighted MR images obtained 1 year apart for every subject. ROIs chosen for analysis were the medulla oblongata and pons. Texture parameters were obtained for these ROIs for every subject, for the 2 sets of images. These parameters were compared longitudinally within groups and transversally between groups.

RESULTS

The comparison between patients and the control group showed a significant differences for the medulla oblongata (t test, P < .05, Bonferroni-corrected) but did not show a statistically significant difference for the pons. Longitudinal comparison of images obtained 1 year apart did not show differences for either patients or for controls, in any of the analyzed structures.

CONCLUSIONS

Gray level co-occurrence matrix-based texture analysis showed statistically significant differences for the medulla oblongata of patients with Friedreich ataxia compared with controls. These results highlight the medulla as an important site of damage in Friedreich ataxia.

Directional fractal signature methods for trabecular bone texture in hand radiographs: data from the Osteoarthritis Initiative

PURPOSE

To develop directional fractal signature methods for the analysis of trabecular bone (TB) texture in hand radiographs. Problems associated with the small size of hand bones and the orientation of fingers were addressed.

METHODS

An augmented variance orientation transform (AVOT) and a quadrant rotating grid (QRG) methods were developed. The methods calculate fractal signatures (FSs) in different directions. Unlike other methods they have the search region adjusted according to the size of bone region of interest (ROI) to be analyzed and they produce FSs defined with respect to any chosen reference direction, i.e., they work for arbitrary orientation of fingers. Five parameters at scales ranging from 2 to 14 pixels (depending on image size and method) were derived from rose plots of Hurst coefficients, i.e., FS in dominating roughness (FSSta), vertical (FSV) and horizontal (FSH) directions, aspect ratio (StrS), and direction signatures (StdS), respectively. The accuracy in measuring surface roughness and isotropy/anisotropy was evaluated using 3600 isotropic and 800 anisotropic fractal surface images of sizes between 20 × 20 and 64 × 64 pixels. The isotropic surfaces had FDs ranging from 2.1 to 2.9 in steps of 0.1, and the anisotropic surfaces had two dominating directions of 30° and 120°. The methods were used to find differences in hand TB textures between 20 matched pairs of subjects with (cases: approximate Kellgren-Lawrence (KL) grade ≥ 2) and without (controls: approximate KL grade <2) radiographic hand osteoarthritis (OA). The OA Initiative public database was used and 20 × 20 pixel bone ROIs were selected on 5th distal and middle phalanges. The performance of the AVOT and QRG methods was compared against a variance orientation transform (VOT) method developed earlier [M. Wolski, P. Podsiadlo, and G. W. Stachowiak, "Directional fractal signature analysis of trabecular bone: evaluation of different methods to detect early osteoarthritis in knee radiographs," Proc. Inst. Mech. Eng., Part H 223, 211-236 (2009)].

RESULTS

The AVOT method correctly quantified the isotropic and anisotropic surfaces for all image sizes and scales. Values of FSSta were significantly different (P < 0.05) between the isotropic surfaces. Using the VOT and QRG methods no differences were found at large scales for the isotropic surfaces that are smaller than 64 × 64 and 48 × 48 pixels, respectively, and at some scales for the anisotropic surfaces with size 48 × 48 pixels. Compared to controls, using the AVOT and QRG methods the authors found that OA TB textures were less rough (P < 0.05) in the dominating and horizontal directions (i.e., lower FSSta and FSH), rougher in the vertical direction (i.e., higher FSV) and less anisotropic (i.e., higher StrS) than controls. No differences were found using the VOT method.

CONCLUSIONS

The AVOT method is well suited for the analysis of bone texture in hand radiographs and it could be potentially useful for early detection and prediction of hand OA.

Three-dimensional textural features of conventional MRI improve diagnostic classification of childhood brain tumours

The aim of this study was to assess the efficacy of three-dimensional texture analysis (3D TA) of conventional MR images for the classification of childhood brain tumours in a quantitative manner. The dataset comprised pre-contrast T1 - and T2-weighted MRI series obtained from 48 children diagnosed with brain tumours (medulloblastoma, pilocytic astrocytoma and ependymoma). 3D and 2D TA were carried out on the images using first-, second- and higher order statistical methods. Six supervised classification algorithms were trained with the most influential 3D and 2D textural features, and their performances in the classification of tumour types, using the two feature sets, were compared. Model validation was carried out using the leave-one-out cross-validation (LOOCV) approach, as well as stratified 10-fold cross-validation, in order to provide additional reassurance. McNemar's test was used to test the statistical significance of any improvements demonstrated by 3D-trained classifiers. Supervised learning models trained with 3D textural features showed improved classification performances to those trained with conventional 2D features. For instance, a neural network classifier showed 12% improvement in area under the receiver operator characteristics curve (AUC) and 19% in overall classification accuracy. These improvements were statistically significant for four of the tested classifiers, as per McNemar's tests. This study shows that 3D textural features extracted from conventional T1 - and T2-weighted images can improve the diagnostic classification of childhood brain tumours. Long-term benefits of accurate, yet non-invasive, diagnostic aids include a reduction in surgical procedures, improvement in surgical and therapy planning, and support of discussions with patients' families. It remains necessary, however, to extend the analysis to a multicentre cohort in order to assess the scalability of the techniques used.

Can Quantitative CT Texture Analysis be Used to Differentiate Fat-poor Renal Angiomyolipoma from Renal Cell Carcinoma on Unenhanced CT Images?

PURPOSE

To determine the accuracy of texture analysis to differentiate fat-poor angiomyolipoma (fp-AML) from renal cell carcinoma (RCC) on unenhanced computed tomography (CT) images.

MATERIALS AND METHODS

In this institutional review board-approved retrospective case-control study, patients with AML and RCC were identified from the pathology database: there were 16 patients with fp-AML (no visible fat at unenhanced CT) and 84 patients with RCC. Axial unenhanced CT images were contoured manually by two independent analysts. Texture analysis was performed for each lesion, and reproducibility was assessed. Texture features related to the gray-level histogram, gray-level co-occurrence, and run-length matrix statistics were evaluated. The most discriminative features were used to generate support vector machine (SVM) classifiers. Diagnostic accuracy of textural features was assessed and 10-fold cross validation was performed. Unenhanced CT images for each patient were independently reviewed by two blinded radiologists who subjectively graded lesion heterogeneity on a five-point scale. Differences in area under the receiver operating characteristic curve (AUC) between subjective heterogeneity ratings and textural features were evaluated by using the DeLong method.

RESULTS

There was lower lesion homogeneity and higher lesion entropy in RCCs (P ≤ .01). A model incorporating several texture features resulted in an AUC of 0.89 ± 0.04. The average SVM accuracy of textural features ranged from 83% to 91% (after 10-fold cross validation). An optimal subjective heterogeneity rating of 2 or higher was identified as a predictor of RCC for both readers, with no significant difference in AUC between readers (P = .06). Each of the three textural-based classifiers was more accurate than either radiologists' subjective heterogeneity ratings for the models incorporating a subset of the top three textural features (difference in AUC between textural features and subjective visual heterogeneity, 0.25; 95% confidence interval: 0.02, 0.47; P = .03).

CONCLUSION

CT texture analysis can be used to accurately differentiate fp-AML from RCC on unenhanced CT images.

Multidimensional texture characterization: on analysis for brain tumor tissues using MRS and MRI

This paper investigates the efficacy of automated pattern recognition methods on magnetic resonance data with the objective of assisting radiologists in the clinical diagnosis of brain tissue tumors. In this paper, the sciences of magnetic resonance imaging (MRI) and magnetic resonance spectroscopy (MRS) are combined to improve the accuracy of the classifier, based on the multidimensional co-occurrence matrices to assess the detection of pathological tissues (tumor and edema), normal tissues (white matter - WM and gray matter - GM), and fluid (cerebrospinal fluid - CSF). The results show the ability of the classifier with iterative training to automatically and simultaneously recover tissue-specific spectral and structural patterns and achieve segmentation of tumor and edema and grading of high and low glioma tumor. Here, extreme learning machine - improved particle swarm optimization (ELM-IPSO) neural network classifier is trained with the feature descriptions in brain magnetic resonance (MR) spectra. This has the characteristics of varying the normal spectral pattern associated with tumor patterns along with imaging features. Validation was performed considering 35 clinical studies. The volumetric features extracted from the vectors of this matrix articulate some important elementary structures, which along with spectroscopic metabolite ratios discriminate the tumor grades and tissue classes. The quantitative 3D analysis reveals significant improvement in terms of global accuracy rate for automatic classification in brain tissues and discriminating pathological tumor tissue from structural healthy brain tissue.

Texture features on T2-weighted magnetic resonance imaging: new potential biomarkers for prostate cancer aggressiveness

To explore contrast (C) and homogeneity (H) gray-level co-occurrence matrix texture features on T2-weighted (T2w) Magnetic Resonance (MR) images and apparent diffusion coefficient (ADC) maps for predicting prostate cancer (PCa) aggressiveness, and to compare them with traditional ADC metrics for differentiating low- from intermediate/high-grade PCas. The local Ethics Committee approved this prospective study of 93 patients (median age, 65 years), who underwent 1.5 T multiparametric endorectal MR imaging before prostatectomy. Clinically significant (volume ≥0.5 ml) peripheral tumours were outlined on histological sections, contoured on T2w and ADC images, and their pathological Gleason Score (pGS) was recorded. C, H, and traditional ADC metrics (mean, median, 10th and 25th percentile) were calculated on the largest lesion slice, and correlated with the pGS through the Spearman correlation coefficient. The area under the receiver operating characteristic curve (AUC) assessed how parameters differentiate pGS = 6 from pGS ≥ 7. The dataset included 49 clinically significant PCas with a balanced distribution of pGS. The Spearman ρ and AUC values on ADC were: -0.489, 0.823 (mean); -0.522, 0.821 (median); -0.569, 0.854 (10th percentile); -0.556, 0.854 (25th percentile); -0.386, 0.871 (C); 0.533, 0.923 (H); while on T2w they were: -0.654, 0.945 (C); 0.645, 0.962 (H). AUC of H on ADC and T2w, and C on T2w were significantly higher than that of the mean ADC (p = 0.05). H and C calculated on T2w images outperform ADC parameters in correlating with pGS and differentiating low- from intermediate/high-risk PCas, supporting the role of T2w MR imaging in assessing PCa biological aggressiveness.

MRI texture analysis predicts p53 status in head and neck squamous cell carcinoma

BACKGROUND AND PURPOSE

Head and neck cancer is common, and understanding the prognosis is an important part of patient management. In addition to the Tumor, Node, Metastasis staging system, tumor biomarkers are becoming more useful in understanding prognosis and directing treatment. We assessed whether MR imaging texture analysis would correctly classify oropharyngeal squamous cell carcinoma according to p53 status.

MATERIALS AND METHODS

A cohort of 16 patients with oropharyngeal squamous cell carcinoma was prospectively evaluated by using standard clinical, histopathologic, and imaging techniques. Tumors were stained for p53 and scored by an anatomic pathologist. Regions of interest on MR imaging were selected by a neuroradiologist and then analyzed by using our 2D fast time-frequency transform tool. The quantified textures were assessed by using the subset-size forward-selection algorithm in the Waikato Environment for Knowledge Analysis. Features found to be significant were used to create a statistical model to predict p53 status. The model was tested by using a Bayesian network classifier with 10-fold stratified cross-validation.

RESULTS

Feature selection identified 7 significant texture variables that were used in a predictive model. The resulting model predicted p53 status with 81.3% accuracy (P < .05). Cross-validation showed a moderate level of agreement (κ = 0.625).

CONCLUSIONS

This study shows that MR imaging texture analysis correctly predicts p53 status in oropharyngeal squamous cell carcinoma with ∼80% accuracy. As our knowledge of and dependence on tumor biomarkers expand, MR imaging texture analysis warrants further study in oropharyngeal squamous cell carcinoma and other head and neck tumors.

Role of the texture features of images in the diagnosis of solitary pulmonary nodules in different sizes

OBJECTIVE

To explore the role of the texture features of images in the diagnosis of solitary pulmonary nodules (SPNs) in different sizes.

MATERIALS AND METHODS

A total of 379 patients with pathologically confirmed SPNs were enrolled in this study. They were divided into three groups based on the SPN sizes: ≤10, 11-20, and >20 mm. Their texture features were segmented and extracted. The differences in the image features between benign and malignant SPNs were compared. The SPNs in these three groups were determined and analyzed with the texture features of images.

RESULTS

These 379 SPNs were successfully segmented using the 2D Otsu threshold method and the self-adaptive threshold segmentation method. The texture features of these SPNs were obtained using the method of grey level co-occurrence matrix (GLCM). Of these 379 patients, 120 had benign SPNs and 259 had malignant SPNs. The entropy, contrast, energy, homogeneity, and correlation were 3.5597±0.6470, 0.5384±0.2561, 0.1921±0.1256, 0.8281±0.0604, and 0.8748±0.0740 in the benign SPNs and 3.8007±0.6235, 0.6088±0.2961, 0.1673±0.1070, 0.7980±0.0555, and 0.8550±0.0869 in the malignant SPNs (all P<0.05). The sensitivity, specificity, and accuracy of the texture features of images were 83.3%, 90.0%, and 86.8%, respectively, for SPNs sized ≤10 mm, and were 86.6%, 88.2%, and 87.1%, respectively, for SPNs sized 
11-20 mm and 94.7%, 91.8%, and 93.9%, respectively, for SPNs sized >20 mm.

CONCLUSIONS

The entropy and contrast of malignant pulmonary nodules have been demonstrated to be higher in comparison to those of benign pulmonary nodules, while the energy, homogeneity correlation of malignant pulmonary nodules are lower than those of benign pulmonary nodules. The texture features of images can reflect the tissue features and have high sensitivity, specificity, and accuracy in differentiating SPNs. The sensitivity and accuracy increase for larger SPNs.

Classification of dynamic contrast enhanced MR images of cervical cancers using texture analysis and support vector machines

Dynamic contrast enhanced MRI (DCE-MRI) provides insight into the vascular properties of tissue. Pharmacokinetic models may be fitted to DCE-MRI uptake patterns, enabling biologically relevant interpretations. The aim of our study was to determine whether treatment outcome for 81 patients with locally advanced cervical cancer could be predicted from parameters of the Brix pharmacokinetic model derived from pre-chemoradiotherapy DCE-MRI. First-order statistical features of the Brix parameters were used. In addition, texture analysis of Brix parameter maps was done by constructing gray level co-occurrence matrices (GLCM) from the maps. Clinical factors and first- and second-order features were used as explanatory variables for support vector machine (SVM) classification, with treatment outcome as response. Classification models were validated using leave-one-out cross-model validation. A random value permutation test was used to evaluate model significance. Features derived from first-order statistics could not discriminate between cured and relapsed patients (specificity 0%-20%, p-values close to unity). However, second-order GLCM features could significantly predict treatment outcome with accuracies (~70%) similar to the clinical factors tumor volume and stage (69%). The results indicate that the spatial relations within the tumor, quantified by texture features, were more suitable for outcome prediction than first-order features.

The role of texture analysis in imaging as an outcome predictor and potential tool in radiotherapy treatment planning

Predicting a tumour's response to radiotherapy prior to the start of treatment could enhance clinical care management by enabling the personalization of treatment plans based on predicted outcome. In recent years, there has been accumulating evidence relating tumour texture to patient survival and response to treatment. Tumour texture could be measured from medical images that provide a non-invasive method of capturing intratumoural heterogeneity and hence could potentially enable a prior assessment of a patient's predicted response to treatment. In this article, work presented in the literature regarding texture analysis in radiotherapy in relation to survival and outcome is discussed. Challenges facing integrating texture analysis in radiotherapy planning are highlighted and recommendations for future directions in research are suggested.

Structural integrity of medial temporal lobes of patients with acute mild traumatic brain injury

Post-traumatic amnesia (PTA) is an acute characteristic of traumatic brain injury (TBI) and the duration of PTA is commonly used to estimate the severity of brain injury. In the context of mild traumatic brain injury (MTBI), PTA is an essential part of the routine clinical assessment. Macroscopic lesions in temporal lobes, especially hippocampal regions, are thought to be connected to memory loss. However, conventional neuroimaging has failed to reveal neuropathological correlates of PTA in MTBI. Texture analysis (TA) is an image analysis technique that quantifies the minor MRI signal changes among image pixels and, therefore, the variations in intensity patterns within the image. The objective of this work was to apply the TA technique to MR images of MTBI patients and control subjects, and to assess the microstructural damage in medial temporal lobes of patients with MTBI with definite PTA. TA was performed for fluid-attenuated inversion recovery (FLAIR) images of 50 MTBI patients and 50 age- and gender-matched controls in the regions of the amygdala, hippocampus, and thalamus. It was hypothesized that 1) there would be statistically significant differences in TA parameters between patients with MTBIs and controls, and 2) the duration of PTA would be related to TA parameters in patients with MTBI. No significant textural differences were observed between patients and controls in the regions of interest (p>0.01). No textural features were observed to correlate with the duration of PTA. Subgroup analyses were conducted on patients with PTA of>1 h, (n=33) and compared the four TA parameters to the age- and gender-matched controls (n=33). The findings were similar. This study did not reveal significant textural changes in medial temporal structures that could be related to the duration of PTA.

Clinical and magnetic resonance imaging (MRI) distinctions between tumefactive demyelination and brain tumors in children

Tumefactive demyelinating lesions can be difficult to distinguish from tumors. Clinical and magnetic resonance imaging features of children with tumefactive demyelination and supratentorial brain tumors were compared. Patients were identified through a 23-site national demyelinating disease study, and from a single-site neuroradiology database. For inclusion, lesions met at least 1 of 3 criteria: maximal cross-sectional diameter >20 mm, local or global cerebral mass effect, or presence of perilesional edema. Thirty-one children with tumefactive demyelination (5 with solitary lesions) were identified: 27 of 189 (14.3%) from the demyelinating disease study and 4 from the database. Thirty-three children with tumors were identified. Children with tumefactive demyelination were more likely to have an abnormal neurologic examination and polyfocal neurologic deficits compared to children with tumors. Tumefactive demyelination was distinguished from tumor by the presence of multiple lesions, absence of cortical involvement, and decrease in lesion size or detection of new lesions on serial imaging.

Progress on the diagnosis and evaluation of brain tumors

Brain tumors are one of the most challenging disorders encountered, and early and accurate diagnosis is essential for the management and treatment of these tumors. In this article, diagnostic modalities including single-photon emission computed tomography, positron emission tomography, magnetic resonance imaging, and optical imaging are reviewed. We mainly focus on the newly emerging, specific imaging probes, and their potential use in animal models and clinical settings.