Tractography

Last updated
Tractography
Tractography animated lateral view.gif
Tractography of human brain
Purposeused to visually represent nerve tracts

In neuroscience, tractography is a 3D modeling technique used to visually represent nerve tracts using data collected by diffusion MRI. [1] It uses special techniques of magnetic resonance imaging (MRI) and computer-based diffusion MRI. The results are presented in two- and three-dimensional images called tractograms. [2]

Contents

In addition to the long tracts that connect the brain to the rest of the body, there are complicated neural circuits formed by short connections among different cortical and subcortical regions. The existence of these tracts and circuits has been revealed by histochemistry and biological techniques on post-mortem specimens. Nerve tracts are not identifiable by direct exam, CT, or MRI scans. This difficulty explains the paucity of their description in neuroanatomy atlases and the poor understanding of their functions.

The most advanced tractography algorithm can produce 90% of the ground truth bundles, but it still contains a substantial amount of invalid results. [3]

MRI technique

DTI of the brachial plexus - see https://doi.org/10.3389/fsurg.2020.00019 for more information Deterministic Tractography of the Adult Brachial Plexus using Diffusion Tensor Imaging.gif
DTI of the brachial plexus - see https://doi.org/10.3389/fsurg.2020.00019 for more information
Tractographic reconstruction of neural connections by diffusion tensor imaging (DTI) DTI-sagittal-fibers.jpg
Tractographic reconstruction of neural connections by diffusion tensor imaging (DTI)
MRI tractography of the human subthalamic nucleus

Tractography is performed using data from diffusion MRI. The free water diffusion is termed "isotropic" diffusion. If the water diffuses in a medium with barriers, the diffusion will be uneven, which is termed anisotropic diffusion. In such a case, the relative mobility of the molecules from the origin has a shape different from a sphere. This shape is often modeled as an ellipsoid, and the technique is then called diffusion tensor imaging [4] . Barriers can be many things: cell membranes, axons, myelin, etc.; but in white matter the principal barrier is the myelin sheath of axons. Bundles of axons provide a barrier to perpendicular diffusion and a path for parallel diffusion along the orientation of the fibers.

Anisotropic diffusion is expected to be increased in areas of high mature axonal order. Conditions where the myelin or the structure of the axon are disrupted, such as trauma, [5] tumors, and inflammation reduce anisotropy, as the barriers are affected by destruction or disorganization.

Anisotropy is measured in several ways. One way is by a ratio called fractional anisotropy (FA). An FA of 0 corresponds to a perfect sphere, whereas 1 is an ideal linear diffusion. Few regions have FA larger than 0.90. The number gives information about how aspherical the diffusion is but says nothing of the direction.

Each anisotropy is linked to an orientation of the predominant axis (predominant direction of the diffusion). Post-processing programs are able to extract this directional information.

This additional information is difficult to represent on 2D grey-scaled images. To overcome this problem, a color code is introduced. Basic colors can tell the observer how the fibers are oriented in a 3D coordinate system, this is termed an "anisotropic map". The software could encode the colors in this way:

The technique is unable to discriminate the "positive" or "negative" direction in the same axis.

Mathematics

Using diffusion tensor MRI, one can measure the apparent diffusion coefficient at each voxel in the image, and after multilinear regression across multiple images, the whole diffusion tensor can be reconstructed. [1]

Suppose there is a fiber tract of interest in the sample. Following the Frenet–Serret formulas, we can formulate the space-path of the fiber tract as a parameterized curve:

where is the tangent vector of the curve. The reconstructed diffusion tensor can be treated as a matrix, and we can compute its eigenvalues and eigenvectors . By equating the eigenvector corresponding to the largest eigenvalue with the direction of the curve:

we can solve for given the data for . This can be done using numerical integration, e.g., using Runge–Kutta, and by interpolating the principal eigenvectors.

See also

Related Research Articles

<span class="mw-page-title-main">Anisotropy</span> In geometry, property of being directionally dependent

Anisotropy is the structural property of non-uniformity in different directions, as opposed to isotropy. An anisotropic object or pattern has properties that differ according to direction of measurement. For example, many materials exhibit very different physical or mechanical properties when measured along different axes, e.g. absorbance, refractive index, conductivity, and tensile strength.

<span class="mw-page-title-main">Principal component analysis</span> Method of data analysis

Principal component analysis (PCA) is a linear dimensionality reduction technique with applications in exploratory data analysis, visualization and data preprocessing.

