Three-Dimensional Selective-Scale Texture Analysis of Computed Tomography Pulmonary Angiograms:

PAHPhysRAD: A Digital Imaging and Communications in Medicine Research Tool for Segmentation and Radiomic Feature Extraction

Introduction

Segmentation and analysis of organs at risks (OARs) and tumor volumes are integral concepts in the development of radiotherapy treatment plans and prediction of patients' treatment outcomes.

Aims

To develop a research tool, PAHPhysRAD, that can be used to semi- and fully automate segmentation of OARs. In addition, the proposed software seeks to extract 3214 radiomic features from tumor volumes and user-specified dose-volume parameters.

Materials and Methods

Developed within MATLAB, PAHPhysRAD provides a comprehensive suite of segmentation tools, including manual, semi-automatic, and automatic options. For semi-autosegmentation, meta AI's Segment Anything Model was incorporated using the bounding box methods. Autosegmentation of OARs and tumor volume are implemented through a module that enables the addition of models in Open Neural Network Exchange format. To validate the radiomic feature extraction module in PAHPhysRAD, radiomic features extracted from gross tumor volume of 15 non-small cell lung carcinoma patients were compared against the features extracted from 3D Slicer™. The dose-volume parameters extraction module was validated using the dose volume data extracted from 28 tangential field-based breast treatment planning datasets. The volume receiving ≥20 Gy (V20) for ipsilateral lung and the mean doses received by the heart and ipsilateral lung, were compared against the parameters extracted from Eclipse.

Results

The Wilcoxon signed-rank test revealed no significant difference between the majority of the radiomic features derived from PAHPhysRAD and 3D Slicer. The average mean lung and heart doses calculated in Eclipse were 5.51 ± 2.28 Gy and 1.64 ± 1.98 Gy, respectively. Similarly, the average mean lung and heart doses calculated in PAHPhysRAD were 5.45 ± 2.89 Gy and 1.67 ± 2.08 Gy, respectively.

Conclusion

The MATLAB-based graphical user interface, PAHPhysRAD, offers a user-friendly platform for viewing and analyzing medical scans with options to extract radiomic features and dose-volume parameters. Its versatility, compatibility, and potential for further development make it an asset in medical image analysis.

Chest CT texture-based radiomics analysis in differentiating COVID-19 from other interstitial pneumonia

PURPOSE

To evaluate the potential role of texture-based radiomics analysis in differentiating Coronavirus Disease-19 (COVID-19) pneumonia from pneumonia of other etiology on Chest CT.

MATERIALS AND METHODS

One hundred and twenty consecutive patients admitted to Emergency Department, from March 8, 2020, to April 25, 2020, with suspicious of COVID-19 that underwent Chest CT, were retrospectively analyzed. All patients presented CT findings indicative for interstitial pneumonia. Sixty patients with positive COVID-19 real-time reverse transcription polymerase chain reaction (RT-PCR) and 60 patients with negative COVID-19 RT-PCR were enrolled. CT texture analysis (CTTA) was manually performed using dedicated software by two radiologists in consensus and textural features on filtered and unfiltered images were extracted as follows: mean intensity, standard deviation (SD), entropy, mean of positive pixels (MPP), skewness, and kurtosis. Nonparametric Mann-Whitney test assessed CTTA ability to differentiate positive from negative COVID-19 patients. Diagnostic criteria were obtained from receiver operating characteristic (ROC) curves.

RESULTS

Unfiltered CTTA showed lower values of mean intensity, MPP, and kurtosis in COVID-19 positive patients compared to negative patients (p = 0.041, 0.004, and 0.002, respectively). On filtered images, fine and medium texture scales were significant differentiators; fine texture scale being most significant where COVID-19 positive patients had lower SD (p = 0.004) and MPP (p = 0.004) compared to COVID-19 negative patients. A combination of the significant texture features could identify the patients with positive COVID-19 from negative COVID-19 with a sensitivity of 60% and specificity of 80% (p = 0.001).

