TY - GEN
T1 - Accelerating global tractography using parallel Markov chain Monte Carlo
AU - Wu, Haiyong
AU - Chen, Geng
AU - Yang, Zhongxue
AU - Shen, Dinggang
AU - Yap, Pew Thian
N1 - Funding Information:
This work was supported in part by an Educational Science in Jiangsu Province grant (D201501125), a UNC BRIC-Radiology start-up fund, and NIH grants (EB006733, EB008374, EB009634, MH088520, AG041721, and MH100217).
PY - 2016
Y1 - 2016
N2 - Global tractography estimates brain connectivity by determining the optimal configuration of signal-generating fiber segments that best describes the measured diffusion-weighted data, promising better stability than local greedy methods with respect to imaging noise. However, global tractography is computationally very demanding and requires computation times that are often prohibitive for clinical applications. We present here a reformulation of the global tractography algorithm for fast parallel implementation amendable to acceleration using multicore CPUs and general-purpose GPUs. Our method is motivated by the key observation that each fiber segment is affected by a limited spatial neighborhood. That is, a fiber segment is influenced only by the fiber segments that are (or can potentially be) connected to its both ends and also by the diffusion-weighted signal in its proximity. This observation makes it possible to parallelize the Markov chain Monte Carlo (MCMC) algorithm used in the global tractography algorithm so that updating of independent fiber segments can be done concurrently. The experiments show that the proposed algorithm can significantly speed up global tractography, while at the same time maintain or improve tractography performance.
AB - Global tractography estimates brain connectivity by determining the optimal configuration of signal-generating fiber segments that best describes the measured diffusion-weighted data, promising better stability than local greedy methods with respect to imaging noise. However, global tractography is computationally very demanding and requires computation times that are often prohibitive for clinical applications. We present here a reformulation of the global tractography algorithm for fast parallel implementation amendable to acceleration using multicore CPUs and general-purpose GPUs. Our method is motivated by the key observation that each fiber segment is affected by a limited spatial neighborhood. That is, a fiber segment is influenced only by the fiber segments that are (or can potentially be) connected to its both ends and also by the diffusion-weighted signal in its proximity. This observation makes it possible to parallelize the Markov chain Monte Carlo (MCMC) algorithm used in the global tractography algorithm so that updating of independent fiber segments can be done concurrently. The experiments show that the proposed algorithm can significantly speed up global tractography, while at the same time maintain or improve tractography performance.
UR - http://www.scopus.com/inward/record.url?scp=84964066800&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84964066800&partnerID=8YFLogxK
U2 - 10.1007/978-3-319-28588-7_11
DO - 10.1007/978-3-319-28588-7_11
M3 - Conference contribution
AN - SCOPUS:84964066800
SN - 9783319285863
T3 - Mathematics and Visualization
SP - 121
EP - 130
BT - Computational Diffusion MRI - MICCAI Workshop, 2015
A2 - Rathi, Yogesh
A2 - Fuster, Andrea
A2 - Ghosh, Aurobrata
A2 - Kaden, Enrico
A2 - Reisert, Marco
PB - Springer Heidelberg
T2 - Workshop on Computational Diffusion MRI, MICCAI 2015
Y2 - 9 October 2015 through 9 October 2015
ER -