The purpose of this study is to develop a MRF-ASL reconstruction algorithm using deep learning (DeepMARS). Compared with the traditional dictionary matching, our DeepMARS achieved higher intra-class correlation (ICC) in B1 (0.971 vs 0.921) and BAT (0.926 vs 0.761), similar ICC in T1 (0.957 vs 0.964) and CBF (0.936 vs 0.948) in the reproducibility test with much shorter calculation time per voxel (0.368 ms vs 2.899 s), suggesting that our DeepMARS may be a better alternative than the conventional MR dictionary matching approach.
The MRF-ASL signal was considered under a single-compartment model [1] and the goal is to estimate four quantitative parameters simultaneously: B1, T1, cerebral blood flow (CBF), and bolus arrive time (BAT). The primarily focus of this study is to compare a novel deep learning (DL) based parametric estimation method to a traditional dictionary matching (DM) [2]:
DM method: First, a MRF dictionary was simulated with 6149000 combinations of B1 values in the range of 0.6~1.1 (in an increment of 0.02), T1 values in the range of 500~5000 ms (in an increment of 20 ms within 2000 ms, and 300 ms beyond 2000 ms), CBF values in the range of 3~150 ml/100g/min (in an increment of 3 ml/100g/min), and BAT values in the range of 300~3000 ms (in an increment of 50 ms). Then, for a given MRF signal time course, its cross-correlations (cc) with all entries in the dictionary were calculated and the highest cc was identified.
DL method: The architecture of our DeepMARS was illustrated in Fig. 1, and five models were trained: Model 1 used dictionary without noise and mean square error (MSE) as the loss function, Model 2 used the dictionary contaminated with Gaussian noise and MSE as loss function, Model 3 used the dictionary with Gaussian noise and mean absolute percentage error (MAPE) as loss function, Model 4 used Model 1 as the pre-trained model and was further trained by the dictionary with Gaussian noise and used MSE as loss function. Model 5 used Model 1 as the pre-trained model and was further trained by the dictionary with Gaussian noise and used MAPE as loss function. The signal to noise ratio of the Gaussian noise was randomly selected from 20 dB to 60 dB. The MSE was defined as: $$$({\parallel True_{Para}-Pred_{Para}\parallel_{2}})^2$$$, and the MAPE was defined as:$$$\sum_{p\in Para}{\left| True_{p}-Pred_{p}\right|}/True_{p}$$$, where $$$Para$$$ was $$$\left[B1, T2, CBF, BAT \right]$$$. The voxel-wise intra-class correlation (ICC) [3] was calculated for each slice to test the reproducibility. The calculation time per voxel of DM and DeepMARS was compared on a 3.2 GHz Inter Core i5 with 26 GB RAM.
1. Pan, S., et al., Multiparametric estimation of brain hemodynamics with MR fingerprinting ASL. Magnetic Resonance in Medicine, 2017. 78(5): p. 1812-1823.
2. Ma, D., et al., Magnetic resonance fingerprinting. Nature, 2013. 495(7440): p. 187-192.
3. McGraw, K. and S.P. Wong, Forming Inferences About Some Intraclass Correlation Coefficients. Vol. 1. 1996. 30-46.