T2w‐MRI signal normalization affects radiomics features reproducibility

Repeatability of radiomic features from brain T1-W MRI after image intensity normalization: Implications for longitudinal studies on structural neurodegeneration

BACKGROUND AND OBJECTIVE

Radiomics extracts quantitative features from magnetic resonance images (MRI) and is especially useful in the presence of subtle pathological changes within human soft tissues. This scenario, however, may not cover the effects of intrinsic, e.g., aging-related, (physiological) neurodegeneration of normal brain tissue. The aim of the work was to study the repeatability of radiomic features extracted from an apparently normal area in longitudinally acquired T1-weighted MR brain images using three different intensity normalization approaches typically used in radiomics: Z-score, WhiteStripe and Nyul.

METHODS

Fifty-nine images of hearing impaired, yet cognitively intact, patients were repeatedly acquired in two different time points within six months. Ninety-one radiomic features were obtained from an area within the pons region, considered to be a healthy brain tissue according to previous analyses and reports. The Intraclass Correlation Coefficient (ICC) and the Concordance Correlation Coefficient (CCC) in the repeatability study were used as metrics.

RESULTS

Features extracted from the MRI normalized with Z-score showed results comparable in both ICC (0.90 (0.82-0.98)) and CCC (0.82 (0.69-0.95)) distribution values, in terms of median and quartiles, with those extracted from the images normalized with WhiteStripe (0.89 (0.84-0.92)) and (0.80 (0.73-0.84)), respectively.

CONCLUSION

Our findings underline the importance of, providing useful guidelines for, the intensity normalization of brain MRI prior to a longitudinal radiomic analysis.

Bronchial wall T2w MRI signal as a new imaging biomarker of severe asthma

OBJECTIVES

Severe asthma patients are prone to severe exacerbations with a need of hospital admission increasing the economic burden on healthcare systems. T2w lung MRI was found to be useful in the assessment of bronchial inflammation. The main goal of this study is to compare quantitative MRI T2 signal bronchial intensity between patients with severe and non-severe asthma.

METHODS

This is an ancillary study of a prospective single-center study (NCT03089346). We assessed the mean T2 intensity MRI signal of the bronchial wall area (BrWall_T2-MIS) in 15 severe and 15 age and sex-matched non-severe asthmatic patients. They also have had pulmonary function tests (PFTs), fractional exhaled nitric oxide (FeNO) and blood eosinophils count (Eos). Comparisons between the two groups were performed using Student's t-test. Correlations were assessed using Pearson coefficients. Reproducibility was assessed using intraclass correlation coefficient and Bland-Altman analysis.

RESULTS

BrWall_T2-MIS was higher in severe than in non-severe asthma patients (74 ± 12 vs 49 ± 14; respectively p < 0.001). BrWall_T2-MIS showed a moderate inverse correlation with PFTs in the whole cohort (r = -0.54, r = -0.44 for FEV1(%pred) and FEV1/FVC respectively, p ≤ 0.01) and in the severe asthma group (r = -0.53, r = -0.44 for FEV1(%pred) and FEV1/FVC respectively, p ≤ 0.01). Eos was moderately correlated with BrWall_T2-MIS in severe asthma group (r = 0.52, p = 0.047). Reproducibility was almost perfect with ICC = 0.99 and mean difference in Bland-Altman analysis of -0.15 [95% CI = -0.48-0.16].

CONCLUSION

Quantification of bronchial wall T2w signal intensity appears to be able to differentiate severe from non-severe asthma and correlates with obstructive PFTs' parameters and inflammatory markers in severe asthma.

CRITICAL RELEVANCE STATEMENT

The development of non-ionizing imaging biomarkers could play an essential role in the management of patients with severe asthma in the current era of biological therapies.

KEY POINTS

Severe asthma exhibits severe exacerbations with a high burden on healthcare systems. T2w bronchial wall signal intensity is related to inflammatory biomarker in severe asthma. T2w MRI may represent a non-invasive tool to follow up severe asthma patients.

Comparison of intensity normalization methods in prostate, brain, and breast cancer multi-parametric magnetic resonance imaging

Objectives

Intensity variation in multi-parametric magnetic resonance imaging (MP-MRI) is a confounding factor in MRI analyses. Previous studies have employed several normalization methods, but there is a lack of consensus on which method results in the most comparable images across vendors and acquisitions. This study used MP-MRI collected from patients with confirmed prostate, brain, or breast cancer to examine common intensity normalization methods to identify which best harmonizes intensity values across cofounds.

Materials and methods

Multiple normalization methods were deployed for intensity comparison between three unique sites, MR vendors, and magnetic field strength. Additionally, we calculated radiomic features before and after intensity normalization to determine how downstream analyses may be affected. Specifically, in the prostate cancer cohort, we tested these methods on T2-weighted imaging (T2WI) and additionally looked at a subset of patients who were scanned with and without the use of an endorectal coil (ERC). In a cohort of glioblastoma (GBM) patients, we tested these methods in T1 pre- and post-contrast enhancement (T1, T1C), fluid attenuated inversion recovery (FLAIR), and apparent diffusion coefficient (ADC) maps. Finally, in the breast cancer cohort, we tested methods on T1-weighted nonfat-suppressed images. All methods were compared using a two one-sided test (TOST) to test for equivalence of mean and standard deviation of intensity distributions.

Results

While each organ had unique results, across every tested comparison, using the Z-score of intensity within a mask of the organ consistently provided an equivalent distribution (all p < 0.001).

Conclusions

Our results suggest that intensity normalization using the Z-score of intensity within prostate, breast, and brain MR images produces the most comparable intensities between sites, MR vendors, magnetic field strength, and prostate endorectal coil usage. Likewise, Z-score normalization provided the highest percentage of radiomic features that were statistically equal across the three organs.

Investigation of radiomic features on MRI images to identify extraprostatic extension in prostate cancer

BACKGROUND AND OBJECTIVE

Detection of extraprostatic extension (EPE) preoperatively is of critical importance in the context of prostate cancer (PCa) management and outcomes. This study aimed to characterize the radiomic features of malignant prostate lesions based on multi-paramagnetic magnetic resonance imaging (mpMRI).

METHODS

We analyzed 20 patients who underwent mpMRI followed by radical prostatectomy. Two experienced radiologists manually segmented the 3D lesions using the T2-weighted (T2WI) and Apparent Diffusion Coefficient (ADC) imaging sequences. A total of 210 radiomic features were extracted from each lesion. We used the Recursive Feature Elimination with Cross-Validation to select key features. Using the selected radiomic features, we developed a Multilayer Perceptron (MLP) neural network to classify the EPE and non-EPE lesions. The pathology results were accepted as gold standard for EPE. We measured the performance of the classifier, calculating the area-under-curve (AUC), sensitivity, and specificity.

RESULTS

A total of 25 lesions were segmented, including 12 lesions with EPE and 13 lesions without EPE, based on the pathology reports. We selected 18 radiomic features (18/210). The MLP classifier using these features provided a good sensitivity (0.75), specificity (0.79), and AUC of 0.82, 95 % CL [0.59 - 0.96] in identifying the EPE lesions.

CONCLUSIONS

This pilot study presents 18 radiomic features derived from T2-weighted and ADC images and demonstrates their potential in the preoperative prediction of EPE in PCa using an MLP model.

Radiogenomics Pilot Study: Association Between Radiomics and Single Nucleotide Polymorphism-Based Microarray Copy Number Variation in Diagnosing Renal Oncocytoma and Chromophobe Renal Cell Carcinoma

RO and ChRCC are kidney tumours with overlapping characteristics, making differentiation between them challenging. The objective of this research is to create a radiogenomics map by correlating radiomic features to molecular phenotypes in ChRCC and RO, using resection as the gold standard. Fourteen patients (6 RO and 8 ChRCC) were included in the prospective study. A total of 1,875 radiomic features were extracted from CT scans, alongside 632 cytobands containing 16,303 genes from the genomic data. Feature selection algorithms applied to the radiomic features resulted in 13 key features. From the genomic data, 24 cytobands highly correlated with histology were selected and cross-correlated with the radiomic features. The analysis identified four radiomic features that were strongly associated with seven genomic features. These findings demonstrate the potential of integrating radiomic and genomic data to enhance the differential diagnosis of RO and ChRCC, paving the way for more precise and non-invasive diagnostic tools in clinical practice.

Image biomarkers and explainable AI: handcrafted features versus deep learned features

Feature extraction and selection from medical data are the basis of radiomics and image biomarker discovery for various architectures, including convolutional neural networks (CNNs). We herein describe the typical radiomics steps and the components of a CNN for both deep feature extraction and end-to-end approaches. We discuss the curse of dimensionality, along with dimensionality reduction techniques. Despite the outstanding performance of deep learning (DL) approaches, the use of handcrafted features instead of deep learned features needs to be considered for each specific study. Dataset size is a key factor: large-scale datasets with low sample diversity could lead to overfitting; limited sample sizes can provide unstable models. The dataset must be representative of all the "facets" of the clinical phenomenon/disease investigated. The access to high-performance computational resources from graphics processing units is another key factor, especially for the training phase of deep architectures. The advantages of multi-institutional federated/collaborative learning are described. When large language models are used, high stability is needed to avoid catastrophic forgetting in complex domain-specific tasks. We highlight that non-DL approaches provide model explainability superior to that provided by DL approaches. To implement explainability, the need for explainable AI arises, also through post hoc mechanisms. RELEVANCE STATEMENT: This work aims to provide the key concepts for processing the imaging features to extract reliable and robust image biomarkers. KEY POINTS: The key concepts for processing the imaging features to extract reliable and robust image biomarkers are provided. The main differences between radiomics and representation learning approaches are highlighted. The advantages and disadvantages of handcrafted versus learned features are given without losing sight of the clinical purpose of artificial intelligence models.

