Physics-Informed Discretization for Reproducible and Robust Radiomic Feature Extraction Using Quantitative MRI
Authors: Zhao Walter , Hu Zheyuan , Kazerooni Anahita Fathi , Körzdörfer Gregor , Nittka Mathias , Davatzikos Christos , Viswanath Satish E. , Wang Xiaofeng , Badve Chaitra , Ma Dan 
Published: October 2023

Objective
Given the limited repeatability and reproducibility of radiomic features derived from weighted magnetic resonance imaging (MRI), there may be significant advantages to using radiomics in conjunction with quantitative MRI. This study introduces a novel physics-informed discretization (PID) method for reproducible radiomic feature extraction and evaluates its performance using quantitative MRI sequences including magnetic resonance fingerprinting (MRF) and apparent diffusion coefficient (ADC) mapping.

Materials and Methods
A multiscanner, scan-rescan dataset comprising whole-brain 3D quantitative (MRF T1, MRF T2, and ADC) and weighted MRI (T1w MPRAGE, T2w SPACE, and T2w FLAIR) from 5 healthy subjects was prospectively acquired. Subjects underwent 2 repeated acquisitions on 3 distinct 3 T scanners each, for a total of 6 scans per subject (30 total scans). First-order statistical (n = 23) and second-order texture (n = 74) radiomic features were extracted from 56 brain tissue regions of interest using the proposed PID method (for quantitative MRI) and conventional fixed bin number (FBN) discretization (for quantitative MRI and weighted MRI). Interscanner radiomic feature reproducibility was measured using the intraclass correlation coefficient (ICC), and the effect of image sequence (eg, MRF T1 vs T1w MPRAGE), as well as image discretization method (ie, PID vs FBN), on radiomic feature reproducibility was assessed using repeated measures analysis of variance. The robustness of PID and FBN discretization to segmentation error was evaluated by simulating segmentation differences in brainstem regions of interest. Radiomic features with ICCs greater than 0.75 following simulated segmentation were determined to be robust to segmentation.

Results
First-order features demonstrated higher reproducibility in quantitative MRI than weighted MRI sequences, with 30% (n = 7/23) features being more reproducible in MRF T1 and MRF T2 than weighted MRI. Gray level co-occurrence matrix (GLCM) texture features extracted from MRF T1 and MRF T2 were significantly more reproducible using PID compared with FBN discretization; for all quantitative MRI sequences, PID yielded the highest number of texture features with excellent reproducibility (ICC > 0.9). Comparing texture reproducibility of quantitative and weighted MRI, a greater proportion of MRF T1 (n = 225/370, 61%) and MRF T2 (n = 150/370, 41%) texture features had excellent reproducibility (ICC > 0.9) compared with T1w MPRAGE (n = 148/370, 40%), ADC (n = 115/370, 32%), T2w SPACE (n = 98/370, 27%), and FLAIR (n = 102/370, 28%). Physics-informed discretization was also more robust than FBN discretization to segmentation error, as 46% (n = 103/222, 46%) of texture features extracted from quantitative MRI using PID were robust to simulated 6 mm segmentation shift compared with 19% (n = 42/222, 19%) of weighted MRI texture features extracted using FBN discretization.

Conclusions
The proposed PID method yields radiomic features extracted from quantitative MRI sequences that are more reproducible and robust than radiomic features extracted from weighted MRI using conventional (FBN) discretization approaches. Quantitative MRI sequences also demonstrated greater scan-rescan robustness and first-order feature reproducibility than weighted MRI.