Optimal imaging parameters for fiber-orientation estimation in diffusion MRI

Abstract

This study uses Monte Carlo simulations to investigate the optimal value of the diffusion weighting factor b for estimating white-matter fiber orientations using diffusion MRI with a standard spherical sampling scheme. We devise an algorithm for determining the optimal echo time, pulse width, and pulse separation in the pulsed-gradient spin-echo sequence for a specific value of b. The Monte Carlo simulations provide an estimate of the optimal value of b for recovering one and two fiber orientations. We show that the optimum is largely independent of the noise level in the measurements and the number of gradient directions and that the optimum depends only weakly on the diffusion anisotropy, the maximum gradient strength, and the spin–spin relaxation time. The optimum depends strongly on the mean diffusivity. In brain tissue, the optima we estimate are in the ranges [0.7, 1.0] × 10⁹ s m⁻² and [2.2, 2.8] × 10⁹ s m⁻² for the one- and two-fiber cases, respectively. The best b for estimating the fractional anisotropy is slightly higher than for estimating fiber directions in the one-fiber case and slightly lower in the two-fiber case. To estimate Tr(D) in the one-fiber case, the optimal setting is higher still. Simulations suggest that a ratio of high to low b measurements of 5 to 1 is a good compromise for measuring fiber directions and size and shape indices.

Introduction

This paper investigates the optimal imaging parameters for estimating white-matter fiber orientations using diffusion MRI. Diffusion MRI measures the density function p of particle displacements over a fixed time t. The anisotropy of p reflects the anisotropy of the microstructure of the material in which the particles move. In human-brain imaging, the particles of interest are usually water molecules. In white matter, water molecules move furthest on average in fiber directions, so the directions of the ridges of p, or the peaks at a fixed radius, provide fiber-orientation estimates.

Diffusion MRI measurements are approximate samples of the Fourier transform of p with wavenumber q, which depends on the length, separation, strength, and direction of gradient pulses in the imaging sequence. Here, we synthesize diffusion MRI measurements from brain tissue with one and two dominant microstructural fiber orientations using simple models of p. We refit the models of p to the synthetic data to recover the fiber orientations. In the one-fiber case, the model is a zero-mean Gaussian density, as in diffusion-tensor MRI (Basser et al., 1994), and in the two-fiber case, a mixture of two zero-mean Gaussians, as in Tuch et al. (2002) and Parker and Alexander (2003). We study the error in the fiber-orientation estimates in the presence of noise as a function of the parameters of the imaging sequence in order to identify optimal settings for brain imaging.

Xing et al. (1997) investigate measurement of the diffusion coefficient using diffusion-weighted MRI. They show numerically that, if the signal decays mono-exponentially as a function of b = t|q|², the signal-to-noise ratio of the diffusion coefficient d, obtained from two measurements with different values b₁ < b₂ of b, is maximum when d(b₂ − b₁) = 1.11. Further, they derive optimal strategies for measuring d with more than two measurements. As the number of measurements increases, the best strategy acquires repeated samples at b₁ and b₂; the optimal number N of measurements at b₂ tends to 3.6 times the number M of measurements at b₁ and the optimal d(b₂ − b₁) increases to an asymptote at 1.28. Jones et al. (1999) extend the work of Xing et al. (1997) to estimation of the full diffusion tensor. They set b₁ = 0 and extend the two-point sampling scheme to three dimensions by acquiring M measurements with b = 0 together with N measurements with fixed |q| > 0, fixed t (and hence fixed b₂) and unique gradient directions qˆ_i, i = 1,…, N, well separated on the unit sphere. They show numerically that, for an approximately isotropic diffusion tensor D, setting b₂ = 1.09 × 3/Tr(D) and M / (N + M) = 0.081 minimizes the sum of the variances of the elements of D. Jones et al. (1999) note that, in practice, higher values of b require longer echo time TE. More transverse relaxation occurs before readout at higher TE, which reduces the signal-to-noise ratio of the measurements. However, the earlier work of Xing et al. (1997) ignores the effects of transverse relaxation. Jones et al. (1999) show that, for a typical measurement sequence, including transverse-relaxation effects in the calculation of the optimal b and M / (N + M), with the spin–spin relaxation constant T₂ = 0.08 s (typical for brain tissue), reduces the optimal b₂ to 0.85 × 3/Tr(D) and increases the optimal M / (N + M) to 0.103. Kingsley and Monahan (2004) correct some formulae in Jones et al. (1999). They use similar simulations, which account for TE increasing with b, to select the b that minimizes the variance of the mean diffusivity for typical diffusivities found in the adult and neonatal human brain with and without the effects of ischemic stroke. They note that diffusion anisotropy reduces the optimal b from Jones et al.'s (1999) optimum and recommend b = 1.0 × 10⁹ s m⁻² for discriminating regions affected by stroke from normal tissue. Armitage and Bastin (2001) extend the work of Xing et al. (1997) by considering samples with a range of values of d; specifically, the range observed in the human brain. Their analysis shows that b₂ − b₁ = 0.9 × 10⁹ s m⁻² minimizes the number of measurements required to measure the whole range of diffusion coefficients within a fixed error threshold. Armitage and Bastin (2001) use Xing et al.'s (1997) two-point sampling strategy, but Brihuega-Moreno et al. (2003) point out that the two-point sampling strategy is optimal only for measuring a single diffusion coefficient. For samples containing a range of values of d, the optimal strategy acquires measurements at more than two values of b. Brihuega-Moreno et al. (2003) derive optimal sampling strategies, specifically for cartilage imaging, with more than two sample points. Armitage and Bastin (2001) and Brihuega-Moreno et al. (2003) do not take into account the reduction in signal from transverse relaxation as b increases.

We model the measurement process to synthesize data for Monte Carlo simulations. We limit investigations to Jones et al.'s (1999) spherical sampling scheme with M measurements at b = 0 and N at fixed b > 0 with well separated qˆ. Unless stated otherwise, we use sets of qˆ_i, i = 1,…, N, that minimize the electrostatic energy of N pairs of equal and opposite points on the sphere with equal charges, which we obtain using the method outlined by Jansons and Alexander (2003). We consider only the pulsed-gradient spin-echo (PGSE) sequence with echo-planar imaging (EPI) readout. We devise an algorithm to choose the TE that maximizes the signal-to-noise ratio of the measurement at b = 0 for a particular b in the N measurements with b > 0. The algorithm accounts for transverse-relaxation effects. We assume that the image resolution, and hence the time required to readout half of k-space, is fixed. Unlike Jones et al. (1999) and Kingsley and Monahan (2004), however, we do not assume that the EPI readout samples the whole of k-space. We assume that the sequence samples exactly half of k-space after TE and samples some variable amount of the other half of k-space within TE in the first part of the echo. The algorithm finds the optimal partial k-space sampling time within TE. The algorithm also determines the optimal duration and separation of the diffusion-weighted gradient pulses. We use Monte Carlo simulations to plot the error in reconstructed directions, for both the one and two-fiber cases, against b for the optimized pulse sequences. We also study the accuracy and stability of indices of size and shape from both one- and two-fiber reconstructions. Specifically, we measure the error in Tr(D) and the fractional anisotropy (Basser and Pierpaoli, 1996) as a function of b.