Optimizing full 3D SPARKLING trajectories for high-resolution T2*-weighted Magnetic Resonance Imaging

Chaithya Giliyar Radhakrishna (1, 2) , Pierre Weiss (3, 4) , Guillaume Daval-Frérot (1, 5, 2) , Aurélien Massire (5) , Alexandre Vignaud (6) , Philippe Ciuciu (1, 2)

1 NEUROSPIN - Service NEUROSPIN
2 MIND - Modèles et inférence pour les données de Neuroimagerie
3 IMT - Institut de Mathématiques de Toulouse UMR5219
4 PRIMO (ITAV)
5 Siemens Healthineers [Saint-Denis]
6 BAOBAB - Building large instruments for neuroimaging: from population imaging to ultra-high magnetic fields

The Spreading Projection Algorithm for Rapid K-space samplING, or SPARKLING, is an optimization-driven method that has been recently introduced for accelerated 2D T2*-w MRI using compressed sensing. It has then been extended to address 3D imaging using either stacks of 2D sampling patterns or a local 3D strategy that optimizes a single sampling trajectory at a time. 2D SPARKLING actually performs variable density sampling (VDS) along a prescribed target density while maximizing sampling efficiency and meeting the gradient-based hardware constraints. However, 3D SPARKLING has remained limited in terms of acceleration factors along the third dimension if one wants to preserve a peaky point spread function (PSF) and thus good image quality. In this paper, in order to achieve higher acceleration factors in 3D imaging while preserving image quality, we propose a new efficient algorithm that performs optimization on full 3D SPARKLING. The proposed implementation based on fast multipole methods (FMM) allows us to design sampling patterns with up to 10^7 k-space samples, thus opening the door to 3D VDS. We compare multi-CPU and GPU implementations and demonstrate that the latter is optimal for 3D imaging in the high-resolution acquisition regime (600µm isotropic). Finally, we show that this novel optimization for full 3D SPARKLING outperforms stacking strategies or 3D twisted projection imaging through retrospective and prospective studies on NIST phantom and in vivo brain scans at 3 Tesla. Overall the proposed method allows for 2.5-3.75x shorter scan times compared to GRAPPA-4 parallel imaging acquisition at 3 Tesla without compromising image quality.

Mots clés: 3D MRI; optimization; compressed sensing; non-Cartesian; acceleration

Domaines: Imagerie médicale; Traitement du signal et de l'image [eess.SP]; Optimisation et contrôle [math.OC]; Calcul parallèle, distribué et partagé [cs.DC]

Chaithya Giliyar Radhakrishna, Pierre Weiss, Guillaume Daval-Frérot, Aurélien Massire, Alexandre Vignaud, et al.. Optimizing full 3D SPARKLING trajectories for high-resolution T2*-weighted Magnetic Resonance Imaging. IEEE Transactions on Medical Imaging, 2022. ⟨hal-03090471v2⟩