Search

Functional magnetic resonance imaging (fMRI) data that are corrupted by temporally colored noise are generally pre-processed (i.e., prewhitened or precolored) prior to functional activation detection. In this paper, we propose likelihood-based hypothesis tests that account for colored noise directly within the framework of functional activation detection. Three likelihood-based tests are proposed: the generalized likelihood ratio (GLR) test, the Wald test, and the Rao test. The fMRI time series is modeled as a linear regression model, where one regressor describes the task-related hemodynamic response, one regressor accounts for a constant baseline and one regressor describes potential drift. The temporal correlation structure of the noise is modeled as an autoregressive (AR) model. The order of the AR model is determined from practical null data sets using Akaike's information criterion (with penalty factor 3) as order selection criterion. The tests proposed are based on exact expressions for the likelihood function of the data. Using Monte Carlo simulation experiments, the performance of the proposed tests is evaluated in terms of detection rate and false alarm rate properties and compared to the current general linear model (GLM) test, which estimates the coloring of the noise in a separate step. Results show that theoretical asymptotic distributions of the GLM, GLR, and Wald test statistics cannot be reliably used for computing thresholds for activation detection from finite length time series. Furthermore, it is shown that, for a fixed false alarm rate, the detection rate of the proposed GLR test statistic is slightly, but (statistically) significantly improved compared to that of the common GLM-based tests. Finally, simulations results reveal that all tests considered show seriously inferior performance if the order of the AR model is not chosen sufficiently high to give an adequate description of the correlation structure of the noise, whereas the effects of (slightly) overmodeling are observed to be less harmful. [ABSTRACT FROM AUTHOR]

Diffusion weighted MRI is an indispensable tool for routine patient screening and diagnostics of pathology. Recently, several deep learning methods have been proposed to quantify diffusion parameters, but poor generalisation to new data prevents broader use of these methods, as they require retraining of the neural network for each new scan protocol. In this work, we present the dtiRIM, a new deep learning method for Diffusion Tensor Imaging (DTI) based on the Recurrent Inference Machines. Thanks to its ability to learn how to solve inverse problems and to use the diffusion tensor model to promote data consistency, the dtiRIM can generalise to variations in the acquisition settings. This enables a single trained network to produce high quality tensor estimates for a variety of cases. We performed extensive validation of our method using simulation and in vivo data, and compared it to the Iterated Weighted Linear Least Squares (IWLLS), the approach of the state-of-the-art MRTrix3 software, and to an implementation of the Maximum Likelihood Estimator (MLE). Our results show that dtiRIM predictions present low dependency on tissue properties, anatomy and scanning parameters, with results comparable to or better than both IWLLS and MLE. Further, we demonstrate that a single dtiRIM model can be used for a diversity of data sets without significant loss in quality, representing, to our knowledge, the first generalisable deep learning based solver for DTI., (Copyright © 2023. Published by Elsevier Inc.)

MR fingerprinting (MRF) is a promising method for quantitative characterization of tissues. Often, voxel-wise measurements are made, assuming a single tissue-type per voxel. Alternatively, the Sparsity Promoting Iterative Joint Non-negative least squares Multi-Component MRF method (SPIJN-MRF) facilitates tissue parameter estimation for identified components as well as partial volume segmentations. The aim of this paper was to evaluate the accuracy and repeatability of the SPIJN-MRF parameter estimations and partial volume segmentations. This was done (1) through numerical simulations based on the BrainWeb phantoms and (2) using in vivo acquired MRF data from 5 subjects that were scanned on the same week-day for 8 consecutive weeks. The partial volume segmentations of the SPIJN-MRF method were compared to those obtained by two conventional methods: SPM12 and FSL. SPIJN-MRF showed higher accuracy in simulations in comparison to FSL- and SPM12-based segmentations: Fuzzy Tanimoto Coefficients (FTC) comparing these segmentations and Brainweb references were higher than 0.95 for SPIJN-MRF in all the tissues and between 0.6 and 0.7 for SPM12 and FSL in white and gray matter and between 0.5 and 0.6 in CSF. For the in vivo MRF data, the estimated relaxation times were in line with literature and minimal variation was observed. Furthermore, the coefficient of variation (CoV) for estimated tissue volumes with SPIJN-MRF were 10.5% for the myelin water, 6.0% for the white matter, 5.6% for the gray matter, 4.6% for the CSF and 1.1% for the total brain volume. CoVs for CSF and total brain volume measured on the scanned data for SPIJN-MRF were in line with those obtained with SPM12 and FSL. The CoVs for white and gray matter volumes were distinctively higher for SPIJN-MRF than those measured with SPM12 and FSL. In conclusion, the use of SPIJN-MRF provides accurate and precise tissue relaxation parameter estimations taking into account intrinsic partial volume effects. It facilitates obtaining tissue fraction maps of prevalent tissues including myelin water which can be relevant for evaluating diseases affecting the white matter., Competing Interests: Declaration of Competing Interest The authors declare that there is no conflict of interest., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

