0782

A Joint Multi-Scale Variational Neural Network for Accelerating Free-breathing Whole-Heart qBOOST-T2 mapping
Niccolo Fuin1, Giorgia Milotta1, Thomas Kuestner1, Aurelien Bustin1, Gastao Cruz1, Rene Botnar1,2, and Claudia Prieto1,2
1Biomedical Engineering Department, School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom, 2Pontificia Universidad Católica de Chile, Santiago, Chile

Synopsis

T2 mapping is a promising technique for the characterization of myocardial inflammation and oedema. We recently proposed a quantitative 3D whole-heart sequence (qBOOST-T2) which provides co-registered 3D high-resolution bright-blood and T2 map volumes from a single free-breathing scan. However, high-resolution qBOOST-T2 requires long scan times of ~10 min. Here we propose a joint Multi-Scale Variational Neural Network (jMS-VNN) to enable the acquisition of 3D high-resolution bright-blood and accurate T2 map volumes in ~3 mins, and their reconstruction in ~30s. The proposed jMS-VNN jointly reconstructs data from multiple contrasts and efficiently apply dictionary-based signal matching for fast T2 map generation.

INTRODUCTION

T2 mapping is a promising technique for the characterization of myocardial inflammation and oedema1. Conventional 2D T2 mapping is performed under several breath-holds with limited spatial resolution and coverage. We recently introduced a novel accelerated quantitative 3D whole-heart sequence2 (qBOOST-T2) which provides co-registered high-resolution 3D bright- blood and T2 map volumes from a single free-breathing scan. qBOOST-T2 relies on image-navigators3 (iNAVs) to achieve 100% respiratory efficiency and predictable scan time. However, high-resolution qBOOST-T2 still requires long acquisition times of ~10 min for 3-fold acceleration (~30 mins for fully sampled acquisition). Here we propose a joint Multi-Scale Variational Neural Network (jMS-VNN) to enable the acquisition of qBOOST-T2 in ~3 mins. This is achieved by extending a novel joint Multi-Scale Variational Neural Network4-5 (jMS-VNN) to jointly reconstruct data from multiple contrasts and to efficiently apply dictionary-based signal matching6 for fast T2 map generation, enabling qBOOST-T2 reconstruction in ~30 seconds. All parameters of this formulation, including the reconstruction parameters, are learned during the offline training procedure. The proposed approach is evaluated for undersampling factor of 9-fold in 4 healthy subjects in comparison to iterative SENSE7 (it-SENSE) and Wavelet-based Compressed Sensing (CS)8 reconstruction.

METHODS

Acquisition: The framework of the 3D qBOOST-T2 sequence is shown in Fig. 1A. Three interleaved bSSFP bright-blood volumes are acquired with a Variable Density Cartesian trajectory with spiral profile order9. A T2-prepared inversion recovery (T2prep-IR) module is applied prior to the first dataset acquisition. T2-preparation is performed prior to the second volume (T2prep), whereas the third volume is acquired with no preparation. iNAVs are aquired at each heartbeat to estimate and correct for foot-head and right-left translational motion. Fourteen healthy subjects were scanned on a 1.5T scanner (Siemens Magnetom Aera). 3D qBOOST-T2 acquisition parameters included: resolution 1x1x2mm3 (reconstructed to 1mm3 isotropic), FOV 320x320x88-120mm3, subject specific diastolic trigger-delay and acquisition window (90-104ms), flip-angle 90o, T2prep 1st-heartbeat 50ms, T2prep 2nd-heartbeat 30ms, TI 110ms. Acquisitions were performed with an undersampling factor of 3x (all subjects) and 9x (4 subjects). Conventional T2prep-based 2D bSSFP T2 mapping sequences were acquired in 3 short axis views in all subjects for comparison purposes: resolution 1.8x1.8x8mm in a ~10sec breath-hold for each slice.
jMS-VNN Architecture: Translational motion corrected undersampled k-space data for the 3 different contrasts and coil sensitivity maps are stacked in the channel dimension and provided as input to the jMS-VNN. In each of the 10 stages of the network, the initial reconstruction is updated according to the joint variational scheme in Fig.1B. A first path enforces data consistency with the motion corrected undersampled k-space data. The second path is a regularization operator that jointly applies two parallel sets of real-valued filter kernels and corresponding activation functions to the magnitude and phase of the complexed-valued images. For the magnitude image, a multi-scale approach is applied with 3 parallel sets of 11x11, 5x5 and 1x1 filters. A joint reconstruction is performed by sharing the parameters (at every stage) along the channel dimension. EPG simulations6, matching the acquisition parameters, are carried out to generate a subject-specific signal evolution dictionary that is given as an additional input to the network. T2 maps are then generated within the network by matching each measured signal evolution to the closest dictionary entry.
MS-VNN Training: All parameters of this formulation, including the filter kernels, activation functions and data term weights, are learned during an offline training procedure. The network was trained on 2D complex-valued images obtained from retrospectively undersampled 9-fold T2prep-IR, T2prep and no-preparation data from 10 subjects (7200 images). During the training procedure, the output of the jMS-VNN is compared to the corresponding 3-fold accelerated reference images (reconstructed with SENSE) and corresponding T2 maps. The root mean-squared error was used as loss function.
Experiments: 9-fold prospectively undersampled data of 4 healthy subjects (not included in the training) were used to evaluate the framework. The proposed jMS-VNN approach was compared to it-SENSE7 and CS8 reconstruction, and to the 3-fold accelerated reference images, in terms of image quality and septal T2 quantification.