CONCLUSIONS

Preliminary evaluation suggests potential role of CTTA in distinguishing COVID-19 pneumonia from other interstitial pneumonia on Chest CT.

Texture Analysis in Brain Tumor MR Imaging

Texture analysis, as well as its broader category radiomics, describes a variety of techniques for image analysis that quantify the variation in surface intensity or patterns, including some that are imperceptible to the human visual system. Cerebral gliomas have been most rigorously studied in brain tumors using MR-based texture analysis (MRTA) to determine the correlation of various clinical measures with MRTA features. Promising results in cerebral gliomas have been shown in the previous MRTA studies in terms of the correlation with the World Health Organization grades, risk stratification in gliomas, and the differentiation of gliomas from other brain tumors. Multiple MRTA studies in gliomas have repeatedly shown high performance of entropy, a measure of the randomness in image intensity values, of either histogram- or gray-level co-occurrence matrix parameters. Similarly, researchers have applied MRTA to other brain tumors, including meningiomas and pediatric posterior fossa tumors.

However, the value of MRTA in the clinical use remains undetermined, probably because previous studies have shown only limited reproducibility of the result in the real world. The low-to-modest generalizability may be attributed to variations in MRTA methods, sampling bias that originates from single-institution studies, and overfitting problems to a limited number of samples.

To enhance the reliability and reproducibility of MRTA studies, researchers have realized the importance of standardizing methods in the field of radiomics. Another advancement is the recent development of a comprehensive assessment system to ensure the quality of a radiomics study. These two-way approaches will secure the validity of upcoming MRTA studies. The clinical use of texture analysis in brain MRI will be accelerated by these continuous efforts.

The role of radiomics in prostate cancer radiotherapy

"Radiomics," as it refers to the extraction and analysis of a large number of advanced quantitative radiological features from medical images using high-throughput methods, is perfectly suited as an engine for effectively sifting through the multiple series of prostate images from before, during, and after radiotherapy (RT). Multiparametric (mp)MRI, planning CT, and cone beam CT (CBCT) routinely acquired throughout RT and the radiomics pipeline are developed for extraction of thousands of variables. Radiomics data are in a format that is appropriate for building descriptive and predictive models relating image features to diagnostic, prognostic, or predictive information. Prediction of Gleason score, the histopathologic cancer grade, has been the mainstay of the radiomic efforts in prostate cancer. While Gleason score (GS) is still the best predictor of treatment outcome, there are other novel applications of quantitative imaging that are tailored to RT. In this review, we summarize the radiomics efforts and discuss several promising concepts such as delta-radiomics and radiogenomics for utilizing image features for assessment of the aggressiveness of prostate cancer and its outcome. We also discuss opportunities for quantitative imaging with the advance of instrumentation in MRI-guided therapies.

Prognostic Value of Texture Analysis Based on Pretreatment DWI-Weighted MRI for Esophageal Squamous Cell Carcinoma Patients Treated With Concurrent Chemo-Radiotherapy

