4222

Human brain shimming performance comparison between matrix-like and stream-function SVD-based multi-coil arrays
Bruno Pinho Meneses1,2 and Alexis Amadon1
1Neurospin/CEA-Saclay, Gif-sur-Yvette, France, 2Université Paris-Saclay, Saclay, France

Synopsis

Last year we presented a novel whole-brain B0-shimming Multi-Coil Array design based on Stream Function Singular Value Decomposition on a cylindrical surface. Here this design is compared to more conventional single-layer matrix arrays in simulations based on a 100-subject B0-map database. While whole-brain shimming performance is naturally improved by our method, the effects of such optimization on slice-by-slice and region-specific shimming are not evident and are therefore explored in this work.

Introduction

Multi-Coil Arrays1-3 (MCA) for $$$B_0$$$ shimming in the human brain have become a popular alternative to Spherical Harmonics (SH) shim coils in the past few years due to their relative easiness to build, better performance on whole-brain shimming and their ability to perform dynamic slice-by-slice4 shimming . These arrays are usually built so that circular coil windings are regularly distributed over a cylindrical surface, creating a matrix of channels, yet not strictly optimized for the human brain anatomy.

We proposed a new method5,6 for designing optimized MCAs based on Singular Value Decomposition (SVD) of optimal Stream Functions7 (SF), providing loops of arbitrary geometry and optimally positioned to mitigate whole-brain inhomogeneity, showing, in simulations, superior performance to that of usual matrix MCAs (M-MCA) in whole-brain shimming while presenting a reduced channel count.

Nevertheless, assessment of slice-by-slice and region-specific shimming performances, of particular interest in high-resolution fMRI or spectroscopy, is also relevant to establish SF-SVD-based designs capable of satisfying several applications. Therefore, three brain Optimized MCA (O-MCA) designs are simulated in this work and compared to M-MCAs with increasing channel count.

Methods

A 100-subject database of three-dimensional $$$\Delta{}B_0$$$ brain fieldmaps was built from 3T acquisitions in a Siemens Magnetom Prisma imager equipped with second-order SH shim coils at 1.7-mm isotropic resolution. The $$$\Delta{}B_0$$$ field intensities were rescaled to 7T, since a shim system for UHF is intended. FSL’s brain extraction tool was used to exclude non-brain voxels.

All fieldmaps are used as target fields, providing 100 SFs8 upon which SVD is applied9; then the optimized loops are extracted from the three first principal modes, each loop being associated with a shimming channel5. Three concentric cylindrical coil formers of equal length 300-mm and radii 140.5-mm, 149.5-mm and 158.5-mm accommodate the optimized loops extracted from first, second and third modes respectively, as shown in Fig. 1b. Single-turn models of each loop are exported to ANSYS®-Maxwell for the calculation of their magnetic field distribution (in a region enclosing all brains) and of their complex impedance. The loops are assumed to be 20-turn windings of copper wire with 1-mm diameter to improve shimming efficiency with limited power supply. Loop current is constrained to 5 A. Each winding is assumed to yield the same normalized magnetic field profile as a single turn.

For comparison purposes, four different M-MCAs of 16, 24, 32 and 48 circular loop channels over cylindrical formers of 140.5-mm radius are designed with similar windings and current constraints (cf. Fig.2a).

For every subject, the shimming performance of each MCA is simulated by computing the electric currents minimizing the $$$L_2$$$ norm of the residual magnetic field in the region of interest, which can be the whole-brain, a thin transverse slice or a slab covering some cerebral lobe depending on the application. Here the frontal and temporal lobes, known to be strongly affected by susceptibility gradients, are addressed with slabs located by hand on 20 randomly-chosen brains of our database. Finally, resulting inhomogeneity (=standard deviation of the residual across all voxels of interest), inhomogeneity reduction and power dissipation are assessed.

Results and Discussion

From Figures 2a and 4, O-MCAs outperform M-MCAs in whole-brain shimming while having much smaller channel count. The 2-layer, 28-channel O-MCA shows slightly better inhomogeneity reduction, 30.9%(6.7), on the database than the 48-channel M-MCA, 28.3%(6.7), and much superior performance than the 23.5%(6.9) achieved by 32-channel M-MCA, a relative improvement of 31.5%. From Figures 3a and 5, still in whole-brain shimming, O-MCAs show improvement in the frontal and temporal lobes compared to M-MCAs, not surprisingly in view of the channel concentration in front of these regions.

In slice-by-slice shimming, however, the regular distribution of channels in M-MCA is more advantageous, as the 32-channel MCA outperforms the 3-layer, 39-channel O-MCA. Nevertheless, inhomogeneity reduction achieved by 2-layer and 3-layer O-MCAs, at 53.4%(6.3) and 56.3%(6.2), respectively, are satisfactory compared to the 57.1%(5.9) and 58.0%(5.9) for 32-ch and 48-ch M-MCA.