Artificial Intelligence for Response Assessment in Neuro Oncology (AI-RANO), part 2: recommendations for standardisation, validation, and good clinical practice

Technological advancements have enabled the extended investigation, development, and application of computational approaches in various domains, including health care. A burgeoning number of diagnostic, predictive, prognostic, and monitoring biomarkers are continuously being explored to improve clinical decision making in neuro-oncology. These advancements describe the increasing incorporation of artificial intelligence (AI) algorithms, including the use of radiomics. However, the broad applicability and clinical translation of AI are restricted by concerns about generalisability, reproducibility, scalability, and validation. This Policy Review intends to serve as the leading resource of recommendations for the standardisation and good clinical practice of AI approaches in health care, particularly in neuro-oncology. To this end, we investigate the repeatability, reproducibility, and stability of AI in response assessment in neuro-oncology in studies on factors affecting such computational approaches, and in publicly available open-source data and computational software tools facilitating these goals. The pathway for standardisation and validation of these approaches is discussed with the view of trustworthy AI enabling the next generation of clinical trials. We conclude with an outlook on the future of AI-enabled neuro-oncology.

A Comprehensive Evaluation of Radiomic Features in Normal Brain Magnetic Resonance Imaging: Investigating Robustness and Region Variations

Background

Despite extensive research on various brain diseases, a few studies have focused on radiomic feature distribution in healthy brain images. The present study applied a novel radiomic framework to investigate the robustness and baseline values of radiomic features in normal brain magnetic resonance imaging (MRIs) regions.

Materials and Methods

Analyses were performed on T1 and T2 images including 276 normal brains and 14 healthy volunteers were scanned with three scanners using the same protocols. The images were divided into 1024 three-dimensional nonoverlap patches with the same pixel size. Seven patches located in the thalamus, putamen, hippocampus and brain stem were selected as volume of interest (VOI). Eighty-five radiomic features were generated. To investigate the variation of features across VOIs, the analysis of variance was performed and coefficient of variation (COV) and intraclass correlation coefficient (ICC) were explored to examine the features repeatability.

Results

Thalamus (right and left) and hippocampus (left) resulted in more stable features (COV ≤ 6%) in T1 and T2 images, respectively. The inter-scanner ICC analysis demonstrated the features of T2 sequences represented more repeatable results and the brain stem and thalamus (both T1 and T2) showed particularly high repeatability (higher ICC values). Robust results (ICC ≥ 0.9) were identified for energy and range features of the first order class and several textures features across different brain regions.

Conclusion

Our results indicated the baselines of the repeatable texture features in healthy brain structural MRI highlighting inter-scanner stability. According to the findings, MRI sequencing and VOI location impact feature robustness and should be considered in brain radiomic studies.

Integrated framework for quantitative T2-weighted MRI analysis following prostate cancer radiotherapy

Purpose

The aim of this study is to develop a framework for quantitative analysis of longitudinal T2-weighted MRIs (T2w) following radiotherapy (RT) for prostate cancer.

Materials and methods

The developed methodology includes: (i) deformable image registration of longitudinal series to pre-RT T2w for automated detection of prostate, peripheral zone (PZ), and gross tumor volume (GTV); and (ii) T2w signal-intensity harmonization based on three reference tissues. The REgistration and HARMonization (REHARM) framework was applied on T2w acquired in a clinical trial consisting of two pre-RT and three post-RT MRI exams. Image registration was assessed by the DICE coefficient between automatic and manual contours, and intensity normalization via inter-patient histogram intersection. Longitudinal consistency was evaluated by the repeatability coefficient and Pearson correlation (r) between the two T2w exams before RT.

Results

T2w from 107 MRI exams (23 patients) were utilized. Following REHARM, the histogram intersections for prostate, PZ and GTV increased from median = 0.43/0.16/0.13 to 0.66/0.44/0.46. The repeatability in T2w intensity estimation was better for the automatic than the manual contours for all three regions of interest (r = 0.9, p < 0.0001, for GTV). The changes in the tissues' T2w values pre- and post-RT became significant, indicating the measurable quantitative signal related to radiation.

Conclusions

The developed methodology allows to automate longitudinal analysis reducing data acquisition-related variation and improving consistency. The quantitative characterization of RT-induced changes in T2w will lead to new understanding of radiation effects enabling prediction modeling of RT response.

Quantifying the reproducibility and longitudinal repeatability of radiomics features in magnetic resonance Image-Guide accelerator Imaging: A phantom study

OBJECTIVE

This study aimed to quantitatively evaluate the inter-platform reproducibility and longitudinal acquisition repeatability of MRI radiomics features in Fluid-Attenuated Inversion Recovery (FLAIR), T2-weighted (T2W), and T1-weighted (T1W) sequences on MR-Linac systems using an American College of Radiology (ACR) phantom.

MATERIALS AND METHODS

This study used two MR-Linac systems (A and B) in different cancer centers. The ACR phantom was scanned on system A daily for 30 consecutive days to evaluate longitudinal repeatability. Additionally, retest data were collected after repositioning the phantom. Inter-platform reproducibility was assessed by conducting scans under identical conditions using system B. Regions of interest were delineated on the T1W sequence from system A and mapped to other sequences via rigid registration. Intra-observer and inter-observer comparisons were conducted. Repeatability and reproducibility were assessed using the intraclass correlation coefficient (ICC) and coefficient of variation (CV). Robust radiomics features were identified based on ICC>0.9 and CV<10 %.

RESULTS

Analysis showed that a higher proportion of radiomics features derived from longitudinal FLAIR sequence (51.65 %) met robustness criteria compared to T2W (48.35 %) and T1W (43.96 %). Additionally, more inter-platform features from the FLAIR sequence (62.64 %) were robust compared to T2W (42.86 %) and T1W (39.56 %). Test-retest and intra-observer repeatability were excellent across all sequences, with a median ICC of 0.99 and CV<5%. However, inter-observer reproducibility was inferior, especially for the T1W sequence.

CONCLUSIONS

Different sequences show variations in repeatability and reproducibility. The FLAIR sequence demonstrated advantages in both longitudinal repeatability and inter-platform reproducibility. Caution is warranted when interpreting data, particularly in longitudinal or multiplatform radiomics studies.

Enhancing the Clinical Utility of Radiomics: Addressing the Challenges of Repeatability and Reproducibility in CT and MRI

Radiomics, which integrates the comprehensive characterization of imaging phenotypes with machine learning algorithms, is increasingly recognized for its potential in the diagnosis and prognosis of oncological conditions. However, the repeatability and reproducibility of radiomic features are critical challenges that hinder their widespread clinical adoption. This review aims to address the paucity of discussion regarding the factors that influence the reproducibility and repeatability of radiomic features and their subsequent impact on the application of radiomic models. We provide a synthesis of the literature on the repeatability and reproducibility of CT/MR-based radiomic features, examining sources of variation, the number of reproducible features, and the availability of individual feature repeatability indices. We differentiate sources of variation into random effects, which are challenging to control but can be quantified through simulation methods such as perturbation, and biases, which arise from scanner variability and inter-reader differences and can significantly affect the generalizability of radiomic model performance in diverse settings. Four suggestions for repeatability and reproducibility studies are suggested: (1) detailed reporting of variation sources, (2) transparent disclosure of calculation parameters, (3) careful selection of suitable reliability indices, and (4) comprehensive reporting of reliability metrics. This review underscores the importance of random effects in feature selection and harmonizing biases between development and clinical application settings to facilitate the successful translation of radiomic models from research to clinical practice.

Magnetic Resonance Image Radiomic Reproducibility: The Impact of Preprocessing on Extracted Features from Gross and High-Risk Clinical Tumor Volumes in Cervical Cancer Patients before Brachytherapy

Background

Radiomic feature reproducibility assessment is critical in radiomics-based image biomarker discovery. This study aims to evaluate the impact of preprocessing parameters on the reproducibility of magnetic resonance image (MRI) radiomic features extracted from gross tumor volume (GTV) and high-risk clinical tumor volume (HR-CTV) in cervical cancer (CC) patients.

Methods

This study included 99 patients with pathologically confirmed cervical cancer who underwent an MRI prior to receiving brachytherapy. The GTV and HR-CTV were delineated on T2-weighted MRI and inputted into 3D Slicer for radiomic analysis. Before feature extraction, all images were preprocessed to a combination of several parameters of Laplacian of Gaussian (1 and 2), resampling (0.5 and 1), and bin width (5, 10, 25, and 50). The reproducibility of radiomic features was analyzed using the intra-class correlation coefficient (ICC).

Results

Almost all shapes and first-order features had ICC values > 0.95. Most second-order texture features were not reproducible (ICC < 0.95) in GTV and HR-CTV. Furthermore, 20% of all neighboring gray-tone difference matrix texture features had ICC > 0.90 in both GTV and HR-CTV.

Conclusion

The results presented here showed that MRI radiomic features are vulnerable to changes in preprocessing, and this issue must be understood and applied before any clinical decision-making. Features with ICC > 0.90 were considered the most reproducible features. Shape and first-order radiomic features were the most reproducible features in both GTV and HR-CTV. Our results also showed that GTV and HR-CTV radiomic features had similar changes against preprocessing sets.

Impact of Preprocessing Parameters in Medical Imaging-Based Radiomic Studies: A Systematic Review

BACKGROUND