Purpose: The purpose of the research was to assess the prognostic value of three-dimensional (3D) texture features based on diffusion-weighted magnetic resonance imaging (DWI) for esophageal squamous cell carcinoma (ESCC) patients undergoing concurrent chemo-radiotherapy (CRT). Methods: We prospectively enrolled 82 patients with ESCC into a cohort study. Two DWI sequences (b = 0 and b = 600 s/mm²) were acquired along with axial T2WI and T1WI before CRT. Two groups of features were examined: (1) clinical and demographic features (e.g., TNM stage, age and sex) and (2) changes in spatial texture characteristics of the apparent diffusion coefficient (ADC), which characterizes gray intensity changes in tumor areas, spatial pattern and distribution, and related changes caused by CRT. Reproducible feature sets without redundancy were statistically filtered and validated. The prognostic values associated with overall survival (OS) for each parameter were studied using Kaplan-Meier and Cox regression models for univariate and multivariate analyses, respectively. Results: Both univariate and multivariate Cox model analyses showed that the energy of intensity histogram texture (IHIST_energy), radiation dose, mean of the contrast in distance 1 of 26 directions (m_contrast_1), extreme difference of the homogeneity in distance 2 of 26 directions (Diff_homogeneity_2), mean of the inverse variance in distance 2 of 26 directions (m_lnversevariance_2), high-intensity small zone emphasis (HISE), and low-intensity large zone emphasis (LILE) were significantly associated with survival. The results showed that 6 texture parameters extracted from the ADC images before treatment could distinguish among high-, medium-, and low-risk groups (log-rank χ² = 9.7; P = 0.00773). The biased C-index value was 0.715 (95% CI: 0.708 to 0.732) based on bootstrapping validation. Conclusions: The ADC 3D texture feature can be used as a useful biomarker to predict the survival of ESCC patients undergoing CRT. Combining ADC 3D texture features with conventional prognostic factors can generate reliable survival prediction models.

Synergistic application of pulmonary 18F-FDG PET/HRCT and computer-based CT analysis with conventional severity measures to refine current risk stratification in idiopathic pulmonary fibrosis (IPF)

INTRODUCTION

To investigate the combined performance of quantitative CT (qCT) following a computer algorithm analysis (IMBIO) and ¹⁸F-FDG PET/CT to assess survival in patients with idiopathic pulmonary fibrosis (IPF).

METHODS

A total of 113 IPF patients (age 70 ± 9 years) prospectively and consecutively underwent ¹⁸F-FDG PET/CT and high-resolution CT (HRCT) at our institution. During a mean follow-up of 29.6 ± 26 months, 44 (48%) patients died. As part of the qCT analysis, pattern evaluation of HRCT (using IMBIO software) included the total extent (percentage) of the following features: normal-appearing lung, hyperlucent lung, parenchymal damage (comprising ground-glass opacification, reticular pattern and honeycombing), and the pulmonary vessels. The maximum (SUV_max) and minimum (SUV_min) standardized uptake value (SUV) for ¹⁸F-FDG uptake in the lungs, and the target-to-background (SUV_max/SUV_min) ratio (TBR) were quantified using routine region-of-interest (ROI) analysis. Pulmonary functional tests (PFTs) were acquired within 14 days of the PET/CT/HRCT scan. Kaplan-Meier (KM) survival analysis was used to identify associations with mortality.

RESULTS

Data from 91 patients were available for comparative analysis. The average ± SD GAP [gender, age, physiology] score was 4.2 ± 1.7 (range 0-8). The average ± SD SUV_max, SUV_min, and TBR were 3.4 ± 1.4, 0.7 ± 0.2, and 5.6 ± 2.8, respectively. In all patients, qCT analysis demonstrated a predominantly reticular lung pattern (14.9 ± 12.4%). KM analysis showed that TBR (p = 0.018) and parenchymal damage assessed by qCT (p = 0.0002) were the best predictors of survival. Adding TBR and qCT to the GAP score significantly increased the ability to differentiate between high and low risk (p < 0.0001).

CONCLUSION

¹⁸F-FDG PET and qCT are independent and synergistic in predicting mortality in patients with IPF.

Magnetic resonance imaging (MRI)-based radiomics for prostate cancer radiotherapy