For region specific shimming, Figures 3a and 5, the 28-channel O-MCA inhomogeneity reduction for frontal and temporal lobes (55.4%(10.9) and 32.8%(7.3)) are very close to that achieved by 32-ch (55.6%(11.6) and 32.8%(7.8)) and 48-ch M-MCA (56.8%(11.2) and 34.3%(7.3)), while the 39-channel O-MCA outperforms them both by a small margin.

Within a design methodology, from Figure 2b, it is observed that power dissipation of the systems tends to increase with the homogeneity improvement, whether in slice-by-slice or whole-brain shimming. The 28-channel O-MCA produces 23.1W average dissipation for global shimming, much superior to the 2.8W dissipation of the 48-ch M-MCA. However, in slice-by-slice shimming, dissipation of the latter surpasses the 42.7W of the former by 18 W.

Conclusion and Perspectives

The SF-SVD-based optimization of MCAs is a powerful method for designing shimming systems with lower channel count, higher whole-brain shimming performance, while maintaining satisfactory results when used for slice-by-slice or focused shimming, in particular in the frontal lobe. The 2-layer, 28-channel system presents a good compromise between performance and channel count.

The present simulations resulted from an SF-SVD design targeting the whole brain. Further improvement of the optimized MCA performance could be attained if application-specific systems are envisioned, designed for addressing specific brain regions.

Acknowledgements

No acknowledgement found.

References

1. C. Juchem, T. W. Nixon, S. McIntyre, V. O. Boer, D. L. Rothman, and R. A. de Graaf. Dynamic multi-coil shimming of the human brain at 7T. Journal of Magnetic Resonance 2011; 212(2):280–288.

2. J. P. Stockmann et al. A 32-channel combined RF and B0 shim array for 3T brain imaging. Magn. Reson. Med. 2016; 75(1):441–451.

3. A. Aghaeifar et al. Dynamic B0 shimming of the human brain at 9.4 T with a 16-channel multi-coil shim setup. Magn. Reson. Med. 2018; 80(4):1714–1725.

4. C. Juchem, S. Umesh Rudrapatna, T. W. Nixon, and R. A. de Graaf. Dynamic multi-coil technique (DYNAMITE) shimming for echo-planar imaging of the human brain at 7 Tesla. NeuroImage 2015; 105:462–472.

5. B. Pinho Meneses, M. Luong, and A. Amadon. Optimized multi-coil array design for human brain shimming at Ultra-High Field. Proc. ISMRM 2019; p.1477.

6. B. Pinho Meneses and A. Amadon. Static-magnetic-field shimming coil system for magnetic resonance imaging. European Patent 072081EP/EPR 2019.

7. G. N. Peeren. Stream function approach for determining optimal surface currents. Journal of Computational Physics 2003; 191(1):305–321.

8. B. Pinho Meneses and A. Amadon. The Dipole Boundary Method: a simple approach to compute stream functions for shim coil design. Proc. ISMRM 2019; p.1476

9. B. Pinho Meneses and A. Amadon. A novel few-channel coil design for human brain shimming based on stream function Singular Value Decomposition. Proc. ISMRM 2019; p.1483.

Figures

Figure 1: a) Illustration of independent loop channels selection from first SVD coil mode for composing one layer of an SVD-based Optimized MCA. b) Optimized MCA layers formed from loop selections out of first, second and third SVD modes of tailored Stream Functions. The layers are composed of 16, 12 and 11 channels, respectively. c) Matrix MCA designs composed of 16, 24, 32 and 48 channels, with circular loop radii of 5.5 cm, 5.5 cm, 4.5 cm and 3.5 cm, respectively.

Figure 2: a) Comparison of B0 inhomogeneity (Hz) and inhomogeneity reduction (%) from our 100-subject database at 7T shimmed by different MCA technologies performing whole-brain or slice-by-slice shimming. The average 65.5-Hz baseline value results from 2nd-order shimming by the scanner. b) Average power dissipation for the above conditions.

Figure 3: a) Comparison of B0 inhomogeneity (Hz) and inhomogeneity reduction (%) in the frontal lobe and temporal lobes of a 20-subject database at 7T, shimmed by different MCA technologies performing whole-brain or region- (slab-) specific shimming. b) Average power dissipation for slab-specific shimming.

Figure 4: Comparison of ΔB0 fieldmaps in selected slices shimmed by a 48-Channel Matrix MCA or a 2-Layer 28-Channel driven for whole-brain and slice-by-slice shimming.

Figure 5: a) Example of brain mask outline (actually the addition of all sagittal slices) and highlighted slab placement for frontal and temporal lobe focused shimming. b) fieldmaps of the frontal and temporal lobes after localized shimming and whole-brain shimming performed by 48-channel Matrix MCA and 2-layer 28-channel optimized MCA. The depicted slices are in the central sagittal plane (left) and inferior axial plane (right). σl = std in entire slab; σs = std in depicted slice.

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