Accelerated Simultaneous Multi-Slice (SMS) fMRI using spiral acquisition and low rank plus sparse (L+S) image reconstruction
Andrii Y Petrov1, Michael Herbst1,2, and V Andrew Stenger1

1Department of Medicine, University of Hawaii, Honolulu, HI, United States, 2Department of Radiology and Medical Physics, University Medical Center Freiburg, Freiburg, Germany

### Synopsis

Sub-second whole brain imaging using Simultaneous Multi-Slice (SMS) imaging is of particular interest in Blood Oxygen Level Dependent (BOLD) functional Magnetic Resonance Imaging (fMRI). Faster acquisitions with higher temporal sampling of the BOLD time course provides several advantages including increased sensitivity in detecting functional activation, the possibility of filtering out physiological noise for improving temporal SNR, and freezing out head motion. The most commonly used strategy to accelerate image acquisition time using SMS involves using parallel imaging methods. We propose to accelerate SMS imaging by under sampling the number of excited slices in the kz-t domain and L+S matrix decomposition method for reconstruction of slice aliased functional images. We present human fMRI results at 3T using 3D spiral sampling with SMS excitation and L+S reconstruction of the aliased slice data.

### Purpose

Sub-second whole brain imaging using Simultaneous Multi-Slice (SMS) imaging is of particular interest in Blood Oxygen Level Dependent (BOLD) functional Magnetic Resonance Imaging (fMRI) 1-3. Faster acquisitions with higher temporal sampling of the BOLD time course provides several advantages including increased sensitivity in detecting functional activation by acquiring a greater number of samples, the possibility of filtering out physiological noise for improving temporal SNR, and freezing out head motion. The most commonly used strategy to accelerate image acquisition time using SMS involves using parallel imaging methods 4-5, where non-acquired k-space data are estimated using receiver array coil sensitivity information. We propose to accelerate SMS imaging by under sampling the number of excited slices in the kz-t domain and L+S matrix decomposition method for reconstruction of slice aliased functional images. We present human fMRI results at 3T using 3D spiral sampling with SMS excitation and L+S reconstruction of the aliased slice data.

### Theory

The fMRI data can be structured as the space-time matrix M, where each column is a temporal frame, and decomposed into low-rank (L) and sparse (S) components by solving the following optimization problem:

$$min \parallel L \parallel_*+\lambda \parallel S \parallel_1, s.t. M=W(L+S)$$

where $\parallel L \parallel _*$ and $\parallel S \parallel_1$ are the nuclear-norm and $l_1$ norms of $L$ and $S$, $\lambda=1/\sqrt{n_1 n_2}$ is the nuclear norm vs $l_1$ norm weight. The decomposition results in matrix $L$ and matrix $S$, where $L$ represents a highly-stuctured information, such as tissue, or slowly-changing periodic information, such as cardic pulsation or respiratory movement. $S$ captures less-structured information, such as new localized dynamic information. $L+S$ results in matrix $M$ that represents background information with overlaid localized dynamic components.

### Methods

We applied the L+S method to fMRI on a healthy adult volunteer on a 3T scanner (Siemens, Tim Trio) with a 12-channel head coil. An SMS 3D spiral sequence was used to acquire fully sampled fMRI data. The imaging parameters were: FOV=22 cm, 64x64 matrix resolution, 16 simultaneous 3 mm thick slices, TR/TE=100/30ms, flip angle 30 degrees and 150 temporal frames. Sixteen z phase encodes were applied such that each slice could be separated during a 3D image reconstruction. The paradigm was a flashing checkerboard consisting of four 30 sec “on” and 30 sec “off” blocks for a total duration of 4 minutes. The data were then retrospectively undersampled by a x4 factor in the kz-t domain by keeping the two centre lines of kz-space and randomly excluding the remaining two planes at every time point (Fig. 1). L+S reconstructions were then performed using a temporal FFT as a sparsifying transform. Statistical analysis of the BOLD based activation was accomplished using a generalized linear model (GLM) after spatial smoothing.