In radiotherapy (RT) of prostate cancer, dose escalation has been shown to reduce biochemical failure. Dose escalation only to determinate prostate tumor habitats has the potential to improve tumor control with less toxicity than when the entire prostate is dose escalated. Other issues in the treatment of the RT patient include the choice of the RT technique (hypo- or standard fractionation) and the use and length of concurrent/adjuvant androgen deprivation therapy (ADT). Up to 50% of high-risk men demonstrate biochemical failure suggesting that additional strategies for defining and treating patients based on improved risk stratification are required. The use of multiparametric MRI (mpMRI) is rapidly gaining momentum in the management of prostate cancer because of its improved diagnostic potential and its ability to combine functional and anatomical information. Currently, the Prostate Imaging, Reporting and Diagnosis System (PIRADS) is the standard of care for region of interest (ROI) identification and risk classification. However, PIRADS was not designed for 3D tumor volume delineation; there is a large degree of subjectivity and PIRADS does not accurately and reproducibly elucidate inter- and intra-lesional spatial heterogeneity. "Radiomics", as it refers to the extraction and analysis of large number of advanced quantitative radiological features from medical images using high throughput methods, is perfectly suited as an engine to effectively sift through the multiple series of prostate mpMRI sequences and quantify regions of interest. The radiomic efforts can be summarized in two main areas: (I) detection/segmentation of the suspicious lesion; and (II) assessment of the aggressiveness of prostate cancer. As related to RT, the goal of the latter is in particular to identify patients at high risk for metastatic disease; and the aim of the former is to identify and segment cancerous lesions and thus provide targets for radiation boost. The article is structured as follows: first, we describe the radiomic approach; and second, we discuss the radiomic pipeline as tailored for RT of prostate cancer. In this process we summarize the current efforts and progress in integrating mpMRI radiomics into the radiotherapeutic management of prostate cancer with emphasis placed on its role in treatment target definition, treatment plan strategizing, and prognostic assessment. The described concepts, methods and tools are not currently applicable to the radiation oncology practice outside of the research setting. More data are required in the form of clinical trials to assess the robustness of radiomics-based predictive models, and to maximize the efficacy of these models.

Benefits of texture analysis of dual energy CT for Computer-Aided pulmonary embolism detection

Pulmonary embolism is an avoidable cause of death if treated immediately but delays in diagnosis and treatment lead to an increased risk. Computer-assisted image analysis of both unenhanced and contrast-enhanced computed tomography (CT) have proven useful for diagnosis of pulmonary embolism. Dual energy CT provides additional information over the standard single energy scan by generating four-dimensional (4D) data, in our case with 11 energy levels in 3D. In this paper a 4D texture analysis method capable of detecting pulmonary embolism in dual energy CT is presented. The method uses wavelet-based visual words together with an automatic geodesic-based region of interest detection algorithm to characterize the texture properties of each lung lobe. Results show an increase in performance with respect to the single energy CT analysis, as well as an accuracy gain compared to preliminary work on a small dataset.

Grey-matter texture abnormalities and reduced hippocampal volume are distinguishing features of schizophrenia

Neurodevelopmental processes are widely believed to underlie schizophrenia. Analysis of brain texture from conventional magnetic resonance imaging (MRI) can detect disturbance in brain cytoarchitecture. We tested the hypothesis that patients with schizophrenia manifest quantitative differences in brain texture that, alongside discrete volumetric changes, may serve as an endophenotypic biomarker. Texture analysis (TA) of grey matter distribution and voxel-based morphometry (VBM) of regional brain volumes were applied to MRI scans of 27 patients with schizophrenia and 24 controls. Texture parameters (uniformity and entropy) were also used as covariates in VBM analyses to test for correspondence with regional brain volume. Linear discriminant analysis tested if texture and volumetric data predicted diagnostic group membership (schizophrenia or control). We found that uniformity and entropy of grey matter differed significantly between individuals with schizophrenia and controls at the fine spatial scale (filter width below 2mm). Within the schizophrenia group, these texture parameters correlated with volumes of the left hippocampus, right amygdala and cerebellum. The best predictor of diagnostic group membership was the combination of fine texture heterogeneity and left hippocampal size. This study highlights the presence of distributed grey-matter abnormalities in schizophrenia, and their relation to focal structural abnormality of the hippocampus. The conjunction of these features has potential as a neuroimaging endophenotype of schizophrenia.

Enhanced visualization of pulmonary perfusion in 4D dual energy CT images

