An Analytical Fiber ODF Reconstruction in 3D Polarized Light Imaging

International audienceThree dimensional polarized light imaging (3D-PLI) utilizes the birefringence in postmortem tissue to map its spatial fiber structure at a submillimeter resolution. We propose an analytical method to compute the fiber orientation distribution function (ODF) from high-resolution vector data provided by 3D-PLI. This strategy enables the bridging of high resolution 3D-PLI to diffusion magnetic resonance imaging with relatively low spatial resolution. First, the fiber ODF is modeled as a sum of K orientations on the unit sphere and expanded with a high order spherical harmonics series. Then, the coefficients of the spherical harmonics are derived directly with the spherical Fourier transform. We quantitatively validate the accuracy of the reconstruction against synthetic data and show that we can recover complex fiber configurations in the human heart at different scales

Insight into the fundamental trade-offs of diffusion MRI from polarization-sensitive optical coherence tomography in ex vivo human brain

In the first study comparing high angular resolution diffusion MRI (dMRI) in the human brain to axonal orientation measurements from polarization-sensitive optical coherence tomography (PSOCT), we compare the accuracy of orientation estimates from various dMRI sampling schemes and reconstruction methods. We find that, if the reconstruction approach is chosen carefully, single-shell dMRI data can yield the same accuracy as multi-shell data, and only moderately lower accuracy than a full Cartesian-grid sampling scheme. Our results suggest that current dMRI reconstruction approaches do not benefit substantially from ultra-high b-values or from very large numbers of diffusion-encoding directions. We also show that accuracy remains stable across dMRI voxel sizes of 1 ​mm or smaller but degrades at 2 ​mm, particularly in areas of complex white-matter architecture. We also show that, as the spatial resolution is reduced, axonal configurations in a dMRI voxel can no longer be modeled as a small set of distinct axon populations, violating an assumption that is sometimes made by dMRI reconstruction techniques. Our findings have implications for in vivo studies and illustrate the value of PSOCT as a source of ground-truth measurements of white-matter organization that does not suffer from the distortions typical of histological techniques.Published versio

Analytical Fiber ODF Reconstruction in 3D Polarized Light Imaging: Performance Assessment

International audienceThree dimensional Polarized Light Imaging (3D-PLI) allows to map the spatial fiber structure of postmortem tissue at a sub-millimeter resolution, thanks to its birefringence property. Different methods have been recently proposed to reconstruct the fiber orientation distribution function (fODF) from high-resolution vector data provided by 3D-PLI. Here, we focus on the analytical fODF computation approach, which uses the spherical harmonics to represent the fODF and analytically computes the spherical harmonics coefficients via the spherical Fourier transform. This work deals with the assessment of the performance of this approach on rich synthetic data which simulates the geometry of the neuronal fibers and on human brain dataset. A computational complexity and robustness to noise analysis demonstrate the interest and great potential of the approach

Towards validation of diffusion MRI tractography: bridging the resolution gap with 3D Polarized Light Imaging

International audienceThree-dimensional Polarized Light Imaging (3D-PLI) is an optical approach presented as a good candidate for validation of diffusion Magnetic Resonance Imaging (dMRI) results such as orientation estimates (fiber Orientation Distribution Functions) and tractography. We developed an anlytical approach to reconstruct fiber ODFs from 3D-PLI datasets. From these fODFs, here we compute brain fiber tracts via dMRI-based probabilistic tractography algorithm. Reconstructed fODFs at different scales proves the ability to bridge the resolution gap between 3D-PLI and dMRI, demonstrating, therefore, a great promise to validate diffusion MRI tractography thanks to multi-scale fiber tracking based on 3D-PLI

Analytical and fast Fiber Orientation Distribution reconstruction in 3D-Polarized Light Imaging

International audienceThree dimensional Polarized Light Imaging (3D-PLI) is an optical technique which allows mapping the spatial fiber architecture of fibrous postmortem tissues, at sub-millimeter resolutions. Here, we propose an analytical and fast approach to compute the fiber orientation distribution (FOD) from high-resolution vector data provided by 3D-PLI. The FOD is modeled as a sum of K orientations/Diracs on the unit sphere, described on a spherical harmonics basis and analytically computed using the spherical Fourier transform. Experiments are performed on rich synthetic data which simulate the geometry of the neuronal fibers and on human brain data. Results indicate the analytical FOD is computationally efficient and very fast, and has high angular precision and angular resolution. Furthermore, investigations on the right occipital lobe illustrate that our strategy of FOD computation enables the bridging of spatial scales from microscopic 3D-PLI information to macro-or mesoscopic dimensions of diffusion Magnetic Resonance Imaging (MRI), while being a means to evaluate prospective resolution limits for diffusion MRI to reconstruct regionspecific white matter tracts. These results demonstrate the interest and great potential of our analytical approach

Quantitative assessment of multi-scale tractography: bridging the resolution gap with 3D-PLI

