Tuesday, June 27, 2017: 12:45 PM
Monday, June 26 & Tuesday, June 27
Divya Varadarajan1, Justin Haldar1
1University of Southern California, Los Angeles, USA
White matter tractography is a powerful tool for medicine and neuroscience, but fiber tracking accuracy is often confounded by crossing fibers. Despite much progress in the development of advanced q-space sampling [1-3] and orientation distribution function (ODF) estimation strategies [4,5], the design of an optimal ODF measurement approach still remains a longstanding unsolved problem. We address this issue by introducing a new model-free theoretical framework for characterizing resolution in ODF estimation. This framework can be used to compare and evaluate different q-space sampling schemes and ODF estimation methods, without requiring modeling assumptions or extensive empirical evaluations. Building off of previous work [4,6,7], our approach relies on a novel theoretical relationship between the estimated ODF and the true ensemble average propagator (EAP), a probability distribution describing local molecular diffusion.
The ODF with constant solid angle correction , is defined as
where u is an orientation unit vector, and fΔ(x) is the EAP.
Extending previous work  and neglecting noise for simplicity (noise is easily included), we can derive that any linearly estimated ODF Ô(u) is related to the EAP by:
where g(u,x) is the "EAP response", and is easily evaluated using similar techniques to . The EAP response can be interpreted as a point spread function (PSF) for ODF estimation that depends only on the q-space sampling scheme and the ODF estimation method. We hypothesize that measures like the main-lobe width of the EAP response will indicate resolution (like with conventional PSFs) and strongly correlate with empirical measures of angular resolution. If true, then it becomes possible to optimize data sampling and ODF estimation based on the EAP response.
To test our hypothesis, we calculated EAP responses for two ODF estimation schemes (FRACT  and SHORE ) and several multi-shell sampling schemes (b=[1000,2000,3000], [2000,3000,4000], [3000,4000,5000] and [4000,5000,6000] s/mm^2). We used a single-shell (b=3000 s/mm^2) acquisition for FRACT and same diffusion orientations as the human connectome project (HCP) 3-shell protocol . Main-lobe width was calculated as the full-width at half maximum at the peak of the EAP response along the axis perpendicular to u.
Empirical simulations and real HCP data were used to assess the correlation between the EAP response characteristics and the empirically-observed ODF resolution. ODF estimation and visualization were implemented using BrainSuite  (http://brainsuite.org/). We quantified empirical ODF resolution by simulating multiple voxels, each containing two diffusion tensors with equal volume fraction and diffusion coefficients in the range 1-3 mm^2/s, and measuring the minimum angle of separation (MAS) between the tensors at which the two distinct orientations are unresolvable.
Fig. 1 shows that both MAS and EAP response main-lobe width were correlated, suggesting that the EAP response predicts angular resolution. Fig. 2 shows the EAP response of FRACT and SHORE for the HCP data. FRACT has a lower main-lobe width, predicting higher angular resolution than SHORE for this data. This prediction is consistent with qualitative examination of the ODFs. Note that our goal in this comparison was to demonstrate the usefulness of the EAP response -- FRACT is not always better than SHORE, and the performance difference is specific to the q-space sampling scheme and the set of estimation parameters that were used.
We have proposed using EAP responses to assess angular resolution in linear ODF measurement. We demonstrated that the proposed metric correlated with empirical angular resolution, and believe that such a theoretical approach will prove useful for improving ODF measurement methods.
Diffusion MRI 2
Modeling and Analysis Methods:
Diffusion MRI Modeling and Analysis 1
Exploratory Modeling and Artifact Removal
Poster Session - Tuesday
Design and Analysis
WHITE MATTER IMAGING - DTI, HARDI, DSI, ETC
1|2Indicates the priority used for review
Would you accept an oral presentation if your abstract is selected for an oral session?
I would be willing to discuss my abstract with members of the press should my abstract be marked newsworthy:
Please indicate below if your study was a "resting state" or "task-activation” study.
By submitting your proposal, you grant permission for the Organization for Human Brain Mapping (OHBM) to distribute the presentation in any format, including video, audio print and electronic text through OHBM OnDemand, social media channels or other electronic media and on the OHBM website.
Healthy subjects only or patients (note that patient studies may also involve healthy subjects):
Internal Review Board (IRB) or Animal Use and Care Committee (AUCC) Approval. Please indicate approval below. Please note: Failure to have IRB or AUCC approval, if applicable will lead to automatic rejection of abstract.
Please indicate which methods were used in your research:
Which processing packages did you use for your study?
Other, Please list
Provide references in author date format
 Tuch, D. S.; Reese, T. G.; Wiegell, M. R.; Makris, N.; Belliveau, J. W. & Wedeen, V. J. (2002), ‘High angular resolution diffusion imaging reveals intravoxel white matter fiber heterogeneity’, Magn. Reson. Med., 48, 577-582
 Descoteaux, M.; Deriche, R.; LeBihan, D.; Mangin, J.-F. & Poupon, C.(2011), ‘Multiple q-shell diffusion propagator imaging’, Med. Image Anal., Elsevier, 15, 603-621
 Wedeen, V. J.; Hagmann, P.; Tseng, W.-Y. I.; Reese, T. G. & Weisskoff, R. M.( 2005), ‘Mapping complex tissue architecture with diffusion spectrum magnetic resonance imaging’, Magn. Reson. Med., 54, 1377-1386
 Tuch, D. S.(2004), ‘Q-Ball imaging’, Magn. Reson. Med., 52, 1358-1372
 Tournier, J.-D.; Calamante, F.; Gadian, D. G. & Connelly, A.(2005), ‘Direct estimation of the fiber orientation density function from diffusion-weighted MRI data using spherical deconvolution’, NeuroImage, 2004, 23, 1176-1185
 Haldar, J. P. & Leahy, R. M.(2013), ‘Linear transforms for Fourier data on the sphere: Application to high angular resolution diffusion MRI of the brain’, NeuroImage, 71, 233-247
 Varadarajan, D. & Haldar, J. P. (2016), ‘A theoretical framework for sampling and reconstructing ensemble average propagators in diffusion MRI’, Proc. Int. Soc. Magn. Reson. Med., 204
 Özarslan, E.; Koay, C.; Shepherd, T. M.; Blackband, S. J. & Basser, P. J. (2009), ‘Simple harmonic oscillator based reconstruction and estimation for three-dimensional q-space’, MRI Proc. Int. Soc. Magn. Reson. Med., 1396
 Van Essen, D. C.; Smith, S. M.; Barch, D. M.; Behrens, T. E.; Yacoub, E.; Ugurbil, K.; & WU-Minn HCP Consortium. (2013), ‘The WU-Minn human connectome project: an overview’, Neuroimage, 80, 62-79.
 Shattuck, D. W. & Leahy, R. M. BrainSuite: An automated cortical surface identification tool Med. Image Anal., 2002, 6, 129-142