Pulmonary embolism (PE) affects up to 600,000 patients and contributes to at least 100,000 deaths every year in the United States alone. Diagnosis of PE can be difficult as most symptoms are unspecific. Computed Tomography (CT) angiography is the reference for diagnosing PE. CT angiography produces grayscale images with darker areas representing any mass filling defects, making the analysis of the images difficult. This article demonstrates a method using the combination of energy levels in Dual Energy CT images to highlight the presence of PE in the lung. The results show that pairing different energy levels from 40 to 140 keV can increase the contrast between well perfused areas and underperfused areas of the lung. In addition, the visualization used in the current study complies with the window/level settings usually employed by radiologists.

Three-dimensional solid texture analysis in biomedical imaging: review and opportunities

Three-dimensional computerized characterization of biomedical solid textures is key to large-scale and high-throughput screening of imaging data. Such data increasingly become available in the clinical and research environments with an ever increasing spatial resolution. In this text we exhaustively analyze the state-of-the-art in 3-D biomedical texture analysis to identify the specific needs of the application domains and extract promising trends in image processing algorithms. The geometrical properties of biomedical textures are studied both in their natural space and on digitized lattices. It is found that most of the tissue types have strong multi-scale directional properties, that are well captured by imaging protocols with high resolutions and spherical spatial transfer functions. The information modeled by the various image processing techniques is analyzed and visualized by displaying their 3-D texture primitives. We demonstrate that non-convolutional approaches are expected to provide best results when the size of structures are inferior to five voxels. For larger structures, it is shown that only multi-scale directional convolutional approaches that are non-separable allow for an unbiased modeling of 3-D biomedical textures. With the increase of high-resolution isotropic imaging protocols in clinical routine and research, these models are expected to best leverage the wealth of 3-D biomedical texture analysis in the future. Future research directions and opportunities are proposed to efficiently model personalized image-based phenotypes of normal biomedical tissue and its alterations. The integration of the clinical and genomic context is expected to better explain the intra class variation of healthy biomedical textures. Using texture synthesis, this provides the exciting opportunity to simulate and visualize texture atlases of normal ageing process and disease progression for enhanced treatment planning and clinical care management.

Abnormalities in fronto-striatal connectivity within language networks relate to differences in grey-matter heterogeneity in Asperger syndrome

Asperger syndrome (AS) is an Autism Spectrum Disorder (ASD) characterised by qualitative impairment in the development of emotional and social skills with relative preservation of general intellectual abilities, including verbal language. People with AS may nevertheless show atypical language, including rate and frequency of speech production. We previously observed that abnormalities in grey matter homogeneity (measured with texture analysis of structural MR images) in AS individuals when compared with controls are also correlated with the volume of caudate nucleus. Here, we tested a prediction that these distributed abnormalities in grey matter compromise the functional integrity of brain networks supporting verbal communication skills. We therefore measured the functional connectivity between caudate nucleus and cortex during a functional neuroimaging study of language generation (verbal fluency), applying psycho-physiological interaction (PPI) methods to test specifically for differences attributable to grey matter heterogeneity in AS participants. Furthermore, we used dynamic causal modelling (DCM) to characterise the causal directionality of these differences in interregional connectivity during word production. Our results revealed a diagnosis-dependent influence of grey matter heterogeneity on the functional connectivity of the caudate nuclei with right insula/inferior frontal gyrus and anterior cingulate, respectively with the left superior frontal gyrus and right precuneus. Moreover, causal modelling of interactions between inferior frontal gyri, caudate and precuneus, revealed a reliance on bottom-up (stimulus-driven) connections in AS participants that contrasted with a dominance of top-down (cognitive control) connections from prefrontal cortex observed in control participants. These results provide detailed support for previously hypothesised central disconnectivity in ASD and specify discrete brain network targets for diagnosis and therapy in ASD.

•

We used MRI techniques to assess the connectivity in language networks in AS.

•

Grey-matter heterogeneity of MR images correlated with volume of caudate in AS.

•

Hence, caudate nuclei were used as seed ROIs in connectivity analyses: PPI, DCM.

•