International audienceThe in vivo validation of diffusion MRI (dMRI)-based tractography has beenshown to be a challenging task [Maier-hein et al.]. Therefore, we have been investigating how 3D Polarized Light Imaging (3D-PLI) could be used as a validation tool for dMRI-based fiber orientation estimation and tractography. PLI is an optical imaging technique that provides us with high-resolution fiber orientation measurements at micrometer scale. For this reason, it has been presented as a good candidate for the afore mentioned validation tasks [Axer et al,2011, Alimi et al, 2019 submitted]. In some previous works [alimi2017,18isbi,18ismrm,19,19submitted] we introduced an approach to close the resolution gap between dMRI and 3D-PLI. The study of the brain network from the topological point of view has seen an increasing interest in the last years [Sizemore et al, 2018, Chung et al, 2017]. In this work, we show how tractograms obtained at different spatial scales using 3D-PLI human brain datasets can bein spected using homology theory to perform a quantitative comparison between them. In particular, we investigate the persistence of the number of connected components in brain networks estimated from data at different resolutions

Doctor of Philosophy

dissertationDiffusion magnetic resonance imaging (dMRI) has become a popular technique to detect brain white matter structure. However, imaging noise, imaging artifacts, and modeling techniques, etc., create many uncertainties, which may generate misleading information for further analysis or applications, such as surgical planning. Therefore, how to analyze, effectively visualize, and reduce these uncertainties become very important research questions. In this dissertation, we present both rank-k decomposition and direct decomposition approaches based on spherical deconvolution to decompose the fiber directions more accurately for high angular resolution diffusion imaging (HARDI) data, which will reduce the uncertainties of the fiber directions. By applying volume rendering techniques to an ensemble of 3D orientation distribution function (ODF) glyphs, which we call SIP functions of diffusion shapes, one can elucidate the complex heteroscedastic structural variation in these local diffusion shapes. Furthermore, we quantify the extent of this variation by measuring the fraction of the volume of these shapes, which is consistent across all noise levels, the certain volume ratio. To better understand the uncertainties in white matter fiber tracks, we propose three metrics to quantify the differences between the results of diffusion tensor magnetic resonance imaging (DT-MRI) fiber tracking algorithms: the area between corresponding fibers of each bundle, the Earth Mover's Distance (EMD) between two fiber bundle volumes, and the current distance between two fiber bundle volumes. Based on these metrics, we discuss an interactive fiber track comparison visualization toolkit we have developed to visualize these uncertainties more efficiently. Physical phantoms, with high repeatability and reproducibility, are also designed with the hope of validating the dMRI techniques. In summary, this dissertation provides a better understanding about uncertainties in diffusion magnetic resonance imaging: where and how much are the uncertainties? How do we reduce these uncertainties? How can we possibly validate our algorithms

Voxel-wise comparisons of cellular microstructure and diffusion-MRI in mouse hippocampus using 3D Bridging of Optically-clear histology with Neuroimaging Data (3D-BOND)

A key challenge in medical imaging is determining a precise correspondence between image properties and tissue microstructure. This comparison is hindered by disparate scales and resolutions between medical imaging and histology. We present a new technique, 3D Bridging of Optically-clear histology with Neuroimaging Data (3D-BOND), for registering medical images with 3D histology to overcome these limitations. Ex vivo 120 × 120 × 200 μm resolution diffusion-MRI (dMRI) data was acquired at 7 T from adult C57Bl/6 mouse hippocampus. Tissue was then optically cleared using CLARITY and stained with cellular markers and confocal microscopy used to produce high-resolution images of the 3D-tissue microstructure. For each sample, a dense array of hippocampal landmarks was used to drive registration between upsampled dMRI data and the corresponding confocal images. The cell population in each MRI voxel was determined within hippocampal subregions and compared to MRI-derived metrics. 3D-BOND provided robust voxel-wise, cellular correlates of dMRI data. CA1 pyramidal and dentate gyrus granular layers had significantly different mean diffusivity (p > 0.001), which was related to microstructural features. Overall, mean and radial diffusivity correlated with cell and axon density and fractional anisotropy with astrocyte density, while apparent fibre density correlated negatively with axon density. Astrocytes, axons and blood vessels correlated to tensor orientation

Investigating the tradeoffs between spatial resolution and diffusion sampling for brain mapping with diffusion tractography: Time well spent?: Spatial vs. Q-Space Sampling for Tractography

Interest in mapping white matter pathways in the brain has peaked with the recognition that altered brain connectivity may contribute to a variety of neurologic and psychiatric diseases. Diffusion tractography has emerged as a popular method for postmortem brain mapping initiatives, including the ex-vivo component of the human connectome project, yet it remains unclear to what extent computer-generated tracks fully reflect the actual underlying anatomy. Of particular concern is the fact that diffusion tractography results vary widely depending on the choice of acquisition protocol. The two major acquisition variables that consume scan time, spatial resolution, and diffusion sampling, can each have profound effects on the resulting tractography. In this analysis we determined the effects of the temporal tradeoff between spatial resolution and diffusion sampling on tractography in the ex-vivo rhesus macaque brain, a close primate model for the human brain. We used the wealth of autoradiography-based connectivity data available for the rhesus macaque brain to assess the anatomic accuracy of six time-matched diffusion acquisition protocols with varying balance between spatial and diffusion sampling. We show that tractography results vary greatly, even when the subject and the total acquisition time are held constant. Further, we found that focusing on either spatial resolution or diffusion sampling at the expense of the other is counterproductive. A balanced consideration of both sampling domains produces the most anatomically accurate and consistent results