# Influence of magnetic field strength on magnetic resonance imaging radiomics features in brain imaging, an in vitro and in vivo study

Samy Ammari, Stéphanie Pitre-Champagnat, Laurent Dercle, Emilie Chouzenoux, Salma Moalla, Sylvain Reuze, Hugues Talbot, Tite Mokoyoko, Joya Hadchiti, Sebastien Diffetocq, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Samy Ammari, Stéphanie Pitre-Champagnat, Laurent Dercle, Emilie Chouzenoux, Salma Moalla, et al. Influence of magnetic field strength on magnetic resonance imaging radiomics features in brain imaging, an in vitro and in vivo study. Frontiers in Oncology, 2021, 10, pp.541663. 10.3389/fonc.2020.541663. hal-03117688

HAL Id: hal-03117688

https://hal.science/hal-03117688

Submitted on 21 Jan 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Influence of Magnetic Field Strength on Magnetic Resonance Imaging Radiomics Features in Brain Imaging, an In Vitro and In Vivo Study

Samy Ammari¹,²*, Stephanie Pitre-Champagnat², Laurent Dercle¹,³,⁴,
Emilie Chouzenoux⁵, Salma Moalla¹, Sylvain Reuze⁶, Hugues Talbot⁵, Tite Mokoyoko¹,
Joya Hadchiti¹, Sebastien Diffetocq¹, Andreas Volk², Mickeal El Haik¹, Sara Lakiss¹,
Corinne Balleyguier¹,², Nathalie Lassau¹,² and Francois Bidault¹,²

OPEN ACCESS

**Edited by:**
Freimut Dankwart Juengling,
Universität Bern, Switzerland

**Reviewed by:**
Prathyush Chirra,
Case Western Reserve University,
United States
Alberto Traverso,
Maastro Clinic, Netherlands

*Correspondence:
Samy Ammari
samy.ammari@gustaveroussy.fr

**Specialty section:**
This article was submitted to
Cancer Imaging and
Image-directed Interventions,
a section of the journal
Frontiers in Oncology

**Received:** 10 March 2020
**Accepted:** 23 November 2020
**Published:** 20 January 2021

**Citation:**
Ammari S, Pitre-Champagnat S, Dercle L, Chouzenoux E, Moalla S, Reuze S, Talbot H, Mokoyoko T, Hadchiti J, Diffetocq S, Volk A, El Haik M, Lakiss S, Balleyguier C, Lassau N and Bidault F (2021) Influence of Magnetic Field Strength on Magnetic Resonance Imaging Radiomics Features in Brain Imaging, an In Vitro and In Vivo Study. *Front. Oncol.* 10:541663.
doi: 10.3389/fonc.2020.541663

¹ Department of Radiology, Gustave Roussy Cancer Campus, Université Paris-Saclay, Villejuif, France, ² BioMaps (UMR1281), Université Paris-Saclay, CNRS, INSERM, CEA, Orsay and Gustave Roussy, Villejuif, France, ³ Immunology of Tumours and Immunotherapy INSERM U1015, Gustave Roussy Cancer Campus, Université Paris Saclay, Villejuif, France, ⁴ Radiology Department, Columbia University Medical Center, New York Presbyterian Hospital, New York, NY, United States,

⁵ Center for Visual Computing, CentraleSupelec, Inria, Université Paris-Saclay, Gif-sur-Yvette, France, ⁶ Department of Radiotherapy - Medical Physics, Gustave Roussy, Université ParisSaclay, Villejuif, France

**Background:** The development and clinical adoption of quantitative imaging biomarkers (radiomics) has established the need for the identification of parameters altering radiomics reproducibility. The aim of this study was to assess the impact of magnetic field strength on magnetic resonance imaging (MRI) radiomics features in neuroradiology clinical practice.

**Methods:** T1 3D SPGR sequence was acquired on two phantoms and 10 healthy volunteers with two clinical MR devices from the same manufacturer using two different magnetic fields (1.5 and 3T). Phantoms varied in terms of gadolinium concentrations and textural heterogeneity. 27 regions of interest were segmented (phantom: 21, volunteers: 6) using the LIFEX software. 34 features were analyzed.

