Abstract. We consider the linear and non-linear enhancement of diffusion weighted magnetic resonance images (DW-MRI) to use contextual information in denoising and inferring fiber crossings. We describe the space of DW-MRI images in a moving frame of reference, attached to fiber fragments which allows for convection-diffusion along the fibers. Because of this approach, our method is naturally able to handle crossings in data. We will perform experiments showing the ability of the enhancement to infer information about crossing structures, even in diffusion tensor images (DTI) which are incapable of representing crossings themselves. We will present a novel non-linear enhancement technique which performs better than linear methods in areas around ventricles, thereby eliminating the need for additional preprocessing steps to segment out the ventricles. We pay special attention to the details of implementation of the various numeric schemes.

The question whether our brain pathways adhere to a geometric grid structure has been a popular topic of debate in the diffusion imaging and neuroscience society. Wedeen et al. (2012a b) proposed that the brain’s white matter is organized like parallel sheets of interwoven pathways. Catani et al. (2012) concluded that this grid pattern is most likely an artifact, resulting from methodological biases that cause the tractography pathways to cross in orthogonal angles. To date, ambiguities in the mathematical conditions for a sheet structure to exist (e.g. its relation to orthogonal angles) combined with the lack of extensive quantitative evidence have prevented wide acceptance of the hypothesis. In this work, we formalize the relevant terminology and recapitulate the condition for a sheet structure to exist. Note that this condition is not related to the presence or absence of orthogonal crossing fibers, and that sheet structure is defined formally as a surface formed by two sets of interwoven pathways intersecting at arbitrary angles within the surface. To quantify the existence of sheet structure, we present a novel framework to compute the sheet probability index (SPI), which reflects the presence of sheet structure in discrete orientation data (e.g. fiber peaks derived from diffusion MRI). With simulation experiments we investigate the effect of spatial resolution, curvature of the fiber pathways, and measurement noise on the ability to detect sheet structure. In real diffusion MRI data experiments we can identify various regions where the data supports sheet structure (high SPI values), but also areas where the data does not support sheet structure (low SPI values) or where no reliable conclusion can be drawn. Several areas with high SPI values were found to be consistent across subjects, across multiple data sets obtained with different scanners, resolutions, and degrees of diffusion weighting, and across various modeling techniques. Under the strong assumption that the diffusion MRI peaks reflect true axons, our results would therefore indicate that pathways do not form sheet structures at every crossing fiber region but instead at well-defined locations in the brain. With this framework, sheet structure location, extent, and orientation could potentially serve as new structural features of brain tissue. The proposed method can be extended to quantify sheet structure in directional data obtained with techniques other than diffusion MRI, which is essential for further validation.

No abstract

The hypothesis that brain pathways form 2D sheet-like structures layered in 3D as “pages of a book” has been a topic of debate in the recent literature. This hypothesis was mainly supported by a qualitative evaluation of “path neighborhoods” reconstructed with diffusion MRI (dMRI) tractography. Notwithstanding the potentially important implications of the sheet structure hypothesis for our understanding of brain structure and development, it is still considered controversial by many for lack of quantitative analysis. A means to quantify sheet structure is therefore necessary to reliably investigate its occurrence in the brain. Previous work has proposed the Lie bracket as a quantitative indicator of sheet structure, which could be computed by reconstructing path neighborhoods from the peak orientations of dMRI orientation density functions. Robust estimation of the Lie bracket, however, is challenging due to high noise levels and missing peak orientations. We propose a novel method to estimate the Lie bracket that does not involve the reconstruction of path neighborhoods with tractography. This method requires the computation of derivatives of the fiber peak orientations, for which we adopt an approach called normalized convolution. With simulations and experimental data we show that the new approach is more robust with respect to missing peaks and noise. We also demonstrate that the method is able to quantify to what extent sheet structure is supported for dMRI data of different species, acquired with different scanners, diffusion weightings, dMRI sampling schemes, and spatial resolutions. The proposed method can also be used with directional data derived from other techniques than dMRI, which will facilitate further validation of the existence of sheet structure.

No abstract