### Results and Discussion

Fig. 2 shows L+S decomposition results for one slice in the vicinity of occipital lobe for the fully sampled data. Activation maps, scaled from t=6 to 10, were overlaid on the reconstructed images. From the time profile among the green line it can be seen that the L component absorbs the static brain background, while S captures localized dynamic components with activation. Fig. 3 shows L+S reconstructions from x4 retrospectively undersampled data of the same slice with overlaid activation maps. Good recovery of activation maps is seen in the S and M components compared with B x4 images computed from undersampled data. Activation patterns in the S and M components are qualitatively consistent to those from the fully sampled data shown in Fig. 2. Finally, Fig. 4 shows L+S reconstruction results from x4 retrospectively undersampled data for multiple slices in the vicinity of motor and visual cortex. It can be seen that L+S was able to produce aliasing artifact corrected images and successfully decompose images into the static component L and the dynamic component S which captures activation. Again, activations patterns in the S and M images are qualitatively consistent with those computed from the fully sampled data seen in Fig. 2.

### Conclusion

In this work we demonstrated that a spiral SMS acquisition combined with L+S reconstruction can accelerate fMRI acquisition by a factor of x4. We have also shown that L+S algorithm is able to produce slice aliasing corrected images and successfully decompose them into L and S components comparable with the those obtained from the fully sampled data. Future work will be to add coil sensitivity encoding for higher acceleration than SMS with parallel imaging alone. We will also add blipped-spiral acquisitors for improved kz-t sampling6.

### Acknowledgements

Work supported by the NIH grants R01DA019912, R01EB011517, and K02DA020569.

### References

1. Feinberg, D.A., Moeller, S., Smith, S.M., et al. Multiplexed echo planar imaging for sub-second whole brain fMRI and fast diffusion imaging. PLoS One . 2010; 5 (12), e15710.

2. Larkman DJ, Hajnal JV, Herlihy AH, Coutts GA, Young IR, Ehnholm G. Use of multicoil arrays for separation of signal from multiple slices simultaneously excited. J Magn Reson Imaging. 2001;13(2):313-7.

3. Setsompop K, Gagoski BA, Polimeni JR, Witzel T, Wedeen VJ, Wald LL. Blipped-controlled aliasing in parallel imaging for simultaneous multislice echo planer imaging with reduced g-factor penalty. Magn Reson Med. 2011. Epub 2011/08/23.

4. Griswold, M.A., Jakob, P.M., Heidemann, R.M., et al. Generalized autocalibrating partially parallel acquisitions (GRAPPA). 612 Magn. Reson. Med. 2002; 47 (6): 1202–1210.

5. Pruessmann KP, Weiger M, Scheidegger MB, Boesiger P. SENSE: Sensitivity encoding for fast MRI. Magnetic Resonance in Medicine. 1999;42:952-62.

6. Zahneisen B, Poser BA, Ernst T, Stenger AV. Simultaneous Multi-Slice fMRI using spiral trajectories. Neuroimage. 2014;92:8-18.

### Figures

Figure 1. Schematic representation of the kz-t 4x under sampling where only two central encoding planes are always fully sampled and the remaining two planes are randomly selected at every time point.

Figure 2. L+S reconstruction from the fully sampled fMRI data of one slice in the vicinity of the visual cortex with overlaid activation maps scaled from 6 to 10. It can be seen that the L captures static highly structured brain information while S detects localized dynamic information.

Figure 3. L+S reconstruction from the x4 retrospectively under sampled fMRI data of one slice in the vicinity of the visual cortex with overlaid activation maps scaled from 6 to 10. It can be seen that good recovery of activation maps is observed in the L and M components compared with B images.

Figure 4. L+S reconstruction from the x4 retrospectively under sampled fMRI data of multiple slices in the vicinity of the motor and visual cortex with overlaid activation maps scaled from 6 to 10. Again, good recovery of activation maps is seen in the S and M components compared with B images.

Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
3258