We propose an imaging sequence for Simultaneous MultiSlice Multi-Shot EPI (SMS-MS EPI) with wave-CAIPI controlled aliasing, to significantly reduce the imaging time and geometric distortion in gradient echo imaging. We extend the MUSSELS low-rank constrained parallel imaging technique to SMS acceleration and exploit the similarities among the EPI shots for improved reconstruction. In simulations, we demonstrate the capability of our sequence to incorporate wave-CAIPI encoding, which allows higher acceleration rates by fully harnessing the three-dimensional encoding capability of multi-channel receive arrays. Using MUSSELS with wave-SMS, whole-brain T2*-weighted images at 1 mm isotropic resolution can be obtained at the total acceleration of Rtotal=24 (RinplanexRSMS=8x3), corresponding to an acquisition with high image quality and geometric fidelity.
Introduction
Wave controlled aliasing in parallel imaging (wave-CAIPI) employs extra sinusoidal gradient modulations during the readout to effectivelyMethod
We designed the imaging sequence for SMS-MS EPI with wave-CAIPI as shown in Figure 1. Multiband RF is used to collect the signals from multiple slices at the same time. The navigator signals are obtained during the first three readouts, and extra gradient modulations for blipped-CAIPI (6) and wave-CAIPI (1) encodings are incorporated into the EPI readout. With the sign of the readout gradient alternating from one EPI-line to the next, the wave-gradients were also made to alternate, with waveforms flipping in time and gradient-polarity directions from one readout to the next.
The forward model for SMS-MS EPI with wave-CAIPI is shown in Figure 2. The image data are modulated by the point-spread function (PSF) from Gy and Gz wave-gradients to generate the wave image, which is spatially shifted along y and slice-collapsed through blipped-CAIPI SMS encoding to obtain the SMS wave data. The SMS wave data is then sub-sampled in k-space through in-plane acceleration and multi-shot acquisition. The aliased SMS-MS-wave image can then be obtained by 2D Fourier transform of the k-space SMS-MS-wave data.
We simulated the k-space SMS-MS-wave data using the forward model in Figure 2. The high-resolution reference images used for the simulation were obtained from several single-band MS EPI acquisition which was previously reconstructed using the MUSSELS. The voxel size is 1mmⅹ1mmⅹ1mm, the FOV is 220mmⅹ208mmⅹ120mm, and the 32-channel head coil was used. Simulations were performed to create a 4-shot SMS-MS acquisition, where the acceleration in each EPI-shot is RinplaneⅹRSMS = 8ⅹ3, with shifted complementary ky sampling between shots and application of blipped-CAIPI to create N/3-FOV shift between slices before the sub-sampling. Data were simulated for acquisition with and without wave-CAIPI encoding and both data modulated by blipped-CAIPI. Wave-encoding with 1 sinusoidal cycle and gradient amplitude of 20 mT/m was used for both y and z directions.
Image reconstructions were performed using generalized parallel imaging reconstruction for the forward model described in Figure 2, with low-rank constraint added as per MUSSELS to mitigate reconstruction error due from shot-to-shot phase variations. Potential reduction in effective resolution of the final reconstructed image due to wave-encoding was also characterized using an intra-voxel phase simulation.
Results and Discussion
The reconstructed images from the simulated data are shown in Figure 3. Results from SMS-MS-EPI data without and with wave-CAIPI encoding are shown, with RMSE of 69 % and 14 %, respectively, compared to the reference images. As can be observed, SMS-MS-EPI without wave-encoding suffers from severe artifacts from the high accelerations, that wave-encoding is shown to nicely mitigate. The PSF estimation of the reconstructed images using the intra-voxel phase simulation is shown in Figure 4. The signals within a voxel experience small de-phasing due to the phase modulations from the wave-gradient as shown in Figure 4a and 4b. The amount of dephasing varies throughout the readout, which causes a small signal modulation along kx. It creates small side-lobes on the PSF of the final reconstructed image, as shown in Figure 4c. For the acquisition parameters used in the simulated reconstruction, the side-lobe amplitude in the PSF remains small at 1.6 %.1. Bilgic B, Gagoski BA, Cauley SF, Fan AP, Polimeni JR, Grant PE, Wald LL, Setsompop K. Wave‐CAIPI for highly accelerated 3D imaging. Magnetic resonance in medicine 2015;73(6):2152-2162.
2. Gagoski BA, Bilgic B, Eichner C, Bhat H, Grant PE, Wald LL, Setsompop K. RARE/turbo spin echo imaging with simultaneous multislice Wave‐CAIPI. Magnetic resonance in medicine 2015;73(3):929-938.
3. Poser B, Bilgic B, Gagoski B, Uludag K, Stenger V, Wald L, Setsompop K. Echo-planar imaging with wave-CAIPI acquisition and reconstruction. In: Proceedings of the 25th Scientific Meeting of ISMRM, Honolulu, 2017 (abstract 1198).
4. Mani M, Jacob M, Kelley D, Magnotta V. Multi-shot multi-channel diffusion data recovery using structured low-rank matrix completion. arXiv preprint arXiv:160207274 2016.
5. Breuer FA, Blaimer M, Heidemann RM, Mueller MF, Griswold MA, Jakob PM. Controlled aliasing in parallel imaging results in higher acceleration (CAIPIR
6. Setsompop K, Gagoski BA, Polimeni JR, Witzel T, Wedeen VJ, Wald LL. Blipped‐controlled aliasing in parallel imaging for simultaneous multislice echo planar imaging with reduced g‐factor penalty. Magnetic resonance in medicine 2012;67(5):1210-1224.
7. Polimeni JR, Setsompop K, Gagoski B, McNab J, Triantafyllou C, Wald L. Rapid multi-shot segmented EPI using the Simultaneous Multi-Slice acquisition method. In: Proceedings of the 20th Annual Meeting of ISMRM, Melbourne, Australia, 2012 (abstract 2222).