In mathematics, a Hermitian matrix is a complex square matrix that is equal to its own conjugate transpose—that is, the element in the i-th row and j-th column is equal to the complex conjugate of the element in the j-th row and i-th column, for all indices i and j:

<span class="mw-page-title-main">Hooke's law</span> Physical law: force needed to deform a spring scales linearly with distance

In physics, Hooke's law is an empirical law which states that the force needed to extend or compress a spring by some distance scales linearly with respect to that distance—that is, Fs = kx, where k is a constant factor characteristic of the spring, and x is small compared to the total possible deformation of the spring. The law is named after 17th-century British physicist Robert Hooke. He first stated the law in 1676 as a Latin anagram. He published the solution of his anagram in 1678 as: ut tensio, sic vis. Hooke states in the 1678 work that he was aware of the law since 1660.

Linear elasticity is a mathematical model of how solid objects deform and become internally stressed by prescribed loading conditions. It is a simplification of the more general nonlinear theory of elasticity and a branch of continuum mechanics.

The diffusion equation is a parabolic partial differential equation. In physics, it describes the macroscopic behavior of many micro-particles in Brownian motion, resulting from the random movements and collisions of the particles. In mathematics, it is related to Markov processes, such as random walks, and applied in many other fields, such as materials science, information theory, and biophysics. The diffusion equation is a special case of the convection–diffusion equation when bulk velocity is zero. It is equivalent to the heat equation under some circumstances.

<span class="mw-page-title-main">Curvilinear coordinates</span> Coordinate system whose directions vary in space

In geometry, curvilinear coordinates are a coordinate system for Euclidean space in which the coordinate lines may be curved. These coordinates may be derived from a set of Cartesian coordinates by using a transformation that is locally invertible at each point. This means that one can convert a point given in a Cartesian coordinate system to its curvilinear coordinates and back. The name curvilinear coordinates, coined by the French mathematician Lamé, derives from the fact that the coordinate surfaces of the curvilinear systems are curved.

<span class="mw-page-title-main">Arcuate fasciculus</span> Neural pathway connecting Brocas area and Wernickes area

In neuroanatomy, the arcuate fasciculus is a bundle of axons that generally connects the Broca's area and the Wernicke's area in the brain. It is an association fiber tract connecting caudal temporal cortex and inferior frontal lobe.

In linear algebra, a generalized eigenvector of an matrix is a vector which satisfies certain criteria which are more relaxed than those for an (ordinary) eigenvector.

In linear algebra, an eigenvector or characteristic vector is a vector that has its direction unchanged by a given linear transformation. More precisely, an eigenvector of a linear transformation is scaled by a constant factor when the linear transformation is applied to it: . It is often important to know these vectors in linear algebra. The corresponding eigenvalue, characteristic value, or characteristic root is the multiplying factor .

<span class="mw-page-title-main">Diffusion MRI</span> Method of utilizing water in magnetic resonance imaging

Diffusion-weighted magnetic resonance imaging is the use of specific MRI sequences as well as software that generates images from the resulting data that uses the diffusion of water molecules to generate contrast in MR images. It allows the mapping of the diffusion process of molecules, mainly water, in biological tissues, in vivo and non-invasively. Molecular diffusion in tissues is not random, but reflects interactions with many obstacles, such as macromolecules, fibers, and membranes. Water molecule diffusion patterns can therefore reveal microscopic details about tissue architecture, either normal or in a diseased state. A special kind of DWI, diffusion tensor imaging (DTI), has been used extensively to map white matter tractography in the brain.

<span class="mw-page-title-main">Maxwell stress tensor</span> Mathematical description in electromagnetism

The Maxwell stress tensor is a symmetric second-order tensor used in classical electromagnetism to represent the interaction between electromagnetic forces and mechanical momentum. In simple situations, such as a point charge moving freely in a homogeneous magnetic field, it is easy to calculate the forces on the charge from the Lorentz force law. When the situation becomes more complicated, this ordinary procedure can become impractically difficult, with equations spanning multiple lines. It is therefore convenient to collect many of these terms in the Maxwell stress tensor, and to use tensor arithmetic to find the answer to the problem at hand.

In mathematics, the structure tensor, also referred to as the second-moment matrix, is a matrix derived from the gradient of a function. It describes the distribution of the gradient in a specified neighborhood around a point and makes the information invariant to the observing coordinates. The structure tensor is often used in image processing and computer vision.