Lately, radiomic studies featuring the development of a signature to use in prediction models in diagnosis or prognosis outcomes have been increasingly published. While the results are shown to be promising, these studies still have many pitfalls and limitations. One of the main issues of these studies is that radiomic features depend on how the images are preprocessed before their computation. Since, in widely known and used software for radiomic features calculation, it is possible to set these preprocessing parameters before the calculation of the radiomic feature, there are ongoing studies assessing the stability and repeatability of radiomic features to find the most suitable preprocessing parameters for every used imaging modality.

MATERIALS AND METHODS

We performed a comprehensive literature search using four electronic databases: PubMed, Cochrane Library, Embase, and Scopus. Mesh terms and free text were modeled in search strategies for databases. The inclusion criteria were studies where preprocessing parameters' influence on feature values and model predictions was addressed. Records lacking information on image acquisition parameters were excluded, and any eligible studies with full-text versions were included in the review process, while conference proceedings and monographs were disregarded. We used the QUADAS-2 (Quality Assessment of Diagnostic Accuracy Studies 2) tool to investigate the risk of bias. We synthesized our data in a table divided by the imaging modalities subgroups.

RESULTS

After applying the inclusion and exclusion criteria, we selected 43 works. This review examines the impact of preprocessing parameters on the reproducibility and reliability of radiomic features extracted from multimodality imaging (CT, MRI, CBCT, and PET/CT). Standardized preprocessing is crucial for consistent radiomic feature extraction. Key preprocessing steps include voxel resampling, normalization, and discretization, which influence feature robustness and reproducibility. In total, 44% of the included works studied the effects of an isotropic voxel resampling, and most studies opted to employ a discretization strategy. From 2021, several studies started selecting the best set of preprocessing parameters based on models' best performance. As for comparison metrics, ICC was the most used in MRI studies in 58% of the screened works.

CONCLUSIONS

From our work, we highlighted the need to harmonize the use of preprocessing parameters and their values, especially in light of future studies of prospective studies, which are still lacking in the current literature.

Machine learning and radiomics analysis by computed tomography in colorectal liver metastases patients for RAS mutational status prediction

PURPOSE

To assess the efficacy of machine learning and radiomics analysis by computed tomography (CT) in presurgical setting, to predict RAS mutational status in colorectal liver metastases.

METHODS

Patient selection in a retrospective study was carried out from January 2018 to May 2021 considering the following inclusion criteria: patients subjected to surgical resection for liver metastases; proven pathological liver metastases; patients subjected to enhanced CT examination in the presurgical setting with a good quality of images; and RAS assessment as standard reference. A total of 851 radiomics features were extracted using the PyRadiomics Python package from the Slicer 3D image computing platform after slice-by-slice segmentation on CT portal phase by two expert radiologists of each individual liver metastasis performed first independently by the individual reader and then in consensus. Balancing technique was performed, and inter- and intraclass correlation coefficients were calculated to assess the between-observer and within-observer reproducibility of features. Receiver operating characteristics (ROC) analysis with the calculation of area under the ROC curve (AUC), sensitivity (SENS), specificity (SPEC), positive predictive value (PPV), negative predictive value (NPV) and accuracy (ACC) were assessed for each parameter. Linear and non-logistic regression model (LRM and NLRM) and different machine learning-based classifiers were considered. Moreover, features selection was performed before and after a normalized procedure using two different methods (3-sigma and z-score).

RESULTS

Seventy-seven liver metastases in 28 patients with a mean age of 60 years (range 40-80 years) were analyzed. The best predictors, at univariate analysis for both normalized procedures, were original_shape_Maximum2DDiameter and wavelet_HLL_glcm_InverseVariance that reached an accuracy of 80%, an AUC ≥ 0.75, a sensitivity ≥ 80% and a specificity ≥ 70% (p value <  < 0.01). However, a multivariate analysis significantly increased the accuracy in RAS prediction when a linear regression model (LRM) was used. The best performance was obtained using a LRM combining linearly 12 robust features after a z-score normalization procedure: AUC of 0.953, accuracy 98%, sensitivity 96%, specificity of 100%, PPV 100% and NPV 96% (p value <  < 0.01). No statistically significant increase was obtained considering the tested machine learning both without normalization and with normalization methods.

CONCLUSIONS

Normalized approach in CT radiomics analysis allows to predict RAS mutational status in colorectal liver metastases patients.

CT and MRI radiomic features of lung cancer (NSCLC): comparison and software consistency

BACKGROUND

Radiomics is a quantitative approach that allows the extraction of mineable data from medical images. Despite the growing clinical interest, radiomics studies are affected by variability stemming from analysis choices. We aimed to investigate the agreement between two open-source radiomics software for both contrast-enhanced computed tomography (CT) and contrast-enhanced magnetic resonance imaging (MRI) of lung cancers and to preliminarily evaluate the existence of radiomic features stable for both techniques.

METHODS

Contrast-enhanced CT and MRI images of 35 patients affected with non-small cell lung cancer (NSCLC) were manually segmented and preprocessed using three different methods. Sixty-six Image Biomarker Standardisation Initiative-compliant features common to the considered platforms, PyRadiomics and LIFEx, were extracted. The correlation among features with the same mathematical definition was analyzed by comparing PyRadiomics and LIFEx (at fixed imaging technique), and MRI with CT results (for the same software).

RESULTS

When assessing the agreement between LIFEx and PyRadiomics across the considered resampling, the maximum statistically significant correlations were observed to be 94% for CT features and 95% for MRI ones. When examining the correlation between features extracted from contrast-enhanced CT and MRI using the same software, higher significant correspondences were identified in 11% of features for both software.

CONCLUSIONS

Considering NSCLC, (i) for both imaging techniques, LIFEx and PyRadiomics agreed on average for 90% of features, with MRI being more affected by resampling and (ii) CT and MRI contained mostly non-redundant information, but there are shape features and, more importantly, texture features that can be singled out by both techniques.

RELEVANCE STATEMENT

Identifying and selecting features that are stable cross-modalities may be one of the strategies to pave the way for radiomics clinical translation.

KEY POINTS

• More than 90% of LIFEx and PyRadiomics features contain the same information. • Ten percent of features (shape, texture) are stable among contrast-enhanced CT and MRI. • Software compliance and cross-modalities stability features are impacted by the resampling method.

Prediction of High-Risk Neuroblastoma Among Neuroblastic Tumors Using Radiomics Features Derived from Magnetic Resonance Imaging: A Pilot Study

PURPOSE

This study aimed to predict high-risk neuroblastoma among neuroblastic tumors using radiomics features extracted from MRI.

MATERIALS AND METHODS

Pediatric patients (age≤18 years) diagnosed with neuroblastic tumors who had pre-treatment MR images available were enrolled from institution A from January 2010 to November 2019 (training set) and institution B from January 2016 to January 2022 (test set). Segmentation was performed with regions of interest manually drawn along tumor margins on the slice with the widest tumor area by two radiologists. First-order and texture features were extracted and intraclass correlation coefficients (ICCs) were calculated. Multivariate logistic regression (MLR) and random forest (RF) models from 10-fold cross-validation were built using these features. The trained MLR and RF models were tested in an external test set.

RESULTS

Thirty-two patients (M:F=23:9, 26.0±26.7 months) were in the training set and 14 patients (M:F=10:4, 33.4±20.4 months) were in the test set with radiomics features (n=930) being extracted. For 10 of the most relevant features selected, intra- and inter-observer variability was moderate to excellent (ICCs 0.633-0.911, 0.695-0.985, respectively). The area under the receiver operating characteristic curve (AUC) was 0.94 (sensitivity 67%, specificity 91%, and accuracy 84%) for the MLR model and the average AUC was 0.83 (sensitivity 44%, specificity 87%, and accuracy 75%) for the RF model from 10-fold cross-validation. In the test set, AUCs of the MLR and RF models were 0.94 and 0.91, respectively.

CONCLUSION

An MRI-based radiomics model can help predict high-risk neuroblastoma among neuroblastic tumors.

Contrast-enhanced CT-based radiomics differentiate anterior mediastinum lymphoma from thymoma without myasthenia gravis and calcification

AIM

To explore the value of a radiomics model based on enhanced computed tomography (CT) in differentiating anterior mediastinal lymphoma (AML) and thymoma without myasthenia gravis (MG) and calcification.

MATERIALS AND METHODS

The present study analysed patients who were diagnosed histologically with AML and thymoma in three independent institutions. All pre-treatment patients underwent enhanced CT. In the training group of patients from institutions 1 (the First Affiliated Hospital of Kunming Medical University) and 3 (the Yunnan Cancer Hospital), two radiologists independently analysed the enhanced CT images and performed manual segmentation of each tumour. Radiomics features were screened using interobserver interclass coefficient (ICC) analysis, feature correlation analysis, and L1 regularisation. The discriminative efficacy of the logistic regression model was evaluated using receiver operating characteristic (ROC) analysis. Validation group of patients from institution 2 (the Second Affiliated Hospital of Zhejiang University School of Medicine) was used to validate the proposed models.

RESULTS

A total of 114 patients were enrolled in this study and 1,743 radiomics features were extracted from the enhanced CT images. After feature screening, the remaining 37 robust radiomics features were used to construct the model. In the training group, the AUC of the model was 0.987 (95% confidence interval [CI]: 0.976-0.999), the sensitivity, specificity, and accuracy were 0.912, 0.946, and 0.924, respectively. In the validation group, the AUC of the model was 0.798 (95% CI: 0.683-0.913), the sensitivity, specificity, and accuracy were 0.760, 0.700, and 0.743, respectively.

CONCLUSION