Purpose: Brain water content provides rich tissue contrast comparable to that of longitudinal relaxation time T 1 , but mapping is usually performed at modest resolution. In particular, the slice thickness in 2D mapping methods is limited. Here, we combine super-resolution reconstruction techniques with a fast water content mapping method to acquire high and isotropic resolution (0.75 mm) water content maps at 3 Tesla., Methods: A high-resolution multi-echo gradient echo image is super-resolution-reconstructed from 3 low-resolution, orthogonal multi-echo gradient echo image acquisitions, followed by water content mapping. The mapping accuracy and SNR of the proposed method are assessed using numerical simulations, phantom studies, and in vivo data acquired from 6 healthy volunteers at 3 Tesla. A high-resolution acquisition with an established mapping method is used as a reference., Results: Whole-brain water content maps with 0.75 mm isotropic resolution are demonstrated. No bias in the water content values was seen following super-resolution reconstruction. In the in vivo experiments, a lower SD of the mean water content values was observed with the proposed method compared to the reference method., Conclusions: Super-resolution reconstruction of multi-echo gradient echo data is demonstrated, enabling whole-brain water content mapping with high and isotropic resolution. The accuracy of the proposed method is shown using phantoms and 6 healthy volunteers and was found to be unchanged compared to the conventional acquisition. The proposed method could increase the sensitivity of water content mapping sufficiently to enable the detection of very small lesions, such as cortical lesions in multiple sclerosis., (© 2022 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals LLC on behalf of International Society for Magnetic Resonance in Medicine.)

In this paper, we propose the use of Recurrent Inference Machines (RIMs) to perform T 1 and T 2 mapping. The RIM is a neural network framework that learns an iterative inference process based on the signal model, similar to conventional statistical methods for quantitative MRI (QMRI), such as the Maximum Likelihood Estimator (MLE). This framework combines the advantages of both data-driven and model-based methods, and, we hypothesize, is a promising tool for QMRI. Previously, RIMs were used to solve linear inverse reconstruction problems. Here, we show that they can also be used to optimize non-linear problems and estimate relaxometry maps with high precision and accuracy. The developed RIM framework is evaluated in terms of accuracy and precision and compared to an MLE method and an implementation of the Residual Neural Network (ResNet). The results show that the RIM improves the quality of estimates compared to the other techniques in Monte Carlo experiments with simulated data, test-retest analysis of a system phantom, and in-vivo scans. Additionally, inference with the RIM is 150 times faster than the MLE, and robustness to (slight) variations of scanning parameters is demonstrated. Hence, the RIM is a promising and flexible method for QMRI. Coupled with an open-source training data generation tool, it presents a compelling alternative to previous methods., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 The Author(s). Published by Elsevier B.V. All rights reserved.)

Purpose: Quantification of the T 2 ∗ relaxation time constant is relevant in various magnetic resonance imaging applications. Mono- or bi-exponential models are typically used to determine these parameters. However, in case of complex, heterogeneous tissues these models could lead to inaccurate results. We compared a model, provided by the fractional-order extension of the Bloch equation with the conventional models., Methods: Axial 3D ultra-short echo time (UTE) scans were acquired using a 3.0 T MRI and a 16-channel surface coil. After image registration, voxel-wise T 2 ∗ was quantified with mono-exponential, bi-exponential and fractional-order fitting. We evaluated all three models repeatability and the bias of their derived parameters by fitting at various noise levels. To investigate the effect of the SNR for the different models, a Monte-Carlo experiment with 1000 repeats was performed for different noise levels for one subject. For a cross-sectional investigation, we used the mean fitted values of the ROIs in five volunteers., Results: Comparing the mono-exponential and the fractional order T 2 ∗ maps, the fractional order fitting method yielded enhanced contrast and an improved delineation of the different tissues. In the case of the bi-exponential method, the long T 2 ∗ component map demonstrated the anatomy clearly with high contrast. Simulations showed a nonzero bias of the parameters for all three mathematical models. ROI based fitting showed that the T 2 ∗ values were different depending on the applied method, and they differed most for the patellar tendon in all subjects., Conclusions: In high SNR cases, the fractional order and bi-exponential models are both performing well with low bias. However, in all observed cases, one of the bi-exponential components has high standard deviation in T 2 ∗ . The bi-exponential model is suitable for T 2 ∗ mapping, but we recommend using the fractional order model for cases of low SNR., (Copyright © 2020 Elsevier Inc. All rights reserved.)