Next Article in Journal
CDK4/6 Inhibitors in Pancreatobiliary Cancers: Opportunities and Challenges
Next Article in Special Issue
Chimp Optimization Algorithm Influenced Type-2 Intuitionistic Fuzzy C-Means Clustering-Based Breast Cancer Detection System
Previous Article in Journal
Prediction of the Malignancy of a Breast Lesion Detected on Breast Ultrasound: Radiomics Applied to Clinical Practice
Previous Article in Special Issue
Region Segmentation of Whole-Slide Images for Analyzing Histological Differentiation of Prostate Adenocarcinoma Using Ensemble EfficientNetB2 U-Net with Transfer Learning Mechanism
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MR Intensity Normalization Methods Impact Sequence Specific Radiomics Prognostic Model Performance in Primary and Recurrent High-Grade Glioma

by
Patrick Salome
1,2,3,4,*,
Francesco Sforazzini
1,2,3,
Gianluca Brugnara
5,
Andreas Kudak
4,6,7,
Matthias Dostal
4,6,7,
Christel Herold-Mende
8,9,
Sabine Heiland
5,
Jürgen Debus
3,4,6,
Amir Abdollahi
1,3,4,6 and
Maximilian Knoll
1,3,4,6,*
1
Clinical Cooperation Unit (CCU) Radiation Oncology, German Cancer Research Centre, INF 280, 69120 Heidelberg, Germany
2
Heidelberg Medical Faculty, Heidelberg University, 69120 Heidelberg, Germany
3
German Cancer Consortium (DKTK) Core Centre Heidelberg, 69120 Heidelberg, Germany
4
Heidelberg Ion-Beam Therapy Centre (HIT), INF 450, 69120 Heidelberg, Germany
5
Department of Neuroradiology, Heidelberg University Hospital, 69120 Heidelberg, Germany
6
Department of Radiation Oncology, Heidelberg University Hospital, INF 400, 69120 Heidelberg, Germany
7
CCU Radiation Therapy, German Cancer Research Centre, INF 280, 69120 Heidelberg, Germany
8
Brain Tumour Group, European Organization for Research and Treatment of Cancer, 1200 Brussels, Belgium
9
Division of Neurosurgical Research, Department of Neurosurgery, Heidelberg University Hospital, 69120 Heidelberg, Germany
*
Authors to whom correspondence should be addressed.
Submission received: 7 January 2023 / Revised: 30 January 2023 / Accepted: 31 January 2023 / Published: 2 February 2023
(This article belongs to the Collection Artificial Intelligence and Machine Learning in Cancer Research)

Abstract

:

Simple Summary

As magnetic resonance (MR) intensities are acquired in arbitrary units, scans from different scanners are not directly comparable; thus, intensity normalization is essential. In this study, we assess the impact of normalization methods on prognostic radiomics models in primary and recurrent high-grade glioma on different MR sequences. Furthermore, we present a methodology that allows for the handling of radiomics performance discrepancy due to MR intensity normalization.

Abstract

Purpose: This study investigates the impact of different intensity normalization (IN) methods on the overall survival (OS) radiomics models’ performance of MR sequences in primary (pHGG) and recurrent high-grade glioma (rHGG). Methods: MR scans acquired before radiotherapy were retrieved from two independent cohorts (rHGG C1: 197, pHGG C2: 141) from multiple scanners (15, 14). The sequences are T1 weighted (w), contrast-enhanced T1w (T1wce), T2w, and T2w-FLAIR. Sequence-specific significant features (SF) associated with OS, extracted from the tumour volume, were derived after applying 15 different IN methods. Survival analyses were conducted using Cox proportional hazard (CPH) and Poisson regression (POI) models. A ranking score was assigned based on the 10-fold cross-validated (CV) concordance index (C-I), mean square error (MSE), and the Akaike information criterion (AICs), to evaluate the methods’ performance. Results: Scatter plots of the 10-CV C-I and MSE against the AIC showed an impact on the survival predictions between the IN methods and MR sequences (C1/C2 C-I range: 0.62–0.71/0.61–0.72, MSE range: 0.20–0.42/0.13–0.22). White stripe showed stable results for T1wce (C1/C2 C-I: 0.71/0.65, MSE: 0.21/0.14). Combat (0.68/0.62, 0.22/0.15) and histogram matching (HM, 0.67/0.64, 0.22/0.15) showed consistent prediction results for T2w models. They were also the top-performing methods for T1w in C2 (Combat: 0.67, 0.13; HM: 0.67, 0.13); however, only HM achieved high predictions in C1 (0.66, 0.22). After eliminating IN impacted SF using Spearman’s rank-order correlation coefficient, a mean decrease in the C-I and MSE of 0.05 and 0.03 was observed in all four sequences. Conclusion: The IN method impacted the predictive power of survival models; thus, performance is sequence-dependent.

1. Introduction

Radiomics, the extraction of features from medical images using data characterization algorithms, is an upcoming field of research expected to yield non-invasive surrogates for important molecular characteristics [1]. The high-dimensional database obtained can be used to create predictive models to help medical professionals make decisions about patient care, such as prognosis, diagnosis, and treatment outcome prediction [1]. Magnetic Resonance Imaging (MRI) has superior soft-tissue contrast, which allows for increased sensitivity and specificity in disease detection [2]. Several MR-based radiomics studies have been performed to determine image biomarkers that can help assess prognosis and improve treatment decisions [3,4,5,6]. Robust radiomics models often require large amounts of data; therefore, MR images are usually collected from multiple centres, sites, and scanners [7]. MR intensities are acquired in scanner-dependent arbitrary units, which leads to scans coming from different scanners and subjects not directly comparable, even when the same scanning protocol is implemented [8,9]. While this intensity variation has no major effects on the clinical diagnosis, it drastically impacts the performance of subsequent MRI preprocessing, such as image registration and segmentation, as well as radiomics feature calculation [10,11,12]. Therefore, intensity normalization methods should be implemented to deal with the intra- and inter-scan image intensity variations. This task has nonetheless been proven non-trivial, as speed, accuracy, and applicability can vary depending on the given data set. Multiple intensity normalization algorithms have been developed over time [13,14,15,16,17,18]. However, even though the image biomarker standardization initiative (IBSI) has defined a more general standardized radiomics image processing workflow, no specific guidelines on the proper choice of intensity normalization methods are currently present [19].
In the study of survival prediction in high-grade glioma patients using MR-based radiomics and deep learning, 23 publications were analyzed. The methods identified are z-score transformation or similar (30%), histogram-matching [13] (27%), the MR brain-specific white strip [20] (18%), tissue-based [21] (3%), and feature-based [22] (1%) [23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47]. Intensity normalization was not reported in the remaining studies (Supplementary Table S1). Furthermore, a single intensity normalization method was applied to all MR sequences. With the rise of published radiomics studies, the issue of reproducibility has become more prominent [48,49]. Therefore, more effort is being put into investigating the effects of different imaging preprocessing steps on the overall robustness and reproducibility of the radiomics models [50,51,52]. Most of these investigations have targeted grey-level discretization, i.e., the clustering of pixels based on intensity values to reduce feature calculation time and noise. Nonetheless, Carre et al. showed that the reproducibility of radiomics features is influenced by grey-level discretization and the normalization method chosen [53]. Noting that radiomics features can be split into four different groups, i.e., shape-based, first-order, second-order [54], and higher-order statistics [55], the authors also report that z-score transformation of first and second-order radiomics features show robust results. In a systematic review of intensity normalization of MRI prior to radiomics feature extraction in glioma datasets, Fatania et al. report that intensity normalization is a key preprocessing step in developing robust radiomics signatures and that few comparative studies of multiple methods exist [56]. Um et al. studied the effect of intensity normalization in radiomics survival model prediction [57]. They found that using histogram matching to normalize MR intensity in primary glioblastoma (GBM) patients improved patient stratification by reducing radiomics feature variability in T1 weighted (w), T1w post-contrast agent (T1wce), and fluid-attenuated inversion recovery (FLAIR) scans. The radiomics signature consisted of a combination of all three radiomics feature groups. Li et al. performed an in vitro and in vivo study by applying seven different normalization methods on T1w images. They demonstrated that the feature-based harmonization method Combat [22] significantly removes scanner effects in brain MR-based radiomics studies [58].
Our study builds upon the studies mentioned above. It aims to investigate the use of different normalization algorithms in multi-scanner brain MRI datasets and, more precisely, the performance of different methods on different sequences and their impact on the survival prediction model’s performance by analyzing the predictive power of the methods’ respective normalized dataset in the prediction of overall survival (OS). Noting that the radiomics survival prediction signatures identified in previous work included features from all radiomics feature groups [23,25,26,27,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,59], our analyses were directly performed on the signature obtained through two multi scanners high-grade glioma (HGG) datasets after applying a stringent feature selection pipeline. The first cohort is a recurrent HGG (rHGG) cohort of 197 patients (C1), and the second is a primary HGG (pHGG) cohort of 141 patients (C2). The MR sequences considered are T1w, T1wce, T2w, and FLAIR. Significant OS-correlated features were first identified through multiple feature reduction and resampling techniques from MR images acquired pre-radiotherapy (RT) and normalized through 15 different normalization approaches. Sequence-specific survival radiomics prediction models were next trained using Cox proportional hazard and Poisson survival regressions and applied to both cohorts. The performances of the intensity normalization algorithms were then compared based on the predictive power of their respective normalized dataset in the predictions of OS. Finally, features affected by the intensity normalization methods were further rejected, and comparisons with models trained with the remaining significant features were performed.