The radiomics model created provided effective information to assist in the selection of clinical strategies, thus reducing unnecessary procedures in patients with AML and guiding direct surgery in patients with thymoma to avoid biopsy.

Achieving imaging and computational reproducibility on multiparametric MRI radiomics features in brain tumor diagnosis: phantom and clinical validation

OBJECTIVES

The Image Biomarker Standardization Initiative has helped improve the computational reproducibility of MRI radiomics features. Nonetheless, the MRI sequences and features with high imaging reproducibility are yet to be established. To determine reproducible multiparametric MRI radiomics features across test-retest, multi-scanner, and computational reproducibility comparisons, and to evaluate their clinical value in brain tumor diagnosis.

METHODS

To assess reproducibility, T1-weighted imaging (T1WI), T2-weighted imaging (T2WI), and diffusion-weighted imaging (DWI) were acquired from three 3-T MRI scanners using standardized phantom, and radiomics features were extracted using two computational algorithms. Reproducible radiomics features were selected when the concordance correlation coefficient value above 0.9 across multiple sessions, scanners, and computational algorithms. Random forest classifiers were trained with reproducible features (n = 117) and validated in a clinical cohort (n = 50) to evaluate whether features with high reproducibility improved the differentiation of glioblastoma from primary central nervous system lymphomas (PCNSLs).

RESULTS

Radiomics features from T2WI demonstrated higher repeatability (65-94%) than those from DWI (38-48%) or T1WI (2-92%). Across test-retest, multi-scanner, and computational comparisons, T2WI provided 41 reproducible features, DWI provided six, and T1WI provided two. The performance of the classification model with reproducible features was higher than that using non-reproducible features in both training set (AUC, 0.916 vs. 0.877) and validation set (AUC, 0.957 vs. 0.869).

CONCLUSION

Radiomics features with high reproducibility across multiple sessions, scanners, and computational algorithms were identified, and they showed higher diagnostic performance than non-reproducible radiomics features in the differentiation of glioblastoma from PCNSL.

CLINICAL RELEVANCE STATEMENT

By identifying the radiomics features showing higher multi-machine reproducibility, our results also demonstrated higher radiomics diagnostic performance in the differentiation of glioblastoma from PCNSL, paving the way for further research designs and clinical application in neuro-oncology.

KEY POINTS

• Highly reproducible radiomics features across multiple sessions, scanners, and computational algorithms were identified using phantom and applied to clinical diagnosis. • Radiomics features from T2-weighted imaging were more reproducible than those from T1-weighted and diffusion-weighted imaging. • Radiomics features with good reproducibility had better diagnostic performance for brain tumors than features with poor reproducibility.

Uncovering prostate cancer aggressiveness signal in T2-weighted MRI through a three-reference tissues normalization technique

Quantitative T2-weighted MRI (T2W) interpretation is impeded by the variability of acquisition-related features, such as field strength, coil type, signal amplification, and pulse sequence parameters. The main purpose of this work is to develop an automated method for prostate T2W intensity normalization. The procedure includes the following: (i) a deep learning-based network utilizing MASK R-CNN for automatic segmentation of three reference tissues: gluteus maximus muscle, femur, and bladder; (ii) fitting a spline function between average intensities in these structures and reference values; and (iii) using the function to transform all T2W intensities. The T2W distributions in the prostate cancer regions of interest (ROIs) and normal appearing prostate tissue (NAT) were compared before and after normalization using Student's t-test. The ROIs' T2W associations with the Gleason Score (GS), Decipher genomic score, and a three-tier prostate cancer risk were evaluated with Spearman's correlation coefficient (rS ). T2W differences in indolent and aggressive prostate cancer lesions were also assessed. The MASK R-CNN was trained with manual contours from 32 patients. The normalization procedure was applied to an independent MRI dataset from 83 patients. T2W differences between ROIs and NAT significantly increased after normalization. T2W intensities in 231 biopsy ROIs were significantly negatively correlated with GS (rS  = -0.21, p = 0.001), Decipher (rS  = -0.193, p = 0.003), and three-tier risk (rS  = -0.235, p < 0.001). The average T2W intensities in the aggressive ROIs were significantly lower than in the indolent ROIs after normalization. In conclusion, the automated triple-reference tissue normalization method significantly improved the discrimination between prostate cancer and normal prostate tissue. In addition, the normalized T2W intensities of cancer exhibited a significant association with tumor aggressiveness. By improving the quantitative utilization of the T2W in the assessment of prostate cancer on MRI, the new normalization method represents an important advance over clinical protocols that do not include sequences for the measurement of T2 relaxation times.

Machine learning-based radiomics analysis in predicting RAS mutational status using magnetic resonance imaging

PURPOSE

To assess the efficacy of radiomics features, obtained by magnetic resonance imaging (MRI) with hepatospecific contrast agent, in pre-surgical setting, to predict RAS mutational status in liver metastases.

METHODS

Patients with MRI in pre-surgical setting were enrolled in a retrospective study. Manual segmentation was made by means 3D Slicer image computing, and 851 radiomics features were extracted as median values using the PyRadiomics Python package. The features were extracted considering the agreement with the Imaging Biomarker Standardization Initiative (IBSI). Balancing was performed through synthesis of samples for the underrepresented classes using the self-adaptive synthetic oversampling (SASYNO) approach. Inter- and intraclass correlation coefficients (ICC) were calculated to assess the between-observer and within-observer reproducibility of all radiomics characteristics. For continuous variables, nonparametric Wilcoxon-Mann-Whitney test was utilized. Benjamini and Hochberg's false discovery rate (FDR) adjustment for multiple testing was used. Receiver operating characteristics (ROC) analysis with the calculation of area under the ROC curve (AUC), sensitivity (SENS), specificity (SPEC), positive predictive value (PPV), negative predictive value (NPV) and accuracy (ACC) were assessed for each parameter. Linear and non-logistic regression model (LRM and NLRM) and different machine learning-based classifiers including decision tree (DT), k-nearest neighbor (KNN) and support vector machine (SVM) were considered. Moreover, features selection were performed before and after a normalized procedure using two different methods (3-sigma and z-score). McNemar test was used to assess differences statistically significant between dichotomic tables. All statistical procedures were done using MATLAB R2021b Statistics and Machine Toolbox (MathWorks, Natick, MA, USA).

RESULTS

Seven normalized radiomics features, extracted from arterial phase, 11 normalized radiomics features, from portal phase, 12 normalized radiomics features from hepatobiliary phase and 12 normalized features from T2-W SPACE sequence were robust predictors of RAS mutational status. The multivariate analysis increased significantly the accuracy in RAS prediction when a LRM was used, combining 12 robust normalized features extracted by VIBE hepatobiliary phase reaching an accuracy of 99%, a sensitivity 97%, a specificity of 100%, a PPV of 100% and a NPV of 98%. No statistically significant increase was obtained, considering the tested classifiers DT, KNN and SVM, both without normalization and with normalization methods.

CONCLUSIONS

Normalized approach in MRI radiomics analysis allows to predict RAS mutational status.

MRI-based radiomic features for identifying recurrent prostate cancer after proton radiation therapy

PURPOSE

Magnetic Resonance Imaging (MRI) evaluation of recurrent prostate cancer (PCa) following proton beam therapy is challenging due to radiation-induced tissue changes. This study aimed to evaluate MRI-based radiomic features so as to identify the recurrent PCa after proton therapy.

METHODS

We retrospectively studied 12 patients with biochemical recurrence (BCR) following proton therapy. Two experienced radiologists identified prostate lesions from multi-parametric MRI (mpMRI) images post-proton therapy and marked control regions of interest (ROIs) on the contralateral side of the prostate gland. A total of 210 radiomic features were extracted from lesions and control regions on the T2-weighted (T2WI) and Apparent Diffusion Coefficient (ADC) image series. Recursive Feature Elimination with Cross-Validation method (RFE-CV) was used for feature selection. A Multilayer Perceptron (MLP) neural network was developed to classify three classes: cancerous, benign, and healthy tissue. The 12-core biopsy results were used as the gold standard for the segmentations. The classifier performance was measured using specificity, sensitivity, the area under receiver operating characteristic curve (AUC), and other statistical indicators.

RESULTS

Based on biopsy results, 10 lesions were identified as PCa recurrence while eight lesions were confirmed to be benign. Ten radiomic features (10/210) were selected to build the multi-class classifier. The radiomics classifier gave an accuracy of 0.83 in identifying cancerous, benign, and healthy tissue with a sensitivity of 0.80 and specificity of 0.85. The model yielded an AUC of 0.87, 95% CI [0.72-1.00] in differentiating cancer from the benign and healthy tissues.

CONCLUSIONS

Our proof-of-concept study demonstrates the potential of using radiomic features as part of the differential diagnosis of PCa on mpMRI following proton therapy. The results need to be validated in a larger cohort.

A Physics-Based Algorithm to Universally Standardize Routinely Obtained Clinical T2-Weighted Images

RATIONALE AND OBJECTIVES

MR images can be challenging for machine learning and other large-scale analyses because most clinical images, for example, T2-weighted (T2w) images, reflect not only the biologically relevant T2 of tissue but also hardware and acquisition parameters that vary from site to site. Quantitative T2 mapping avoids these confounds because it quantitatively isolates the biological parameter of interest, thus representing a universal standardization across sites. However, efforts to incorporate quantitative mapping sequences into routine clinical practice have seen slow adoption. Here we show, for the first time, that the routine T2w complex raw dataset can be successfully regarded as a quantitative mapping sequence that can be reconstructed with classical optimization methods and physics-based constraints.

MATERIALS AND METHODS