Grey-matter heterogeneity differently tuned caudate connectivity in AS, controls.

•

DCM of language circuitry featured bottom-up models in AS and top-down in controls.

Measurements of heterogeneity in gliomas on computed tomography relationship to tumour grade

To undertake a preliminary study that uses CT texture analysis (CTTA) to quantify heterogeneity in gliomas on contrast-enhanced CT and to assess the relationship between tumour heterogeneity and grade. Retrospective analysis of contrast enhanced CT images was performed in 44 patients with histologically proven cerebral glioma between 2007 and 2010. 11 tumours were low grade (Grade I = 3; Grade II, = 8) and 33 high grade (Grade III = 10, Grade IV = 23). CTTA assessment of tumour heterogeneity was performed using a proprietary software algorithm (TexRAD) that selectively filters and extracts textures at different anatomical scales between filter values 1.0 (fine detail) and 2.5 (coarse features). Heterogeneity was quantified as standard deviation (SD) with or without filtration. Tumour heterogeneity, size and attenuation were correlated with tumour grade. For each parameter, receiver operating characteristics characterised the diagnostic performance for discrimination of high grade from low grade glioma and of grade III tumours from grade IV. Further the CTTA was compared to the radiological diagnosis. Tumour heterogeneity correlated significantly with grade (SD without filtration rs = 0.664, p < 0.001, SD with coarse filtration (rs = 0.714, p < 0.001). Tumour size and attenuation showed only moderate correlations with tumour grade (rs = 0.426, p = 0.004 and rs = 0.447, p = 0.002 respectively). Coarse texture was the best discriminator between high and low grade tumours (AUC 0.832, p < 0.0001) and between grade III and grade IV gliomas (AUC = 0.878 p = 0.0001). Compared to the radiological diagnosis, CTTA further characterised the indetermined cases. By quantifying tumour heterogeneity, CTTA has the potential to provide a marker of tumour grade for patients with cerebral glioma. By differentiating between high and low grade tumours, CTTA could possibly assist clinical management.

Heterogeneity of focal breast lesions and surrounding tissue assessed by mammographic texture analysis: preliminary evidence of an association with tumor invasion and estrogen receptor status

AIM

This pilot study investigates whether heterogeneity in focal breast lesions and surrounding tissue assessed on mammography is potentially related to cancer invasion and hormone receptor status.

MATERIALS AND METHODS

Texture analysis (TA) assessed the heterogeneity of focal lesions and their surrounding tissues in digitized mammograms from 11 patients randomly selected from an imaging archive [ductal carcinoma in situ (DCIS) only, n = 4; invasive carcinoma (IC) with DCIS, n = 3; IC only, n = 4]. TA utilized band-pass image filtration to highlight image features at different spatial frequencies (filter values: 1.0-2.5) from fine to coarse texture. The distribution of features in the derived images was quantified using uniformity.

RESULTS

Significant differences in uniformity were observed between patient groups for all filter values. With medium scale filtration (filter value = 1.5) pure DCIS was more uniform (median = 0.281) than either DCIS with IC (median = 0.246, p = 0.0102) or IC (median = 0.249, p = 0.0021). Lesions with high levels of estrogen receptor expression were more uniform, most notably with coarse filtration (filter values 2.0 and 2.5, r(s) = 0.812, p = 0.002). Comparison of uniformity values in focal lesions and surrounding tissue showed significant differences between DCIS with or without IC versus IC (p = 0.0009).

CONCLUSION

This pilot study shows the potential for computer-based assessments of heterogeneity within focal mammographic lesions and surrounding tissue to identify adverse pathological features in mammographic lesions. The technique warrants further investigation as a possible adjunct to existing computer aided diagnosis systems.

Tumour heterogeneity in oesophageal cancer assessed by CT texture analysis: preliminary evidence of an association with tumour metabolism, stage, and survival

AIM

To undertake a pilot study assessing whether tumour heterogeneity evaluated using computed tomography texture analysis (CTTA) has the potential to provide a marker of tumour aggression and prognosis in oesophageal cancer.

MATERIALS AND METHODS

In 21 patients, unenhanced CT images of the primary oesophageal lesion obtained using positron-emission tomography (PET)-CT examinations underwent CTTA. CTTA was carried out using a software algorithm that selectively filters and extracts textures at different anatomical scales between filter values 1.0 (fine detail) and 2.5 (coarse features) with quantification as entropy and uniformity (measures image heterogeneity). Texture parameters were correlated with average tumour 2-[(18)F]-fluoro-2-deoxy-d-glucose (FDG) uptake [standardized uptake values (SUV(mean) and SUV(max))] and clinical staging as determined by endoscopic ultrasound (nodal involvement) and PET-CT (distant metastases). The relationship between tumour stage, FDG uptake, and texture with survival was assessed using Kaplan-Meier analysis.

RESULTS

Tumour heterogeneity correlated with SUV(max) and SUV(mean). The closest correlations were found for SUV(mean) measured as uniformity and entropy with coarse filtration (r=-0.754, p<0.0001; and r=0.748, p=0.0001 respectively). Heterogeneity was also significantly greater in patients with clinical stage III or IV for filter values between 1.0 and 2.0 (maximum difference at filter value 1.5: entropy: p=0.027; uniformity p=0.032). The median (range) survival was 21 (4-34) months. Tumour heterogeneity assessed by CTTA (coarse uniformity) was an independent predictor of survival [odds ratio (OR)=4.45 (95% CI: 1.08, 18.37); p=0.039].

CONCLUSION

CTTA assessment of tumour heterogeneity has the potential to identify oesophageal cancers with adverse biological features and provide a prognostic indicator of survival.

Feasibility of texture analysis for the assessment of biochemical changes in meniscal tissue on T1 maps calculated from delayed gadolinium-enhanced magnetic resonance imaging of cartilage data: comparison with conventional relaxation time measurements

OBJECTIVES

To (1) establish the feasibility of texture analysis for the in vivo assessment of biochemical changes in meniscal tissue on delayed gadolinium-enhanced magnetic resonance imaging of cartilage (dGEMRIC), and (2) compare textural with conventional T1 relaxation time measurements calculated from dGEMRIC data ("T1(Gd) relaxation times").

MATERIALS AND METHODS

We enrolled 10 asymptomatic volunteers (7 men and 3 women; mean age, 27.2 +/- 4.5 years), without a history of meniscus damage, in our study. MRI of the right knee was performed at 3.0 T. An isotropic, 3-dimensional (3D), double-echo steady-state sequences was used for morphologic evaluation, and a dual flip angle 3D gradient echo sequence was used for T1(Gd) mapping. All MRI scans were performed 90 minutes after injection of 0.2 mmol/kg of Gd-diethylenetriamine pentaacetic acid (DTPA), and subsequently, during application of a compressive force (50% of the body weight) in the axial direction. Regions of interest, covering the central portions of the posterior horn of the medial meniscus, were defined on 3 adjacent sagittal sections. Based on the relaxation time maps, mean T1(Gd), as well as the T1(Gd) texture features derived from the co-occurrence matrix (COC: Angular Second Moment, Entropy, Inverse Difference Moment) and wavelet transform (WAV: WavEnLL, WavEnHL, WavEnHH, WavEnLH), were calculated. Paired t tests were used to assess differences between baseline and compression, and intraclass correlation coefficients (ICC) were calculated to establish the intrarater reliability of the measurements.

RESULTS

Mean T1(Gd) (-67.3 ms, P = 0.011), Angular Second Moment (-0.0002, P = 0.009), Entropy (+0.033, P = 0.025), WavEnLL (+1011.16, P = 0.002), WavEnHL (+18.64, P = 0.012), and WavEnLH (+72.74, P = 0.035) differed significantly between baseline and compression. Intrarater reliability was substantial for mean T1(Gd) relaxation times (ICC = 0.99-1.0), and also for T1(Gd) co-occurrence matrix (ICC = 0.63-0.92) and WAV (ICC = 0.86-0.98) features.

CONCLUSIONS

Texture features extracted from T1 maps calculated from dGEMRIC data are feasible for the in vivo assessment of biochemical changes in the menisci, such as might be induced by mechanical loading. Thus, T1(Gd) texture features complement conventional relaxation time measurements. Further studies are necessary to determine whether the mechanical compression, or a prolonged Gd-DTPA uptake, or both, are responsible for the observed decrease in mean T1(Gd) relaxation times in the menisci.

Texture-based classification of focal liver lesions on MRI at 3.0 Tesla: a feasibility study in cysts and hemangiomas

PURPOSE

To determine the feasibility of texture analysis for the classification of liver cysts and hemangiomas, on nonenhanced, zero-fill interpolated T1- and T2-weighted MR images.

MATERIALS AND METHODS

Forty-five patients (26 women and 19 men; mean age, 58.1 +/- 16.9 years) with liver cysts or hemangiomas were enrolled in the study. After exclusion of images with artifacts, T1-weighted images of 42 patients, and T2-weighted images of 39 patients, obtained at 3.0 Tesla (T), were available for further analysis. Texture features derived from the gray-level histogram, co-occurrence and run-length matrix, gradient, autoregressive model, and wavelet transform were calculated. Fisher, probability of classification error and average correlation (POE+ACC), and mutual information coefficients were used to extract subsets of optimized texture features. Linear discriminant analysis (LDA) in combination with k nearest neighbor (k-NN) classification, and k-means clustering, were used for lesion classification.

RESULTS

LDA/k-NN produced misclassification rates of 16-18% on T1-weighted, and 12-18% on T2-weighted images. K-means clustering yielded misclassification rates of 15-23% on T1-weighted, and 15-25% on T2-weighted images.

CONCLUSION

Texture-based classification of liver cysts and hemangiomas is feasible on zero-fill interpolated MR images obtained at 3.0T. Further studies are warranted to investigate the value of texture-based classification of other liver lesions, such as hepatocellular and cholangiocellular carcinoma, on MRI.

Texture analysis of non-small cell lung cancer on unenhanced computed tomography: initial evidence for a relationship with tumour glucose metabolism and stage

The aim was to undertake an initial study of the relationship between texture features in computed tomography (CT) images of non-small cell lung cancer (NSCLC) and tumour glucose metabolism and stage. This retrospective pilot study comprised 17 patients with 18 pathologically confirmed NSCLC. Non-contrast-enhanced CT images of the primary pulmonary lesions underwent texture analysis in 2 stages as follows: (a) image filtration using Laplacian of Gaussian filter to differentially highlight fine to coarse textures, followed by (b) texture quantification using mean grey intensity (MGI), entropy (E) and uniformity (U) parameters. Texture parameters were compared with tumour fluorodeoxyglucose (FDG) uptake (standardised uptake value (SUV)) and stage as determined by the clinical report of the CT and FDG-positron emission tomography imaging. Tumour SUVs ranged between 2.8 and 10.4. The number of NSCLC with tumour stages I, II, III and IV were 4, 4, 4 and 6, respectively. Coarse texture features correlated with tumour SUV (E: r = 0.51, p = 0.03; U: r = −0.52, p = 0.03), whereas fine texture features correlated with tumour stage (MGI: r_s = 0.71, p = 0.001; E: r_s = 0.55, p = 0.02; U: r_s = −0.49, p = 0.04). Fine texture predicted tumour stage with a kappa of 0.7, demonstrating 100% sensitivity and 87.5% specificity for detecting tumours above stage II ( p = 0.0001). This study provides initial evidence for a relationship between texture features in NSCLC on non-contrast-enhanced CT and tumour metabolism and stage. Texture analysis warrants further investigation as a potential method for obtaining prognostic information for patients with NSCLC undergoing CT.