GPU accelerated Cartesian GRAPPA reconstruction using CUDA

被引:6
作者
Inam, Omair [1 ]
Qureshi, Mahmood [1 ]
Laraib, Zoia [1 ]
Akram, Hamza [1 ]
Omer, Hammad [1 ]
机构
[1] COMSATS Univ Islamabad, Dept Elect & Comp Engn, Med Image Proc Res Grp MIPRG, Islamabad, Pakistan
关键词
GRAPPA; Parallel MRI; Graphic processing units; CUDA; Computational complexity; REAL-TIME MRI; IMPLEMENTATION; RESOLUTION;
D O I
10.1016/j.jmr.2022.107175
中图分类号
Q5 [生物化学];
学科分类号
071010 ; 081704 ;
摘要
Background and Objective: GRAPPA (Generalized Auto-calibrating Partially Parallel Acquisition) is an advanced parallel MRI reconstruction method (pMRI) that enables under-sampled data acquisition with multiple receiver coils to reduce the MRI scan time and reconstructs artifact free image from the acquired under-sampled data. However, the reduction in MRI scan time comes at the expense of long reconstruction time. It is because the GRAPPA reconstruction time shows exponential growth with increasing number of receiver coils. Consequently, the conventional CPU platforms may not adhere to the requirements of fast data processing for MR image reconstruction.& nbsp;Methods: Graphics Processing Units (GPUs) have recently emerged as a viable commodity hardware to reduce the reconstruction time of pMRI methods. This paper presents a novel GPU based implementation of GRAPPA using custom built CUDA kernels, to meet the rising demands of fast MRI processing. The proposed framework exploits intrinsic parallelism in the calibration and synthesis phases of GRAPPA reconstruction process, aiming to achieve high speed MR image reconstruction for various GRAPPA configuration settings using different number of receiver coils, auto-calibration signals (ACS), sizes of GRAPPA kernel and acceleration factors. In-vivo experiments (using 8, 12 and 30 receiver coils) are performed to compare the performance of the proposed GPU accelerated GRAPPA with the CPU based GRAPPA extensions and GPU counterpart.& nbsp;Results: The results indicate that the proposed method achieves up to asymptotic to 47.8x, asymptotic to 17x & nbsp;and asymptotic to 3.8x & nbsp;speed up gains over multicore CPU (single thread), multicore CPU (8 thread) and Gadgetron (GPU based GRAPPA) respectively, without compromising the reconstruction accuracy.& nbsp;Conclusions: The proposed method reduces the GRAPPA reconstruction time by employing the calibration phase (GRAPPA weights estimation) and synthesis phase (interpolation) on GPU. Our study shows that the proposed GPU based parallel framework for GRAPPA reconstruction provides a solution for highspeed image reconstruction while maintaining the quality of the reconstructed images. (C)& nbsp;2022 Elsevier Inc. All rights reserved.
引用
收藏
页数:17
相关论文
共 44 条
[1]  
[Anonymous], 2011, CUDA by example: an introduction to general-purpose GPU programming
[2]  
Barney B., 2009, POSIX THREADS PROGRA, V5, P46
[3]   Comparison of reconstruction accuracy and efficiency among autocalibrating data-driven parallel imaging methods [J].
Brau, Anja C. S. ;
Beatty, Philip J. ;
Skare, Stefan ;
Bammer, Roland .
MAGNETIC RESONANCE IN MEDICINE, 2008, 59 (02) :382-395
[4]   Real-time MRI Guidance of Cardiac Interventions [J].
Campbell-Washburn, Adrienne E. ;
Tavallaei, Mohammad A. ;
Pop, Mihaela ;
Grant, Elena K. ;
Chubb, Henry ;
Rhode, Kawal ;
Wright, Graham A. .
JOURNAL OF MAGNETIC RESONANCE IMAGING, 2017, 46 (04) :935-950
[5]   Magnetic resonance imaging in children: common problems and possible solutions for lung and airways imaging [J].
Ciet, Pierluigi ;
Tiddens, Harm A. W. M. ;
Wielopolski, Piotr A. ;
Wild, Jim M. ;
Lee, Edward Y. ;
Morana, Giovanni ;
Lequin, Maarten H. .
PEDIATRIC RADIOLOGY, 2015, 45 (13) :1901-1915
[6]   Nomenclature for real-time magnetic resonance imaging [J].
Dietz, Bryson ;
Fallone, B. Gino ;
Wachowicz, Keith .
MAGNETIC RESONANCE IN MEDICINE, 2019, 81 (03) :1483-1484
[7]   XD-GRASP: Golden-Angle Radial MRI with Reconstruction of Extra Motion-State Dimensions Using Compressed Sensing [J].
Feng, Li ;
Axel, Leon ;
Chandarana, Hersh ;
Block, Kai Tobias ;
Sodickson, Daniel K. ;
Otazo, Ricardo .
MAGNETIC RESONANCE IN MEDICINE, 2016, 75 (02) :775-788
[8]   Channel Reduction in Massive Array Parallel MRI [J].
Feng, Shuo ;
Ji, Jim .
2009 ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY, VOLS 1-20, 2009, :4045-4048
[9]   Efficient large-array k-domain parallel MRI using channel-by-channel array reduction [J].
Feng, Shuo ;
Zhu, Yudong ;
Ji, Jim .
MAGNETIC RESONANCE IMAGING, 2011, 29 (02) :209-215
[10]   An in-depth analysis of the impact of processor affinity on network performance [J].
Foong, A ;
Fung, J ;
Newell, D .
2004 12TH IEEE INTERNATIONAL CONFERENCE ON NETWORKS, VOLS 1 AND 2 , PROCEEDINGS: UNITY IN DIVERSITY, 2004, :244-250