While previous constrained reconstruction methods are unable to reconstruct a T2 map based on this data, the expanding-constrained alternating minimization for parameter mapping (e-CAMP), which employs stepwise initialization, a linearized version of the exponential model and a phase conjugacy constraint, is demonstrated to provide useful quantitative maps directly from a vendor T2w single image data.

RESULTS

This paper introduces the method and demonstrates its performance using simulations, retrospectively undersampled brain images, and prospectively acquired T2w images taken on both phantom and brain.

CONCLUSION

Because T2w scans are included in nearly every protocol, this approach could open the door to creating large, standardized datasets without requiring widespread changes in clinical protocols.

A Critical Analysis of the Robustness of Radiomics to Variations in Segmentation Methods in 18F-PSMA-1007 PET Images of Patients Affected by Prostate Cancer

BACKGROUND

Radiomics shows promising results in supporting the clinical decision process, and much effort has been put into its standardization, thus leading to the Imaging Biomarker Standardization Initiative (IBSI), that established how radiomics features should be computed. However, radiomics still lacks standardization and many factors, such as segmentation methods, limit study reproducibility and robustness.

AIM

We investigated the impact that three different segmentation methods (manual, thresholding and region growing) have on radiomics features extracted from 18F-PSMA-1007 Positron Emission Tomography (PET) images of 78 patients (43 Low Risk, 35 High Risk). Segmentation was repeated for each patient, thus leading to three datasets of segmentations. Then, feature extraction was performed for each dataset, and 1781 features (107 original, 930 Laplacian of Gaussian (LoG) features, 744 wavelet features) were extracted. Feature robustness and reproducibility were assessed through the intra class correlation coefficient (ICC) to measure agreement between the three segmentation methods. To assess the impact that the three methods had on machine learning models, feature selection was performed through a hybrid descriptive-inferential method, and selected features were given as input to three classifiers, K-Nearest Neighbors (KNN), Support Vector Machines (SVM), Linear Discriminant Analysis (LDA), Random Forest (RF), AdaBoost and Neural Networks (NN), whose performance in discriminating between low-risk and high-risk patients have been validated through 30 times repeated five-fold cross validation.

CONCLUSIONS

Our study showed that segmentation methods influence radiomics features and that Shape features were the least reproducible (average ICC: 0.27), while GLCM features the most reproducible. Moreover, feature reproducibility changed depending on segmentation type, resulting in 51.18% of LoG features exhibiting excellent reproducibility (range average ICC: 0.68-0.87) and 47.85% of wavelet features exhibiting poor reproducibility that varied between wavelet sub-bands (range average ICC: 0.34-0.80) and resulted in the LLL band showing the highest average ICC (0.80). Finally, model performance showed that region growing led to the highest accuracy (74.49%), improved sensitivity (84.38%) and AUC (79.20%) in contrast with manual segmentation.

Clinico-biological-radiomics (CBR) based machine learning for improving the diagnostic accuracy of FDG-PET false-positive lymph nodes in lung cancer

BACKGROUND

The main problem of positron emission tomography/computed tomography (PET/CT) for lymph node (LN) staging is the high false positive rate (FPR). Thus, we aimed to explore a clinico-biological-radiomics (CBR) model via machine learning (ML) to reduce FPR and improve the accuracy for predicting the hypermetabolic mediastinal-hilar LNs status in lung cancer than conventional PET/CT.

METHODS

A total of 260 lung cancer patients with hypermetabolic mediastinal-hilar LNs (SUVmax ≥ 2.5) were retrospectively reviewed. Patients were treated with surgery with systematic LN resection and pathologically divided into the LN negative (LN-) and positive (LN +) groups, and randomly assigned into the training (n = 182) and test (n = 78) sets. Preoperative CBR dataset containing 1738 multi-scale features was constructed for all patients. Prediction models for hypermetabolic LNs status were developed using the features selected by the supervised ML algorithms, and evaluated using the classical diagnostic indicators. Then, a nomogram was developed based on the model with the highest area under the curve (AUC) and the lowest FPR, and validated by the calibration plots.

RESULTS

In total, 109 LN- and 151 LN + patients were enrolled in this study. 6 independent prediction models were developed to differentiate LN- from LN + patients using the selected features from clinico-biological-image dataset, radiomics dataset, and their combined CBR dataset, respectively. The DeLong test showed that the CBR Model containing all-scale features held the highest predictive efficiency and the lowest FPR among all of established models (p < 0.05) in both the training and test sets (AUCs of 0.90 and 0.89, FPRs of 12.82% and 6.45%, respectively) (p < 0.05). The quantitative nomogram based on CBR Model was validated to have a good consistency with actual observations.

CONCLUSION

This study presents an integrated CBR nomogram that can further reduce the FPR and improve the accuracy of hypermetabolic mediastinal-hilar LNs evaluation than conventional PET/CT in lung cancer, thereby greatly reducing the risk of overestimation and assisting for precision treatment.

Magnetic resonance imaging preprocessing and radiomic features for classification of autosomal dominant polycystic kidney disease genotype

Purpose

Our study aims to investigate the impact of preprocessing on magnetic resonance imaging (MRI) radiomic features extracted from the noncystic kidney parenchyma of patients with autosomal dominant polycystic kidney disease (ADPKD) in the task of classifying PKD1 versus PKD2 genotypes, which differ with regard to cyst burden and disease outcome.

Approach

The effect of preprocessing on radiomic features was investigated using a single T2-weighted fat saturated (T2W-FS) MRI scan from PKD1 and PKD2 subjects (29 kidneys in total) from the Consortium for Radiologic Imaging Studies of Polycystic Kidney Disease study. Radiomic feature reproducibility using the intraclass correlation coefficient (ICC) was computed across MRI normalizations (z -score, reference-tissue, and original image), gray-level discretization, and upsampling and downsampling pixel schemes. A second dataset for genotype classification from 136 subjects T2W-FS MRI images previously enrolled in the HALT Progression of Polycystic Kidney Disease study was matched for age, gender, and Mayo imaging classification class. Genotype classification was performed using a logistic regression classifier and radiomic features extracted from (1) the noncystic kidney parenchyma and (2) the entire kidney. The area under the receiver operating characteristic curve (AUC) was used to evaluate the classification performance across preprocessing methods.

Results

Radiomic features extracted from the noncystic kidney parenchyma were sensitive to preprocessing parameters, with varying reproducibility depending on the parameter. The percentage of features with good-to-excellent ICC scores ranged from 14% to 58%. AUC values ranged between 0.47 to 0.68 and 0.56 to 0.73 for the noncystic kidney parenchyma and entire kidney, respectively.

Conclusions

Reproducibility of radiomic features extracted from the noncystic kidney parenchyma was dependent on the preprocessing parameters used, and the effect on genotype classification was sensitive to preprocessing parameters. The results suggest that texture features may be indicative of genotype expression in ADPKD.

Biological and clinical significance of radiomics features obtained from magnetic resonance imaging preceding pre-carbon ion radiotherapy in prostate cancer based on radiometabolomics

Introduction

We aimed to investigate the feasibility of metabolomics to explain the underlying biological implications of radiomics features obtained from magnetic resonance imaging (MRI) preceding carbon ion radiotherapy (CIRT) in patients with prostate cancer and to further explore the clinical significance of radiomics features on the prognosis of patients, based on their biochemical recurrence (BCR) status.

Methods

Metabolomic results obtained using high-performance liquid chromatography coupled with tandem mass spectrometry of urine samples, combined with pre-RT radiomic features extracted from MRI images, were evaluated to investigate their biological significance. Receiver operating characteristic (ROC) curve analysis was subsequently conducted to examine the correlation between these biological implications and clinical BCR status. Statistical and metabolic pathway analyses were performed using MetaboAnalyst and R software.

Results

Correlation analysis revealed that methionine alteration extent was significantly related to four radiomic features (Contrast, Difference Variance, Small Dependence High Gray Level Emphasis, and Mean Absolute Deviation), which were significantly correlated with BCR status. The area under the curve (AUC) for BCR prediction of these four radiomic features ranged from 0.704 to 0.769, suggesting that the higher the value of these four radiomic features, the greater the decrease in methionine levels after CIRT and the lower the probability of BCR. Pre-CIRT MRI radiomic features were associated with CIRT-suppressed metabolites.

Discussion

These radiomic features can be used to predict the alteration in the amplitude of methionine after CIRT and the BCR status, which may contribute to the optimization of the CIRT strategy and deepen the understanding of PCa.

Towards the Definition of Radiomic Features and Clinical Indices to Enhance the Diagnosis of Clinically Significant Cancers in PI-RADS 4 and 5 Lesions

Prostate cancer (PC) is the most frequently diagnosed cancer among adult men, and its incidence is increasing worldwide [...].

An image-based modeling framework for predicting spatiotemporal brain cancer biology within individual patients

Imaging is central to the clinical surveillance of brain tumors yet it provides limited insight into a tumor's underlying biology. Machine learning and other mathematical modeling approaches can leverage paired magnetic resonance images and image-localized tissue samples to predict almost any characteristic of a tumor. Image-based modeling takes advantage of the spatial resolution of routine clinical scans and can be applied to measure biological differences within a tumor, changes over time, as well as the variance between patients. This approach is non-invasive and circumvents the intrinsic challenges of inter- and intratumoral heterogeneity that have historically hindered the complete assessment of tumor biology and treatment responsiveness. It can also reveal tumor characteristics that may guide both surgical and medical decision-making in real-time. Here we describe a general framework for the acquisition of image-localized biopsies and the construction of spatiotemporal radiomics models, as well as case examples of how this approach may be used to address clinically relevant questions.

