Fractional anisotropy

A Diffusion Ellipsoid is completely represented by the Diffusion Tensor, D. FA is calculated from the eigenvalues (${\displaystyle \lambda _{1},\lambda _{2},\lambda _{3}}$) of the diffusion tensor.[1] The eigenvectors ${\displaystyle \epsilon }$ give the directions in which the ellipsoid has major axes, and the corresponding eigenvalues ${\displaystyle \lambda }$ give the magnitude of the peak in that direction.

Diffusion Tensor Schematic

${\displaystyle {\text{FA}}={\sqrt {\frac {3}{2}}}{\frac {\sqrt {(\lambda _{1}-{\hat {\lambda }})^{2}+(\lambda _{2}-{\hat {\lambda }})^{2}+(\lambda _{3}-{\hat {\lambda }})^{2}}}{\sqrt {\lambda _{1}^{2}+\lambda _{2}^{2}+\lambda _{3}^{2}}}}}$

with ${\displaystyle {\hat {\lambda }}=(\lambda _{1}+\lambda _{2}+\lambda _{3})/3}$ being the mean value of the eigenvalues.

An equivalent formula for FA is

${\displaystyle {\text{FA}}={\sqrt {\frac {1}{2}}}{\frac {\sqrt {(\lambda _{1}-\lambda _{2})^{2}+(\lambda _{2}-\lambda _{3})^{2}+(\lambda _{3}-\lambda _{1})^{2}}}{\sqrt {\lambda _{1}^{2}+\lambda _{2}^{2}+\lambda _{3}^{2}}}},}$

which is further equivalent to:[2]

${\displaystyle {\text{FA}}={\sqrt {{\frac {1}{2}}\left(3-{\frac {1}{{\text{trace}}({\text{R}}^{2})}}\right)}}}$

where R is the "normalized" diffusion tensor:

${\displaystyle {\text{R}}={\frac {\text{D}}{{\text{trace}}({\text{D}})}}}$

Note that if all the eigenvalues are equal, which happens for isotropic (spherical) diffusion, as in free water, the FA is 0. The FA can reach a maximum value of 1 (this rarely happens in real data), in which case D has only one nonzero eigenvalue and the ellipsoid reduces to a line in the direction of that eigenvector. This means that the diffusion is confined to that direction alone.

This can be visualized with an ellipsoid, which is defined by the eigenvectors and eigenvalues of D. The FA of a sphere is 0 since the diffusion is isotropic, and there is equal probability of diffusion in all directions. The eigenvectors and eigenvalues of the Diffusion Tensor give a complete representation of the diffusion process. FA quantifies the pointedness of the ellipsoid, but does not give information about which direction it is pointing to.

Note that the FA of most liquids, including water, is 0 unless the diffusion process is being constrained by structures such as network of fibers. The measured FA may depend on the effective length scale of the diffusion measurement. If the diffusion process is not constrained on the scale being measured (the constraints are too far apart) or the constraints switch direction on a smaller scale than the measured one, then the measured FA will be attenuated. For example, the brain can be thought of as a fluid permeated by many fibers (nerve axons). However, in most parts the fibers go in all directions, and thus although they constrain the diffusion the FA is 0. In some regions, such as the corpus callosum the fibers are aligned over a large enough scale (on the order of a mm) for their directions to mostly agree within the resolution element of a magnetic resonance image, and it is these regions that stand out in an FA image. Liquid crystals can also exhibit anisotropic diffusion because the needle or plate-like shapes of their molecules affect how they slide over one another. When the FA is 0 the tensor nature of D is often ignored, and it is called the diffusion constant.

• 
  FA value of 0.7698, the DT matrix is diagonal([10 2 2])
• 
  FA value of 0, the DT matrix is diagonal([4 4 4])
• 
  FA value of 0.6030, the DT matrix is diagonal([4 4 2])

One drawback of the Diffusion Tensor model is that it can account only for Gaussian diffusion processes, which has been found to be inadequate in accurately representing the true diffusion process in the human brain. Due to this, higher order models using spherical harmonics and Orientation Distribution Functions (ODF) have been used to define newer and richer estimates of the anisotropy, called Generalized Fractional Anisotropy. GFA computations use samples of the ODF to evaluate the anisotropy in diffusion. They can also be easily calculated by using the Spherical Harmonic coefficients of the ODF model.[3]