**Results:** In the phantom dataset, 10 (67%) out of 15 radiomics features were significantly different when measured at 1.5T or 3T (student's t-test, p < 0.05). Gray levels resampling, and pixel size also influence part of texture features. These findings were validated in healthy volunteers.

**Conclusions:** According to daily used protocols for clinical examinations, radiomic features extracted on 1.5T should not be used interchangeably with 3T when evaluating texture features. Such confounding factor should be adjusted when adapting the results of a study to a different platform, or when designing a multicentric trial.

**Keywords:** tissue features, heterogeneous phantom, homogeneous phantom, magnetic fields, texture, magnetic resonance imaging
## HIGHLIGHTS

1. Radiomic features at 1.5T are not interchangeable with 3T when evaluating tumor texture

2. Field strength should be taken into account in the interpretation of the texture indices

3. Signal to noise ratio should be taken into account in the interpretation of the texture indices

## INTRODUCTION

Radiomics is a fast growing discipline, which is undergoing growing interest in computational medical imaging (1). This field of medical study aims at extracting a large amount of quantitative features from medical images using data-characterization algorithms. This research field is faced with multiple challenges (2). Radiomics is used in oncology to analyze features that are invisible to the naked eye and that may be associated with gene expression, tumor histology, treatment response and patient outcome (3).

MRI has several advantages and disadvantages for radiomics analysis (4-7). Among imaging modalities, it offers the best soft tissue contrast. Conversely, differences in MRI parameters (field strength, gradient characteristics), image acquisition protocols (8, 9), sequences, pixel size (10), and the signal-to-noise ratio (SNR) (11, 12) might impact radiomics features that are sensitive to image quality.

The impact of MRI acquisition and processing on radiomics reproducibility is scarcely reported. As in nuclear medicine, the validation of a biologically relevant and reproducible clinical biomarker based on radiomics implies standardization of protocols across several centers (13). For example, 3T compared with 1.5T MRI not only provides higher SNR, allows increased image resolution, and modifies relaxation times T1 and T2 but also induces some artifacts. Thus, there is a clear need to evaluate the influence of field strength and related settings like image resolution (pixel size, field of view [FOV], and matrix) on radiomic features.

Recent articles using radiomics as biomarkers consider that the major challenge is that grayscale MRI intensities, contrary to X-ray CT, are not standardized and are highly dependent on manufacturer, sequence type and acquisition parameters (14, 15).

To address this problem, authors focused on image pre-processing techniques that effectively minimize MR intensity inhomogeneity in a tissue region (16-18), spatial resampling (17-20) and brain arch extraction prior to image intensity normalisation (21, 22).

Although several studies have shown variability in texture analysis as a function of MRI acquisition parameters and gray

level discretization steps, none of them have evaluated the combined impact of magnetic field strength, matrix size, pixel size, intensity normalization, and gray level discretization pre-processing methods on MRI radiomic feature values (9, 10, 23).

This study was designed to evaluate the impact of the field strength (1.5T vs. 3T) on radiomic features. Two clinical 1.5 and 3T MR devices from the same manufacturer were used to image the same phantoms, mimicking homogeneous and heterogeneous tissues. We also imaged healthy volunteers with the same sequences acquisition and image processing parameters.

It is believed that this step is crucial before proposing recommendations to standardize brain MRI pre-processing techniques, which is in turn essential for guaranteeing reliable radiomics-based models with big data and artificial intelligence.

## MATERIAL AND METHODS

This prospective study was approved by the local institution review board.

### Homogeneous Phantoms

Our homogeneous phantom was designed to mimic cerebrospinal fluid and opacified blood vessels. It was defined with different gadolinium concentrations: eight 30-ml tubes were filled with demineralized water mixed with increasing gadolinium chelate concentration (Gadoteric Acid; Dotarem®, Guerbet): 0.25, 0.5, 0.75, 1, 1.25, 1.5, 1.75, and 2 mmol/l (Figures 1A, B). Those eight homogeneous tubes were designated C1 to C8.

### Heterogeneous Phantoms

Our heterogeneous phantom, designed to mimic the brain white matter, was composed of agarose gel coating polystyrene beads. The proton density of agarose has similar characteristics and relaxation times to biological tissues (10). We defined six heterogeneous tubes and two homogeneous tubes. The 30-ml tubes were filled with polystyrene beads of different diameters (1, 2 and 3 mm) in an agarose gel solution (2%), either pure or mixed with 0.25 mmol/l Gadolinium chelate (Figures 1C, D) (9). The tubes were displayed in the tube slots of the Eurospin phantom (A Eurospin II-(TO5) phantom; Diagnostic Sonar) glass cylinder, filled with 1% copper sulfate (24). The eight tubes were designated T1 to T8.

### Healthy Volunteers

Clinical MRI sequences were performed in 10 healthy volunteers aged 21-26 years (six men, four women). Two MRI acquisitions were performed on each MRI device within a 40-min time interval.

### MRI Devices and Protocols

The influence of field strength was tested on two different MRI devices from the same manufacturer (General Electric): an Optima MR450w 1.5T superconducting magnet MRI installed in 2016 with a 70 cm tunnel, 32 channels, 50 cm FOV (Z axis), and gradients 40 mT SR 200 mT/m/s, and a Discovery MR750w 3T superconducting magnet MRI installed in 2012, with a 70 cm tunnel, 32 channels, 50 cm FOV (Z axis), and gradients 44 mT/m SR 200 T/m/s.

<footnote>**Abbreviations:** 3D SPGR, 3D rapid gradient echo sequence; EORTC, European Organization for Research and Treatment of Cancer; TR, Repeat time; MTX, Image matrix; CoMat, The co-occurrence matrix; GRLM, The gray-level run length matrix; GZLM, The gray-level zone length matrix; SBW, Sampling bandwidth.</footnote>
<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | (A) Homogeneous phantoms: eight 30-ml tubes filled with demineralized water mixed with increasing gadolinium chelate concentration (Gadoteric Acid; Dotarem®, Guerbet): 0.25, 0.5, 0.75, 1, 1.25, 1.5, 1.75 and 2 mmol/l and a central demineralized water tube. (B) T1 weighted MRI image of the homogeneous phantom (256x256 matrix, 18 cm FOV). (C, D) T1 weighted MRI image of the heterogeneous phantom (256;256 matrix, 24 cm FOV): in the two central columns are six 30-ml tubes filled with polystyrene beads of different diameters (1, 2 and 3 mm) in an agarose gel solution (2%), either pure or mixed with 0.25 mmol/l gadolinium chelate, and two 30-ml tubes filled with agarose gel solution (2%) and two different concentration of gadolinium chelate.</figcaption></figure>

We used 'head and neck' coils with 32 channels with a 35 cm diameter adapted to the frequency of each MRI.

## Acquisition

MRI imaging acquired was based on a T1-weighted 3D rapid gradient echo sequence (3D SPGR). This sequence is used in clinical imaging for rapid volumetric imaging and can be acquired before or after contrast agent injection. The same sequences parameters (matrix, FOV, plane) were used for both MRIs. Parameters were chosen as recommended by the European Organization for Research and Treatment of Cancer (EORTC) to explore brain tumors: repeat time (TR) 6.1 ms; echo time (TE) inphase 1.2-2.1 ms; NEX1; thickness 1 mm in contiguous sections and bandwidth at 31 with 165 slices; 2-5-min acquisition depending on the FOV and matrices. The temperature was maintained between 19 and 21°C in each MRI during acquisition. For the healthy volunteers, the same 3DT1 SPGR sequence was used covering the whole encephalon.

## Acquisition Parameters: Field of View and Matrix Size

For the phantom study, five couples of matrices and FOV, determining various pixel size, were applied on both machines (Table 1). The couples were chosen to be suitable with clinical acquisition. In healthy volunteers only one acquisition per device was performed, the FOV was fixed at 24 cm with an image matrix of 256 x 256 on both machines.

## Texture Analysis

Segmentations and texture features calculations were performed with the freely available LIFEX software package (http://www.lifexsoft.org) (25). A 6 cm³ VOI was placed in each homogeneous and heterogeneous phantom's tube. VOI position along Z axis of tubes was controlled according to the reference markers located on the phantom surface. For healthy volunteers, VOI were displayed in six bilateral areas: corpus callosum, central gray nuclei, and white matter of the centrum semiovale.
<table><caption>**TABLE 1** | Matrices and FOVs studied at 1.5 and 3T with the corresponding pixel size.</caption><thead><tr><td>MTX (pixels)</td><td>Field of view (cm)</td><td>Pixel size (mm)</td></tr></thead><tbody><tr><td>256 x 256</td><td>24</td><td>0.938</td></tr><tr><td>256 x 256</td><td>18</td><td>0.703</td></tr><tr><td>256 x 256</td><td>12</td><td>0.468</td></tr><tr><td>256 x 128</td><td>24</td><td>0.937 x 1.875</td></tr><tr><td>128 x 128</td><td>24</td><td>1.875</td></tr></tbody></table>

In each VOI, 38 texture features (indices) were extracted from first order and second-order statistics, arising from the analysis of the intensity histogram and the calculation of three texture matrices: the co-occurrence matrix (CM), the gray-level run length matrix (GRLM) and the gray-level zone length matrix (GZLM) (26). The texture analysis software needs to resample the 16-bit gray levels of MRI images in order to calculate these features. Voxel intensities were resampled in three different ways, using 256, 128 and 64 discrete values.

## Calculation of Signal-to-Noise Ratio

Mean SNR values were calculated on both MRIs by measuring the ratio between the mean signal intensity in each pure agarose sample and the standard deviation (SD) of the background noise selected in the frequency encoding direction. Calculation of SNR in healthy volunteers was also performed by measuring the mean signal intensity in each white substance and the SD of the background noise selected in the frequency encoding direction.

## Statistical Methods

Two-by-two correlation texture features were calculated. The statistical significant differences between 1.5 and 3T were determined by the Student's t-test using paired data. In the event of a lack of normal distribution (N < 30), results were also presented as boxplots to demonstrate the trend of texture values on 1.5T vs. 3T. The association between acquisition protocols (magnetic field, matrix, FOV), textural changes in the phantom (homogeneous vs. heterogeneous and polystyrene size), and imaging features, was calculated with Spearman's rho correlation coefficients with P-values corrected for multiple tests. Statistical analyses were performed using SPSS v24.0. The datasets generated and analyzed are available from the corresponding author.

# RESULTS

## Phantom Study

### Influence of Software Resampling Step

The analysis software needs to resample the images gray scale by sub-sampling the number gray levels, in order to be able to calculate texture features. Thus, we firstly explored the influence of this sub-sampling step before exploring the influence of magnetic field strength. We did this through the two-by-two coefficient correlation calculation between the texture features (Figure 2). A colored (red/blue) hierarchical clustering is presented, showing the importance of the correlation between coefficients (Figure 2) for each level of resampling. A strong

correlation was observed for a large number of texture features at the 256 gray level (only 9 features were totally independent with a Spearman's rho <0,70). The correlation value decreased with the sub-sampling in gray levels resolution. The gray level sub-sampling influence on texture features was not significant for parameters extracted from the histogram. In contrast, the gray levels sub-sampling had a significant impact on parameters extracted from matrices.

## Influence of Field Strength on Homogeneous Phantom

The majority of texture features values were significantly different between the two magnetic fields (1.5T vs. 3T). For example the mean value is presented for comparison between the two fields strenght (Figure 3). The entropy mean value was higher on 3T versus 1.5T by a factor of four (Figure 3).

## Influence of Field Strength on Heterogeneous Phantoms

The majority of texture features were significantly different between the two magnetic fields (1.5T vs. 3T) according to the Student's paired t-test (Table 2). Concerning features from histogram analysis, only Kurtosis, Entropy, and Energy did not significantly differ between 1.5 and 3T. For the matrix-based texture features, only LZHGE shows no significant difference.

## Influence of Pixel Size

The pixels size, dependant on both matrix size and FOV size, altered the radiomics output in homogeneous phantom according to three behaviors. First of all, pixel size altered the values of the texture features on both the 1.5 and 3T magnetic field, for 15 indices: stdvalue, skewness, Homogeneity, Contrast, Correlation, Entropy, Dissimilarity, ZLNU, GLNU, RLNU, Coarseness, SZE, LZE, GLNU_1, and ZP. Figure 4A illustrates the "Correlation" named feature. Second, pixel size altered the value of texture features only on 1.5T, but not on 3T: Kurtosis, Entropy H, Energy H, Energy, LRE, RP. Figure 4B illustrates "Energy" named feature. Note that in contrast "Contrast_1" named feature value was isolatedly modified only on 3T. Thirdly, pixel size did not significantly alter the value of texture features for: minvalue, meanvalue, maxvalu, HGRE, SRHGE, LRHGE, HGZE, SZHGE, and LZHE.

The impact of matrix size and FOV as well as pixel size, on the radiomics output was also studied with the heterogeneous phantoms. Also three behaviors were observed. First, pixel size altered the values of the texture features on the two magnetic fields for: RLNU, Coarseness, SZLGE, GLNU and ZLNU. Figure 5A shows the example of "GLNU" named feature. Second, pixel size altered the value of texture features only on 1.5T, but not on 3T for: Correlation, LGRE, SRLGE, GLNU, and LGZE. Figure 5B shows an example with the "SRLGE" named feature. Third, pixel size did not significantly alter the value of texture features for the other 29 texture parameters.

Ability of texture features to identify a difference between phantom tubes, in both field strength separately, regardless of the pixel size, is observed through the absolute Spearman's rho correlation coefficient (Table 3). Eight imaging features
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**FIGURE 2** | Graphs showing influence of software resampling step on phantom study. Each graph represents the two-by-two coefficient correlation calculation between the texture parameters at one resampling level. The difference in resulting patterns shows the influence of the software resampling step on texture parameters calculation, particularly on parameters extracted from co occurrence matrix. (A) 256 gray levels resampling. (B) 128 gray levels resampling. (C) 64 gray levels resampling.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**FIGURE 3** | Difference of phantoms textures features values at 1.5 and 3T. Example of the mean value. Each colored dot groups shows the mean value calculated from repeated phantom MRI acquisition (one dot per MRI acquisition) in different homogeneous tubes with different concentration of Gadolinium chelate (on color per tube), respectively at 1.5 and 3T.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**FIGURE 4** | *In vitro* parameters value according to magnetic field strength and pixel size: For correlation (A) and energy (B). Texture features extracted to the homogeneous phantoms. Each colored dot groups and related box plot shows the value of correlation (A) and energy (B) calculated from repeated phantom MRI acquisition (one dot per MRI acquisition) with different pixel size (one color per matrix/FOV couple), respectively, at 1.5 and 3T.</figcaption></figure>

identified the (visually obvious) difference between homogeneous and heterogeneous phantoms, with an absolute Spearman's rho correlation coefficient above 0.5. The following texture features identified the difference between the different heterogeneous phantoms (mainly different according to their polystyrene beads size and spatial distribution) in both field strength separately:

dissimilarity, LZHGE, entropy, homogeneity, SRE, SZE, LRE, GLNU, RLNU, ZLNU, KurtosisH, LZE, entropyH, coarseness, SRLGE, SRLGE, SZLGE, and LGRE. Two texture features (coarseness, RLNU) identified the difference between homogeneous tubes (mainly different by their Gadolinium chelate concentration). As long as the pixel size remains sufficiently tiny
<table><caption>**TABLE 2** | Texture features presenting significant or non-significant differences at 1.5 and 3T based on the heterogeneous phantoms examination (student's paired t-test).</caption><tr><td>Significant difference in texture features value at 1.5 and 3T (p-values &lt; 10<sup>-5</sup>)</td><td>LGRE, Homogeneity, SZE, ZP, Contrast, LGZE, Correlation, RLNU Entropy, SRLGE, LRHGE, HGRE, mean value, SRHGE, LRE, GLNU, HGZE, GLNU, SRE, SZHGE, SZLGE, LRLGE, max value, Dissimilarity Contrast, RP, std value, LZE, ZLNU, skewness, H min, value Coarseness, LZLGE</td></tr><tr><td>Non-significant difference in texture features value at 1.5 and 3Tp-values &gt; 0.005</td><td>Energy H, Entropy H, Kurtosis H, LZHGE, Energy H</td></tr></table>

small (1 mm or less in our study), 26 features were sensitive to texture alteration.

## Effect of Field Strength on MRI Texture Features Values in Healthy Volunteers

Among the 38 parameters, significant differences were observed with 15 texture features measured in healthy volunteers between 1.5 and 3T MRIs for the same anatomical region (Figure 6). With constant fields (1.5T vs. 1.5T and 3T vs. 3T), those 15 texture features values appeared identical when measured in symmetrical anatomical structures (eg corpus callosum or caudate nucleus).

## DISCUSSION

Our study demonstrates that field strength (1.5T vs. 3T) influences numerous texture features values for both phantom and human exploration. Our study also demonstrated that gray-level resampling and pixel size influence some texture features. Those results are of importance for clinical study design and for patient examination and follow-up. In addition to that, our study identifies texture parameters that are able to differentiate phantom textures in our setup. It is important to mention that these relevant textures parameters are significantly influenced by field strength.

An increasing number of radiomics or deep-learning studies use imaging based on MRI acquisitions (24–28). Confounding parameters that can alter radiomic output have been largely addressed for CT scans but less for MRIs (27, 29–31). Our study aims to complete knowledge in the field of radiomics applied to neuroradiology MRI acquisitions.

The presented study identifies the importance to consider the impact of magnetic field strength on radiomics features values in clinical practice. SNR is well-known to increase when the field strength increases. We hypothesize that our result may be partially due to this difference in SNR between the two MRI machines. Even given that field strength is a major factor

<figure><img src="image_10.png" /><figcaption>**FIGURE 5** | *In vitro* parameters value according to magnetic field strength and pixel size: for GLNU_1 (A) and SRLGE (B). Texture features extracted to the heterogeneous phantoms. Each colored dot groups and related box plot shows the value of GLNU_1 (A) and SRLGE (B) calculated from repeated phantom MRI acquisition (one dot per MRI acquisition) with different pixel size (one color per matrix/FOV couple), respectively, at 1.5 and 3T.</figcaption></figure>
<table><caption>TABLE 3 | Tabulation of textures features and their ability to dissociate the different phantom's tubes at both field strengths separately, regardless to pixel size.</caption><thead><tr><th>Heterogeneous vs. Homogeneous media</th><th></th><th>Size of polystyrene</th><th></th><th>Gadolinium concentration</th><th></th></tr></thead><tbody><tr><td>RLNU</td><td>0.796</td><td>Contrast</td><td>0.801</td><td>Coarseness</td><td>0.577</td></tr><tr><td>Coarseness</td><td>0.795</td><td>Contrast_1</td><td>0.795</td><td>RLNU</td><td>0.513</td></tr><tr><td>minValue</td><td>0.683</td><td>Dissimilarity</td><td>0.794</td><td>minValue</td><td>0.482</td></tr><tr><td>GLNU 1</td><td>0.617</td><td>LZHGE</td><td>0.789</td><td>GLNU</td><td>0.414</td></tr><tr><td>SZLGE</td><td>0.552</td><td>stdValue</td><td>0.786</td><td>SZLGE</td><td>0.371</td></tr><tr><td>SRLGE</td><td>0.534</td><td>Entropy</td><td>0.783</td><td>SRLGE</td><td>0.354</td></tr><tr><td>Entropy</td><td>0.645</td><td>Homogeneity</td><td>0.773</td><td>GLNU</td><td>0.347</td></tr><tr><td>Homogeneity</td><td>0.622</td><td>ZP</td><td>0.77</td><td>EntropyH</td><td>0.33</td></tr><tr><td>ZLNU</td><td>0.491</td><td>SRE</td><td>0.768</td><td>meanValue</td><td>0.328</td></tr><tr><td>EntropyH</td><td>0.484</td><td>RP</td><td>0.762</td><td>LGZE</td><td>0.325</td></tr><tr><td>LGRE</td><td>0.481</td><td>Energy</td><td>0.756</td><td>LGRE</td><td>0.315</td></tr><tr><td>LGZE</td><td>0.477</td><td>SZE</td><td>0.756</td><td>EnergyH</td><td>0.308</td></tr><tr><td>GLNU</td><td>0.468</td><td>LRE</td><td>0.741</td><td>KurtosisH</td><td>0.302</td></tr><tr><td>Contrast</td><td>0.458</td><td>GLNU</td><td>0.738</td><td>SkewnessH</td><td>0.293</td></tr><tr><td>EnergyH</td><td>0.439</td><td>RLNU</td><td>0.732</td><td>Correlation</td><td>0.262</td></tr><tr><td>Contrast</td><td>0.429</td><td>minValue</td><td>0.727</td><td>maxValue</td><td>0.246</td></tr><tr><td>Dissimilarity</td><td>0.423</td><td>ZLNU</td><td>0.698</td><td>Contrast</td><td>0.246</td></tr><tr><td>stdValue</td><td>0.421</td><td>KurtosisH</td><td>0.677</td><td>ZLNU</td><td>0.237</td></tr><tr><td>Energy</td><td>0.418</td><td>LZE</td><td>0.646</td><td>LRLGE</td><td>0.228</td></tr><tr><td>meanValue</td><td>0.409</td><td>EntropyH</td><td>0.624</td><td>LRHGE</td><td>0.198</td></tr><tr><td>KurtosisH</td><td>0.405</td><td>Coarseness</td><td>0.617</td><td>Entropy</td><td>0.195</td></tr><tr><td>RP</td><td>0.378</td><td>SRLGE</td><td>0.583</td><td>Energy</td><td>0.182</td></tr><tr><td>SRE</td><td>0.373</td><td>EnergyH</td><td>0.554</td><td>Contrast</td><td>0.181</td></tr><tr><td>Correlation</td><td>0.363</td><td>SZLGE</td><td>0.554</td><td>Dissimilarity</td><td>0.18</td></tr><tr><td>SZE</td><td>0.346</td><td>LGRE</td><td>0.524</td><td>LZLGE</td><td>0.173</td></tr><tr><td>LRE</td><td>0.329</td><td>LGZE</td><td>0.46</td><td>Homogeneity</td><td>0.165</td></tr><tr><td>ZP</td><td>0.329</td><td>LRHGE</td><td>0.373</td><td>stdValue</td><td>0.155</td></tr><tr><td>LRHGE</td><td>0.316</td><td>LZLGE</td><td>0.321</td><td>HGZE</td><td>0.15</td></tr><tr><td>LRLGE</td><td>0.3</td><td>meanValue</td><td>0.275</td><td>HGRE</td><td>0.145</td></tr><tr><td>SkewnessH</td><td>0.229</td><td>SkewnessH</td><td>0.267</td><td>RP</td><td>0.141</td></tr><tr><td>maxValue</td><td>0.228</td><td>LRLGE</td><td>0.244</td><td>SZHGE</td><td>0.137</td></tr><tr><td>LZHGE</td><td>0.21</td><td>maxValue</td><td>0.112</td><td>SRE</td><td>0.132</td></tr><tr><td>HGZE</td><td>0.185</td><td>HGRE</td><td>0.101</td><td>SRHGE</td><td>0.132</td></tr><tr><td>HGRE</td><td>0.174</td><td>GLNU</td><td>0.092</td><td>SZE</td><td>0.132</td></tr><tr><td>SRHGE</td><td>0.14</td><td>HGZE</td><td>0.087</td><td>LRE</td><td>0.111</td></tr><tr><td>SZHGE</td><td>0.132</td><td>SRHGE</td><td>0.041</td><td>ZP</td><td>0.104</td></tr><tr><td>LZLGE</td><td>0.097</td><td>SZHGE</td><td>0.037</td><td>LZE</td><td>0.077</td></tr><tr><td>LZE</td><td>0.096</td><td>Correlation</td><td>0.017</td><td>LZHGE</td><td>0.01</td></tr></tbody></table>

Given are the names of the feature variables and the absolute Spearman's rho correlation coefficient between imaging features and the type of phantom tube. First column: ability to dissociate homogeneous versus heterogeneous tubes. Second column: ability to dissociate heterogeneous tubes according to their polystyrene beads sizes and their spatial distribution. Third column: ability to dissociate homogeneous tubes according to their Gadolinium chelate concentration.

influencing signal intensity, further factors modifying signal have to be considered, as discussed in a recent literature review. Mayerhofer et al. previously demonstrated that MRI texture features are influenced by change in SNR (5, 9). The difference in SNR could be explained by magnetic field strength, but SNR is also linked to the entire signal acquisition system (coils, electronic device etc.). The European research project COST B11 (31, 32), aimed at developing quantitative methods for MRI imaging features extraction (between 1998 and 2002), demonstrated that signal and subsequently texture features were dependent on parameters such as configuration of the transmitting and/or receiving RF coil (antenna), and the number of active segments in the coil due to change in tilt angle.

Our results showed that most texture features are also sensitive to the variations in matrix and FOV and consequently spatial resolution, even if the clinical ranges in matrix size and FOV for a dedicated organ may be thinner than those investigated in literature physics studies. Of note, we acquired images at different pixel sizes whereas most studies

changed the pixel size during post-processing. Our study showed that as long as the pixel size remains sufficiently small (1 mm in our experience), 26 features are sensitive to texture alteration. This is in accordance with Jirák study, which also compared MRI features measured on phantoms with nodular patterns in a multicenter settings and showed that texture classification was influenced by low resolution causing large errors (10).

For a defined magnetic field strength, our study showed that part of the texture features differentiated the heterogeneous phantoms from the homogeneous phantoms on all sequences, irrespective of the pixel size, and part of radiomics features differentiated phantoms with large and small spherical objects (polystyrene beads) scattered in agarose gel. These results are in agreement with literature (8, 10, 28, 29, 33, 34).

In our study, gray level resampling has a significant impact on radiomic parameters. This observation is in accordance with literature. Collewet et al. resampled gray scale using four methods: original gray levels, same maximum for all images, same average for all images, and limited dynamics. Results also
<figure><img src="image_11.png" /><figcaption>**FIGURE 6** | *In vivo* parameters value difference according to magnetic field strength: (A) entropy and (B) homogeneity. Each box plot represents the values dispersion of the 10 healthy volunteers for one of the six anatomical regions.</figcaption></figure>

showed the influence of the normalization method on the texture classification accuracy.

While our study focuses on field strength, pixel size and gray level resampling, previous studies have identified further MRI sequence parameters of interest, including repeat time (TR), echo time (TE) and receiver bandwidth (BW). Mayerhoefer et al. evaluated how texture features varied according to acquisition parameters (e.g., TR, TE, sampling bandwidth [SBW], and spatial resolution). Interestingly, variations in TR, TE and SBW had little effect on pattern discrimination results, as long as a high-spatial resolution was observed. In addition, both studies (Mayerhoefer's and ours) showed higher stability with the texture features derived from the co-occurrence matrix over first order texture features for characterizing patterns close to the resolution limits in varying conditions of TR/TE or pixel size.

Our approach to evaluate measures of radiomic features on standardized phantoms modeled to the organ under investigation is in line with the study by Bianchini et al. who explicitely designed and used an organ specific phantom for the optimization of radiomic studies of the female pelvis (35).

The question of the stability of our phantoms during the delay between the different acquisition/machine can be raised. Jirák et al. evaluated the long-term stability of the agarose phantoms, the optimal choice of texture parameters, and compared different MRI magnetic fields (3T, 4T, and 7T) (8, 10). They demonstrated that phantoms are stable over 12 months. In our study, all sequences were acquired on the same day within the same hour of examination, consequently we did not expect to incur phantom stability issue in our setup.

Optimal methodological guidelines are currently being defined to optimize data analysis strategies (4, 36–38) in the field of radiomics. Our results are critical since they demonstrate that standardization strategies for the use of radiomics in MRI should focus on addressing the challenge of heterogeneity in protocols since there is a dramatical difference between

radiomics features extracted from the same patient or the same phantom at 1.5 or 3T.

Our study demonstrates that field strength had a strong influence on texture feature values, also spatial resolution and gray scale resampling. Quantification of the impact of the various parameters of imaging features is of major interest. Our findings are clinically relevant as image acquisition was performed according to daily used protocols for clinical examinations. Those confounding factors need to be adjusted when designing a multicentric trial and when adapting the results of a study to different platforms.

## DATA AVAILABILITY STATEMENT

The datasets generated for this study are available on request to the corresponding author.

## ETHICS STATEMENT

This study was approved by the Scientific Review Committee of Gustave Roussy.

## AUTHOR CONTRIBUTIONS

All authors contributed to the article and approved the submitted version.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2020.541663/full#supplementary-material
REFERENCES

1. Savadjiev P, Chong J, Dohan A, Vakalopoulou M, Reinhold C, Paragios N, et al. Demystification of AI-driven medical image interpretation: past, present and future. *Eur Radiol* (2019) 29:1616-24. doi: 10.1007/s00330-018-5674-x

2. Rizzo S, Botta F, Raimondi S, Origgi D, Fanciullo C, Morganti AG, et al. Radiomics: the facts and the challenges of image analysis. *Eur Radiol Exp* (2018) 2:36. doi: 10.1186/s41747-018-0068-z

3. Sun R, Orlhac F, Robert C, Reuzé S, Schernberg A, Buvat I, et al. In Regard to Mattonen et al. *Int J Radiat Oncol Biol Phys* (2016) 95:1544-5. doi: 10.1016/j.ijrobp.2016.03.038

4. Lambin P, Leijenaar RTH, Deist TM, Peerlings J, de Jong EEC, van Timmeren J, et al. Radiomics: the bridge between medical imaging and personalized medicine. *Nat Rev Clin Oncol* (2017) 14:749-62. doi: 10.1038/nrclonin.2017.141

5. Limkin EJ, Sun R, Dercle L, Zacharaki EI, Robert C, Reuzé S, et al. Promises and challenges for the implementation of computational medical imaging (radiomics) in oncology. *Ann Oncol* (2017) 28:1191-206. doi: 10.1093/annonc/mdx034

6. Parmar C, Leijenaar RT, Grossmann P, Rios Velazquez E, Bussink J, Rietveld D, et al. Radiomic feature clusters and prognostic signatures specific for Lung and Head & Neck cancer. *Sci Rep* (2015) 5:11044. doi: 10.1038/srep11044

7. Sun R, Limkin EJ, Dercle L, Reuzé S, Zacharaki EI, Chargari C, et al. [Computational medical imaging (radiomics) and potential for immuno-oncology]. *Cancer Radiother* (2017) 21:648-54. doi: 10.1016/j.canrad.2017.07.035

8. Lerski RA, Schad LR, Luypaert R, Amorison A, Muller RN, Mascaro L, et al. Multicentre magnetic resonance texture analysis trial using reticulated foam test objects. *Magn Reson Imaging* (1999) 17:1025-31. doi: 10.1016/S0730-725X(99)00034-X

9. Mayerhoefer ME, Szomolanyi P, Jirak D, Materka A, Trattnig S. Effects of MRI acquisition parameter variations and protocol heterogeneity on the results of texture analysis and pattern discrimination: an application-oriented study. *Med Phys* (2009) 36:1236-43. doi: 10.1118/1.3081408

10. Jirák D, Dezortová M, Hájek M. Phantoms for texture analysis of MR images. *Long-term Multi-center Study Med Phys* (2004) 31:616-22. doi: 10.1118/1.1646231

11. Hajek M, Dezortova M, Materka A, Lerski R. *Texture analysis for magnetic resonance imaging*. Prague: Med4publishing (2006).

12. Orlhac F, Nioche C, Soussan M, Buvat I. Understanding Changes in Tumor Texture Indices in PET: A Comparison Between Visual Assessment and Index Values in Simulated and Patient Data. *J Nucl Med* (2017) 58:387-92. doi: 10.2967/jnumed.116.181859

13. Shiri I, Rahmim A, Ghaffarian P, Geramifar P, Abdollahi H, Bitarafan-Rajabi A. The impact of image reconstruction settings on 18F-FDG PET radiomic features: multi-scanner phantom and patient studies. *Eur Radiol* (2017) 27:4498-509. doi: 10.1007/s00330-017-4859-z

14. Carré A, Klausner G, Edjlali M, Lerousseau M, Briend-Diop J, Sun R, et al. Standardization of brain MR images across machines and protocols: bridging the gap for MRI-based radiomics. *Sci Rep* (2020) 10:12340. doi: 10.1038/s41598-020-69298-z

15. Simmons A, Tofts PS, Barker GJ, Arridge SR. Sources of intensity nonuniformity in spin echo images at 1.5 T. *Magn Reson Med* (1994) 32:121-8. doi: 10.1002/mrm.1910320117

16. Antunes J, Viswanath S, Rusu M, Valls L, Hoimes C, Avril N, et al. Radiomics Analysis on FLT-PET/MRI for Characterization of Early Treatment Response in Renal Cell Carcinoma: A Proof-of-Concept Study. *Transl Oncol* (2016) 9:155-62. doi: 10.1016/j.tranon.2016.01.008

17. Chirra P, Leo P, Yim M, Bloch BN, Rastinehad AR, Purysko A, et al. Multisite evaluation of radiomic feature reproducibility and discriminability for identifying peripheral zone prostate tumors on MRI. *J Med Imaging* (Bellingham) (2019) 6:024502. doi: 10.1117/1.JMI.6.2.04502

18. Moradmand H, Aghamiri SMR, Ghaderi R. Impact of image preprocessing methods on reproducibility of radiomic features in multimodal magnetic resonance imaging in glioblastoma. *J Appl Clin Med Phys* (2020) 21:179-90. doi: 10.1002/acm2.12795

19. Bologna M, Corino V, Mainardi L. Technical Note: Virtual phantom analyses for preprocessing evaluation and detection of a robust feature set for MRI-radiomics of the brain. *Med Phys* (2019) 46:5116-23. doi: 10.1002/ mp.13834

20. Shafiq-Ul-Hassan M, Zhang GG, Latifi K, Ullah G, Hunt DC, Balagunathan Y, et al. Intrinsic dependencies of CT radiomic features on voxel size and number of gray levels. *Med Phys* (2017) 44:1050-62. doi: 10.1002/mp.12123

21. Rathore S, Akbari H, Rozycki M, Abdullah KG, Nasrallah MP, Binder ZA, et al. Radiomic MRI signature reveals three distinct subtypes of glioblastoma with different clinical and molecular characteristics, offering prognostic value beyond IDH1. *Sci Rep* (2018) 8:5087. doi: 10.1038/s41598-018-22739-2

22. Zinn PO, Singh SK, Kotrotsou A, Hassan I, Thomas G, Luedi MM, et al. A Coclinical Radiogenomic Validation Study: Conserved Magnetic Resonance Radiomic Appearance of Periostin-Expressing Glioblastoma in Patients and Xenograft Models. *Precision Medicine and Imaging* (2018) 24:6288-99. doi: 10.1158/1078-0432.CCR-17-3420

23. Yip SS, Aerts HJ. Applications and limitations of radiomics. *Phys Med Biol* (2016) 61:R150-66. doi: 10.1088/0031-9155/61/13/R150

24. Lerski RA, de Certaines JD. Performance assessment and quality control in MRI by Eurospin test objects and protocols. *Magn Reson Imaging* (1993) 11:817-33. doi: 10.1016/0730-725X(93)90199-N

25. Buvat I, Orlhac F, Soussan M. Tumor Texture Analysis in PET: Where Do We Stand? *J Nucl Med* (2015) 56:1642-4. doi: 10.2967/jnumed.115.163469

26. Orlhac F, Soussan M, Maisonobe JA, Garcia CA, Vanderlinden B, Buvat I. Tumor texture analysis in 18F-FDG PET: relationships between texture parameters, histogram indices, standardized uptake values, metabolic volumes, and total lesion glycolysis. *J Nucl Med* (2014) 55:414-22. doi: 10.2967/jnumed.113.129858

27. Shukla G, Alexander GS, Bakas S, Nikam R, Talekar K, Palmer JD, et al. Advanced magnetic resonance imaging in glioblastoma: a review. *Chin Clin Oncol* (2017) 6:40. doi: 10.21037/ccco.2017.06.28

28. Xi YB, Guo F, Xu ZL, Li C, Wei W, Tian P, et al. Radiomics signature: A potential biomarker for the prediction of MGMT promoter methylation in glioblastoma. *J Magn Reson Imaging* (2018) 47:1380-7. doi: 10.1002/jmri.25860

29. Dercle L, Ammari S, Bateson M, Durand PB, Haspinger E, Massard C, et al. Limits of radiomic-based entropy as a surrogate of tumor heterogeneity: ROI-area, acquisition protocol and tissue site exert substantial influence. *Sci Rep* (2017) 7:7952. doi: 10.1038/s41598-017-08310-5

30. Lao J, Chen Y, Li ZC, Li Q, Zhang J, Liu J, et al. A Deep Learning-Based Radiomics Model for Prediction of Survival in Glioblastoma Multiforme. *Sci Rep* (2017) 7:10353. doi: 10.1038/s41598-017-10649-8

31. Leger S, Zwanenburg A, Pilz K, Lohaus F, Linge A, Zöphel K, et al. A comparative study of machine learning methods for time-to-event survival data for radiomics risk modelling. *Sci Rep* (2017) 7:13206. doi: 10.1038/s41598-017-13448-3

32. *Quantification of magnetic resonance image texture*. European Commission (1998).

33. 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. doi: 10.1016/j.mri.2003.09.001

34. Waugh SA, Lerski RA, Bidaut L, Thompson AM. The influence of field strength and different clinical breast MRI protocols on the outcome of texture analysis using foam phantoms. *Med Phys* (2011) 38:5058-66. doi: 10.1118/1.3622605

35. Bianchini L, Botta F, Origgi D, Rizzo S, Mariani M, Summers P, et al. PETER PHAN: An MRI phantom for the optimisation of radiomic studies of the female pelvis. *Phys Med* (2020) 71:71-81. doi: 10.1016/j.ejmp.2020.02.003

36. Mongan J, Moy L, Kahn CEJr. Checklist for Artificial Intelligence in Medical Imaging (CLAIM): A guide for authors and reviewers. *Radiol Soc North Am* (2020) 2:200029. doi: 10.1148/ryai.2020200029

37. Zwanenburg A. Radiomics in nuclear medicine: robustness, reproducibility, standardization, and how to avoid data analysis traps and replication crisis. *Eur J Nucl Med Mol Imaging* (2019) 46:2638-55. doi: 10.1007/s00259-019-04391-8

38. Zwanenburg A, Leger S, Vallières M, Löck S. Image biomarker standardisation initiative—feature definitions. (2016). eprint, 2019.
**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Ammari, Pitre-Champagnat, Dercle, Chouzenoux, Moalla, Reuze, Talbot, Mokoyoko, Hadchiti, Difetocq, Volk, El Haik, Lakiss, Balleyguier, Lassau and

Bidault. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.