2. Materials and Methods

2.1. Datasets

This study analyses multiparametric pre-RT MR sequences from two independent data cohorts. Patients’ eligibility was based on the availability of clinical information of at least two MR sequences taken no longer than 30 days before RT and of the RT-DICOM data, specifically the DICOM structure set (SS). The first cohort (C1) consisted of 197 patients with pathologically confirmed rHGG collected retrospectively from 15 different MR scanners at the Heidelberg Ion-Beam Therapy Centre (HIT) and University Clinic Heidelberg (UKHD) from 2009 to 2018. All 197 patients received particle irradiation. The second cohort (C2) consists of 141 pHGG patients collected retrospectively from 14 different MR scanners at the UKHD from 2011 to 2016. All 141 patients received standard photon RT. Patients between cohorts were matched for the frequency in gender, tumour grade, and MR sequence. OS was calculated as the number of days between the start of the re-RT (C1) or RT(C2) and death. MR scans were acquired post-surgical tumour resection prior to radiotherapy treatment (RT). Conventional multislice (2D) acquired in the axial, sagittal, or coronal plane and 3D scans are present. The MR sequences found in the cohorts are the widely used sequences for brain tumour imaging [60] in clinical routines and trials [61,62]. However, the classes considered in this study are T1w, T1wce, T2w, and FLAIR. The four considered sequences were identified at different rates in both cohorts. The in-plane resolution ranged from 0.45 × 0.45 to 1.40 × 1.40 mm in the discovery cohort and 0.33 × 0.33 to 2 × 2 mm in the test set. Slice thickness ranged from 0.9 to 5 mm in all MR scans. A summary of both cohorts is shown in Table 1. An overview of the MR scanners and protocols found are reported in Supplementary Tables S2 and S3.

2.2. MRI Preprocessing Workflow

DICOM dataset curation and MR image classification were performed using pyCuRT and MR-Class [63]. All images were first reoriented to a common orientation. T1w images were corrected for signal inhomogeneities using the N4 bias field correction algorithm [64]. Brain extraction with the HD-BET tool was next performed [65]. When available, 3D MR sequences were mainly selected. Motion correction and volumetric image reconstruction were performed when 2D transversal, sagittal, and coronal MR scans were present. Reconstruction of the low-resolution 2D slices to a high-resolution 3D MR was performed using NiftyMic [66]. Next, cross-sectional linear co-registrations with six degrees of freedom (DOF) of the present MR images were performed on the T1wce using advanced normalization tools (ANTs) [67]. Furthermore, cross-sectional linear co-registrations with 6 DOF of the T1wce were performed on the RT planning CT. This registration was solely to generate the MR to CT transformation matrix, used to bring the target volume (TV) segmentations extracted from the DICOM SS objects to the MR space. Next, intensity normalization was performed. The different intensity normalization methods implemented in the comparison study are described in the next section. All images and segmentations were then resampled to a matrix size of 2 × 2 mm and a slice thickness of 2 mm using a cubic spline and linear interpolation, respectively. As for image discretization in an attempt to neutralize the impact of grey-level discretization on the overall result, five different bin counts were implemented, resulting in five sets of features per normalization algorithm. A bin count discretization approach was implemented since it was more frequently seen in HGG radiomics survival prediction studies. The image preprocessing diagram is shown in Figure 1.

2.3. Intensity Normalization Methods

Intensity normalization was performed with the help of the intensity normalization package by Reinhold et al. [68] and the FMRIB’s Automated Segmentation Tool (FAST) [69]. The intensity-normalization methods considered are: Fuzzy C-Means (FCM) [21] (9 different masks combinations), kernel density estimation (KDE), Gaussian mixture models (GMM) [70], Nyul’s and Udupa’s histogram matching-based abbreviated in this study as HM [13], white-strips (WS) [20], z-score normalization, and the feature-based batch adjustment method, i.e., Combat [22], resulting in 15 different MRI normalized datasets. A brief description of the methods is given in this section. For a broader description, we refer to the original normalization method papers and the manuscript by Reinhold et al. [68].

2.3.1. Standard Score

The standard score, also known as the z-score, represents the distance of a raw score from the mean measured in standard deviations. In MR brain image normalization, given that B is the brain mask in image I, the z-score calculates the mean μ and standard deviation σ of the intensities inside the brain image (excluding the background) as follows: μ = 1 B b B     I b ,   σ = b B     I b μ 2 B 1 with the normalized image being I norm x = I x μ σ . A disadvantage of this method is that the high intensities in the images are usually attenuated, risking a loss of information.

2.3.2. Fuzzy Clustering

Clustering is a method for analyzing data that aims to discover structures or groups in a data set. Fuzzy clustering allows a piece of data to be part of more than one cluster. In a fuzzy c-means algorithm, a data point is assigned a membership function, with 0 being the farthest from a cluster’s centre and one being the closest to a cluster’s centre, with the data point theoretically being able to belong to all clusters. Used as a normalization technique in brain MRI, the fuzzy c-means algorithm uses the segmentation of specific brain tissue to normalize the image to the mean intensity of the tissue. If the tissue mean is μ = 1 T t T     I t , then the normalized image is I norm x = I x μ , where x is the image voxels, and T is the tissue mask. The segmentations of the brain tissue masks, i.e., white matter (wm), grey matter (gm), and cerebrospinal fluid (csf), were performed using FSL’s FAST. In conjunction with the most common intensity value (mode) in a particular image, nine different mask combinations were implemented to generate nine fuzzy c-means normalized datasets. The masks are csf, gm, wm, csf-gm, wm-csf, wm-gm, csf-mode, wm-mode, and gm-mode. The normalization with two brain tissue masks is performed as: With µ 1 = 1 T 1 t T 1     I t and µ 2  = 1 T 2 t T 2     I t the normalized image is derived as Inorm(x) = I x a b a   with a = min( µ 1,   µ 2) and b = max( µ 1,   µ 2) The normalization with a brain tissue mask and the mode is performed as: as Inorm(x) = I x diff   with diff = µT − mode(B) with T as the tissue mask and B as the brain mask.

2.3.3. Kernel Density Estimation

A density estimator aims to find a function for the probability distribution from which a dataset is generated. The kernel density estimation (KDE) is an empirical calculation in a parametrized form. The formula for calculating the KDE for the probability distribution function is p x = 1 N M L h i = 1 N M L   K x x i h , where N, M, and L are the sizes of the images, K is the kernel (normalized to one), and h is the bandwidth parameter that scales the kernel. This method provides a smoother version of the histogram, making it easier to find the maxima π , which is used to normalize the entire image as I norm x =   c     I x π , where c is a positive, real constant. For the MR brain images, the KED finds the peak of the white matter histogram and translates it to a standard value.

2.3.4. Mixture Models

A mixture model assumes that a data set comprises subsets whose individual distributions are the respective probability distributions in the overall data set. A specific mixture model is the Gaussian mixture model, where the subsets are considered to be generated from a finite number of Gaussian distributions with undefined parameters. The method used in our study fits three Gaussian distributions to the histogram of the brain mask and normalizes the white matter mean to a standard value.

2.3.5. Landmark-Based Histogram Matching

The landmark-based histogram matching method by Nyúl et al. deforms the input image intensity histogram to match a reference histogram. The reference histogram is commonly obtained by averaging histograms in a data set and setting the landmarks of interest. Each input image histogram is then matched to the reference through linear interpolation based on the defined landmarks, usually quantiles.

2.3.6. White Stripe Normalization

The white stripe normalization approach by Shinohara et al. normalizes an image based on the normal-appearing white matter (NAWM) [20]. The NAWM values are obtained through a smoothening of the image histogram, followed by selecting the largest peak μ. The so-called white stripe contains intensity values up to 10% around μ. The white stripe can be defined as Ω T = I x   |   F 1 F μ   τ <   I x <   F 1 F μ +   τ , where F(x) is the cumulative distribution function of the image I and τ = 5%. If σ is the standard deviation in the white stripe, the normalized image is Inorm(x) = I x   μ σ .

2.3.7. Combat

Combat is a feature-based method originally developed for microarray expression data [22]. However, it has also been applied in imaging data and radiomics studies in recent years [71,72,73]. It eliminates batch effects through a known batch covariate using parametric or non-parametric empirical Bayes frameworks. In this study, an empirical Bayes Combat method was applied through the Surrogate Variable Analysis (sva) package (v 3.20.0) to eliminate batch effects due to the MR scanner. Adjustment of the following covariates was performed: age, tumour grade, and gender.

2.4. Comparison Study Design