RESULTS

The average scan time was 3min±30sec with 100% respiratory scan efficiency. The average time for image reconstruction and T2 map generation was ~16 mins for CS and ~30s for jMS-VNN reconstructions. 3D co-registered T2prep-IR, T2prep and no preparation volumes for one representative healthy subject are presented in Fig. 2. Corresponding 3D T2 maps obtained with it-SENSE, CS and jMS-VNN are presented in Fig. 3. The proposed jMS-VNN framework presents higher image quality than it-SENSE7, and CS8 for both bright- and T2 map images. Reformatted short-axis T2 maps obtained with the proposed approach and reference standard 2D T2 maps are shown in Fig. 4A and quantitatively compared in Fig. 4B . A representative T2 map volume obtained with the proposed approach is shown in Fig. 5.

DISCUSSION

The combination of a highly accelerated qBOOST-T2 with a novel jMS-VNN reconstruction framework allows the acquisition of 3D co-registered high-resolution bright-blood volumes and T2 maps for comprehensive assessment of cardiovascular disease in a clinically feasible scan time of ~3min and their reconstruction in ~30s. Preliminary results showed good agreement in T2 quantification between the proposed approach and reference standard 2D T2 mapping.

Acknowledgements

This work was supported by EPSRC (EP/P001009, EP/P007619, EP/P032311/1) and Wellcome EPSRC Centre for Medical Engineering (NS/A000049/1).

References

1. Montant P, et al. MR imaging assessment of myocardial edema with T2 mapping. Diagn Interv Imag. 2015;96:885–890.

2. Milotta G, et al. 3D Whole-heart free-breathing qBOOST-T2 mapping. Magn Reson Med. 2019, Oct 21. doi:10.1002/mrm.28039.

3. Henningsson M, et al. Whole-heart coronary MR angiography with 2D self-navigated image reconstruction. Magn Reson Med. 2012;67:437–445.

4. Hammernik K, et al. Learning a variational network for reconstruction of accelerated MRI data. Magn Reson Med. 2018;79:3055–3071.

5. Fuin N, et al. A Variational Neural Network for Accelerating Free-breathing Whole-Heart Coronary MR Angiography. ISMRM 2019, p 0478, Montreal, Canada.

6. Weigel M, Extended phase graphs: Dephasing, RF pulses, and echoes - pure and simple. J Magn Reson Imag. 2015;41:266–295.

7. Pruessmann KP, et al. Advances in sensitivity encoding with arbitrary k-space trajectories. Magn Reson Med. 2001;46:638 – 651.

8. Lustig M, et al. Sparse MRI: The application of compressed sensing for rapid MR imaging. Magn Reson Med. 2007;58(6):1182–1195.

9. Prieto C, et al. Highly efficient respiratory motion compensated free-breathing coronary MRA using golden-step Cartesian acquisition. J Magn Res Imag. 2015;41:738–746.

Figures

Figure 1 - Acquisition: Data is acquired using an undersampled qBOOST-T2 sequence. EPG simulations are carried out to generate a signal evolution dictionary that is given as an additional input to the network. jMS-VNN Architecture: It consists in a data-consistency proximal term and in two sets of filter kernels and activation functions that are applied to the magnitude and phase of the complex-valued images. For the magnitude image a multi-scale approach is applied. T2 maps are generated within the network by matching each measured signal evolution to the closest dictionary entry.

Figure 2 - High-resolution T2prep-IR (first row), T2prep (second row), no preparation (third row) images and corresponding T2 maps (fourth column) for one representative healthy subject are presented. Undersampled acquisitions with acceleration factors of 9x and 3x (reference) are shown. Images were reconstructed using it-SENSE7, CS8, and the proposed jMS-VNN reconstruction framework. jMS-VNN provides higher image quality than it-SENSE and CS for both bright blood and T2 map images, achieving similar image quality to the 3x accelerated reference (different acquisition).

Figure 3 - T2 maps for one representative healthy subject acquired with the proposed qBOOST-T2 sequence with 9-fold and 3-fold (reference) acceleration. 9-fold accelerated T2 map images were reconstructed using it-SENSE7, CS8 and the proposed jMS-VNN reconstruction framework. The average time for volume reconstruction and T2 map generation was ~10min for it-SENSE, ~12 mins for CS and ~30s for jMS-VNN reconstructions.

Figure 4 - A: Comparison between 2D short‐axis standard T2 maps and short‐axis reformatted 3D qBOOST‐T2 maps for two healthy volunteers. Apical, mid, and basal slices are shown. B: Bland Altman plot comparing the proposed qBOOST‐T2 sequence with the conventional 2D bSSFP T2 mapping technique for four healthy volunteers. Good agreement is observed: A slight T2 overestimation is obtained with qBOOST‐T2mapping (bias = 4.77 ms), however, T2 quantification is within the 95% interval.

Figure 5 (animated) - T2 map volume for one representative healthy subject acquired with the a 9-fold accelerated qBOOST-T2 sequence and reconstructed with the proposed jMS-VNN reconstruction framework.

Proc. Intl. Soc. Mag. Reson. Med. 28 (2020)
0782