Oriented tensor reconstruction: tracing neural pathways from diffusion tensor MRI

In this paper we develop a new technique for tracing anatomical fibers from 3D tensor fields. The technique extracts salient tensor features using a local regularization technique that allows the algorithm to cross noisy regions and bridge gaps in the data. We applied the method to human brain DT-MRI data and recovered identifiable anatomical structures that correspond to the white matter brain-fiber pathways. The images in this paper are derived from a dataset having 121x88x60 resolution. We were able to recover fibers with less than the voxel size resolution by applying the regularization technique, i.e., using a priori assumptions about fiber smoothness. The regularization procedure is done through a moving least squares filter directly incorporated in the tracing algorithm

Probing white-matter microstructure with higher-order diffusion tensors and susceptibility tensor MRI.

Diffusion MRI has become an invaluable tool for studying white matter microstructure and brain connectivity. The emergence of quantitative susceptibility mapping and susceptibility tensor imaging (STI) has provided another unique tool for assessing the structure of white matter. In the highly ordered white matter structure, diffusion MRI measures hindered water mobility induced by various tissue and cell membranes, while susceptibility sensitizes to the molecular composition and axonal arrangement. Integrating these two methods may produce new insights into the complex physiology of white matter. In this study, we investigated the relationship between diffusion and magnetic susceptibility in the white matter. Experiments were conducted on phantoms and human brains in vivo. Diffusion properties were quantified with the diffusion tensor model and also with the higher order tensor model based on the cumulant expansion. Frequency shift and susceptibility tensor were measured with quantitative susceptibility mapping and susceptibility tensor imaging. These diffusion and susceptibility quantities were compared and correlated in regions of single fiber bundles and regions of multiple fiber orientations. Relationships were established with similarities and differences identified. It is believed that diffusion MRI and susceptibility MRI provide complementary information of the microstructure of white matter. Together, they allow a more complete assessment of healthy and diseased brains

Random noise in Diffusion Tensor Imaging, its Destructive Impact and Some Corrections

The empirical origin of random noise is described, its influence on DTI variables is illustrated by a review of numerical and in vivo studies supplemented by new simulations investigating high noise levels. A stochastic model of noise propagation is presented to structure noise impact in DTI. Finally, basics of voxelwise and spatial denoising procedures are presented. Recent denoising procedures are reviewed and consequences of the stochastic model for convenient denoising strategies are discussed

PDEs for tensor image processing

Methods based on partial differential equations (PDEs) belong to those image processing techniques that can be extended in a particularly elegant way to tensor fields. In this survey paper the most important PDEs for discontinuity-preserving denoising of tensor fields are reviewed such that the underlying design principles becomes evident. We consider isotropic and anisotropic diffusion filters and their corresponding variational methods, mean curvature motion, and selfsnakes. These filters preserve positive semidefiniteness of any positive semidefinite initial tensor field. Finally we discuss geodesic active contours for segmenting tensor fields. Experiments are presented that illustrate the behaviour of all these methods

Curvature-driven PDE methods for matrix-valued images

Matrix-valued data sets arise in a number of applications including diffusion tensor magnetic resonance imaging (DT-MRI) and physical measurements of anisotropic behaviour. Consequently, there arises the need to filter and segment such tensor fields. In order to detect edgelike structures in tensor fields, we first generalise Di Zenzo\u27s concept of a structure tensor for vector-valued images to tensor-valued data. This structure tensor allows us to extend scalar-valued mean curvature motion and self-snakes to the tensor setting. We present both two-dimensional and three-dimensional formulations, and we prove that these filters maintain positive semidefiniteness if the initial matrix data are positive semidefinite. We give an interpretation of tensorial mean curvature motion as a process for which the corresponding curve evolution of each generalised level line is the gradient descent of its total length. Moreover, we propose a geodesic active contour model for segmenting tensor fields and interpret it as a minimiser of a suitable energy functional with a metric induced by the tensor image. Since tensorial active contours incorporate information from all channels, they give a contour representation that is highly robust under noise. Experiments on three-dimensional DT-MRI data and an indefinite tensor field from fluid dynamics show that the proposed methods inherit the essential properties of their scalar-valued counterparts

Spatial Smoothing for Diffusion Tensor Imaging with low Signal to Noise Ratios

Though low signal to noise ratio (SNR) experiments in DTI give key information about tracking and anisotropy, e.g. by measurements with very small voxel sizes, due to the complicated impact of thermal noise such experiments are up to now seldom analysed. In this paper Monte Carlo simulations are presented which investigate the random fields of noise for different DTI variables in low SNR situations. Based on this study a strategy for spatial smoothing, which demands essentially uniform noise, is derived. To construct a convenient filter the weights of the nonlinear Aurich chain are adapted to DTI. This edge preserving three dimensional filter is then validated in different variants via a quasi realistic model and is applied to very new data with isotropic voxels of the size 1x1x1 mm3 which correspond to a spatial mean SNR of approximately 3

A Convex Semi-Definite Positive Framework for DTI Estimation and Regularization

International audienceIn this paper we introduce a novel variational method for joint estimation and regularization of diffusion tensor fields from noisy raw data. To this end, we use the classic quadratic data fidelity term derived from the Stejskal-Tanner equation with a new smoothness term leading to a convex objective function. The regularization term is based on the assumption that the signal can be reconstructed using a weighted average of observations on a local neighborhood. The weights measure the similarity between tensors and are computed directly from the diffusion images. We preserve the positive semi-definiteness constraint using a projected gradient descent. Experimental validation and comparisons with a similar method using synthetic data with known noise model, as well as classification of tensors towards understanding the anatomy of human skeletal muscle demonstrate the potential of our method

A generic approach to diffusion filtering of matrix-fields

Diffusion tensor magnetic resonance imaging (DT-MRI), is a image acquisition method, that provides matrix-valued data, so-called matrix fields. Hence image processing tools for the filtering and analysis of these data types are in demand. In this artricle we propose a generic framework that allows us to find the matrix-valued counterparts of the Perona-Malik PDEs with various diffusivity functions. To this end we extend the notion of derivatives and associated differential operators to matrix fields of symmetric matrices by adopting an operator-algebraic point of view. In order to solve these novel matrix-valued PDEs successfully we develop truly matrix-valued analogs to numerical solution schemes of the scalar setting. Numerical experiments performed on both synthetic and real world data substantiate the effectiveness of our novel matrix-valued Perona-Malik diffusion filters

DTI denoising for data with low signal to noise ratios

Low signal to noise ratio (SNR) experiments in diffusion tensor imaging (DTI) give key information about tracking and anisotropy, e. g., by measurements with small voxel sizes or with high b values. However, due to the complicated and dominating impact of thermal noise such data are still seldom analysed. In this paper Monte Carlo simulations are presented which investigate the distributions of noise for different DTI variables in low SNR situations. Based on this study a strategy for the application of spatial smoothing is derived. Optimal prerequisites for spatial filters are unbiased, bell shaped distributions with uniform variance, but, only few variables have a statistics close to that. To construct a convenient filter a chain of nonlinear Gaussian filters is adapted to peculiarities of DTI and a bias correction is introduced. This edge preserving three dimensional filter is then validated via a quasi realistic model. Further, it is shown that for small sample sizes the filter is as effective as a maximum likelihood estimator and produces reliable results down to a local SNR of approximately 1. The filter is finally applied to very recent data with isotropic voxels of the size 1×1×1mm^3 which corresponds to a spatially mean SNR of 2.5. This application demonstrates the statistical robustness of the filter method. Though the Rician noise model is only approximately realized in the data, the gain of information by spatial smoothing is considerable