Characterization of high-grade prostate cancer at multiparametric MRI using a radiomic-based computer-aided diagnosis system as standalone and second reader

PURPOSE

The purpose of this study was to develop and test across various scanners a zone-specific region-of-interest (ROI)-based computer-aided diagnosis system (CAD) aimed at characterizing, on MRI, International Society of Urological Pathology (ISUP) grade≥2 prostate cancers.

MATERIALS AND METHODS

ROI-based quantitative models were selected in multi-vendor training (265 pre-prostatectomy MRIs) and pre-test (112 pre-biopsy MRIs) datasets. The best peripheral and transition zone models were combined and retrospectively assessed in internal (158 pre-biopsy MRIs) and external (104 pre-biopsy MRIs) test datasets. Two radiologists (R1/R2) retrospectively delineated the lesions targeted at biopsy in test datasets. The CAD area under the receiver operating characteristic curve (AUC) for characterizing ISUP≥2 cancers was compared to that of the Prostate Imaging-Reporting and Data System version2 (PI-RADSv2) score prospectively assigned to targeted lesions.

RESULTS

The best models used the 25th apparent diffusion coefficient (ADC) percentile in transition zone and the 2nd ADC percentile and normalized wash-in rate in peripheral zone. The PI-RADSv2 AUCs were 82% (95% confidence interval [CI]: 74-87) and 86% (95% CI: 81-91) in the internal and external test datasets respectively. They were not different from the CAD AUCs obtained with R1 and R2 delineations, in the internal (82% [95% CI: 76-89], P = 0.95 and 85% [95% CI: 78-91], P = 0.55) and external (82% [95% CI: 74-91], P = 0.41 and 86% [95% CI:78-95], P = 0.98) test datasets. The CAD yielded sensitivities of 86-89% and 90-91%, and specificities of 64-65% and 69-75% in the internal and external test datasets respectively.

CONCLUSION

The CAD performance for characterizing ISUP grade≥2 prostate cancers on MRI is not different from that of PI-RADSv2 score across two test datasets.

Bayesian statistical modeling to predict observer-specific optimal windowing parameters in magnetic resonance imaging

Magnetic resonance (MR) images require a process known as windowing for optimizing the display conditions. However, the conventional windowing process often fails to achieve the preferred display conditions for observers due to various factors. This study proposes a novel framework for predicting the preferred windowing parameters for each observer using Bayesian statistical modeling. MR images obtained from 1000 patients were divided into training and test sets at a 7:3 ratio. The image intensity and windowing parameters were standardized using previously reported methods. Bayesian statistical modeling was utilized to predict the windowing parameters preferred by three MR imaging (MRI) operators. The performance of the proposed framework was evaluated by assessing the mean relative error (MRE), mean absolute error (MAE), and Pearson's correlation coefficient (ρ) of the test set. In addition, the naive method, which presumes that the average value of the windowing parameters for each acquisition sequence and body region in the training set is optimal, was also used for comparison. Three MRI operators and three radiologists conducted visual assessments. The mean MRE, MAE, and ρ values for the window level and width (WL/WW) in the proposed framework were 12.6 and 13.9, 42.9 and 85.4, and 0.98 and 0.98, respectively. These results outperformed those obtained using the naive method. The visual assessments revealed no significant differences between the original and predicted display conditions, indicating that the proposed framework accurately predicts individualized windowing parameters with the additional advantages of robustness and ease of use. Thus, the proposed framework can effectively predict the windowing parameters preferred by each observer.

Associations between Vertebral Localized Contrast Changes and Adjacent Annular Fissures in Patients with Low Back Pain: A Radiomics Approach

Low back pain (LBP) is multifactorial and associated with various spinal tissue changes, including intervertebral disc fissures, vertebral pathology, and damaged endplates. However, current radiological markers lack specificity and individualized diagnostic capability, and the interactions between the various markers are not fully clear. Radiomics, a data-driven analysis of radiological images, offers a promising approach to improve evaluation and deepen the understanding of spinal changes related to LBP. This study investigated possible associations between vertebral changes and annular fissures using radiomics. A dataset of 61 LBP patients who underwent conventional magnetic resonance imaging followed by discography was analyzed. Radiomics features were extracted from segmented vertebrae and carefully reduced to identify the most relevant features associated with annular fissures. The results revealed three important texture features that display concentrated high-intensity gray levels, extensive regions with elevated gray levels, and localized areas with reduced gray levels within the vertebrae. These features highlight patterns within vertebrae that conventional classification systems cannot reflect on distinguishing between vertebrae adjacent to an intervertebral disc with or without an annular fissure. As such, the present study reveals associations that contribute to the understanding of pathophysiology and may provide improved diagnostics of LBP.

Stability of Multi-Parametric Prostate MRI Radiomic Features to Variations in Segmentation

Stability analysis remains a fundamental step in developing a successful imaging biomarker to personalize oncological strategies. This study proposes an in silico contour generation method for simulating segmentation variations to identify stable radiomic features. Ground-truth annotation provided for the whole prostate gland on the multi-parametric MRI sequences (T2w, ADC, and SUB-DCE) were perturbed to mimic segmentation differences observed among human annotators. In total, we generated 15 synthetic contours for a given image-segmentation pair. One thousand two hundred twenty-four unfiltered/filtered radiomic features were extracted applying Pyradiomics, followed by stability assessment using ICC(1,1). Stable features identified in the internal population were then compared with an external population to discover and report robust features. Finally, we also investigated the impact of a wide range of filtering strategies on the stability of features. The percentage of unfiltered (filtered) features that remained robust subjected to segmentation variations were T2w-36% (81%), ADC-36% (94%), and SUB-43% (93%). Our findings suggest that segmentation variations can significantly impact radiomic feature stability but can be mitigated by including pre-filtering strategies as part of the feature extraction pipeline.

RESOLVE-based radiomics in cervical cancer: improved image quality means better feature reproducibility?

AIM

To compare the reproducibility of apparent diffusion coefficient (ADC)-based radiomic features between readout-segmented echo-planar diffusion-weighted imaging (RESOLVE) and single-shot echo-planar diffusion-weighted imaging (SS-EPI DWI) in cervical cancer.

MATERIALS AND METHODS

The RESOLVE and SS-EPI DWI images of 36 patients with histopathologically confirmed cervical cancer were collected retrospectively. Two observers independently delineated the whole tumour on RESOLVE and SS-EPI DWI, and then copied them to the corresponding ADC maps. Shape, first-order, and texture features were extracted from ADC maps in the original and filtered (Laplacian of Gaussian [LoG] and wavelet) images. Thereafter, 1,316 features were generated in each RESOLVE and SS-EPI DWI, respectively. The reproducibility of radiomic features was assessed using intraclass correlation coefficient (ICC).

RESULTS

In the original images, RESOLVE showed 92.86%, 66.67%, and 86.67% of features with excellent reproducibility in shape, first-order, and texture features, while SS-EPI DWI showed 85.71%, 72.22%, and 60% of features with excellent reproducibility, respectively. In the LoG and wavelet filtered images, RESOLVE had 56.77% and 65.32% of features with excellent reproducibility and SS-EPI DWI had 44.95% and 61.96% of features with excellent reproducibility, respectively.

CONCLUSION

Compared with SS-EPI DWI, the feature reproducibility of RESOLVE was better in cervical cancer, especially for texture features. The filtered images cannot improve the feature reproducibility compared with the original images for both SS-EPI DWI and RESOLVE.

Radiomics-Based Machine Learning Model for Predicting Overall and Progression-Free Survival in Rare Cancer: A Case Study for Primary CNS Lymphoma Patients

Primary Central Nervous System Lymphoma (PCNSL) is an aggressive neoplasm with a poor prognosis. Although therapeutic progresses have significantly improved Overall Survival (OS), a number of patients do not respond to HD–MTX-based chemotherapy (15–25%) or experience relapse (25–50%) after an initial response. The reasons underlying this poor response to therapy are unknown. Thus, there is an urgent need to develop improved predictive models for PCNSL. In this study, we investigated whether radiomics features can improve outcome prediction in patients with PCNSL. A total of 80 patients diagnosed with PCNSL were enrolled. A patient sub-group, with complete Magnetic Resonance Imaging (MRI) series, were selected for the stratification analysis. Following radiomics feature extraction and selection, different Machine Learning (ML) models were tested for OS and Progression-free Survival (PFS) prediction. To assess the stability of the selected features, images from 23 patients scanned at three different time points were used to compute the Interclass Correlation Coefficient (ICC) and to evaluate the reproducibility of each feature for both original and normalized images. Features extracted from Z-score normalized images were significantly more stable than those extracted from non-normalized images with an improvement of about 38% on average (p-value < 10−12). The area under the ROC curve (AUC) showed that radiomics-based prediction overcame prediction based on current clinical prognostic factors with an improvement of 23% for OS and 50% for PFS, respectively. These results indicate that radiomics features extracted from normalized MR images can improve prognosis stratification of PCNSL patients and pave the way for further study on its potential role to drive treatment choice.

Impact of Wavelet Kernels on Predictive Capability of Radiomic Features: A Case Study on COVID-19 Chest X-ray Images

