GPU-Accelerated Reconstruction of T2 Maps in Magnetic Resonance Imaging

Open access


The main tissue parameters targeted by MR tomography include, among others, relaxation times T1 and T2. This paper focuses on the computation of the relaxation time T2 measured with the Spin Echo method, where the sensing coil of the tomograph provides a multi-echo signal. The maxima of these echoes must be interleaved with an exponential function, and the T2 relaxation can be determined directly from the exponential waveform. As this procedure needs to be repeated for each pixel of the scanned tissue, the processing of large images then becomes very intensive. For example, given the common resolution of 256×256 with 20 slices and five echoes at different times TE, it is necessary to reconstruct 1.3∙106 exponential functions. At present, such computation performed on a regular PC may last even several minutes. This paper introduces the results provided by accelerated computation based on parallelization and carried out with a graphics card. By using the simple method of linear regression, we obtain a processing time of less than 36 ms. Another effective option consists in the Levenberg-Marquardt algorithm, which enables us to reconstruct the same image in 96 ms. This period is at least 900 times shorter than that achievable with professional software. In this context, the paper also comprises an analysis of the results provided by the above-discussed techniques.

If the inline PDF is not rendering correctly, you can download the PDF file here.

  • [1] Dvořák P. Bartušek K. Smékal Z. (2014). Unsupervised pathological area extraction using 3D T2 and FLAIR MR images. Measurement Science Review 14 (2) 357-364.

  • [2] Mikulka J. Burget R. Říha K. Gescheidtová E. (2013). Segmentation of brain tumor parts in magnetic resonance images. In 36th International Conference on Telecommunications and Signal Processing 2-4 July 2013. IEEE 565-568.

  • [3] Hnilicová P. Bittšanský M. Dobrota D. (2014). Optimization of brain T2 mapping using standard CPMG sequence in a clinical scanner. Measurement Science Review 14 (2) 117-125.

  • [4] Hahn E. (1950) Spin echoes. Physical Review 80 (4) 580-594.

  • [5] Jiménez J. Ruiz de Miras J. (2012). Fast boxcounting algorithm on GPU. Computer Methods and Programs in Biomedicine 108 (3) 1229-1242.

  • [6] Gogola D. Krafčík A. Štrbák O. Frollo I. (2013). Magnetic resonance imaging of surgical implants made from weak magnetic materials. Measurement Science Review 13 (4) 165-168.

  • [7] Eklund A. Andersson M. Knutsson H. (2012). fMRI analysis on the GPU. Computer Methods and Programs in Biomedicine 105 (2) 145-161.

  • [8] Zhu F. Rodriguez-Gonzalez D. Carpenter T. Atkinson M. Wardlaw J. (2012). Parallel perfusion imaging processing using GPGPU. Computer Methods and Programs in Biomedicine 108 (3) 1012-1021.

  • [9] Broxvall M. Emilsson K. Thunberg P. (2012). Fast GPU based adaptive filtering of 4D echocardiography. IEEE Transactions on Medical Imaging 31 (6) 1165-1172.

  • [10] Prassni J.S. Ropinski T. Hinrichs K. (2010). Uncertainty-aware guided volume segmentation. IEEE Transactions on Visualization and Computer Graphics 16 (6) 1358-1365.

  • [11] Shams R. Sadeghi P. Kennedy R.A. Hartley R.I. (2010). A survey of medical image registration on multicore and the GPU. IEEE Signal Processing Magazine 27 (2) 50-60.

  • [12] Fluck O. Vetter C. Wein W. Kamen A. Preim B. Westermann R. (2011). A survey of medical image registration on graphics hardware. Computer Methods and Programs in Biomedicine 104 (3) e45-e57.

  • [13] Pratx G. Xing L. (2011) GPU computing in medical physics: A review. Medical Physics 38 (5) 2011 2685-2697.

  • [14] Eklund A. et al. (2013) Medical image processing on the GPU - Past present and future. Medical Image Analysis 17 (8) 1073-1094.

  • [15] Giri S. Chung Y.C. Merchant A. Mihai G. Rajagopalan S. Raman S.V. Simonetti O.P. (2009). T2 quantification for improved detection of myocardial edema. Journal of Cardiovascular Magnetic Resonance 11 (1) 56.

  • [16] Jones C.K. (2003). T2 decay curve acquisition and analysis in MRI noise considerations short T2 and B1 field encoding. PhD thesis University of British Columbia.

  • [17] Yoo Y. (2003). Robust and fast T2 decay analysis for measuring myelin water in MRI. PhD thesis University of British Columbia.

  • [18] Fraile R. García-Ortega E. (2005). Fitting an exponential distribution. Journal of Applied Meteorology and Climatology 44 (10) 1620-1625.

  • [19] Rektorys K. (1995). Přehled užite matematiky II. Praha: Prometheus.

  • [20] Prodanov D. (2009). MRI Processor v. 1.1.6 22. ri_processor:start.

  • [21] Starcuk J. Starcuk Z. Kozlowski P. (2008). Marevisi v. 8.2. More information at: http://www.nrccnrc. mri_data_software.html.

  • [22] Press W.H. Teukolsky S.A. Vetterling W.T. Flannery B.P. (2007). Numerical Recipes: The Art of Scientific Computing. Cambridge University Press.

  • [23] Lee V.W. et al. (2010). Debunking the 100X GPU vs. CPU myth: An evaluation of throughput computing on CPU and GPU. In ISCA’10 19-23 June 2010 Saint- Malo France

Journal information
Impact Factor

IMPACT FACTOR 2018: 1.122
5-year IMPACT FACTOR: 1.157

CiteScore 2018: 1.39

SCImago Journal Rank (SJR) 2018: 0.325
Source Normalized Impact per Paper (SNIP) 2018: 0.881

Cited By
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 169 69 1
PDF Downloads 82 52 2