A polynomial regression-based approach to estimate relaxation rate maps suitable for multiparametric segmentation of clinical brain MRI studies in multiple sclerosis

Background and Objective: Relaxation parameter maps (RPMs) calculated from spin-echo data have provided a basis for the segmentation of normal brain tissues and white matter lesions in multiple sclerosis (MS) MRI studies. However, Conventional Spin-Echo (CSE) sequences, once the core of clinical MRI studies, have been largely replaced by faster ones, which do not allow the calculation a-posteriori of RPMs from clinical studies. Aim of the study was to develop and validate a method to estimate RPMs (pseudo-RPMs) from routine clinical MRI protocols (including 3D-Gradient Echo T1w, FLAIR and fast-T2w sequences), suitable for fully automatic multiparametric segmentation of normal-appearing and pathological brain tissues in MS. Methods: The proposed method processes spatially normalized clinical MRI studies through a multistep pipeline, to collect a set of data points of matched signal intensities (from MRI studies) and relaxation parameters (from a CSE-derived digital template and an MS lesion database), which are then fitted by a multiple and multivariate 4-th degree polynomial regression, providing pseudo-RPMs. The method was applied to a dataset of 59 clinical MRI studies providing pseudo-RPMs that were segmented through a method originally developed for the CSE-derived RPMs. Results of the segmentation in 12 studies were used to iteratively optimize method parameters. Accuracy of segmentation of normal-appearing brain tissues from the pseudo-RPMs was assessed by comparing their age-related changes, as measured in 47 clinical studies, against those measured acquired using CSE sequences in a comparable dataset of 47 patients. Lesion segmentation was validated against manual segmentation carried out by three neuroradiologists. Results: Age-related changes of normal-appearing brain tissue volumes measured using the pseudo-RPMs substantially overlapped those measured using the RPMs obtained from CSE sequences, and segmentation of MS lesions showed a moderate-high spatial overlap with manual segmentation, comparable to that achieved by the widely used Lesion Segmentation Tool on FLAIR images, with a greater volumetric agreement. Conclusions: The proposed approach allows calculation from clinical studies of pseudo-RPMs, which are equivalent to those obtainable from CSE sequences, avoiding the need for the acquisition of additional, dedicated sequences for segmentation purposes.