RT info:eu-repo/semantics/article
T1 A clinically viable vendor-independent and device-agnostic solution for accelerated cardiac MRI reconstruction
A1 Martín González, Elena
A1 Moya Sáez, Elisa
A1 Menchón Lara, Rosa María
A1 Royuela del Val, Javier
A1 Palencia de Lara, César
A1 Rodríguez Cayetano, Manuel
A1 Simmross Wattenberg, Federico Jesús
A1 Alberola López, Carlos
K1 Magnetic resonance imaging
K1 Imagen por resonancia magnética
K1 Graphics processing unit
K1 Unidad de procesamiento gráfico
K1 OpenCLIPER
AB Background and objective: Recent research has reported methods that reconstruct cardiac MR images acquired with acceleration factors as high as 15 in Cartesian coordinates. However, the computational cost of these techniques is quite high, taking about 40 min of CPU time in a typical current machine. This delay between acquisition and final result can completely rule out the use of MRI in clinical environments in favor of other techniques, such as CT. In spite of this, reconstruction methods reported elsewhere can be parallelized to a high degree, a fact that makes them suitable for GPU-type computing devices. This paper contributes a vendor-independent, device-agnostic implementation of such a method to reconstruct 2D motion-compensated, compressed-sensing MRI sequences in clinically viable times. Methods: By leveraging our OpenCLIPER framework, the proposed system works in any computing device (CPU, GPU, DSP, FPGA, etc.), as long as an OpenCL implementation is available, and development is significantly simplified versus a pure OpenCL implementation. In OpenCLIPER, the problem is partitioned in independent black boxes which may be connected as needed, while device initialization and maintenance is handled automatically. Parallel implementations of both a groupwise FFD-based registration method, as well as a multicoil extension of the NESTA algorithm have been carried out as processes of OpenCLIPER. Our platform also includes significant development and debugging aids. HIP code and precompiled libraries can be integrated seamlessly as well since OpenCLIPER makes data objects shareable between OpenCL and HIP. This also opens an opportunity to include CUDA source code (via HIP) in prospective developments. Results: The proposed solution can reconstruct a whole 12–14 slice CINE volume acquired in 19–32 coils and 20 phases, with an acceleration factor of ranging 4–8, in a few seconds, with results comparable to another popular platform (BART). If motion compensation is included, reconstruction time is in the order of one minute. Conclusions: We have obtained clinically-viable times in GPUs from different vendors, with delays in some platforms that do not have correspondence with its price in the market. We also contribute a parallel groupwise registration subsystem for motion estimation/compensation and a parallel multicoil NESTA subsystem for -norm problem solving.
PB Elsevier
SN 0169-2607
YR 2021
FD 2021
LK https://uvadoc.uva.es/handle/10324/48168
UL https://uvadoc.uva.es/handle/10324/48168
LA eng
NO Computer Methods and Programs in Biomedicine, 2021, vol. 207, 106143
NO Producción Científica
DS UVaDOC
RD 24-nov-2024