After MR image preprocessing, radiomics features were calculated automatically from the gross tumour volume (GTV) segmentations extracted from the DICOM RT structure set and the original image, as well as from derived images from each normalized/discretized dataset using Pyradiomics (v 3.0) [74]. The derived images were retrieved from first Wavelet filtering, which yielded eight decompositions per level, each representing a combination of either a high or a low pass filter in each of the three dimensions, and then by applying a Laplacian of Gaussian filter with spatial scaling factors (SSFs) of 2, 3, and 4 mm. The total yielded features were around 1200 per MR sequence. The different feature classes and corresponding feature numbers can be seen in Table 2.
A Spearman rank-order correlation coefficient was used on the total number of features to exclude redundant features (rs > 0.80). Three feature selection methods, including a univariate analysis under Cox proportional hazard (CPH) models (p < 0.05), a random forest (RF)-based method, and lasso regression, were applied on 1000 random subsamples of C1 and C2 (10% left out) separately to identify features correlated to OS. Sequence-specific significant features identified at least 950 times were selected, and survival analyses were conducted using CPH [75] and Poisson survival regressions (POI) models [76]. A ranking score was next assigned to each normalization approach based on the converted standardized z score of the CPH averaged 10-fold cross-validated (CV) concordance index (C-I), the POI averaged 10-fold CV mean square error (mse), and the respective Akaike information criterion (AIC) of the OS prediction models. Lastly, after identifying the top-ranked methods for the different MR sequences, correlation heatmaps between the different normalization approaches for each significant feature forming the sequence-specific radiomics signature were plotted. Stable features that showed a high correlation (rs > 0.80) between at least 12 intensity normalization methods were further used to train CPH and POI models again, and the effects on the model predictions were studied. Finally, the performance of the feature-based method Combat was assessed in combination with the top-ranked image-based normalization method for each sequence in both cohorts. A flowchart of the study design is shown in Figure 2.

3. Results

3.1. Performance Assessment of the Intensity Normalization Method-Specific Survival Prediction Models for the Different MR Sequence

Scatter plots of the CPH averaged (over the five bin counts investigated) C-index and POI averaged mse, plotted against the respective AIC, for the 15 different intensity normalization-specific OS models derived from cohorts C1 and C2 are shown in Figure 3 and Figure 4. The OS model derived from the non-normalized (nn) dataset is also plotted.
Due to result variations and for a better interpretation of the performance of the different intensity methods, a ranking score was assigned to each normalization approach based on the converted standardized z score of the CPH averaged 10-fold CV C-I, the POI averaged 10-fold CV mean mse and the respective AIC of the OS prediction models. Table 3 summarizes and ranks the performance scores of the intensity normalization methods for each of the four MR sequences considered in cohorts C1 and C2.
The white stripe method is ranked first for T1wce in both cohorts (C1/C2 10-fold CV C-I: 0.71/0.65, AIC: 1033/547, 10-CV mse: 0.21/0.14, AIC: 410/252). For T1w, the feature-based batch adjustment method, i.e., Combat, had the best performance in C1 (0.68, 964, 0.22), while z-score transformation in C2 (0.65, 494, 0.15, 239). Nevertheless, the HM method was ranked second for both cohorts (C1/C2, 0.66/0.64, 970/494, 0.21/0.15, 389/2371). Furthermore, the top two ranked methods for T2w in both cohorts were Combat (C1/C2 0.62/0.67, 661/417, 0.22/0.13, 292/199) and the HM method (C1/C2 0.65/0.67, 667/415, 0.22/0.13, 294/200). As for T2w-FLAIR, the Fuzzy C-Means algorithm showed the best performance in C1 and C2, however, with different masks. For C1, the mask combination of wm and mode (0.67, 907, 0.21, 366) had the best performance, while the mask combination of wm and csf (0.72, 508, 0.15, 230) showed the best results for C2. Nevertheless, the former was ranked second in C2 (0.72, 517, 0.18, 235). Performance metrics of the remaining models in both cohorts are summarized in Supplementary Tables S4 and S5.

3.2. Significant Feature Correlation between the Normalized Datasets

Pairwise correlation tests were conducted to determine which features are impacted by the intensity normalization methods. Heatmaps displaying the Spearman correlation between the significant features and normalization methods for the various bin counts were generated. An example of T1wce significant features from C1 and bin count 32 is shown in Figure 5. The remaining heatmaps can be seen in Supplementary Figure S1.
Although all features were found to be significantly correlated with OS, some feature distributions varied when different intensity normalization methods were applied. To assess the performance of the top-ranked image normalization method before and after the elimination of the intensity normalization impacted significant feature for cohort C1 and C2 for each MR sequence, the 10-CV C-I and mse of the CPH and POI models with only the stable features that have a high correlation (rs > 0.8) between at least 12 methods are reported (Table 4). Figure 6 displays boxplots showing the differences in C-I and MSE before and after the elimination of the intensity normalization impacted significant features for each modality in both cohorts.

3.3. Performance Comparison of the Feature-Based and Top-Ranked Image-Based Normalisation Methods

Table 5 summarizes the performance of the top-ranked image normalization method separate and in combination with the feature-based method Combat for cohorts C1 and C2. Since Combat ranked first for the T1w models from C1 and T2w models from C2, the second-ranked method, i.e., the HM method, was the image-based intensity normalization method for these two datasets.

4. Discussion

This study evaluated the impact of MRI intensity normalization algorithms on MR-based radiomics survival prediction models in primary and recurrent high-grade glioma. The sequences considered are T1w, T1wce, T2w, and T2w-FLAIR. Performance assessment of the intensity normalization method-specific CPH and POI survival prediction models showed an impact on the survival predictions between the different intensity normalization methods and the different MR sequences. Therefore, it can be concluded that the MR intensity normalization approach directly impacts the overall power of the radiomics-based MR predictive models. Moreover, considering the variability of the acquired results for the different MR sequences, it can be seen that the intensity normalization algorithm performance is correlated with the MR sequence and that the problem cannot be simplified to one intensity normalization method.
Due to these variations and for a better interpretation of the results, the ranking score was developed. The WS method showed promising results in T1wce models as it was ranked first in two independent multi-scanner datasets. Combat and the HM method showed consistent prediction results between the two cohorts for T2w models. These two methods were the top-performing methods for T1w in C1; however, only HM achieved high predictions in C2. Combat performance in C2 might be due to the higher number of batches and the number of images per batch, as 22% of T1w images in C2 were missing, making batch effect removal more challenging. As for T2w-FLAIR, the FCM showed favorable results in both cohorts; however, with different mask combinations, including the wm and csf or wm and mode. A tighter intensity range is observed in T2w-FLAIR than the other sequences, as csf signals are attenuated. These results might indicate that a mask-based normalization approach might be more favourable when dealing with images with tighter intensity ranges. The application of both an image-based and feature-based normalization method had little impact on the performance of the CPH and POI models. Exceptions were observed in the dataset where Combat was ranked first, i.e., T1w in C1 and T2w in C2.
As CPH models were part of the radiomics signature building pipeline, POI models were also trained to assess whether model performances were biased to CPH models. Comparably to CPH models, the impact of the intensity normalization methods was also observed in the POI models. Furthermore, the performance of both models was similarly affected after the elimination of the intensity normalization impacted significant features. A mean decrease in the 10-CV C-I and 10-CV MSE of 0.05 and 0.03 was observed in all four sequences across both cohorts. The use of a correlation coefficient between different normalization methods as a feature robustness check leads to a trade-off between model stability and the risk of eliminating important imaging biomarkers.
This study included two independent HGG cohorts collected from a single university hospital, UKHD. However, since the data cohorts included data between 2008 and 2019, 19 different scanners from three vendors with a 0.5 to 3.0 Tesla range were identified. Intensity normalization improved OS prediction in radiomics survival models as the non-normalized datasets generally ranked low in both cohorts. Therefore, the need for normalization is based on the number of scanners and image protocols, not just the number of centers. However, an exemption is seen in the T1w dataset in C1. This exception might be because many images in the T1w dataset from C1 were reconstructed using NiftyMic (as mostly 2D MR scans were present) and therefore preprocessed before applying the intensity normalization methods [66].
Since multiple MR scanners were found in both cohorts, where some have been withdrawn from clinical practice, phantoms could not be applied to assess the impact of the intensity normalization methods. Therefore, the hard endpoint OS was used in this study as a possible appropriate surrogate.
In literature, multiple intensity normalization methods have been reported in HGG radiomics studies, where all implemented the same method across all MR sequences [56]. However, as demonstrated in this study, the performance of the different methods varies. This study shows that the variations are big and that if reproducibility of the radiomics model is to be possible, the method of intensity normalization should be reported. Another way is to eliminate features impacted by the different normalization methods. When unstable features are impacted, the performance of the individual MR sequence prediction models is reduced, a necessary trade-off for stable radiomics models. However, combining multiple stable radiomics signatures from multiple MR sequences or modalities might mitigate that reduction and lead to high survival prediction models.
All in all, the main strengths of our study are as follows. First, we found that different intensity normalization methods produce varied results across different sequences. Therefore, to generate a good predictive model in the context of MRI radiomics, it may be necessary to apply different normalization methods to different sequences. Secondly, the methodology presented in this study, where features impacted by normalization methods are screened using Spearman pairwise correlations, can result in more generalizable models and facilitate the reproducibility of MRI-based radiomics models.
The following limitations exist in this work. The application of different preprocessing methods makes it generally hard to assess the impact of different normalization methods seamlessly. The changes in the radiomics values are as much affected by other preprocessing methods as image discretization or delineating the region of interest, which suggests that the application of intensity normalization alone may not be enough. In this study, we attempted to limit the effect of intensity discretization by applying five different bin counts and reporting the average score. Nevertheless, as demonstrated by many radiomics robustness studies, the overall performance and reproducibility of the radiomics models are indeed affected by choice of discretization approach [50,51,52]. Nevertheless, similarly to using correlation coefficient heatmaps between the different normalization methods to determine stable radiomics features, the same can be implemented across different bin counts or widths. Furthermore, all GTVs were segmented following institutional guidelines for RT treatment. Nonetheless, delineation variabilities are known to impact radiomics features, and the impact of intensity normalization and ROI segmentation should also be evaluated. As automatic tumour segmentation networks become more robust and popular, these inter-observer variabilities would be reduced, thus eliminating another layer of variability. However, as segmentation networks are also impacted by the intensity normalization method, future work will evaluate the performance of different normalization methods on automatic segmentation networks. Another limitation of our approach is the need for a sufficient amount of observation for the survival models since these models rely on having enough data to accurately reflect the underlying relationship between the predictors and the response variable. A lack of observations can lead to unreliable or unstable model estimates and increase the risk of overfitting.
Moreover, differences in the performance of the different IN methods across both cohorts can be possibly due to the alterations in the structure of intra-tumour heterogeneity, which differ between pHGG and rHGG, as well as the difference in the treatment of rHGG in comparison to pHGG, since the treatment of rHGG is not standardized as for pHGG, i.e., incorporating surgery, postoperative adjuvant RT, and adjuvant chemotherapy [77]. In addition, heterogeneity of cohorts, such as in MGMT methylation, IDH1/2 mutation, and 1p/19q deletion, can also lead to survival prediction differences [78]. More detailed studies are needed to assess the impact in more stringently defined cohorts.

