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. 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, [4] 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 properties when measured along different axes: physical or mechanical properties.

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

Principal component analysis (PCA) is a popular technique for analyzing large datasets containing a high number of dimensions/features per observation, increasing the interpretability of data while preserving the maximum amount of information, and enabling the visualization of multidimensional data. Formally, PCA is a statistical technique for reducing the dimensionality of a dataset. This is accomplished by linearly transforming the data into a new coordinate system where the variation in the data can be described with fewer dimensions than the initial data. Many studies use the first two principal components in order to plot the data in two dimensions and to visually identify clusters of closely related data points. Principal component analysis has applications in many fields such as population genetics, microbiome studies, and atmospheric science.

In linear algebra, a square matrix  is called diagonalizable or non-defective if it is similar to a diagonal matrix. That is, if there exists an invertible matrix  and a diagonal matrix such that . This is equivalent to . This property exists for any linear map: for a finite-dimensional vector space , a linear map  is called diagonalizable if there exists an ordered basis of  consisting of eigenvectors of . These definitions are equivalent: if  has a matrix representation as above, then the column vectors of  form a basis consisting of eigenvectors of , and the diagonal entries of  are the corresponding eigenvalues of ; with respect to this eigenvector basis,  is represented by .

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

Fourier optics is the study of classical optics using Fourier transforms (FTs), in which the waveform being considered is regarded as made up of a combination, or superposition, of plane waves. It has some parallels to the Huygens–Fresnel principle, in which the wavefront is regarded as being made up of a combination of spherical wavefronts whose sum is the wavefront being studied. A key difference is that Fourier optics considers the plane waves to be natural modes of the propagation medium, as opposed to Huygens–Fresnel, where the spherical waves originate in the physical medium.

A directional derivative is a concept in multivariable calculus that measures the rate at which a function changes in a particular direction at a given point.

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">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 mathematics, the discrete Laplace operator is an analog of the continuous Laplace operator, defined so that it has meaning on a graph or a discrete grid. For the case of a finite-dimensional graph, the discrete Laplace operator is more commonly called the Laplacian matrix.

In linear algebra, it is often important to know which vectors have their directions unchanged by a linear transformation. An eigenvector or characteristic vector is such a vector. Thus an eigenvector of a linear transformation is scaled by a constant factor when the linear transformation is applied to it: . 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.

In the mathematics of evolving systems, the concept of a center manifold was originally developed to determine stability of degenerate equilibria. Subsequently, the concept of center manifolds was realised to be fundamental to mathematical modelling.

The expander mixing lemma intuitively states that the edges of certain -regular graphs are evenly distributed throughout the graph. In particular, the number of edges between two vertex subsets and is always close to the expected number of edges between them in a random -regular graph, namely .

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 respect 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.

<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.

Group actions are central to Riemannian geometry and defining orbits. The orbits of computational anatomy consist of anatomical shapes and medical images; the anatomical shapes are submanifolds of differential geometry consisting of points, curves, surfaces and subvolumes,. This generalized the ideas of the more familiar orbits of linear algebra which are linear vector spaces. Medical images are scalar and tensor images from medical imaging. The group actions are used to define models of human shape which accommodate variation. These orbits are deformable templates as originally formulated more abstractly in pattern theory.

<span class="mw-page-title-main">MRI sequence</span>

An MRI 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–32. doi:10.1002/1522-2594(200010)44:4<625::AID-MRM17>3.0.CO;2-O. PMID   11025519.
  2. Catani, Marco; Thiebaut De Schotten, Michel (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. Wade, Ryckie G.; Tanner, Steven F.; Teh, Irvin; Ridgway, John P.; Shelley, David; Chaka, Brian; Rankine, James J.; Andersson, Gustav; Wiberg, Mikael; Bourke, Grainne (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.