In linear algebra, eigendecomposition is the factorization of a matrix into a canonical form, whereby the matrix is represented in terms of its eigenvalues and eigenvectors. Only diagonalizable matrices can be factorized in this way. When the matrix being factorized is a normal or real symmetric matrix, the decomposition is called "spectral decomposition", derived from the spectral theorem.

In the fields of computer vision and image analysis, the Harris affine region detector belongs to the category of feature detection. Feature detection is a preprocessing step of several algorithms that rely on identifying characteristic points or interest points so to make correspondences between images, recognize textures, categorize objects or build panoramas.

Fractional anisotropy (FA) is a scalar value between zero and one that describes the degree of anisotropy of a diffusion process. A value of zero means that diffusion is isotropic, i.e. it is unrestricted (or equally restricted) in all directions. A value of one means that diffusion occurs only along one axis and is fully restricted along all other directions. FA is a measure often used in diffusion imaging where it is thought to reflect fiber density, axonal diameter, and myelination in white matter. The FA is an extension of the concept of eccentricity of conic sections in 3 dimensions, normalized to the unit range.

The following outline is provided as an overview of and topical guide to brain mapping:

<span class="mw-page-title-main">Symmetry in quantum mechanics</span> Properties underlying modern physics

Symmetries in quantum mechanics describe features of spacetime and particles which are unchanged under some transformation, in the context of quantum mechanics, relativistic quantum mechanics and quantum field theory, and with applications in the mathematical formulation of the standard model and condensed matter physics. In general, symmetry in physics, invariance, and conservation laws, are fundamentally important constraints for formulating physical theories and models. In practice, they are powerful methods for solving problems and predicting what can happen. While conservation laws do not always give the answer to the problem directly, they form the correct constraints and the first steps to solving a multitude of problems. In application, understanding symmetries can also provide insights on the eigenstates that can be expected. For example, the existence of degenerate states can be inferred by the presence of non commuting symmetry operators or that the non degenerate states are also eigenvectors of symmetry operators.

<span class="mw-page-title-main">High-definition fiber tracking</span>

High definition fiber tracking (HDFT) is a tractography technique where data from MRI scanners is processed through computer algorithms to reveal the detailed wiring of the brain and to pinpoint fiber tracts. Each tract contains millions of neuronal connections. HDFT is based on data acquired from diffusion spectrum imaging and processed by generalized q-sampling imaging. The technique makes it possible to virtually dissect 40 major fiber tracts in the brain. The HDFT scan is consistent with brain anatomy unlike diffusion tensor imaging (DTI). Thus, the use of HDFT is essential in pinpointing damaged neural connections.

<span class="mw-page-title-main">MRI pulse sequence</span> A pulse sequence during a medical test

An MRI pulse sequence in magnetic resonance imaging (MRI) is a particular setting of pulse sequences and pulsed field gradients, resulting in a particular image appearance.

References

  1. 1 2 Basser PJ, Pajevic S, Pierpaoli C, Duda J, Aldroubi A (October 2000). "In vivo fiber tractography using DT-MRI data". Magnetic Resonance in Medicine. 44 (4): 625–632. doi:10.1002/1522-2594(200010)44:4<625::AID-MRM17>3.0.CO;2-O. PMID   11025519.
  2. Catani M, Thiebaut De Schotten M (2012). Atlas of Human Brain Connections. doi:10.1093/med/9780199541164.001.0001. ISBN   978-0-19-954116-4.[ page needed ]
  3. Maier-Hein KH, Neher PF, Houde JC, Côté MA, Garyfallidis E, Zhong J, et al. (November 2017). "The challenge of mapping the human connectome based on diffusion tractography". Nature Communications. 8 (1): 1349. Bibcode:2017NatCo...8.1349M. doi:10.1038/s41467-017-01285-x. PMC   5677006 . PMID   29116093.
  4. Basser PJ, Mattiello J, LeBihan D (January 1994). "MR diffusion tensor spectroscopy and imaging". Biophysical Journal. 66 (1): 259–267. Bibcode:1994BpJ....66..259B. doi:10.1016/S0006-3495(94)80775-1. PMC   1275686 . PMID   8130344.
  5. Wade RG, Tanner SF, Teh I, Ridgway JP, Shelley D, Chaka B, et al. (16 April 2020). "Diffusion Tensor Imaging for Diagnosing Root Avulsions in Traumatic Adult Brachial Plexus Injuries: A Proof-of-Concept Study". Frontiers in Surgery. 7: 19. doi: 10.3389/fsurg.2020.00019 . PMC   7177010 . PMID   32373625.