Fast computation of myelin maps from MRI T2 relaxation data using multicore CPU and graphics card parallelization. Issue 3 (27th February 2014)

• Record Type: Journal Article
• Authors: Yoo, Youngjin; Prasloski, Thomas; Vavasour, Irene; MacKay, Alexander; Traboulsee, Anthony L.; Li, David K.B.; Tam, Roger C.
• Is Part Of: Journal of magnetic resonance imaging. Volume 41:Issue 3(2015)
• Page Start: 700
• Page End: 707
• Publication Date: 2014-02-27
• Subjects: Magnetic resonance imaging -- Periodicals; 616
• DOI: 10.1002/jmri.24604
• Languages: English
• ISSNs: 1053-1807

Abstract

<abstract abstract-type="main"> <title> <x xml:space="preserve">Abstract</x> </title> <sec id="jmri24604-sec-0001" sec-type="section"> <title>Purpose</title> <p>To develop a fast algorithm for computing myelin maps from multiecho T<sub>2</sub> relaxation data using parallel computation with multicore CPUs and graphics processing units (GPUs).</p> </sec> <sec id="jmri24604-sec-0002" sec-type="section"> <title>Materials and Methods</title> <p>Using an existing MATLAB (MathWorks, Natick, MA) implementation with basic (nonalgorithm‐specific) parallelism as a guide, we developed a new version to perform the same computations but using C++ to optimize the hybrid utilization of multicore CPUs and GPUs, based on experimentation to determine which algorithmic components would benefit from CPU versus GPU parallelization. Using 32‐echo T<sub>2</sub> data of dimensions 256 × 256 × 7 from 17 multiple sclerosis patients and 18 healthy subjects, we compared the two methods in terms of speed, myelin values, and the ability to distinguish between the two patient groups using Student's <italic>t</italic>‐tests.</p> </sec> <sec id="jmri24604-sec-0003" sec-type="section"> <title>Results</title> <p>The new method was faster than the MATLAB implementation by 4.13 times for computing a single map and 14.36 times for batch‐processing 10 scans. The two methods produced very similar myelin values, with small and explainable differences that did not impact the ability to distinguish the two patient groups.</p> </sec> <sec id="jmri24604-sec-0004" sec-type="section"> <title>Conclusion</title> <p>The proposed hybrid multicore approach represents a more efficient alternative to MATLAB, especially for large‐scale batch processing. <bold>J. Magn. Reson. Imaging 2015;41:700–707.</bold> © <bold>2014 Wiley Periodicals, Inc.</bold></p> </sec> </abstract>