
ORIGINAL ARTICLE 



Year : 2007  Volume
: 32
 Issue : 1  Page : 3442 

Principles and limitations of NMR diffusion measurements
Jan Hrabe, Gurjinder Kaur, David N Guilfoyle
Center for Advanced Brain Imaging, Nathan S. Kline Institute Orangeburg, New York, USA
Correspondence Address: Jan Hrabe Center for Advanced Brain Imaging, Nathan S. Kline Institute, 140 Old Orangeburg Road, Orangeburg, NY 10962 USA
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/09716203.31148
Abstract   
Diffusion spectroscopy, imaging and particularly diffusion tensor imaging have become popular thanks to their numerous clinical and research applications which span from brain stroke evaluation to fiber tracking. With a few exceptions, these methods are rooted in the classic StejskalTanner formula for the diffusionattenuated signal, usually obtained by solving the BlochTorrey partial differential equations. Here we derive the StejskalTanner formula in the simplest possible manner, avoiding integrals and differential equations. This approach makes it easy to understand the origin of the diffusion signal attenuation, the effects of various diffusion sequence parameters, and also the numerous important pitfalls, which are discussed in the last section.
Keywords: b value, diffusion time, random walk, tortuosity, volume fraction
How to cite this article: Hrabe J, Kaur G, Guilfoyle DN. Principles and limitations of NMR diffusion measurements. J Med Phys 2007;32:3442 
Introduction   
Most readers are probably aware of the diffusion tensor imaging (DTI) method,[1] which is ubiquitous in the contemporary medical MR literature. However, DTI is a relatively recent advance in diffusion experiments. NMR was used to measure diffusion long before MR imaging was first proposed in the early 1970s. The undesirable influence of selfdiffusion on spin echo amplitudes was first recognized by Hahn as early as 1950.[2] Carr and Purcell[3] further extended Hahn's idea by using multiple echoes as a means of minimizing the diffusion effect. Shortly after these momentous discoveries, Torrey[4] extended the basic MR formalism of Bloch differential equations by including additional terms to accommodate the diffusion effect.
Stejskal and Tanner in their seminal paper[5] introduced pulsed gradients into the basic spin echo sequence, which resulted in much improved sensitivity to diffusion in comparison to the steady state gradients used previously. They solved the BlochTorrey partial differential equations for a symmetric pair of pulsed gradients and obtained the wellknown StejskalTanner formula
where
In Eqs. (1) and (2), S is the signal strength in a pulse sequence with a pair of balanced diffusionsensitizing gradients of strength G , each of a duration δ and with a delay t between them. S_{ 0} is the signal strength in an identical experiment but without the diffusion gradient pair. When it can be safely assumed that δ << t, the expression for b (usually called b value) simplifies to
Tanner and Stejskal were also the first to propose the idea of measuring restricted diffusion of water molecules[6] by varying the delay τbetween the gradient pulses. They discovered that the attenuation of signal from water molecules restricted by cellular walls in their diffusion movement was less pronounced than the signal from freely moving molecules and used this method to estimate the diameter of yeast cells.
Restricted diffusion plays a central role in DTI. Geometrically complex biological environments, such as the nervous tissue, can be characterized by two aggregate parameters,α and λ.[7] The volume fraction α is simply a percentage of the total tissue volume where molecules can diffuse. For example, in a healthy brain, the extracellular volume fraction is about 20% and the intracellular about 65%. The tortuosity parameter, λ, is more complicated.[8] It describes geometrical hindrance of an environment relative to an obstaclefree medium. If the obstacles exhibit some directional preference, the hindrance becomes anisotropic, that is, it depends on direction. For example, the molecules diffuse more readily along the white matter fibers than across them. In a macroscopically homogeneous and anisotropic environment, the tortuosity takes the shape of a symmetric tensor of the second order, which can be represented by a 3 x 3 matrix with six independent values. The tortuosity tensor combines with the scalar free diffusion coefficient into a tensor of apparent diffusion D . However, despite all this complexity, Eqs. (1) and (2) or (3) are still used to calculate signal attenuation due to diffusion along any single direction of the diffusion gradient. The only real difference from the homogeneous and isotropic case is that the experiment must be repeated using at least six noncollinear directions of the diffusion gradient to obtain six independent components of the apparent diffusion tensor. It is therefore fair to conclude that the pulse field gradient method proposed by Stejskal and Tanner is at the heart of most modern DTI experiments.
Unfortunately, there is an important caveat. A single imaging voxel in a nervous tissue often includes compartments with different diffusion properties both on a cellular level (e.g., intracellular and extracellular compartments) and on a larger scale (e.g., several crossing axonal bundles). Spins can also migrate from one compartment to another during the measurement. These complications violate the assumption of a macroscopically homogeneous and anisotropic environment with a single compartment. They may result in a nonGaussian behavior of the diffusion signal for which the second order diffusion tensor no longer represents an adequate description.
DTI has become a very popular MR imaging modality and is developing into an important tool for noninvasive study and characterization of the brain white matter. It has been applied to the study of many neurological brain disorders such as schizophrenia, cocaine addiction, HIV infection, alcoholism, geriatric depression and Alzheimer's disease. An overview of theoretical issues surrounding the DTI technique can be found in.[9],[10],[11] A review of DTI applications in neuroscience is presented in.[12]
It is our view that good understanding of the StejskalTanner formulae (Eqs. (1), (2) and (3)) is essential for sensible design or interpretation of any MR diffusion experiment. Unfortunately, the derivation in a completely general case is nontrivial as it involves either solving the BlochTorrey partial differential equations[4],[5] or integration in a complex plane.[13],[14] The goal of this report is to present a detailed derivation of Eqs. (1) and (3) using the simplest possible model of onedimensional (1D) diffusion. While this approach lacks the complete generality afforded by the BlochTorrey description, it does preserve all important aspects of the physics involved and it does not require familiarity with higher mathematics beyond the binomial theorem and l'H τpital's rule.
The simplest diffusion model
This section reviews a connection between Gaussian diffusion and random walks under very simple circumstances. A comprehensive treatment can be found, e.g., in.[15] We shall consider a 1D diffusion of a particle, e.g., a hydrogen atom, along the x axis. At time t = 0, the atom will be located at a position x = 0 [Figure  1]. We will allow it freedom of movement of a very limited kind: during every constant time interval Δt , it may move either left or right, with equal probabilities p_{ l} = p_{ r} =½ . The amount of movement is restricted to a jump of a constant size Δx . It is obvious that the atom cannot get anywhere else than to positions x = 0, ±Δx , ± 2Δx , .... We want to know the probabilities that it will be found at any of these positions at times t = 0, Δt , 2Δt , ....
Let us monitor the atom's progress step by step. [Figure  1] shows that the atom is at position x = 0 at time t = 0 with certainty P (0, 0) = 1. Proceeding to time t = Δt , there is one possible pathway to the left and one to the right, a total of two possibilities. The next time step is only slightly more complicated. At time t = 2Δt , the atom can be at one of three positions. A single pathway leads to each of the two peripheral positions but there are two different pathways (that is, two possibilities) leading to the center position. The number of pathways to the center position is equal to the sum of pathways for the previous positions from where the center position could be reached, that is, 1+1 = 2. The total number of all possible pathways is four.
The pattern of possible pathways is built up further in every time step. Counting the possible pathways, we end up with a number pattern called the Pascal triangle. Interestingly, these numbers also represent multiplication factors in the binomial theorem:
A general case at time t = nΔt is depicted in the lower part of [Figure  1]. There are (n/k ) ways to get to a position placed k steps from the leftmost position. The total number of all pathways grows by a factor of 2 at each time step. We can verify this observation using the binomial theorem, Eq. (4), with a = b = 1:
The probability P ( k , n ) of finding an atom at a position placed k steps from the extreme left edge after n time steps is calculated as a number of pathways leading to it, divided by the total number of all possible pathways:
This is a special case of the binomial distribution (also called Bernoulli distribution). The reason why the simple leftright jump pattern represents a fairly adequate 1D diffusion model rests in the central limit theorem of the probability theory, which mandates that the binomial distribution asymptotically approaches the Gaussian distribution if n becomes very large. The resemblance can already be seen for n = 30 in [Figure  2].
Because the atoms move to the left and to the right with the same probability, they are equally likely to be on either side of the central position x = 0 at any time step. In other words, the mean value of the position is always
The mean position is therefore not very helpful in trying to estimate how far the atoms really moved because the positive and negative signs of the x coordinate simply cancel out in the averaging. Instead, we can calculate the mean square distance from the origin, which is always positive. For our simple diffusion model, it can be obtained by mathematical induction: n = 1 at time t = Δt and the distance from origin is either x = +Δx or x = Δx . The average square of distance is then simply
Moving on to time t =( n 1) t , we can calculate the effect of the n th time step:
It is clear that the mean square distance grows by the amount of Δx^{ 2} in every time step:
Index n represents progress of time and the mean square distance thus increases linearly with time. To make this relationship more explicit, we define a diffusion constant D as
and rewrite the equation for mean square distance in a more familiar form
This is the Einstein's famous formula for 1D diffusion in a "discrete form". However, it is equally valid for a "continuous" 1D diffusion. A transition to continuous time and space variables can be made by taking smaller and smaller time steps Δt and smaller and smaller spatial steps Δx while taking care to preserve the value of the diffusion constant. This is easily arranged by always choosing
The time interval Δt then approaches zero and n goes to infinity in such a way that
t = nΔt
still shows the correct diffusion time. In this way, a generally valid result for Gaussian diffusion is finally obtained:
The NMR diffusion signal
The principle of NMR diffusion measurement is depicted in [Figure  3]. Let us begin with a qualitative description of diffusion on the NMR signal. First, an RF pulse turns all the equilibrium magnetization M_{ 0} into the transverse plane, perpendicular to the main static magnetic field B_{ 0} . The magnetization vector will rotate around it at angular frequency ω=d φ/dt given by the Larmor equation:
where γ is the gyromagnetic ratio (γ =2π x 42 . 576 rad s^{ 1} T^{ 1} for a hydrogen proton).
After the excitation, a short and strong gradient G is applied along the x axis, changing the constant main field B_{ 0} to a spatially variable field B ( x ) = B_{ 0} + Gx . Larmor frequencies therefore become different at different places along the x axis. When the gradient is switched off again, some phase differences will have accumulated between the spins at different positions. At this point, we just wait for a period equal to the diffusion time t = nΔt . Then an opposite but otherwise identical gradient, G , is applied. If the atoms did not change their positions during the diffusion time, all the phase differences would be perfectly reversed and the magnetization would be fully restored (apart from the neglected relaxation effects). This is the situation exemplified by the center atom in [Figure  3]. However, if the spins move, the second gradient finds them at different locations than the first one and the phases will be reversed "incorrectly". The result is a phase dispersion in the measured sample and loss of signal when all the spins are eventually summed up to form the magnetization vector. Faster diffusion (larger D ) means that the spins have bigger chance to travel farther and therefore experience larger magnetic field changes due to diffusion gradients. This causes larger spread in the phases and therefore results in a smaller signal.
We can now proceed to a quantitative description and calculate the NMR signal using the simple diffusion model introduced earlier. The diffusion starts at time t = 0 from the position x = 0. At time t = t = nΔt , the atom has undergone a total of k steps to the right and k_ steps to the left:
k + k_= n.
Its position measured in Δx steps from x = 0 is therefore
kx = k  k _ = k  ( n  k )=2 k  n.
The probability of finding it at x = kx Δx is given by Eq. (5), using the steps to the right (index k ) to measure the position.
The atom located at x = x_{ 0} when the first gradient is switched on has its phase φ changed by
where we took advantage of the Larmor Eq. (11) and assumed that the gradient time δ was very short so that the atom did not significantly move during the gradient application. The effect of the main field (γB_{ 0}δ) was left out because it is the same for any atom regardless of its position, its movement or the gradient strengths. During the diffusion time t, the atom moves to a new position x_{ 0} + kxΔx . The second (negative) gradient then alters the phase again, this time by
Assuming the phase was zero in the beginning ( φ = 0 at t = 0), it is now, at t = τ,
Clearly, the result does not depend on the original position x_{ 0} and we can therefore safely return to the assumption
x_{ 0} = 0.
To compute the NMR signal, it is necessary to sum up all the phaseshifted magnetization vectors from all possible positions using the corresponding probabilities as weighting factors. If the equilibrium magnetization is M_{ 0} , the available signal S at t = τ is
where we substituted kx = 2 k  n for the phase function parameter.
This sum is best tackled using the Euler theorem
exp( iα) = cos(α)+ i sin(α).
Abbreviating
α = γGΔxδ
and employing again Eq. (4), we derive
The NMR signal undergoes various stages of amplification, filtering and other transformations. It is therefore not possible to measure M_{ 0} in absolute terms. We can remedy this unfortunate drawback of NMR by measuring the same sample twice: once without the diffusion gradients, to obtain unattenuated signal S_{ 0} , and once with them, to obtain signal S . The M_{ 0} term in the signal will stay the same but the attenuation term will disappear from S_{ 0} . We then calculate the ratio of the signals with and without the diffusion gradients:
The last step is to eliminate the Dx and express the result in terms of experimentally accessible variables. We substitute
where we introduced the socalled b value,
a quantity which depends on the spectrometer hardware and the pulse program controlling it but does not depend on the diffusion constant. Note that the expression for the b value becomes more complicated if various NMR sequence intricacies are taken into account. Most importantly, we assumed that the diffusion gradients are switched on for a negligible period of time in comparison with the diffusion time. However, Eq. (13) does capture the most important featuressquare dependency on the gradient moment Gδ and linear dependency on the diffusion time t.
To summarize, the NMR signal from nuclei following our simple diffusion model is attenuated due to phase randomization as
This is a "discrete" form of the StejskalTanner equation. Eq. (14) offers the possibility to measure the diffusion constant. This can be done by acquiring signals with and without diffusion gradients and calculating D from Eq. (14). Better still, one can obtain the signal many times with different b values and obtain D by a fitting procedure.
Continuous space and time
Almost everything is prepared for the continuous case in Eq. (14). It remains to make the discretization steps finer and finer to achieve, in the limit, a continuous change in space and time. For some constant diffusion time τ = nΔt , this is accomplished by simultaneously sending the Δt to zero and n to infinity, while keeping the D constant. The space discretization step Δx is always adjusted by Eq. (9) and is therefore also gradually diminishing towards zero.
In Eq. (14), this step refinement procedure simply means keeping b and D constant and calculating a limit for n→☆. The calculation will employ l'H τpital's rule for limits of the ^{0}/_{ 0} type.
We first simplify
and concentrate on the limit in the exponent:
The original limit can therefore be evaluated as
and used to generalize Eq. (14) for continuous space and time by substituting Eq. (15) for its right side:
This is the StejskalTanner formula, labeled as Eq. (1) in the introduction.
Practical application of Eq. (16) usually means taking its natural logarithm
Diffusion constant D is then extracted either directly from two experiments performed with and without the diffusion gradients or preferably by linear fitting of a series of signals acquired with different b values.
Even though we derived the StejskalTanner formula using a very simple 1D diffusion model, it is valid more generally and correctly describes NMR diffusion experiment in a homogeneous 3D environment as long as the diffusion behavior is Gaussian. Because this is not always the case, it is prudent to refer to the diffusion coefficient obtained from the StejskalTanner formula as "apparent diffusion coefficient".
In a 3D case, it is the diffusion gradient vector G = ( G_{ x} ,G_{ y} ,G_{ z} ) that determines the direction along which diffusion is measured. Diffusion perpendicular to the gradient vector does not alter the phase of magnetization in any way and is therefore invisible. Essentially then, NMR experiments measure a 1D diffusion along the diffusion gradient vector.
In an isotropic environment, application of a diffusion gradient along an arbitrary direction leads to the same signal attenuation predicted by the StejskalTanner formula. In an anisotropic environment, at least six independent 1D diffusion measurements are required to fully assess the six independent components of the diffusion tensor. Each of these measurements still obeys the StejskalTanner formula except that the diffusion coefficient D must be replaced by a scalar quantity
is a unit vector in the direction of the gradient, D is the diffusion tensor, and the dot signifies scalar multiplication. For example, in a coordinate system rotated to align u with the x axis, D˜ becomes the D_{ xx} component of the diffusion tensor D . Eq. (13) for the b value remains valid regardless of the gradient direction or the coordinate system rotation, as long as the gradient pulses are much shorter than the diffusion time.
Diffusion measurement pitfalls
In this last section, we endeavor to briefly point out some of the problems associated with NMR diffusion experiments. Particular emphasis is placed on brain diffusion measurements which have become very popular.
Subject motion
Because diffusion sequences are so sensitive to even minute movement of spins, macroscopic sample motion is the biggest enemy of diffusion experiments.[16] Multishot sequences are virtually impossible to use in a living subject because even a very small motion affects the signal phase, making it difficult to combine the individual echoes into a consistent dataset. The usual answer comes in the form of singleshot echo planar imaging (EPI) which acquires the entire image in about 100 ms. Although the EPI method ensures much better phase consistency with respect to the diffusion gradients, it brings a slew of its own problems: low signal to noise ratio, low resolution and numerous distortions. Not even EPI sequences are completely immune to macroscopic motion, especially when the motion has variable velocity in the form of linear acceleration or rotation.[17]
Imaging gradients
Although the diffusion gradient pair has dominant effect on the randomly walking spins, any other gradient present during the experiment also acts upon the magnetization phase and thus alters the diffusion signal. Besides the diffusion gradients, a typical MRI pulse sequence has a considerable number of other gradients used for slice selection, phase encoding, and readout. They can all act as additional and unwanted diffusion gradients. Their effect on the b value can be estimated, given an accurate pulse sequence timing diagram.[18] This is important whenever a precise quantitative diffusion measurement is to be undertaken.
Local gradients
Local magnetic field gradients due to inhomogeneities of magnetic susceptibility not only enhance the T_{ 2} * decay but also act as small spatially variable diffusion gradients. Unlike the pulse sequence gradients, the local gradients are usually not known in any detail and their precise effect therefore cannot be calculated. However, spin echo sequences decrease this type of diffusion weighting.[3] It can be shown that the b value associated with the local parasitic gradients diminishes with a number of spin echoes squared.[19] Furthermore, alternation of the diffusion gradient polarity in a multiecho sequence minimizes the effect of the background local gradients.[20]
Dependence on the diffusion time
The result of a diffusion experiment very much depends on the diffusion time τ.[13] For example, if the diffusion time was extremely short, most molecules would not collide with any cellular walls or other obstacles and the process would closely approximate free diffusion. Longer diffusion time gives the molecules more opportunity to explore the complex cellular environment and sample its geometric tortuosity, as well as to move between the intracellular and extracellular compartments. Dependency on the diffusion time means that one cannot reasonably mix results obtained with different diffusion times. This also implies that the diffusion time should always be reported because data interpretation is not possible without it. Unfortunately, limited maximum gradient strength and the need to avoid eddy currents often result in relatively long diffusion gradients distributed throughout the sequence. It is impossible to precisely define the diffusion time of such sequences. Depending on the exact sequence design, such experiments can become essentially qualitative in nature and difficult to compare with those obtained by another sequence. The results acquired with such sequences would be truly quantitative only if the diffusion was strictly Gaussian, that is, if the diffusion constant (or tensor) did not depend on the diffusion time for the whole range of the diffusion times involved (e.g., in an agar gel phantom).
Tissue compartments
The most common molecule used in NMR diffusion measurements is, by a large margin, water. Because water is ubiquitous in a living tissue, the detected signal typically comprises a mixture of components from the extracellular and intracellular compartments, generally of different diffusion properties.[21] There is a lively debate about their relative contributions under various circumstances, e.g., during ischemia.[22] The situation is considerably complicated by the exchange of water between these compartments.[23] The intracellular space may be compartmentalized even further due to the presence of organelles such as mitochondria.
Inhomogeneous voxels
The voxel size in most diffusion imaging experiments, particularly in human imaging, is quite large in comparison to many brain features. However, the StejskalTanner formula assumes constant D in the entire measured region. If the voxel contains both grey matter and white matter or if fiber tracts of different prevailing directions cross inside it, the assumption of homogeneity is violated. Similarly to the extracellular and intracellular contributions to the signal, the diffusion signal becomes a mixture of two or more compartments and the result will depend both on their individual properties and the rates of exchange between them.[24] However, the voxel partial volume effects can also be useful, aiding the diffusionbased tissue segmentation.[25] The degree to which diffusion is nonGaussian can be quantified with diffusional kurtosis imaging.[26] For example, if the voxel environment comprises a number of separate compartments, each with its own Gaussian diffusion, kurtosis imaging can establish the degree of diffusion coefficient variability inside this voxel.
Apportioning the DTI imaging time
Many brain regions, most notably the white matter tracts, are significantly anisotropic. As we discussed earlier, the StejskalTanner formulae can still be used in this case but the experiment has to be repeated with different gradient directions. There is a considerable body of literature on the optimum number and arrangement of the gradient vectors used.[9],[10],[11] However, one should keep in mind that the optimum choice depends to a large degree on the subject under study and thus varies from one type of experiment to another. For example, larger anisotropy may require more sampling directions. In general, given a finite available imaging time, the tradeoff is made by spending it on more signal averages, more gradient directions, more different b values, or more different diffusion times. There is no simple and universal answer to this puzzle.
Temperature and viscosity
Diffusion depends not only on the size of the diffusing molecule (its Stokes hydrodynamic radius) and on geometric tortuosity of the environment but also on temperature and viscosity. Viscosity in turn also depends on temperature. The overall effect is a faster diffusion at higher temperatures.[27],[28] This dependency has to be taken into account, e.g., when comparing human inviv o data with invitr o samples examined at room temperature.
Larger molecules
It is perfectly feasible to measure diffusion of other substances than water. Standard MR spectroscopic techniques can be used, with addition of diffusion gradients.[28] The attenuation of signal peaks corresponding to the individual metabolites still obeys the StejskalTanner formulae. However, smaller gyromagnetic ratios of these substances are reflected in smaller b values (see Eq. (13)). The only remedy is to use larger diffusion gradients. Furthermore, because efficient signal detection requires that the exponent bD does not become too small, larger and slower metabolites require even larger b values to compensate for their lower diffusion coefficient and are therefore particularly affected.
Scanner hardware
Diffusion sequences are very demanding on the scanner. The gradient set is required to perform very accurately and achieve perfect symmetry of the diffusion gradient pairs.[16] Equally important is their spatial linearity. Because the b value depends on a square of the gradient, a nonlinearity of 5% over the subject head will translate into 10% error in diffusion coefficient. Suppression of eddy currents is also essential. They frequently produce observable artefacts quite capable of ruining the diffusion experiment. Careful pulse sequence design[29] and various postprocessing corrections[30],[31] certainly help in this respect but the most important factor is the quality of the scanner hardware and its maintenance. Good performance on standard clinical sequences in no way guarantees that the same system can be successfully used for diffusion experiments.
The list presented above may create a somewhat pessimistic impression regarding the ability of the NMR diffusion method to obtain truly quantitative results, particularly in living nervous tissue. However, despite all the problems, much valuable and relevant information has been extracted by their employment. Clearly, good understanding of the StejskalTanner formulae is crucial to avoid suboptimal experimental designs and misinterpretations of results.
Acknowledgments   
We gratefully acknowledge generous help of Dr. S. Hrabetova who provided numerous valuable comments. This work was in part supported by National Institutes of Health Grant NS045797 to J. Hrabe.
References   
1.  Basser PJ, Mattiello J, LeBihan D. MR diffusion tensor spectroscopy and imaging. Biophys J 1994;66:25967. [PUBMED] [FULLTEXT] 
2.  Hahn E. Spin echoes. Phys Rev 1950;80:58094. 
3.  Carr E, Purcell E. Effects of diffusion on free precession in nuclear magnetic resonance experiments. Phys Rev 1954;94:6308. 
4.  Torrey HC. Bloch equations with diffusion terms. Phys Rev 1956;104:5636. 
5.  Stejskal EO, Tanner JE. Spin diffusion measurements: Spin echoes in the presence of a timedependent field gradient. J Chem Phys 1965;42:28892. 
6.  Tanner J, Stejskal E. Restricted selfdiffusion of protons in colloidal systems by the pulsedgradient, spin echo method. J Chem Phys 1968;49:176877. 
7.  Nicholson C, Phillips JM. Ion diffusion modified by tortuosity and volume fraction in the extracellular microenvironment of the rat cerebellum. J Physiol 1981;321:22557. [PUBMED] [FULLTEXT] 
8.  Hrabe J, Hrabetova S, Segeth K. A model of effective diffusion and tortuosity in the extracellular space of the brain. Biophys J 2004;87:160617. [PUBMED] [FULLTEXT] 
9.  Kingsley PB. Introduction to diffusion tensor imaging mathematics: Part I. Tensors, rotations and eigenvectors. Concepts Magn Reson Part A 2006;28:10122. 
10.  Kingsley PB. Introduction to diffusion tensor imaging mathematics: Part II. Anisotropy, diffusionweighting factors and gradient encoding schemes. Concepts Magn Reson Part A 2006;28A:12354. 
11.  Kingsley PB. Introduction to diffusion tensor imaging mathematics: Part III. Tensor calculation, noise, simulations and optimization. Concepts Magn Reson Part A 2006;28A:15579. 
12.  Lim KO, Helpern JA. Neuropsychiatric applications of DTIa review. NMR Biomed 2002;15:58793. [PUBMED] [FULLTEXT] 
13.  Callaghan PT. Principles of nuclear magnetic resonance microscopy. Oxford University Press: New York; 1991. 
14.  Haacke ME, Brown RW, Thompson MR, Venkatesan R. Magnetic resonance imaging: Physical principles and sequence design. John Wiley and Sons: 1999. 
15.  Chandrasekhar S. Stochastic problems in physics and astronomy. Rev Modern Physics 1943;15:189. 
16.  LeBihan D, Poupon C, Amadon A, Lethimonnier F. Artifacts and pitfalls in diffusion MRI. J Magn Reson Imaging 2006;24:47888. 
17.  Norris DG. Implications of bulk motion for diffusionweighted imaging experiments: Effects, mechanisms and solutions. J Magn Reson Imaging 2001;13:48695. [PUBMED] [FULLTEXT] 
18.  Mattiello J, Basser P, LeBihan D. Analytical expressions for the b matrix in NMR diffusion imaging and spectroscopy. J Magn Reson A 1994;108:13141. 
19.  Meiboom S, Gill D. Modified spinecho method for measuring nuclear relaxation times. Rev Sci Instruments 1958;29:68891. 
20.  Karlicek RF Jr, Lowe IJ. A modified pulsed gradient technique for measuring diffusion in the presence of large background gradients. J Magn Reson 1980;37:7591. 
21.  Kroenke CD, Ackerman JJ, Neil JJ. Magnetic resonance measurement of tetramethylammonium diffusion in rat brain: Comparison of magnetic resonance and ionophoresis in vivo diffusion measurements. Magn Reson Med 2003;50:71726. [PUBMED] [FULLTEXT] 
22.  Nicolay K, Braun KP, Graaf RA, Dijkhuizen RM, Kruiskamp MJ. Diffusion NMR spectroscopy. NMR Biomed 2001;14:94111. [PUBMED] [FULLTEXT] 
23.  Quirk JD, Bretthorst GL, Duong TQ, Snyder AZ, Springer CS, Ackerman JJ, et al . Equilibrium water exchange between the intraand extracellular spaces of mammalian brain. Magn Reson Med 2003;50:4939. 
24.  Assaf Y, Freidlin RZ, Rohde GK, Basser PJ. New modeling and experimental framework to characterize hindered and restricted water diffusion in brain white matter. Magn Reson Med 2004;52:96578. [PUBMED] [FULLTEXT] 
25.  Ardekani BA, Bappal A, D'Angelo D, Ashtari M, Lencz T, Szeszko PR, et al . Brain morphometry using diffusionweighted magnetic resonance imaging: Application to schizophrenia. Neuroreport 2005;16:14559. [PUBMED] [FULLTEXT] 
26.  Jensen JH, Helpern JA, Ramani A, Lu H, Kaczynski K. Diffusional kurtosis imaging: the quantification of nongaussian water diffusion by means of magnetic resonance imaging. Magn Reson Med 2005;53:143240. [PUBMED] [FULLTEXT] 
27.  LeBihan D. Diffusion and perfusion magnetic resonance imaging: Applications to functional MRI. Raven Press: New York; 1995. 
28.  de Graaf RA. In vivo NMR spectroscopy: Principles and techniques. John Wiley and Sons: Chichester; 2002. 
29.  Reese TG, Heid O, Weisskoff RM, Wedeen VJ. Reduction of eddy currentinduced distortion in diffusion MRI using a twicerefocused spin echo. Magn Reson Med 2003;49:17782. [PUBMED] [FULLTEXT] 
30.  Bodammer N, Kaufmann J, Kanowski M, Tempelmann C. Eddy current correction in diffusionweighted imaging using pairs of images acquired with opposite diffusion gradient polarity. Magn Reson Med 2004;51:18893. [PUBMED] [FULLTEXT] 
31.  Zhuang J, Hrabe J, Kangarlu A, Xu D, Bansal R, Branch CA, et al . Correction of eddycurrent distortions in diffusion tensor images using the known directions and strengths of diffusion gradients. J Magn Reson Imaging 2006;24:118893. [PUBMED] [FULLTEXT] 
Figures
[Figure  1], [Figure  2], [Figure  3]
This article has been cited by  1 
Diffusion tensor imaging in the human spinal cord: Development, limitations, and clinical applications 

 Bosma, R.L., Stroman, P.W.   Critical Reviews in Biomedical Engineering. 2012; 40(1): 120   [Pubmed]   2 
Évolution de l’imagerie par résonance magnétique 

 I. IdyPeretti   IRBM. 2009; 30(2): 53   [Pubmed]  [DOI]   3 
Advanced imaging in acute stroke managementPart II: Magnetic resonance imaging 

 Saini, M., Butcher, K.   Neurology India. 2009; 57(5): 550558   [Pubmed]   4 
Release profiles of encapsulated actives from colloidosomes sintered for various durations 

 Yow, H.N., Routh, A.F.   Langmuir. 2009; 25(1): 159166   [Pubmed]  