5. Conclusions

Variations in the results for the different MR sequences showed that the intensity normalization method performance is sequence-dependent and directly impacts the predictive power of glioma survival models. Therefore, the documentation of the adapted normalization approach is highly recommended and necessary to enable the reproducibility of the MRI-based radiomics model. The methodology presented in this study can be further implemented in different entities to determine the stable radiomics features for signature building. Future work includes the study of additional sequences and anatomy sites. One major limitation of this study is the difficulty in solely evaluating the effect of a specific preprocessing method due to the application of various preprocessing techniques. Therefore, future work will involve assessing the various preprocessing methods using the methodology proposed in this study.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/cancers15030965/s1, Table S1: Intensity normalization algorithms applied in MR-based radiomics and deep learning-based survival prediction studies in high-grade glioma patients. Table S2: MR scanner models found in the cohorts; Table S3: MR image protocols found in the cohorts; Table S4: Model performance metrics for each MR sequence and normalization method for cohort 1; Table S5: Model performance metrics for each MR sequence and normalization method for cohort 2; Figure S1: Correlation heatmaps between the 16 different normalization methods considered and the reference non-normalized dataset for each MR sequence in cohorts C1 and C2 discretized with bin counts of 16, 32, 48, 64, and 128.

Author Contributions

Conceptualization, P.S., F.S., A.A. and M.K.; Data curation, P.S. and F.S.; Formal analysis, P.S., F.S. and M.K.; Funding acquisition, J.D. and A.A.; Investigation, G.B., C.H.-M. and S.H.; Methodology, P.S. and M.K.; Resources, A.K. and M.D.; Software, P.S. and F.S.; Supervision, J.D., A.A. and M.K.; Visualization, P.S.; Writing—original draft, P.S.; Writing—review and editing, F.S., G.B., C.H.-M., S.H., A.A. and M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie Innovative Training Network grant agreement (MSCA-ITN-PREDICT, No 766276), the collaborative research center of the German Research Foundation (DFG, Unite, SFB-1389, Project 404521405), Zentrum für Personalisierte Medizin (ZPM-Network BW, Project PROMI), the intramural funds of the National Center for Tumor Diseases (NCT) and the German Cancer Consortium (DKTK) Radiation Oncology programs.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of the Medical Faculty of Heidelberg University (approval number S-540/2010, date of last updated approval: 20 July 2020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the S-540 study.

Data Availability Statement

C1 and C2 are available from the corresponding author on reasonable request.

Conflicts of Interest

P.S. No relevant relationships. F.S. No relevant relationships. A.K. No relevant relationships. N.B. No relevant relationships. J.D. Grants/contracts from/with Viewray, CRI—The Clinical Research Institute, Accuray International Sarl, RaySearch Laboratories, Vision RT, Merck Serono, Astellas Pharma, AstraZeneca, Siemens Healthcare, Solution Akademie, Ergomed PLC Surrey Research Park, Quintiles, Pharmaceutical Research Associates, Boehringer Ingelheim Pharma & CoKG, PTW-Freiburg Dr. Pychlau, Nanobiotix, Accuray, Varian; participation on a data safety monitoring board or advisory board for Merck Serono. A.A. Predict MarieCurie innovative training network (ITN), in frame of Horizon 2020 from the European Union, Marie Skłodowska-Curie grant agreement No 766276. M.K. No relevant relationships.

References

  1. Lambin, P.; Rios-Velazquez, E.; Leijenaar, R.; Carvalho, S.; van Stiphout, R.G.; Granton, P.; Zegers, C.M.; Gillies, R.; Boellard, R.; Dekker, A.; et al. Radiomics: Extracting more information from medical images using advanced feature analysis. Eur. J. Cancer 2012, 48, 441–446. [Google Scholar] [CrossRef] [PubMed]
  2. Liang, Z.-P.; Lauterbur, P.C. Principles of Magnetic Resonance Imaging: A Signal Processing Perspective; SPIE Optical Engineering Press: Bellingham, WA, USA, 2000. [Google Scholar]
  3. Kickingereder, P.; Bonekamp, D.; Nowosielski, M.; Kratz, A.; Sill, M.; Burth, S.; Wick, A.; Eidel, O.; Schlemmer, H.-P.; Radbruch, A.; et al. Radiogenomics of glioblastoma: Machine learning–based classification of molecular characteristics by using multiparametric and multiregional MR imaging features. Radiology 2016, 281, 907–918. [Google Scholar] [CrossRef] [PubMed]
  4. Bonekamp, D.; Kohl, S.; Wiesenfarth, M.; Schelb, P.; Radtke, J.P.; Götz, M.; Kickingereder, P.; Yaqubi, K.; Hitthaler, B.; Gählert, N.; et al. Radiomic machine learning for characterization of prostate lesions with MRI: Comparison to ADC values. Radiology 2018, 289, 128–137. [Google Scholar] [CrossRef] [PubMed]
  5. Vallières, M.; Freeman, C.R.; Skamene, S.R.; El Naqa, I. A radiomics model from joint FDG-PET and MRI texture features for the prediction of lung metastases in soft-tissue sarcomas of the extremities. Phys. Med. Biol. 2015, 60, 5471. [Google Scholar] [CrossRef]
  6. Tian, Q.; Yan, L.-F.; Zhang, X.; Zhang, X.; Hu, Y.-C.; Han, Y.; Liu, Z.-C.; Nan, H.-Y.; Sun, Q.; Sun, Y.-Z.; et al. Radiomics strategy for glioma grading using texture features from multiparametric MRI. J. Magn. Reson. Imaging 2018, 48, 1518–1528. [Google Scholar] [CrossRef]
  7. L’Heureux, A.; Grolinger, K.; Elyamany, H.F.; Capretz, M.A. Machine learning with big data: Challenges and approaches. IEEE Access 2017, 5, 7776–7797. [Google Scholar] [CrossRef]
  8. Stonnington, C.M.; Tan, G.; Klöppel, S.; Chu, C.; Draganski, B.; Jack, C.R., Jr.; Chen, K.; Ashburner, J.; Frackowiak, R.S. Interpreting scan data acquired from multiple scanners: A study with Alzheimer’s disease. Neuroimage 2008, 39, 1180–1185. [Google Scholar] [CrossRef]
  9. Preboske, G.M.; Gunter, J.L.; Ward, C.P.; Jack, C.R., Jr. Common MRI acquisition non-idealities significantly impact the output of the boundary shift integral method of measuring brain atrophy on serial MRI. Neuroimage 2006, 30, 1196–1202. [Google Scholar] [CrossRef]
  10. Collewet, G.; Strzelecki, M.; Mariette, F. Influence of MRI acquisition protocols and image intensity normalization methods on texture classification. Magn. Reson. Imaging 2004, 22, 81–91. [Google Scholar] [CrossRef]
  11. Alam, F.; Sami, U.; Aziz, U.; Fawad, Q. Medical image registration: Classification, applications and issues. J. Postgrad. Med. Inst. 2018, 32, 300–3007. [Google Scholar]
  12. Chu, R.; Hurwitz, S.; Tauhid, S.; Bakshi, R. Automated segmentation of cerebral deep gray matter from MRI scans: Effect of field strength on sensitivity and reliability. BMC Neurol. 2017, 17, 172. [Google Scholar] [CrossRef]
  13. Nyúl, L.G.; Udupa, J.K. On standardizing the MR image intensity scale. Magn. Reson. Med. 1999, 42, 1072–1081. [Google Scholar] [CrossRef]
  14. Shah, M.; Xiao, Y.; Subbanna, N.; Francis, S.; Arnold, D.L.; Collins, D.L.; Arbel, T. Evaluating intensity normalization on MRIs of human brain with multiple sclerosis. Med. Image Anal. 2011, 15, 267–282. [Google Scholar] [CrossRef]
  15. Jäger, F.; Deuerling-Zheng, Y.; Frericks, B.; Wacker, F.; Hornegger, J. A new method for MRI intensity standardization with application to lesion detection in the brain. Proc. Vis. Model. Vis. 2006, 2006, 276–296. [Google Scholar]
  16. Hellier, P. Consistent intensity correction of MR images. In Proceedings of the 2003 International Conference on Image Processing, Barcelona, Spain, 14–17 September 2003; Volume 1, p. I–1109. [Google Scholar]
  17. Loizou, C.P.; Pantziaris, M.; Seimenis, I.; Pattichis, C.S. Brain MR image normalization in texture analysis of multiple sclerosis. In Proceedings of the 2009 9th International Conference on Information Technology and Applications in Biomedicine, Larnaka, Cyprus, 4–7 November 2009; pp. 1–5. [Google Scholar]
  18. Bergeest, J.-P.; Jäger, F. A comparison of five methods for signal intensity standardization in MRI. In Bildverarbeitung für die Medizin 2008; Springer: Berlin/Heidelberg, Germany, 2008; pp. 36–40. [Google Scholar]
  19. Zwanenburg, A.; Vallières, M.; Abdalah, M.A.; Aerts, H.J.; Andrearczyk, V.; Apte, A.; Ashrafinia, S.; Bakas, S.; Beukinga, R.J.; Boellaard, R.; et al. The image biomarker standardization initiative: Standardized quantitative radiomics for high-throughput image-based phenotyping. Radiology 2020, 295, 328–338. [Google Scholar] [CrossRef]
  20. Shinohara, R.T.; Sweeney, E.M.; Goldsmith, J.; Shiee, N.; Mateen, F.J.; Calabresi, P.A.; Jarso, S.; Pham, D.L.; Reich, D.S.; Crainiceanu, C.M.; et al. Statistical normalization techniques for magnetic resonance imaging. NeuroImage Clin. 2014, 6, 9–19. [Google Scholar] [CrossRef]
  21. Bezdek, J.C.; Ehrlich, R.; Full, W. FCM: The fuzzy c-means clustering algorithm. Comput. Geosci. 1984, 10, 191–203. [Google Scholar] [CrossRef]
  22. Johnson, W.E.; Li, C.; Rabinovic, A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 2007, 8, 118–127. [Google Scholar] [CrossRef]
  23. Ruan, Z.; Mei, N.; Lu, Y.; Xiong, J.; Li, X.; Zheng, W.; Liu, L.; Yin, B. A Comparative and Summative Study of Radiomics-based Overall Survival Prediction in Glioblastoma Patients. J. Comput. Assist. Tomogr. 2022, 46, 470–479. [Google Scholar] [CrossRef]
  24. Lao, J.; Chen, Y.; Li, Z.-C.; Li, Q.; Zhang, J.; Liu, J.; Zhai, G. A Deep Learning-Based Radiomics Model for Prediction of Survival in Glioblastoma Multiforme. Sci. Rep. 2017, 7, 10353. [Google Scholar] [CrossRef]
  25. Li, Q.; Bai, H.; Chen, Y.; Sun, Q.; Liu, L.; Zhou, S.; Wang, G.; Liang, C.; Li, Z.-C. A Fully-Automatic Multiparametric Radiomics Model: Towards Reproducible and Prognostic Imaging Signature for Prediction of Overall Survival in Glioblastoma Multiforme. Sci. Rep. 2017, 7, 14331. [Google Scholar] [CrossRef] [PubMed]
  26. Liu, X.; Li, Y.; Qian, Z.; Sun, Z.; Xu, K.; Wang, K.; Liu, S.; Fan, X.; Li, S.; Zhang, Z.; et al. A radiomic signature as a non-invasive predictor of progression-free survival in patients with lower-grade gliomas. NeuroImage Clin. 2018, 20, 1070–1077. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, X.; Lu, H.; Tian, Q.; Feng, N.; Yin, L.; Xu, X.; Du, P.; Liu, Y. A radiomics nomogram based on multiparametric MRI might stratify glioblastoma patients according to survival. Eur. Radiol. 2019, 29, 5528–5538. [Google Scholar] [CrossRef] [PubMed]
  28. Yang, Y.; Han, Y.; Hu, X.; Wang, W.; Cui, G.; Guo, L.; Zhang, X. An Improvement of Survival Stratification in Glioblastoma Patients via Combining Subregional Radiomics Signatures. Front. Neurosci. 2021, 15, 683452. Available online: https://www.frontiersin.org/article/10.3389/fnins.2021.683452 (accessed on 18 May 2022). [CrossRef] [PubMed]
  29. Li, G.; Li, L.; Li, Y.; Qian, Z.; Wu, F.; He, Y.; Jiang, H.; Li, R.; Wang, D.; Zhai, Y.; et al. An MRI radiomics approach to predict survival and tumour-infiltrating macrophages in gliomas. Brain 2022, 145, 1151–1161. [Google Scholar] [CrossRef]
  30. Wang, J.; Zheng, X.; Zhang, J.; Xue, H.; Wang, L.; Jing, R.; Chen, S.; Che, F.; Heng, X.; Li, G.; et al. An MRI-based radiomics signature as a pretreatment non-invasive predictor of overall survival and chemotherapeutic benefits in lower-grade gliomas. Eur. Radiol. 2021, 31, 1785–1794. [Google Scholar] [CrossRef]
  31. Li, Z.; Liu, P.; An, T.; Yang, H.; Zhang, W.; Wang, J. Construction of a prognostic immune signature for lower grade glioma that can be recognized by MRI radiomics features to predict survival in LGG patients. Transl. Oncol. 2021, 14, 101065. [Google Scholar] [CrossRef]
  32. Chaddad, A.; Daniel, P.; Zhang, M.; Rathore, S.; Sargos, P.; Desrosiers, C.; Niazi, T. Deep radiomic signature with immune cell markers predicts the survival of glioma patients. Neurocomputing 2022, 469, 366–375. [Google Scholar] [CrossRef]
  33. Han, W.; Qin, L.; Bay, C.; Chen, X.; Yu, K.-H.; Miskin, N.; Li, A.; Xu, X.; Young, G. Deep Transfer Learning and Radiomics Feature Prediction of Survival of Patients with High-Grade Gliomas. Am. J. Neuroradiol. 2020, 41, 40–48. [Google Scholar] [CrossRef]
  34. Shboul, Z.A.; Alam, M.; Vidyaratne, L.; Pei, L.; Elbakary, M.I.; Iftekharuddin, K.M. Feature-Guided Deep Radiomics for Glioblastoma Patient Survival Prediction. Front. Neurosci. 2019, 13, 966. Available online: https://www.frontiersin.org/article/10.3389/fnins.2019.00966 (accessed on 18 May 2022). [CrossRef]
  35. Tan, Y.; Mu, W.; Wang, X.; Yang, G.; Gillies, R.J.; Zhang, H. Improving survival prediction of high-grade glioma via machine learning techniques based on MRI radiomic, genetic and clinical risk factors. Eur. J. Radiol. 2019, 120, 108609. [Google Scholar] [CrossRef]
  36. Choi, Y.S.; Ahn, S.S.; Chang, J.H.; Kang, S.-G.; Kim, E.H.; Kim, S.H.; Jain, R.; Lee, S.-K. Machine learning and radiomic phenotyping of lower grade gliomas: Improving survival prediction. Eur. Radiol. 2020, 30, 3834–3842. [Google Scholar] [CrossRef]
  37. Chaddad, A.; Daniel, P.; Desrosiers, C.; Toews, M.; Abdulkarim, B. Novel Radiomic Features Based on Joint Intensity Matrices for Predicting Glioblastoma Patient Survival Time. IEEE J. Biomed. Health Inform. 2019, 23, 795–804. [Google Scholar] [CrossRef]
  38. Bakas, S.; Shukla, G.; Akbari, H.; Erus, G.; Sotiras, A.; Rathore, S.; Sako, C.; Ha, S.M.; Rozycki, M.; Shinohara, R.T.; et al. Overall survival prediction in glioblastoma patients using structural magnetic resonance imaging (MRI): Advanced radiomic features may compensate for lack of advanced MRI modalities. JMI 2020, 7, 031505. [Google Scholar] [CrossRef]
  39. Baid, U.; Rane, S.U.; Talbar, S.; Gupta, S.; Thakur, M.H.; Moiyadi, A.; Mahajan, A. Overall Survival Prediction in Glioblastoma With Radiomic Features Using Machine Learning. Front. Comput. Neurosci. 2020, 14, 61. Available online: https://www.frontiersin.org/article/10.3389/fncom.2020.00061 (accessed on 18 May 2022). [CrossRef]
  40. Chaddad, A.; Sabri, S.; Niazi, T.; Abdulkarim, B. Prediction of survival with multi-scale radiomic analysis in glioblastoma patients. Med. Biol. Eng. Comput. 2018, 56, 2287–2300. [Google Scholar] [CrossRef]
  41. Tixier, F.; Um, H.; Bermudez, D.; Iyer, A.; Apte, A.; Graham, M.S.; Nevel, K.S.; Deasy, J.O.; Young, R.J.; Veeraraghavan, H. Preoperative MRI-radiomics features improve prediction of survival in glioblastoma patients over MGMT methylation status alone. Oncotarget 2019, 10, 660–672. [Google Scholar] [CrossRef]
  42. Han, K.; Ren, M.; Wick, W.; Abrey, L.; Das, A.; Jin, J.; Reardon, D.A. Progression-free survival as a surrogate endpoint for overall survival in glioblastoma: A literature-based meta-analysis from 91 trials. Neuro-Oncol. 2014, 16, 696–706. [Google Scholar] [CrossRef]
  43. Yan, J.; Zhang, B.; Zhang, S.; Cheng, J.; Liu, X.; Wang, W.; Dong, Y.; Zhang, L.; Mo, X.; Chen, Q.; et al. Quantitative MRI-based radiomics for noninvasively predicting molecular subtypes and survival in glioma patients. npj Precis. Oncol. 2021, 5, 72. [Google Scholar] [CrossRef]
  44. Prasanna, P.; Patel, J.; Partovi, S.; Madabhushi, A.; Tiwari, P. Radiomic features from the peritumoural brain parenchyma on treatment-naïve multiparametric MR imaging predict long versus short-term survival in glioblastoma multiforme: Preliminary findings. Eur. Radiol. 2017, 27, 4188–4197. [Google Scholar] [CrossRef]
  45. Bae, S.; Choi, Y.S.; Ahn, S.S.; Chang, J.H.; Kang, S.-G.; Kim, E.H.; Kim, S.H.; Lee, S.-K. Radiomic MRI phenotyping of glioblastoma: Improving survival prediction. Radiology 2018, 289, 797–806. [Google Scholar] [CrossRef] [PubMed]
  46. Kickingereder, P.; Burth, S.; Wick, A.; Götz, M.; Eidel, O.; Schlemmer, H.-P.; Maier-Hein, K.H.; Wick, W.; Bendszus, M.; Radbruch, A.; et al. Radiomic profiling of glioblastoma: Identifying an imaging predictor of patient survival with improved performance over established clinical and radiologic risk models. Radiology 2016, 280, 880–889. [Google Scholar] [CrossRef] [PubMed]
  47. Choi, Y.; Nam, Y.; Jang, J.; Shin, N.-Y.; Lee, Y.S.; Ahn, K.-J.; Kim, B.; Park, J.-S.; Jeon, S.; Hong, Y.G. Radiomics may increase the prognostic value for survival in glioblastoma patients when combined with conventional clinical and genetic prognostic models. Eur. Radiol. 2021, 31, 2084–2093. [Google Scholar] [CrossRef] [PubMed]
  48. Ibrahim, A.; Primakov, S.; Beuque, M.; Woodruff, H.; Halilaj, I.; Wu, G.; Refaee, T.; Granzier, R.; Widaatalla, Y.; Hustinx, R.; et al. Radiomics for precision medicine: Current challenges, future prospects, and the proposal of a new framework. Methods 2020, 188, 20–29. [Google Scholar] [CrossRef] [PubMed]
  49. Rizzo, S.; Botta, F.; Raimondi, S.; Origgi, D.; Fanciullo, C.; Morganti, A.G.; Bellomi, M. Radiomics: The facts and the challenges of image analysis. Eur. Radiol. Exp. 2018, 2, 36. [Google Scholar] [CrossRef]
  50. Molina, D.; Pérez-Beteta, J.; Martínez-González, A.; Martino, J.; Velásquez, C.; Arana, E.; Pérez-García, V.M. Influence of gray level and space discretization on brain tumour heterogeneity measures obtained from magnetic resonance images. Comput. Biol. Med. 2016, 78, 49–57. [Google Scholar] [CrossRef]
  51. Bologna, M.; Corino, V.; Mainardi, L. Virtual phantom analyses for preprocessing evaluation and detection of a robust feature set for MRI-radiomics of the brain. Med. Phys. 2019, 46, 5116–5123. [Google Scholar] [CrossRef]
  52. Duron, L.; Balvay, D.; Perre, S.V.; Bouchouicha, A.; Savatovsky, J.; Sadik, J.-C.; Thomassin-Naggara, I.; Fournier, L.; Lecler, A. Gray-level discretization impacts reproducible MRI radiomics texture features. PLoS ONE 2019, 14, e0213459. [Google Scholar] [CrossRef]
  53. Carré, A.; Klausner, G.; Edjlali, M.; Lerousseau, M.; Briend-Diop, J.; Sun, R.; Ammari, S.; Reuzé, S.; Andres, E.A.; Estienne, T.; et al. Standardization of brain MR images across machines and protocols: Bridging the gap for MRI-based radiomics. Sci. Rep. 2020, 10, 12340. [Google Scholar] [CrossRef]
  54. Haralick, R.M.; Shanmugam, K.; Dinstein, I. Textural Features for Image Classification. IEEE Trans. Syst. Man Cybern. 1973, SMC-3, 610–621. [Google Scholar] [CrossRef]
  55. Scapicchio, C.; Gabelloni, M.; Barucci, A.; Cioni, D.; Saba, L.; Neri, E. A deep look into radiomics. Radiol. Med. 2021, 126, 1296–1311. [Google Scholar] [CrossRef]
  56. Fatania, K.; Mohamud, F.; Clark, A.; Nix, M.; Short, S.C.; O’Connor, J.; Scarsbrook, A.F.; Currie, S. Intensity standardization of MRI prior to radiomic feature extraction for artificial intelligence research in glioma—A systematic review. Eur. Radiol. 2022, 32, 7014–7025. [Google Scholar] [CrossRef]
  57. Um, H.; Tixier, F.; Bermudez, D.; Deasy, J.O.; Young, R.J.; Veeraraghavan, H. Impact of image preprocessing on the scanner dependence of multiparametric MRI radiomic features and covariate shift in multi-institutional glioblastoma datasets. Phys. Med. Biol. 2019, 64, 165011. [Google Scholar] [CrossRef]
  58. Li, Y.; Ammari, S.; Balleyguier, C.; Lassau, N.; Chouzenoux, E. Impact of Preprocessing and Harmonization Methods on the Removal of Scanner Effects in Brain MRI Radiomic Features. Cancers 2021, 13, 3000. [Google Scholar] [CrossRef]
  59. Lu, Y.; Patel, M.; Natarajan, K.; Ughratdar, I.; Sanghera, P.; Jena, R.; Watts, C.; Sawlani, V. Machine learning-based radiomic, clinical and semantic feature analysis for predicting overall survival and MGMT promoter methylation status in patients with glioblastoma. Magnetic resonance imaging 2020, 74, 161–170. [Google Scholar] [CrossRef]
  60. Ellingson, B.M.; Bendszus, M.; Boxerman, J.; Barboriak, D.; Erickson, B.J.; Smits, M.; Nelson, S.J.; Gerstner, E.; Alexander, B.; Goldmacher, G.; et al. Consensus recommendations for a standardized Brain Tumour Imaging Protocol in clinical trials. Neuro-Oncol. 2015, 17, 1188–1198. [Google Scholar] [CrossRef]
  61. Combs, S.E.; Burkholder, I.; Edler, L.; Rieken, S.; Habermehl, D.; Jäkel, O.; Haberer, T.; Haselmann, R.; Unterberg, A.; Wick, W.; et al. Randomised phase I/II study to evaluate carbon ion radiotherapy versus fractionated stereotactic radiotherapy in patients with recurrent or progressive gliomas: The CINDERELLA trial. BMC Cancer 2010, 10, 533. [Google Scholar] [CrossRef]
  62. Combs, S.E.; Kieser, M.; Rieken, S.; Habermehl, D.; Jäkel, O.; Haberer, T.; Nikoghosyan, A.; Haselmann, R.; Unterberg, A.; Wick, W.; et al. Randomized phase II study evaluating a carbon ion boost applied after combined radiochemotherapy with temozolomide versus a proton boost after radiochemotherapy with temozolomide in patients with primary glioblastoma: The CLEOPATRA Trial. BMC Cancer 2010, 10, 478. [Google Scholar] [CrossRef]
  63. Sforazzini, F.; Salome, P.; Kudak, A.; Ulrich, M.; Bougatf, N.; Debus, J.; Knoll, M.; Abdollahi, A. pyCuRT: An Automated Data Curation Workflow for Radiotherapy Big Data Analysis using Pythons’ NyPipe. Int. J. Radiat. Oncol. Biol. Phys. 2020, 108, e772. [Google Scholar] [CrossRef]
  64. Tustison, N.J.; Avants, B.B.; Cook, P.A.; Zheng, Y.; Egan, A.; Yushkevich, P.A.; Gee, J.C. N4ITK: Improved N3 bias correction. IEEE Trans. Med. Imaging 2010, 29, 1310–1320. [Google Scholar] [CrossRef]
  65. Isensee, F.; Schell, M.; Pflueger, I.; Brugnara, G.; Bonekamp, D.; Neuberger, U.; Wick, A.; Schlemmer, H.-P.; Heiland, S.; Wick, W.; et al. Automated brain extraction of multisequence MRI using artificial neural networks. Hum. Brain Mapp. 2019, 40, 4952–4964. [Google Scholar] [CrossRef] [PubMed]
  66. Ebner, M.; Wang, G.; Li, W.; Aertsen, M.; Patel, P.A.; Aughwane, R.; Melbourne, A.; Doel, T.; Dymarkowski, S.; De Coppi, P.; et al. An automated framework for localization, segmentation and super-resolution reconstruction of fetal brain MRI. NeuroImage 2020, 206, 116324. [Google Scholar] [CrossRef] [PubMed]
  67. Avants, B.B.; Tustison, N.; Song, G. Advanced normalization tools (ANTS). Insight J. 2009, 2, 1–35. [Google Scholar]
  68. Reinhold, J.C.; Dewey, B.E.; Carass, A.; Prince, J.L. Evaluating the impact of intensity normalization on MR image synthesis. Proc. SPIE Int. Soc. Opt. Eng. 2019, 10949, 109493H. [Google Scholar] [PubMed]
  69. Zhang, Y.; Brady, M.; Smith, S. Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans. Med. Imaging 2001, 20, 45–57. [Google Scholar] [CrossRef]
  70. Reynolds, D.A. Gaussian Mixture Models. In Encyclopedia of Biometrics; Springer: Berlin/Heidelberg, Germany, 2009; Volume 741, pp. 659–663. [Google Scholar]
  71. Beer, J.C.; Tustison, N.J.; Cook, P.A.; Davatzikos, C.; Sheline, Y.I.; Shinohara, R.T.; Linn, K.A. Alzheimer’s Disease Neuroimaging Initiative Longitudinal Combat: A method for harmonizing longitudinal multi-scanner imaging data. Neuroimage 2020, 220, 117129. [Google Scholar] [CrossRef]
  72. Da-Ano, R.; Masson, I.; Lucia, F.; Doré, M.; Robin, P.; Alfieri, J.; Rousseau, C.; Mervoyer, A.; Reinhold, C.; Castelli, J.; et al. Performance comparison of modified Combat for harmonization of radiomic features for multicentre studies. Sci. Rep. 2020, 10, 10248. [Google Scholar] [CrossRef]
  73. Orlhac, F.; Lecler, A.; Savatovski, J.; Goya-Outi, J.; Nioche, C.; Charbonneau, F.; Ayache, N.; Frouin, F.; Duron, L.; Buvat, I. How can we Combat multicentre variability in MR radiomics? Validation of a correction procedure. Eur. Radiol. 2021, 31, 2272–2280. [Google Scholar] [CrossRef]
  74. Van Griethuysen, J.J.; Fedorov, A.; Parmar, C.; Hosny, A.; Aucoin, N.; Narayan, V.; Beets-Tan, R.G.; Fillion-Robin, J.-C.; Pieper, S.; Aerts, H.J. Computational radiomics system to decode the radiographic phenotype. Cancer Res. 2017, 77, e104–e107. [Google Scholar] [CrossRef]
  75. Lin, D.Y.; Wei, L.-J. The robust inference for the Cox proportional hazards model. J. Am. Stat. Assoc. 1989, 84, 1074–1078. [Google Scholar] [CrossRef]
  76. Frome, E.L. The analysis of rates using Poisson regression models. Biometrics 1983, 39, 665–674. [Google Scholar] [CrossRef]
  77. Campos, B.; Olsen, L.R.; Urup, T.; Poulsen, H.S. A comprehensive profile of recurrent glioblastoma. Oncogene 2016, 35, 5819–5825. [Google Scholar] [CrossRef]
  78. Boots-Sprenger, S.H.E.; Sijben, A.; Rijntjes, J.; Tops, B.B.J.; Idema, A.J.; Rivera, A.L.; Bleeker, F.E.; Gijtenbeek, A.M.; Diefes, K.; Heathcock, L.; et al. Significance of complete 1p/19q co-deletion, IDH1 mutation and MGMT promoter methylation in gliomas: Use with caution. Mod. Pathol. 2013, 26, 922–929. [Google Scholar] [CrossRef]
Figure 1. MR preprocessing diagram was applied to both cohorts. RT SS GTV represents the gross tumour volume segmentation extracted from the DICOM RT structure set. T2w-FL: T2w-FLAIR.
Figure 1. MR preprocessing diagram was applied to both cohorts. RT SS GTV represents the gross tumour volume segmentation extracted from the DICOM RT structure set. T2w-FL: T2w-FLAIR.
Cancers 15 00965 g001
Figure 2. Study workflow—following MR image preprocessing on cohorts C1 and C2, features were extracted from each normalized dataset, intensity normalization method-specific radiomics signatures were derived, and Cox proportional hazards and Poisson regression models were trained.
Figure 2. Study workflow—following MR image preprocessing on cohorts C1 and C2, features were extracted from each normalized dataset, intensity normalization method-specific radiomics signatures were derived, and Cox proportional hazards and Poisson regression models were trained.
Cancers 15 00965 g002
Figure 3. Scatter plots of the averaged (over the five bin counts considered) C-index vs AIC obtained by the CPH models for all four sequences in the study. Upper panel: cohort C1, Lower panel: cohort C2. csf: cerebrospinal fluid, wm: white matter, gm: grey matter, md: mode, gmm: Gaussian mixture models, kde: kernel density estimation, hm: Nyúl/Udupa histogram matching, ws: white stripe, nn: no normalization.
Figure 3. Scatter plots of the averaged (over the five bin counts considered) C-index vs AIC obtained by the CPH models for all four sequences in the study. Upper panel: cohort C1, Lower panel: cohort C2. csf: cerebrospinal fluid, wm: white matter, gm: grey matter, md: mode, gmm: Gaussian mixture models, kde: kernel density estimation, hm: Nyúl/Udupa histogram matching, ws: white stripe, nn: no normalization.
Cancers 15 00965 g003
Figure 4. Scatter plots of the averaged (over the five bin counts considered) mse vs AIC obtained by the POI models for all four sequences in the study. Upper panel: cohort C1, Lowe panel: cohort C2. csf: cerebrospinal fluid, wm: white matter, gm: grey matter, md: mode, gmm: Gaussian mixture models, kde: kernel density estimation, hm: Nyúl/Udupa histogram matching, ws: white stripe, nn: no normalization. Table 2 summarizes and ranks the performance scores of the intensity normalization methods for each of the four MR sequences considered in both cohorts.
Figure 4. Scatter plots of the averaged (over the five bin counts considered) mse vs AIC obtained by the POI models for all four sequences in the study. Upper panel: cohort C1, Lowe panel: cohort C2. csf: cerebrospinal fluid, wm: white matter, gm: grey matter, md: mode, gmm: Gaussian mixture models, kde: kernel density estimation, hm: Nyúl/Udupa histogram matching, ws: white stripe, nn: no normalization. Table 2 summarizes and ranks the performance scores of the intensity normalization methods for each of the four MR sequences considered in both cohorts.
Cancers 15 00965 g004
Figure 5. Significant feature correlation heatmaps between the 15 different normalization methods considered and the reference non-normalized dataset for T1wce images from cohort C1 discretized with a bin count of 32. Features with a high correlation (rs > 0.8) between at least 12 methods were further selected for modelling. Wv: Wavelet filter transformation, or: original image, log 3: Laplacian of gaussian transformation with a sigma of 3 mm.
Figure 5. Significant feature correlation heatmaps between the 15 different normalization methods considered and the reference non-normalized dataset for T1wce images from cohort C1 discretized with a bin count of 32. Features with a high correlation (rs > 0.8) between at least 12 methods were further selected for modelling. Wv: Wavelet filter transformation, or: original image, log 3: Laplacian of gaussian transformation with a sigma of 3 mm.
Cancers 15 00965 g005
Figure 6. Box plots of the top-ranked image normalization method evaluation metrics C-I and MSE before and after the elimination of the intensity normalization impacted significant features for cohorts C1 and C2 for each MR sequence. The average (across all bin counts) 10-CV C-I/MSE with 95% confidence intervals are plotted. Performance of both models was similarly affected after the elimination of the intensity normalization impacted significant features, with a mean decrease in the 10-CV C-I and 10-CV MSE of 0.05 and 0.03 in all four sequences across both cohorts C-I: Concordance-index, mse: mean squared error, T2w-FL: T2w-FLAIR.
Figure 6. Box plots of the top-ranked image normalization method evaluation metrics C-I and MSE before and after the elimination of the intensity normalization impacted significant features for cohorts C1 and C2 for each MR sequence. The average (across all bin counts) 10-CV C-I/MSE with 95% confidence intervals are plotted. Performance of both models was similarly affected after the elimination of the intensity normalization impacted significant features, with a mean decrease in the 10-CV C-I and 10-CV MSE of 0.05 and 0.03 in all four sequences across both cohorts C-I: Concordance-index, mse: mean squared error, T2w-FL: T2w-FLAIR.
Cancers 15 00965 g006
Table 1. Descriptions of the data cohorts C1 and C2 considered in this study.
Table 1. Descriptions of the data cohorts C1 and C2 considered in this study.
C1C2
n%n%
Patients197100141100
Gender
Male120618661
Female77395539
Age
<5084644733
50–69105537352
≥708172115
Tumour grade
III71363424
IV126646546
MR sequence
T1wce197100141100
T1w1869413596
T2w-FLAIR1688511883
T2w1417110071
Table 2. The number of shape, first and second-order statistics derived per sequence and calculated on both the original and derived images.
Table 2. The number of shape, first and second-order statistics derived per sequence and calculated on both the original and derived images.
ClassNo. Features
First-order statistics19
Shape-based (3D)16
Second-order statistics
Gray Level Co-occurrence Matrix24
Gray Level Run Length Matrix16
Gray Level Size Zone Matrix16
Neighboring Gray Tone Difference Matrix5
Gray Level Dependence Matrix14
Table 3. Ranking with scores of the intensity normalizations of the MR sequences for cohorts C1 and C2. T2w-FL: T2w-FLAIR, IN: Intensity normalization method.
Table 3. Ranking with scores of the intensity normalizations of the MR sequences for cohorts C1 and C2. T2w-FL: T2w-FLAIR, IN: Intensity normalization method.
C1T1wceT1wT2wT2w-FL
INScoreINScoreINScoreINScore
1ws0.71Combat0.13hm0.27wm-md0.02
2kde−0.13hm−0.28Combat−0.03wm-gm−0.11
3csf-gm−0.20csf-md−0.90z-score−0.28kde−0.13
4z-score−0.48nn−1.00gmm−0.38gm-md−0.23
5wm-gm−0.85z-score−1.14csf-gm−0.61gm−0.24
6csf−0.97csf-gm−1.58kde−0.71wm−0.42
7hm−1.04wm-csf−1.65nn−0.76csf-gm−0.46
8gmm−1.11wm−1.85csf−0.78Combat−0.77
9gm−1.13kde−1.88wm-md−0.80csf-md−0.77
10wm−1.24wm-md−1.95gm-md−0.96hm−0.80
11wm-md−1.67gm-md−2.05csf-md−1.09wm-csf−1.01
12csf-md−1.71ws−2.15ws−1.18gmm−1.02
13gm-md−1.72csf−2.16wm−1.22ws−1.29
14wm-csf−2.16gm−2.23wm-gm−1.72csf−1.75
15Combat−2.25wm-gm−2.37gm−1.79z-score−2.21
16nn−2.27gmm−2.48wm-csf−2.01nn−2.65
C2
1ws1.00z-score0.64Combat0.07wm-csf0.66
2csf−0.54hm−0.11hm−0.09wm-md−0.32
3hm−0.73csf−0.34gm-md−0.21gmm−0.56
4z-score−0.76gmm−0.35wm-csf−0.24kde−0.63
5gm−0.77csf-md−0.81gmm−0.41csf-gm−0.71
6wm−0.87kde−0.93wm-md−0.78wm−0.72
7wm-gm−0.87gm-md−0.97gm−1.00gm−0.76
8csf-gm−0.96csf-gm−0.97csf-md−1.12hm−0.81
9kde−0.98ws−1.04ws−1.13gm-md−0.90
10wm-csf−1.07gm−1.18z-score−1.21csf-md−1.05
11gmm−1.10Combat−1.20csf−1.31nn−1.25
12wm-md−1.13nn−1.41kde−1.36csf−1.35
13Combat−1.19wm-csf−1.43wm−1.52Combat−1.42
14gm-md−1.28wm-md−1.64nn−1.60ws−1.50
15csf-md−1.39wm-gm−2.01wm-gm−1.69wm-gm−1.59
16nn−1.82wm−2.11csf-gm−1.81z-score−2.11
Table 4. Performance of the top-ranked image normalization method separate before and after the elimination of the intensity normalization impacted significant features for cohort C1 and C2 for each MR sequence. The average (across all bin counts) 10-CV C-I/MSE with 95% confidence intervals are reported.
Table 4. Performance of the top-ranked image normalization method separate before and after the elimination of the intensity normalization impacted significant features for cohort C1 and C2 for each MR sequence. The average (across all bin counts) 10-CV C-I/MSE with 95% confidence intervals are reported.
C1C2
BeforeAfterBeforeAfter
T1wce0.71 [0.69 0.74]/
0.21 [0.19 0.23]
0.65 [0.63 0.69]/
0.23 [0.21 0.25]
0.65 [0.62 0.67]/
0.15 [0.13 0.17]
0.62 [0.60 0.65]/
0.19 [0.17 0.21]
T1w0.68 [0.64 0.70]/
0.22 [0.20 0.25]
0.63 [0.61 0.67]/
0.24 [0.22 0.26]
0.65 [0.61 0.69]/
0.15 [0.12 0.18]
0.62 [0.58 0.65]/
0.18 [0.15 0.20]
T2w0.65 [0.62 0.67]/
0.22 [0.19 0.25]
0.63 [0.60 0.67]/
0.25 [0.22 0.28]
0.67 [0.64 0.69]/
0.13 [0.11 0.17]
0.60 [0.58 0.65]/
0.16 [0.14 0.20]
T2w-FL0.67 [0.64 0.69]/
0.20 [0.18 0.23]
0.62 [0.59 0.67]/
0.23 [0.21 0.25]
0.72 [0.65 0.76]/
0.18 [0.15 0.21]
0.66 [0.64 0.69]/
0.20 [0.17 0.22]
Table 5. Performance of the top-ranked image normalization method separate and in combination with the feature-based method Combat for cohort C1 and C2 for each MR sequence. The average (across all bin counts) 10-CV C-I/MSE with 95% confidence intervals are reported.
Table 5. Performance of the top-ranked image normalization method separate and in combination with the feature-based method Combat for cohort C1 and C2 for each MR sequence. The average (across all bin counts) 10-CV C-I/MSE with 95% confidence intervals are reported.
C1 C2
CombatI. Norm.CombinedCombatI. Norm.Combined
T1wce0.68 [0.66 0.70]/
0.21 [0.19 0.23]
0.71 [0.690.74]/
0.21 [0.19 0.23]
0.68 [0.66 0.69]/
0.21 [0.19 0.23]
0.64 [0.62 0.68]/
0.15 [0.13 0.17]
0.65 [0.62 0.67]/
0.15 [0.13 0.17]
0.63 [0.61 0.66]/
0.17 [0.15 0.19]
T1w0.68 [0.64 0.70]/
0.22 [0.20 0.24]
0.66 [0.64 0.68]/
0.22 [0.19 0.24]
0.62 [0.59 0.64]/
0.23 [0.20 0.26]
0.62 [0.60 0.66]/
0.15 [0.12 0.17]
0.65 [0.61 0.69]/
0.15 [0.12 0.18]
0.62 [0.59 0.65]/
0.15 [0.11 0.16]
T2w0.62 [0.59 0.64]/
0.23 [0.21 0.23]
0.65 [0.62 0.67]/
0.22 [0.19 0.25]
0.61 [0.58 0.63]/
0.25 [0.23 0.27]
0.67 [0.64 0.69]/
0.13 [0.11 0.17]
0.67 [0.64 0.69]/
0.13 [0.11 0.15]
0.62 [0.59 0.65]/
0.15 [0.13 0.19]
T2w-FL0.67 [0.64 0.69]/
0.21 [0.19 0.24]
0.67 [0.64 0.69]/
0.20 [0.18 0.23]
0.64 [0.61 0.66]/
0.24 [0.22 0.26]
0.70 [0.67 0.72]/
0.16 [ 0.14 0.19]
0.72 [0.65 0.76]/
0.14 [0.12 0.17]
0.68 [0.65 0.70]/
0.17 [0.15 0.21]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Salome, P.; Sforazzini, F.; Brugnara, G.; Kudak, A.; Dostal, M.; Herold-Mende, C.; Heiland, S.; Debus, J.; Abdollahi, A.; Knoll, M. MR Intensity Normalization Methods Impact Sequence Specific Radiomics Prognostic Model Performance in Primary and Recurrent High-Grade Glioma. Cancers 2023, 15, 965. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers15030965

AMA Style

Salome P, Sforazzini F, Brugnara G, Kudak A, Dostal M, Herold-Mende C, Heiland S, Debus J, Abdollahi A, Knoll M. MR Intensity Normalization Methods Impact Sequence Specific Radiomics Prognostic Model Performance in Primary and Recurrent High-Grade Glioma. Cancers. 2023; 15(3):965. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers15030965

Chicago/Turabian Style

Salome, Patrick, Francesco Sforazzini, Gianluca Brugnara, Andreas Kudak, Matthias Dostal, Christel Herold-Mende, Sabine Heiland, Jürgen Debus, Amir Abdollahi, and Maximilian Knoll. 2023. "MR Intensity Normalization Methods Impact Sequence Specific Radiomics Prognostic Model Performance in Primary and Recurrent High-Grade Glioma" Cancers 15, no. 3: 965. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers15030965

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop