Numerical Simulations Comparing Pore Imaging Methods Based on Diffusion-Weighted MR Imaging
Yasar Goedecke1 and Jürgen Finsterbusch1

1Systems Neuroscience, University Medical Center Hamburg-Eppendorf, Hamburg, Germany

Synopsis

In a conventional diffusion-weighted MRI experiment, the signal amplitude depends on the squared magnitude of the Fourier transformation of the pore or cell geometry, i.e. the underlying cell or pore geometry cannot be reconstructed. Several approaches have been proposed that determine the otherwise missing phase information and, thus, can image the pore or cell geometry directly. Here, the performance of these methods is compared with respect to their applicability in practice, e.g. considering the impact of the noise level, mixtures of pore sizes, orientations, and shapes, and gradient pulse durations and diffusion times achievable on standard MRI systems.

Introduction

In a conventional diffusion-weighted MRI experiment, the signal amplitude depends on the squared magnitude of the Fourier transformation of the pore or cell geometry [1]. From such data, the pore or cell geometry cannot be reconstructed because the phase information is missing [1]. Recently, several approaches have been proposed that offer access to the phase and, thus, allow to reconstruct the pore or cell geometry directly [2-4]. Here, the performance of these methods is compared with respect to their applicability in practice, in particular in biological systems. considering the impact of the noise level, mixtures of pore sizes, orientations, and shapes, a freely diffusing compartment, and gradient pulse durations and diffusion times achievable on conventional MRI systems.

Methods

The pore imaging methods considered are the method of Laun et al. based on a short and a long diffusion-weighting gradient pulse [2] (”asymmetric” method), the method of Shemesh et al. involving the signal ratio of a conventional and a double diffusion encoding experiment [3] (“ratio” method), and the method of Kuder et al. that estimates the phase information iteratively from a double diffusion encoding experiment [4] (“phase iteration” method). The pulse sequences required for these methods are presented in Fig. 1.

Numerical simulations were performed using an IDL algorithm (version 7.1, ITT Visual Information Solutions, Boulder, Colorado, USA) of a random-walk model without or within confining geometries with diffuse reflection at the boundaries (see, e.g., [5]). Only two dimensions were considered for simplicity but all results can be expected to hold for three dimensions as well. Up to 107 spins were simulated with a temporal resolution of 10 µs using a diffusion coefficient of 2.0×10‑3 mm2 s‑1. Different pore shapes (square, rectangle, circle), orientations, and sizes (radius and edge length between 5 µm and 20 µm) were considered. Simulations were performed for a radial grid, regridded to a cartesian grid, and Fourier transformed to obtain the pore image. For some simulations, Gaussian noise was added prior to the Fourier transformation to investigate the performance in the presence of a finite signal-to-noise ratio (SNR).

For most simulations ideal timing parameters were used, however, some simulations were performed with timing parameters compatible with state-of-the-art gradient systems for animal and whole-body MR systems (500 mT m‑1 and 80 mT m‑1, respectively) or total durations of the diffusion weighting limited to 200 ms. For the “ratio” method, a thresholding approach was used that reduces noise amplification by setting the ratio to 0 for small values of the denominator (zero crossings of the denominator are also zero crossings of the nominator).

Results and Discussion

Pore images obtained for different noise levels are presented in Fig. 2. While the “phase iteration” method amplifies noise due to the iterative calculation of the phase, the “ratio” method suffers from blurring making the “asymmetric” method the most reliable one.

Results for mixtures of different pore orientations, sizes, and shapes are presented in Fig. 3. Again the “asymmetric” method provides the best performance resolving the orientations, sizes, and shapes while the other methods fail to reconstruct the underlying geometries correctly.

Images obtained for realistic timing parameters and gradient amplitudes are presented in Fig. 4. With finite gradient amplitudes, the pore size is underestimated, in particular for the weaker gradient amplitudes available on whole-body MR systems (about 13-14 µm for a pore size of 20 µm; Fig. 4c). This is because for finite gradient pulses the center-of-mass of the spin trajectory defines the effective encoding position [6] which moves towards the pore’s center for longer gradient pulses. Considering a finite time for the diffusion weighting does not have a significant impact on the “ratio” and “phase iteration” method but significantly affects the “asymmetric” method that yields a featureless image approaching the pore’s autocorrelation function (Fig. 4d and e) as for a conventional experiment. Since in practice the duration of the long gradient pulse is limited by T2 relaxation, improved gradient hardware cannot solve this problem but larger diffusion coefficients, smaller pores, or longer relaxation times are beneficial.

Adding a free compartment (Fig. 5) as can be relevant in biological systems does seem to have a significant impact on either of the methods.

Conclusion

The method of Laun et al. ("asymmetric" method) is the most robust one regarding the noise level and, unlike the other two methods, is able to resolve the pore geometry in mixtures which is important for more complex samples like biological tissues. However, in contrast to the other methods, it faces limitations for the timing parameters that are required in practice.

Acknowledgements

No acknowledgement found.

References

[1] Cory DG, Garroway AN. Measurement of translational displacement probabilities by NMR: an indicator of compartmentation. Magn. Reson. Med 1990; 14: 435–444

[2] Laun FB, Kuder TA, Semmler W, Stieltjes B. Determination of the defining boundary in nuclear magnetic resonance diffusion experiments. Phys. Rev. Lett. 2011; 107: 048102.

[3] Shemesh N, Westin CF, Cohen Y. Magnetic resonance imaging by synergistic diffusion-diffraction patterns. Phys. Rev. Lett. 2012; 108: 058103.

[4] Kuder TA, Laun FB. NMR-based diffusion pore imaging by double wave vector measurements. Magn. Reson. Med. 2013; 70: 836–841.

[5] Finsterbusch J, Numerical simulations of short-mixing-time double-wave-vector diffusion-weighting experiments with multiple concatenations on whole-body MR systems, J. Magn. Res. 2010; 207: 274–282.

[6] Mitra PP, Halperin BI, Effects of finite gradient-pulse width in pulsed-field-gradient diffusion measurements. J. Magn. Reson. A 1995; 113: 94–101.

Figures

Fig. 1: Basic pulse sequences considered in the present study: (a) Stejskal-Tanner diffusion weighting required for the “ratio” method (denominator), (b) “asymmetric” diffusion weighting with a prolonged rephasing pulse as used for the “asymmetric” method, and (c) double-wave-vector or double-diffusion-encoding sequence required for the “ratio” and the “phase iteration” method.

Fig. 2: Pore images obtained for simulations with different levels of Gaussian noise added. The “phase iteration” method shows the highest noise sensitivity. Similar noise effects can be avoided for the “ratio” method with the thresholding approach chosen but is at the expense of significant blurring.

Fig. 3: Pore images obtained for simulations with a mixture of (a) three pore orientations for a rectangular pore shape, (b) square pores with three different sizes (5 μm, 10 μm, and 20 μm), and (c) a square and a hollow cylinder.

Fig. 4: Pore images obtained with (a) ideal timing parameters and timing parameters compatible with a maximum gradient amplitude of (b, d) 500 mT m−1 and (c, e) 80 mT m−1 and (d, e) an echo time of 200 ms.

Fig. 5: Pore images obtained for simulations with realistic parameters (maximum gradient amplitude 80 mT m−1, echo time 200 ms) and a freely diffusing compartment with a volume fraction of 0% and 60%.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
1993