Radiomic analysis allows for the detection of imaging biomarkers supporting decision-making processes in clinical environments, from diagnosis to prognosis. Frequently, the original set of radiomic features is augmented by considering high-level features, such as wavelet transforms. However, several wavelets families (so called kernels) are able to generate different multi-resolution representations of the original image, and which of them produces more salient images is not yet clear. In this study, an in-depth analysis is performed by comparing different wavelet kernels and by evaluating their impact on predictive capabilities of radiomic models. A dataset composed of 1589 chest X-ray images was used for COVID-19 prognosis prediction as a case study. Random forest, support vector machine, and XGBoost were trained (on a subset of 1103 images) after a rigorous feature selection strategy to build-up the predictive models. Next, to evaluate the models generalization capability on unseen data, a test phase was performed (on a subset of 486 images). The experimental findings showed that Bior1.5, Coif1, Haar, and Sym2 kernels guarantee better and similar performance for all three machine learning models considered. Support vector machine and random forest showed comparable performance, and they were better than XGBoost. Additionally, random forest proved to be the most stable model, ensuring an appropriate balance between sensitivity and specificity.

Robustness of Radiomics in Pre-Surgical Computer Tomography of Non-Small-Cell Lung Cancer

Background: Radiomic features are increasingly used in CT of NSCLC. However, their robustness with respect to segmentation variability has not yet been demonstrated. The aim of this study was to assess radiomic features agreement across three kinds of segmentation. Methods: We retrospectively included 48 patients suffering from NSCLC who underwent pre-surgery CT. Two expert radiologists in consensus manually delineated three 3D-ROIs on each patient. To assess robustness for each feature, the intra-class correlation coefficient (ICC) across segmentations was evaluated. The ‘sensitivity’ of ICC upon some parameters affecting features computation (such as bin-width for first-order features and pixel-distances for second-order features) was also evaluated. Moreover, an assessment with respect to interpolator and isotropic resolution was also performed. Results: Our results indicate that ‘shape’ features tend to have excellent agreement (ICC > 0.9) across segmentations; moreover, they have approximately zero sensitivity to other parameters. ‘First-order’ features are in general sensitive to parameters variation; however, a few of them showed excellent agreement and low sensitivity (below 0.1) with respect to bin-width and pixel-distance. Similarly, a few second-order features showed excellent agreement and low sensitivity. Conclusions: Our results suggest that a limited number of radiomic features can achieve a high level of reproducibility in CT of NSCLC.

A Framework of Analysis to Facilitate the Harmonization of Multicenter Radiomic Features in Prostate Cancer

Pooling radiomic features coming from different centers in a statistical framework is challenging due to the variability in scanner models, acquisition protocols, and reconstruction settings. To remove technical variability, commonly called batch effects, different statistical harmonization strategies have been widely used in genomics but less considered in radiomics. The aim of this work was to develop a framework of analysis to facilitate the harmonization of multicenter radiomic features extracted from prostate T2-weighted magnetic resonance imaging (MRI) and to improve the power of radiomics for prostate cancer (PCa) management in order to develop robust non-invasive biomarkers translating into clinical practice. To remove technical variability and correct for batch effects, we investigated four different statistical methods (ComBat, SVA, Arsynseq, and mixed effect). The proposed approaches were evaluated using a dataset of 210 prostate cancer (PCa) patients from two centers. The impacts of the different statistical approaches were evaluated by principal component analysis and classification methods (LogitBoost, random forest, K-nearest neighbors, and decision tree). The ComBat method outperformed all other methods by achieving 70% accuracy and 78% AUC with the random forest method to automatically classify patients affected by PCa. The proposed statistical framework enabled us to define and develop a standardized pipeline of analysis to harmonize multicenter T2W radiomic features, yielding great promise to support PCa clinical practice.

Effect of magnetic resonance imaging pre-processing on the performance of model-based prostate tumor probability mapping

Objective. Multi-parametric magnetic resonance imaging (mpMRI) has become an important tool for the detection of prostate cancer in the past two decades. Despite the high sensitivity of MRI for tissue characterization, it often suffers from a lack of specificity. Several well-established pre-processing tools are publicly available for improving image quality and removing both intra- and inter-patient variability in order to increase the diagnostic accuracy of MRI. To date, most of these pre-processing tools have largely been assessed individually. In this study we present a systematic evaluation of a multi-step mpMRI pre-processing pipeline to automate tumor localization within the prostate using a previously trained model.Approach. The study was conducted on 31 treatment-naïve prostate cancer patients with a PI-RADS-v2 compliant mpMRI examination. Multiple methods were compared for each pre-processing step: (1) bias field correction, (2) normalization, and (3) deformable multi-modal registration. Optimal parameter values were estimated for each step on the basis of relevant individual metrics. Tumor localization was then carried out via a model-based approach that takes both mpMRI and prior clinical knowledge features as input. A sequential optimization approach was adopted for determining the optimal parameters and techniques in each step of the pipeline.Main results. The application of bias field correction alone increased the accuracy of tumor localization (area under the curve (AUC) = 0.77;p-value = 0.004) over unprocessed data (AUC = 0.74). Adding normalization to the pre-processing pipeline further improved diagnostic accuracy of the model to an AUC of 0.85 (p-value = 0.000 12). Multi-modal registration of apparent diffusion coefficient images to T2-weighted images improved the alignment of tumor locations in all but one patient, resulting in a slight decrease in accuracy (AUC = 0.84;p-value = 0.30).Significance. Overall, our findings suggest that the combined effect of multiple pre-processing steps with optimal values has the ability to improve the quantitative classification of prostate cancer using mpMRI. Clinical trials: NCT03378856 and NCT03367702.

Harmonization Strategies in Multicenter MRI-Based Radiomics

Radiomics analysis is a powerful tool aiming to provide diagnostic and prognostic patient information directly from images that are decoded into handcrafted features, comprising descriptors of shape, size and textural patterns. Although radiomics is gaining momentum since it holds great promise for accelerating digital diagnostics, it is susceptible to bias and variation due to numerous inter-patient factors (e.g., patient age and gender) as well as inter-scanner ones (different protocol acquisition depending on the scanner center). A variety of image and feature based harmonization methods has been developed to compensate for these effects; however, to the best of our knowledge, none of these techniques has been established as the most effective in the analysis pipeline so far. To this end, this review provides an overview of the challenges in optimizing radiomics analysis, and a concise summary of the most relevant harmonization techniques, aiming to provide a thorough guide to the radiomics harmonization process.

Phantom Study on the Robustness of MR Radiomics Features: Comparing the Applicability of 3D Printed and Biological Phantoms

The objectives of our study were to (a) evaluate the feasibility of using 3D printed phantoms in magnetic resonance imaging (MR) in assessing the robustness and repeatability of radiomic parameters and (b) to compare the results obtained from the 3D printed phantoms to metrics obtained in biological phantoms. To this end, three different 3D phantoms were printed: a Hilbert cube (5 × 5 × 5 cm³) and two cubic quick response (QR) code phantoms (a large phantom (large QR) (5 × 5 × 4 cm³) and a small phantom (small QR) (4 × 4 × 3 cm³)). All 3D printed and biological phantoms (kiwis, tomatoes, and onions) were scanned thrice on clinical 1.5 T and 3 T MR with 1 mm and 2 mm isotropic resolution. Subsequent analyses included analyses of several radiomics indices (RI), their repeatability and reliability were calculated using the coefficient of variation (CV), the relative percentage difference (RPD), and the interclass coefficient (ICC) parameters. Additionally, the readability of QR codes obtained from the MR images was examined with several mobile phones and algorithms. The best repeatability (CV ≤ 10%) is reported for the acquisition protocols with the highest spatial resolution. In general, the repeatability and reliability of RI were better in data obtained at 1.5 T (CV = 1.9) than at 3 T (CV = 2.11). Furthermore, we report good agreements between results obtained for the 3D phantoms and biological phantoms. Finally, analyses of the read-out rate of the QR code revealed better texture analyses for images with a spatial resolution of 1 mm than 2 mm. In conclusion, 3D printing techniques offer a unique solution to create textures for analyzing the reliability of radiomic data from MR scans.

Radiomics software comparison using digital phantom and patient data: IBSI-compliance does not guarantee concordance of feature values

Introduction Radiomics is a promising imaging-based tool which could enhance clinical observation and identify representative features. To avoid different interpretations, the Image Biomarker Standardisation Initiative (IBSI) imposed conditions for harmonisation. This study evaluates IBSI-compliant radiomics applications against a known benchmark and clinical datasets for agreements. Materials and methods The applications compared were RadiomiX Research Toolbox, LIFEx v7.0.0, and syngo.via Frontier Radiomics v1.2.5 (based on PyRadiomics v2.1). Basic assessment included comparing feature names and their formulas. The IBSI digital phantom was used for evaluation against reference values. For agreement evaluation (including same software but different versions), two clinical datasets were used: 27 contrast-enhanced computed tomography (CECT) of colorectal liver metastases and 39 magnetic resonance imaging (MRI) of breast cancer, including intravoxel incoherent motion and dynamic contrast-enhanced (DCE) MRI. The lower 95% confidence interval of intraclass correlation coefficient was used to define excellent agreement (>0.9). Results The three radiomics applications share 41 (3 shape, 8 intensity, 30 texture) out of 172, 84 and 110 features for RadiomiX, LIFEx and syngo.via, respectively, as well as wavelet filtering. The naming convention is, however, different between them. Syngo.via had excellent agreement with the IBSI benchmark, while LIFEx and RadiomiX showed slightly worse agreement. Excellent reproducibility was achieved for shape features only, while intensity and texture features varied considerably with the imaging type. For intensity, excellent agreement ranged from 46% for the DCE maps to 100% for CECT, while this lowered to 44% and 73% for texture features, respectively. Wavelet features produced the greatest variation between applications, with an excellent agreement for only 3% to 11% features. Conclusion Even with IBSI-compliance, the reproducibility of features between radiomics applications is not guaranteed. To evaluate variation, quality assurance of radiomics applications should be performed and repeated when updating to a new version or adding a new modality.

