Spatio-Temporal dMRI Acquisition Design: Reducing the Number of Samples

Synopsis Acquisition time is a major limitation in recovering brain white matter microstructure with diffusion magnetic resonance imaging. Finding a sampling scheme that maximizes signal quality and satisfies given time constraints is NP-hard. Therefore, we propose a heuristic method based on genetic algorithm that finds sub-optimal solutions in reasonable time. Our diffusion model is defined in the qτ-space, so that it captures both spacial and temporal phenomena. The experiments on synthetic data and in-vivo diffusion images of the C57Bl6 wild-type mouse corpus callosum reveal superiority of our approach over random sampling and even distribution in the qτ-space. Introduction Brain white matter (WM) microstructure recovery with diffusion Magnetic Resonance Imaging (dMRI) requires lengthy acquisition which is unattainable in clinical practice. Dense scanning schemes studied by researchers [1-5] typically take few hours of imaging time, whereas human subjects can tolerate a little more than one hour [6, 7]. Nonetheless, recent in vivo studies of the WM microstructure [7-9] call for more fine-grained investigation of both space-and time-dependent diffusion. In this work, we aim at bridging the gap between growing demands on spatio-temporal (qτ) probing of dMRI signal [10] and acquisition time limitations. To this end, we propose an acquisition design that reduces the number of samples under adjustable quality loss. Most of the current acquisition schemes assume the fixed τ case, focusing on a dense sampling of the q-space instead [3-5]. However, a pronounced time-dependence in dMRI was recently reported by De Santis et al. [9], Burcaw et al. [11], and Novikov et al. [12]. Their results incline towards paying more attention to temporal phenomena in dMRI signal by incorporating multiple τ variants into acquisition schemes. Methods The main goal of our study is to find a qτ-indexed sampling scheme that best preserves the dMRI signal while satisfying given acquisition time limits [10,13]. We formulate the acquisition design task as an optimization problem. Furthermore, we want our approach to be applicable for real data. To this end, we discretize the spatio-temporal search space by performing a state-of-the-art dense pre-acquisition of dMRI signal. The problem thus boils down to selecting an optimal subset of Diffusion Weighted Images (DWIs), which is NP-hard [13]. Taking into account that the time complexity of our problem grows exponentially with the size of domain, such that global optima cannot be found deterministically within few hours or even few days, we apply a stochastic search engine instead. We use Standard Genetic Algorithm (SGA) [14] for this purpose due to its fast convergence rate, ability to avoid local optima, and the fact that it is based on the mathematically profound Markov Chain model [15]. Experiments For evaluation of our approach, we used both synthetic diffusion data and in vivo dMRI images of the C57Bl6 wild-type mouse. The dense pre-acquisition of signals covered 40 shells, each of which comprised 20 directions and one b 0-image, i.e. 40 × 20 = 800 DWIs plus 40 non-weighted images. We used combinations of 5 separation times Δ ∈ {10.8, 13.1, 15.4, 17.7, 20.0} [ms] and 8 gradient strengths G ∈ {50, 100, 150, 200, 250, 300, 350, 400} [mT/m]. The gradient duration δ = 5 ms remained constant throughout the experiments. We considered four variants of time limits expressed as budget sizes n max = {100, 200, 300, 400} out of 800 DWIs. We compared our method with two alternative sampling schemes. One of them, called random, used the uniform random distribution of qτ samples in the index space {1, …, N}. In the second one, referred to as even, we picked each i-th sample for i = ⌊kN / n max ⌋ and k = 1,. , n max. Discussion As Figure 1 shows, our method outperformed the other two in all analyzed cases, assuring lowest mean squared errors (MSEs) and standard deviations (STDs). We verified statistical significance of the results with the two-sample Student's t level α = 10 − 5