Radiomics: a primer on high-throughput image phenotyping

Radiomics is a high-throughput approach to image phenotyping. It uses computer algorithms to extract and analyze a large number of quantitative features from radiological images. These radiomic features collectively describe unique patterns that can serve as digital fingerprints of disease. They may also capture imaging characteristics that are difficult or impossible to characterize by the human eye. The rapid development of this field is motivated by systems biology, facilitated by data analytics, and powered by artificial intelligence. Here, as part of Abdominal Radiology's special issue on Quantitative Imaging, we provide an introduction to the field of radiomics. The technique is formally introduced as an advanced application of data analytics, with illustrating examples in abdominal radiology. Artificial intelligence is then presented as the main driving force of radiomics, and common techniques are defined and briefly compared. The complete step-by-step process of radiomic phenotyping is then broken down into five key phases. Potential pitfalls of each phase are highlighted, and recommendations are provided to reduce sources of variation, non-reproducibility, and error associated with radiomics.

Effect of Gray Value Discretization and Image Filtration on Texture Features of the Pancreas Derived from Magnetic Resonance Imaging at 3T

Radiomics of pancreas magnetic resonance (MR) images is positioned well to play an important role in the management of diseases characterized by diffuse involvement of the pancreas. The effect of image pre-processing configurations on these images has been sparsely investigated. Fifteen individuals with definite chronic pancreatitis (an exemplar diffuse disease of the pancreas) and 15 healthy individuals were included in this age- and sex-matched case-control study. MR images of the pancreas were acquired using a single 3T scanner. A total of 93 first-order and second-order texture features of the pancreas were compared between the study groups, by subjecting MR images of the pancreas to 7 image pre-processing configurations related to gray level discretization and image filtration. The studied parameters of intensity discretization did not vary in terms of their effect on the number of significant first-order texture features. The number of statistically significant first-order texture features varied after filtering (7 with the use of logarithm filter and 3 with the use of Laplacian of Gaussian filter with 5 mm σ). Intensity discretization generally affected the number of significant second-order texture features more markedly than filtering. The use of fixed bin number of 16 yielded 42 significant second-order texture features, fixed bin number of 128–38 features, fixed bin width of 6–24 features, and fixed bin width of 42–26 features. The specific parameters of filtration and intensity discretization had differing effects on radiomics signature of the pancreas. Relative discretization with fixed bin number of 16 and use of logarithm filter hold promise as pre-processing configurations of choice in future radiomics studies in diffuse diseases of the pancreas.

An Investigation on Radiomics Feature Handling for HNSCC Staging Classification

The incidence of Head and Neck Squamous Cell Carcinoma (HNSCC) has been growing in the last few decades. Its diagnosis is usually performed through clinical evaluation and analyzing radiological images, then confirmed by histopathological examination, an invasive and time-consuming operation. The recent advances in the artificial intelligence field are leading to interesting results in the early diagnosis, personalized treatment and monitoring of HNSCC only by analyzing radiological images, without performing a tissue biopsy. The large amount of radiological images and the increasing interest in radiomics approaches can help to develop machine learning (ML) methods to support diagnosis. In this work, we propose an ML method based on the use of radiomics features, extracted from CT and PET images, to classify the disease in terms of pN-Stage, pT-Stage and Overall Stage. After the extraction of radiomics features, a selection step is performed to remove dataset redundancy. Finally, ML methods are employed to complete the classification task. Our pipeline is applied on the “Head-Neck-PET-CT” TCIA open-source dataset, considering a cohort of 201 patients from four different institutions. An AUC of 97%, 83% and 93% in terms of pN-Stage, pT-Stage and Overall Stage classification, respectively, is achieved. The obtained results are promising, showing the potential efficiency of the use of radiomics approaches in staging classification.

Clinical and radiomics prediction of complete response in rectal cancer pre-chemoradiotherapy

Background and purpose

Patients with rectal cancer could avoid major surgery if they achieve clinical complete response (cCR) post neoadjuvant treatment. Therefore, prediction of treatment outcomes before treatment has become necessary to select the best neo-adjuvant treatment option. This study investigates clinical and radiomics variables' ability to predict cCR in patients pre chemoradiotherapy.

Materials and methods

Using the OnCoRe database, we recruited a matched cohort of 304 patients (152 with cCR; 152 without cCR) deriving training (N = 200) and validation (N = 104) sets. We collected pre-treatment MR (magnetic resonance) images, demographics and blood parameters (haemoglobin, neutrophil, lymphocyte, alkaline phosphate and albumin). We segmented the gross tumour volume on T2 Weighted MR Images and extracted 1430 stable radiomics features per patient. We used principal component analysis (PCA) and receiver operating characteristic area under the curve (ROC AUC) to reduce dimensionality and evaluate the models produced.

Results

Using Logistic regression analysis, PCA-derived combined model (radiomics plus clinical variables) gave a ROC AUC of 0.76 (95% CI: 0.69-0.83) in the training set and 0.68 (95% CI 0.57-0.79) in the validation set. The clinical only model achieved an AUC of 0.73 (95% CI 0.66-0.80) and 0.62 (95% CI 0.51-0.74) in the training and validation set, respectively. The radiomics model had an AUC of 0.68 (95% CI 0.61-0.75) and 0.66 (95% CI 0.56-0.77) in the training and validation sets.

Conclusion

The predictive characteristics of both clinical and radiomics variables for clinical complete response remain modest but radiomics predictability is improved with addition of clinical variables.

Robustness Analysis of DCE-MRI-Derived Radiomic Features in Breast Masses: Assessing Quantization Levels and Segmentation Agreement

Machine learning models based on radiomic features allow us to obtain biomarkers that are capable of modeling the disease and that are able to support the clinical routine. Recent studies have shown that it is fundamental that the computed features are robust and reproducible. Although several initiatives to standardize the definition and extraction process of biomarkers are ongoing, there is a lack of comprehensive guidelines. Therefore, no standardized procedures are available for ROI selection, feature extraction, and processing, with the risk of undermining the effective use of radiomic models in clinical routine. In this study, we aim to assess the impact that the different segmentation methods and the quantization level (defined by means of the number of bins used in the feature-extraction phase) may have on the robustness of the radiomic features. In particular, the robustness of texture features extracted by PyRadiomics, and belonging to five categories—GLCM, GLRLM, GLSZM, GLDM, and NGTDM—was evaluated using the intra-class correlation coefficient (ICC) and mean differences between segmentation raters. In addition to the robustness of each single feature, an overall index for each feature category was quantified. The analysis showed that the level of quantization (i.e., the `bincount’ parameter) plays a key role in defining robust features: in fact, in our study focused on a dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) dataset of 111 breast masses, sets with cardinality varying between 34 and 43 robust features were obtained with `binCount’ values equal to 256 and 32, respectively. Moreover, both manual segmentation methods demonstrated good reliability and agreement, while automated segmentation achieved lower ICC values. Considering the dependence on the quantization level, taking into account only the intersection subset among all the values of `binCount’ could be the best selection strategy. Among radiomic feature categories, GLCM, GLRLM, and GLDM showed the best overall robustness with varying segmentation methods.

The stability of oncologic MRI radiomic features and the potential role of deep learning: a review

The use of MRI radiomic models for the diagnosis, prognosis and treatment response prediction of tumors has been increasingly reported in literature. However, its widespread adoption in clinics is hampered by issues related to features stability. In the MRI radiomic workflow, the main factors that affect radiomic features computation can be found in the image acquisition and reconstruction phase, in the image pre-processing steps, and in the segmentation of the region of interest on which radiomic indices are extracted. Deep Neural Networks (DNNs), having shown their potentiality in the medical image processing and analysis field, can be seen as an attractive strategy to partially overcome the issues related to radiomic stability and mitigate their impact. In fact, DNN approaches can be prospectively integrated in the MRI radiomic workflow to improve image quality, obtain accurate and reproducible segmentations and generate standardized images. In this review, DNN methods that can be included in the image processing steps of the radiomic workflow are described and discussed, in the light of a detailed analysis of the literature in the context of MRI radiomic reliability.

A tri-modal tissue-equivalent anthropomorphic phantom for PET, CT and multi-parametric MRI radiomics

PURPOSE

Radiomics has emerged as an advanced image processing methodology to define quantitative imaging biomarkers for prognosis and prediction of treatment response and outcome. The development of quantitative imaging biomarkers requires careful analysis to define their accuracy, stability and reproducibility through phantom measurements. Few efforts were devoted to develop realistic anthropomorphic phantoms. In this work, we developed a multimodality image phantom suitable for PET, CT and multiparametric MRI imaging.

METHODS

A tissue-equivalent gel-based mixture was designed and tested for compatibility with different imaging modalities. Calibration measurements allowed to assess gel composition to simulate PET, CT and MRI contrasts of oncological lesions. The characterized gel mixture was used to create realistic synthetic lesions (e.g. lesions with irregular shape and non-uniform image contrast), to be inserted in a standard anthropomorphic phantom. In order to show phantom usefulness, issues related to accuracy, stability and reproducibility of radiomic biomarkers were addressed as proofs-of-concept.

RESULTS

The procedure for gel preparation was straightforward and the characterized gel mixture allowed to mime simultaneously oncological lesion contrast in CT, PET and MRI imaging. Proofs-of-concept studies suggested that phantom measurements can be customized for specific clinical situations and radiomic protocols.

CONCLUSIONS

We developed a strategy to manufacture an anthropomorphic, tissue-equivalent, multimodal phantom to be customized on specific radiomics protocols, for addressing specific methodological issues both in mono and multicentric studies.