Automated High-Order Shimming for Neuroimaging Studies

B0 inhomogeneity presents a significant challenge in MRI and MR spectroscopy, particularly at high-field strengths, leading to image distortion, signal loss, and spectral broadening. Existing high-order shimming methods can alleviate these issues but often require time-consuming and subjective manual selection of regions of interest (ROIs). To address this, we proposed an automated high-order shimming (autoHOS) method, incorporating deep-learning-based brain extraction and image-based high-order shimming. This approach performs automated real-time brain extraction to define the ROI of the field map to be used in the shimming algorithm. The shimming performance of autoHOS was assessed through in vivo echo-planar imaging (EPI) and spectroscopic studies at both 3T and 7T field strengths. AutoHOS outperforms linear shimming and manual high-order shimming, enhancing both the image and spectral quality by reducing the EPI image distortion and narrowing the MRS spectral lineshapes. Therefore, autoHOS demonstrated a significant improvement in correcting B0 inhomogeneity while eliminating the need for additional user interaction.


Introduction
Magnetic field (B 0 ) inhomogeneity in magnetic resonance imaging (MRI) can cause unwanted signal loss, image distortion, and spectral line broadening.Robust and automated shimming procedures are fundamental to the success of MRI techniques sensitive to B 0 inhomogeneity, such as magnetic resonance spectroscopy (MRS) [1,2] and echo-planar imaging (EPI) [3].
B 0 shimming is a procedure that the MRI scanner performs to minimize the B 0 variation, thereby enhancing image quality.Specifically, based on the measured B 0 field maps, i.e., the B 0 field distribution within the region of interest (ROI), MRI hardware shim coils generate additional B 0 fields to compensate for these inhomogeneities.There are mainly two classes of automatic shimming techniques: projection-based shimming and imagebased shimming [1,4].Projection-based shimming methods, such as the Fast Automatic Shimming Technique by Mapping Along Projections (FASTMAP) [5][6][7], utilize fields measured along six 1D column projections to calculate the first-and second-order spherical harmonic (SH) shim currents.While projection-based shimming has been proven to be fast and efficient for small-volume applications such as single-voxel MRS, its selected column projections may not capture localized field inhomogeneities [7].Image-based shimming methods acquire 3D field maps over a ROI and hence are more suitable for arbitrarily shaped volumes such as human brains [8,9].
High-order shimming methods employ higher-order spherical harmonics to model the B 0 inhomogeneity, allowing for a more accurate B 0 inhomogeneity correction and improved image quality.The introduction of high-order shims confers significant benefits to both brain and body imaging, especially at high-field strengths of 3T and above, by mitigating increased field inhomogeneity.For example, the High Order Shimming (HOS) software on the GE Signa 7.0T scanner calculates all shimming currents up to third order to minimize the B 0 inhomogeneity within a selected ROI of the acquired gradient-echo-based field maps [8].However, to date, most high-order shimming methods need a manually defined shim ROI if different from the imaging volume, which introduces intra-and inter-operator variability.Various automated brain extraction algorithms have emerged to address such variabilities, removing non-brain tissues from MR images, including the B 0 field map [10,11].In addition, skull-stripping by brain extraction can reduce the field magnitude variations caused by the presence of extracranial lipid signals [1,12].However, using the extracted brain as the shimming ROI is typically not recommended since brain extraction is time-consuming and may not be robust when applied to lower resolution data [1].Recently, the development of deep-learning-based image segmentation tools now enables fast and robust brain extraction [13,14].Among the existing (semi-)automated brain extraction algorithms, the deep-learning-based HD-BET stands out, as it can generate brain masks from diverse MRI scans, even in the presence of pathology, with a high degree of accuracy [10].Based on a 3D U-net architecture, HD-BET outperforms conventional and 2D deep-learning (DL) brain extraction tools, albeit at a higher computational cost [13,14].However, modern GPUs enable the runtime of HD-BET to be reduced to less than one minute [10].
In this work, we propose an automated High Order Shimming (autoHOS) prototype to conduct objective and automated high-order shimming based on automated brain extraction.Our method employs a combination of automated brain extraction [10] and image-based high-order shimming [8] to automatically define the shim ROI and compute the shim currents.We assessed our method using in vivo echo-planar imaging and spectroscopic studies performed on 3T and 7T scanners.Our findings indicate that autoHOS markedly improves the B 0 homogeneity, leading to superior image and spectral quality.

Automated HOS Pipeline and Implementation
The autoHOS prototype is developed in Python 3.8 and Tcl/Tk and depends on HD-BET and GE (GE HealthCare, Waukesha, WI, USA) HOS software.Figure 1A shows the flowchart of automated HOS.First, a 3D field map covering the whole brain (e.g., 28.0 × 28.0 × 21.0 cm 3 box) with 128 × 128 × 64 matrix size is acquired using a 3D gradient-echo pulse sequence.The magnitude images of the 3D field map are transferred to Volume Recon Engine (VRE), the image reconstruction architecture of GE MRI scanners, for automated brain extraction using a GPU.Alternatively, the magnitude images can be transferred to a remote computer with a GPU if one is not available on the VRE.Brain masks are generated by HD-BET in GPU mode and applied to magnitude images of the field map.The obtained skull-stripped magnitude images are then transferred back to the MRI host computer for least-squares calculation of the updated shim currents.By default, the autoHOS procedure will automatically adjust the high order shims over the whole brain (Figure 1B).However, it also supports selection of smaller ROIs based on voxel placement for single-voxel MRS applications (Figure S1).

Participants and Experimental Protocol
Eleven participants were recruited into this study after receiving informed written consent using an Institutional Review Board (IRB) protocol.The subjects had an average age of 51 ± 22 years and an average weight of 68 ± 13 kg.The image quality of single-shot gradient echo EPI images was used to benchmark the performance of the linear shimming, HOS, and autoHOS at both 3T and 7T.In addition, the spectral linewidth of 31 P MRS/MRSI at 7T was used to compare the shimming performance of HOS and autoHOS.

Image Analysis
Three-dimensional T2-weighted anatomical images were downsampled to match the spatial resolution of the T2*-weighted single-shot gradient echo EPI images using AFNI and served as an anatomical reference assumed to be free of image distortion [15].All T2*-weighted EPI images and B 0 field maps were aligned to these T2-weighted reference images using rigid-body transformation through the ANTs registration tool [16].Brain extraction of all images was performed using HD-BET [10].The edges of the skull-stripped images were obtained using AFNI's 3dEdge command.To evaluate the efficiency of different shimming strategies, the standard deviation of B 0 values was estimated both in each slice and across the whole field map, serving as an indicator of field homogeneity, respectively.The structural similarity index (SSIM) and peak signal-to-noise ratio (PSNR) were calculated using scikit-image to compare the image quality of the EPI images to the reference scans acquired with different shimming strategies [17].
The 31 P MRS/MRSI data were processed and analyzed using in-house Python scripts.For nonlocalized 31 P MRS, the obtained FIDs were zero-filled to 8K and apodized with a 5 Hz line-broadening before Fourier transformation.Nmrglue was used to conduct peak picking of the phase-corrected absorptive spectrum and quantify the width of the picked peaks as full width at half maximum (FWHM) [18].For 2D 31 P MRSI, a 20 Hz Gaussian line-broadening and a Hamming filter were applied to both the spectral and spatial dimensions of the k-space data prior to Fourier transformation.Because the 31 P MRSI spectrum cannot be phased due to its prolonged dead time, AMARES was used to quantify the linewidth as damping factors [19].
The Wilcoxon paired-sum test was used to compare the performance of various shimming methods.

Comparison of B 0 Homogeneities by Shimming Techniques
The standard deviation (SD) of the field map acquired after shimming within the brain region was used to assess the B 0 homogeneity (Figure 2) for each of the methods.Overall, the HOS and autoHOS displayed a similar performance at both 3T and 7T.Among the seven subjects studied, both high order shimming methods consistently demonstrated reduced smaller B 0 standard deviation values across slices as compared to linear shimming, with most of these reductions reaching statistical significance (p < 0.05) as highlighted in red on the sagittal view in Figure 2K-N.At 3T, HOS's reduction in B 0 inhomogeneity was significant only in tissues bordering the sinuses and ear canals, while autoHOS showed significant reductions across most slices.In contrast, at 7T, while both autoHOS and HOS displayed a comparable number of significantly improved slices (Figure 2M,N), when comparing the global SD(B 0 ) of HOS and autoHOS with linear shimming, notable distinctions were observed (Figure 2H,J), with both HOS and autoHOS significantly outperforming linear shimming (p < 0.001).Furthermore, autoHOS was superior to HOS at both 3T (p < 0.001) and 7T (p < 0.01).

Comparison of EPI Distortion by Shimming Techniques
In addition to assessing the field map, we evaluated the performance of various shimming methods by comparing the quality of the single-shot gradient-echo EPI image acquired using the three shimming approaches.Specifically, we calculated the PSNR and SSIM between the slices of the acquired T2*-weighted gradient-echo EPI images and the corresponding slices of the reference anatomical T2-weighted CUBE images.
At 3T, EPI image quality acquired with different shimming methods was compared using PSNR (Figure 3A-F).Both HOS and autoHOS exhibited modest increases in PSNR values across slices over the brain (Figure 3E), and on a global scale, both methods showed significant superiority over linear shimming (p < 0.001, Figure 3F).

Comparison of EPI Distortion by Shimming Techniques
In addition to assessing the field map, we evaluated the performance of various shimming methods by comparing the quality of the single-shot gradient-echo EPI image acquired using the three shimming approaches.Specifically, we calculated the PSNR and SSIM between the slices of the acquired T2*-weighted gradient-echo EPI images and the corresponding slices of the reference anatomical T2-weighted CUBE images.
At 3T, EPI image quality acquired with different shimming methods was compared using PSNR (Figure 3A-F).Both HOS and autoHOS exhibited modest increases in PSNR values across slices over the brain (Figure 3E), and on a global scale, both methods showed significant superiority over linear shimming (p < 0.001, Figure 3F).

Comparison of EPI Distortion by Shimming Techniques
In addition to assessing the field map, we evaluated the performance of various shimming methods by comparing the quality of the single-shot gradient-echo EPI image acquired using the three shimming approaches.Specifically, we calculated the PSNR and SSIM between the slices of the acquired T2*-weighted gradient-echo EPI images and the corresponding slices of the reference anatomical T2-weighted CUBE images.
At 3T, EPI image quality acquired with different shimming methods was compared using PSNR (Figure 3A-F).Both HOS and autoHOS exhibited modest increases in PSNR values across slices over the brain (Figure 3E), and on a global scale, both methods showed significant superiority over linear shimming (p < 0.001, Figure 3F).The SSIM of the edge images was used to gauge the underlying distortions in the acquired EPI images as compared to the corresponding slice acquired using a standard T2 acquisition (Figure 3G-L).Similar to PSNR comparison, both HOS and autoHOS displayed slight improvements in SSIM across slices (Figure 3K).Notably, although both HOS and autoHOS techniques significantly improved the SSIM values (p < 0.01 for HOS, p < 0.001 for autoHOS), autoHOS outperformed HOS (p < 0.001, Figure 3L).
Following our analyses at 3T, similar assessments were conducted at 7T.The overall findings closely mirrored those observed at the lower field strength.Overall, both HOS and autoHOS demonstrated statistically significant improvements over linear shimming (p < 0.001).However, between these two high-order shimming methods, no noticeable differences were observed (Figure S2A), except for the SSIM index for the outer brain contour, where autoHOS modestly but significantly outperformed HOS (Figure S2B).In a slice-by-slice comparison, EPI images obtained using HOS and autoHOS exhibited higher PSNR and SSIM values than those acquired with linear shimming across the entire brain (Figure 4E,K).Importantly, slices including regions close to air cavities, such as the frontal lobes and ear canals, exhibited significant improvements with both HOS and autoHOS (p < 0.05, highlighted in red in Figure 4F,L; also, see representative examples in Figure S3).
the EPI images and the corresponding outer edge of the brain in the EPI images, respectively.Panels (E) and (K) illustrate the PSNR and SSIM of the outer brain contour, respectively.Panels (F) and (L) depict comparisons of global PSNR and SSIM of outer brain contour, respectively.(** p < 0.01, *** p < 0.001).
The SSIM of the edge images was used to gauge the underlying distortions in the acquired EPI images as compared to the corresponding slice acquired using a standard T2 acquisition (Figure 3G-L).Similar to PSNR comparison, both HOS and autoHOS displayed slight improvements in SSIM across slices (Figure 3K).Notably, although both HOS and autoHOS techniques significantly improved the SSIM values (p < 0.01 for HOS, p < 0.001 for autoHOS), autoHOS outperformed HOS (p < 0.001, Figure 3L).
Following our analyses at 3T, similar assessments were conducted at 7T.The overall findings closely mirrored those observed at the lower field strength.Overall, both HOS and autoHOS demonstrated statistically significant improvements over linear shimming (p < 0.001).However, between these two high-order shimming methods, no noticeable differences were observed (Figure S2A), except for the SSIM index for the outer brain contour, where autoHOS modestly but significantly outperformed HOS (Figure S2B).In a slice-by-slice comparison, EPI images obtained using HOS and autoHOS exhibited higher PSNR and SSIM values than those acquired with linear shimming across the entire brain (Figure 4E,K).Importantly, slices including regions close to air cavities, such as the frontal lobes and ear canals, exhibited significant improvements with both HOS and autoHOS (p < 0.05, highlighted in red in Figure 4F,L; also, see representative examples in Figure S3).

Superiority of autoHOS in MRS Spectral Linewidths
The MRS spectral linewidth serves as a sensitive indicator of B0 homogeneity within the selected MRS voxel.We chose the spectral peak width of bioenergetic peaks that are distinctly separated in 31 P MRS, including ATPs and phosphocreatine (PCr), to assess the shimming performance.The comparison between MRS lineshapes acquired with HOS and autoHOS underscored the distinctive advantage of autoHOS in MRS.This advantage is exemplified by the narrower linewidth of representative cases of both nonlocalized 31 P MRS (Figure 5A) and 31 P MRSI (Figure S4) at 7T.In 22 nonlocalized 31 P MRS sessions, autoHOS significantly reduced the peak width of PCr, αATP, and γATP (Figure 5B).

Superiority of autoHOS in MRS Spectral Linewidths
The MRS spectral linewidth serves as a sensitive indicator of B 0 homogeneity within the selected MRS voxel.We chose the spectral peak width of bioenergetic peaks that are distinctly separated in 31 P MRS, including ATPs and phosphocreatine (PCr), to assess the shimming performance.The comparison between MRS lineshapes acquired with HOS and autoHOS underscored the distinctive advantage of autoHOS in MRS.This advantage is exemplified by the narrower linewidth of representative cases of both nonlocalized 31 P MRS (Figure 5A) and 31 P MRSI (Figure S4) at 7T.In 22 nonlocalized 31 P MRS sessions, autoHOS significantly reduced the peak width of PCr, αATP, and γATP (Figure 5B).Interestingly, even though the linewidths for βATP were narrower in spectra acquired with autoHOS, the difference did not achieve statistical significance.Interestingly, even though the linewidths for βATP were narrower in spectra acquired with autoHOS, the difference did not achieve statistical significance.

Discussion
B0 inhomogeneity has always been a challenge in MRI, particularly at higher magnetic field strengths like 7T.It introduces artifacts, reduces the SNR, and hampers the accuracy of quantitative measurements.The significance of the issue marks the importance of developing advanced shimming methods such as autoHOS.While both HOS and au-toHOS significantly reduce B0 inhomogeneity, autoHOS outperformed HOS on a broader scale.First, while not quantified, autoHOS is inherently faster and more objective than HOS because it eliminates the need for manual ROI selection via GUI interaction.In addition, the automated tool provides similar or better results as compared to the manual approach for high order shimming and is not biased based on the technologist running the scan.Compared to projection-based shimming techniques, such as FASTMAP and its derivatives [5][6][7]20], autoHOS does not rely on small spherical or cubic ROIs.Thus, it is more suitable for applications involving larger and often irregularly-shaped ROIs, such as EPI images, MRSI, and non-localized MRS.However, the limited number of field map points within a smaller spectroscopic voxel may affect the effectiveness of image-based shimming.In contrast, the 1D nature of projection-based shimming allows extended encoding times, which enhances the sensitivity to frequency variations and maximizes the accuracy of the magnetic field maps within such a small voxel.Although projection-based shimming has proven be to more efficient than full 3D B0 field mapping for single-voxel MRS [1,7], the new rapid 3D B0 field mapping techniques on which autoHOS is based are comparably fast and can easily incorporate third-or higher-order shimming terms.For instance, on a GE 7.0T SIGNA scanner, acquiring the B0 field map by a 3D gradient-echo pulse sequence with a matrix of 128 × 128 × 64 takes only 43 s.Immediately following the acquisition of the B0 field map, the brain extraction takes less than 30 s, thanks to the SIGNA VRE equipped with the NVIDIA T4 Tensor Core GPU.As a result, the entire au-toHOS procedure can be completed within 2 min.

High-Order Shimming Techniques Improve B0 Homogeneity Similarly at 3T and 7T
For B0 field maps at both 3T and 7T, the standard deviation of B0 values (SD(B0)) indicates that high-order shimming techniques, both manual HOS and autoHOS, significantly reduce B0 inhomogeneity across the entire brain at both 3T and 7T.As depicted in Figure 2G-J, the magnitude of SD(B0) of 7T is considerably larger than that at 3T,

Discussion
B 0 inhomogeneity has always been a challenge in MRI, particularly at higher magnetic field strengths like 7T.It introduces artifacts, reduces the SNR, and hampers the accuracy of quantitative measurements.The significance of the issue marks the importance of developing advanced shimming methods such as autoHOS.While both HOS and autoHOS significantly reduce B 0 inhomogeneity, autoHOS outperformed HOS on a broader scale.First, while not quantified, autoHOS is inherently faster and more objective than HOS because it eliminates the need for manual ROI selection via GUI interaction.In addition, the automated tool provides similar or better results as compared to the manual approach for high order shimming and is not biased based on the technologist running the scan.Compared to projection-based shimming techniques, such as FASTMAP and its derivatives [5][6][7]20], autoHOS does not rely on small spherical or cubic ROIs.Thus, it is more suitable for applications involving larger and often irregularly-shaped ROIs, such as EPI images, MRSI, and non-localized MRS.However, the limited number of field map points within a smaller spectroscopic voxel may affect the effectiveness of image-based shimming.In contrast, the 1D nature of projection-based shimming allows extended encoding times, which enhances the sensitivity to frequency variations and maximizes the accuracy of the magnetic field maps within such a small voxel.Although projection-based shimming has proven be to more efficient than full 3D B 0 field mapping for single-voxel MRS [1,7], the new rapid 3D B 0 field mapping techniques on which autoHOS is based are comparably fast and can easily incorporate third-or higher-order shimming terms.For instance, on a GE 7.0T SIGNA scanner, acquiring the B 0 field map by a 3D gradient-echo pulse sequence with a matrix of 128 × 128 × 64 takes only 43 s.Immediately following the acquisition of the B 0 field map, the brain extraction takes less than 30 s, thanks to the SIGNA VRE equipped with the NVIDIA T4 Tensor Core GPU.As a result, the entire autoHOS procedure can be completed within 2 min.

High-Order Shimming Techniques Improve B 0 Homogeneity Similarly at 3T and 7T
For B 0 field maps at both 3T and 7T, the standard deviation of B 0 values (SD(B 0 )) indicates that high-order shimming techniques, both manual HOS and autoHOS, significantly reduce B 0 inhomogeneity across the entire brain at both 3T and 7T.As depicted in Figure 2G-J, the magnitude of SD(B 0 ) of 7T is considerably larger than that at 3T, suggesting the increased inhomogeneity of ultra-high field strength 7T.Globally, HOS outperforms linear shimming significantly at both 3T and 7T (Figure 2H).A slice-by-slice comparison reveals that HOS substantially outperforms linear shimming in several slices adjacent to air cavities (Figure 2K), suggesting that the inclusion of the extracranial lipid pixels particularly affect the low-resolution field map (64 × 64) at 3T and introduce errors.While for 7T, HOS demonstrates significance in most slices (Figure 2M), suggesting the high-resolution field map (128 × 128) may alleviate the issues introduced by extracranial lipid pixels.In contrast, autoHOS, which is based on skull-stripped B 0 field maps, shows significant improvement across most slices at both 3T and 7T (Figure 2L,N).

HOS and AutoHOS Reduces EPI Distortions
For EPI images, high-order shimming techniques have been shown to reduce image distortions caused by B 0 inhomogeneity.At both 3T and 7T, although autoHOS is evidently superior to linear shimming, it does not show considerable improvement compared to HOS when gauging image quality with PSNR.This could be attributed to the brain extraction process, which does not remove all pixels with possible phase offsets [1], although the erroneous extracranial lipid pixels are already stripped.Remarkably, autoHOS produces a modest yet significantly improved brain outer edge compared to HOS at both 3T and 7T (Figure 3L and Figure S2B).This improvement can enhance fMRI studies by better co-registering fMRI images with structural MRI and facilitating region-of-interest analyses.

AutoHOS Significantly Improves Nonlocalized 31 P MRS
For nonlocalized 31 P MRS at 7T, the benefits of autoHOS become particularly evident.AutoHOS substantially narrows the lineshapes, suggesting enhanced global B 0 homogeneity, which is further corroborated by the notably improved B 0 maps.Such improvement is crucial as B 0 variation is the primary hurdle in achieving accurate spectral quantification.Importantly, spectral linewidth data points of autoHOS show tighter clustering across ATP and PCr peaks, indicating superior reproducibility, a benefit that probably can also be attributed to the elimination of operator variability.An intriguing observation is that although autoHOS visibly reduces the spectral width of βATP peak (Figure 5B), this narrowing does not achieve statistical significance like γATP and PCr.This reason behind this is that when B 0 -induced line broadening (e.g., ~25 Hz as shown in Figure 5A) affects all 31 P peaks, the relative percentage change is inherently smaller for broadened peaks like βATP than it is for the sharper peaks.

Limitations
Our study still has several limitations.First, the relatively small sample size (n = 7) potentially limits the statistical power of our findings related to B 0 field maps and EPI images.Second, due to constraints on scanner access and subject availability, the 22 sessions of nonlocalized MRS were acquired from six subjects over a span of 4 months-11 sessions were before and 11 after the deployment of autoHOS on our 7T scanner.Ideally, we would have compared the three shimming techniques using the same subjects within the same sessions, following the approach we took for assessing shimming performance with B 0 field map and EPI images.
It is worth noting that although we presented an example where autoHOS narrows lineshapes in a representative MRSI voxel (Figure S4), this comparison was made between a 31 P MRSI acquired with autoHOS and one acquired with HOS, where the entire brain was manually selected.However, it is common to select the ROI as defined by the voxel of the single-voxel or multi-voxel MRS scan.In such cases, there may be no difference between autoHOS and HOS since brain extraction is unnecessary for manually selecting a region smaller than the entire brain.Similarly, autoHOS does not impact local neuroimaging studies focusing on specific regions that do not include the skull.For this reason, the current autoHOS tool retains the ability to define the shimming ROI of interest based on the MRS prescription volume.

Future Work
In the future, we plan to improve the field map for shim calculation by utilizing brain masks generated as part of the automated prescription algorithm on the scanner.This would eliminate the most time-consuming aspect of the current solution.Ultimately, our goal is to integrate this into the pre-scan process, similar to coil calibration and center frequency adjustment.Furthermore, since autoHOS is deployed with Docker and is inherently cross-platform, it can be easily ported to MRI systems manufactured by vendors other than GE.

Conclusions
To address the important challenges of B 0 inhomogeneity, especially at high field strengths like 3T and 7T, we proposed and rigorously assessed the autoHOS method against existing linear shimming and manual HOS techniques on GE 3T and 7T scanners.Our findings underscored that both HOS and autoHOS effectively mitigate B 0 inhomogeneities at 3T and 7T.Importantly, autoHOS not only has the advantage of eliminating manual ROI selection but also shows broader superiority, particularly in reducing EPI image distortions and narrowing MRS/MRSI spectral linewidths.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/tomography9060168/s1, Figure S1 Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Figure 1 .
Figure 1.Flowchart of automated HOS: (A) First, a 3D gradient-echo pulse sequence is acquired to generate a B0 field map.Subsequentially, the DL-based brain extraction software operates on the image reconstruction computer or a remote server to produce a skull-stripped field map.By default, the brain-extracted field map is used for least-squares optimization of shim currents without ROI selection.(B) The top image shows an ROI (red) drawn manually for the brain region to be used for shim current optimization, while the bottom image shows the brain region automatically extracted also used for the same purpose.

Figure 1 .
Figure 1.Flowchart of automated HOS: (A) First, a 3D gradient-echo pulse sequence is acquired to generate a B 0 field map.Subsequentially, the DL-based brain extraction software operates on the image reconstruction computer or a remote server to produce a skull-stripped field map.By default, the brain-extracted field map is used for least-squares optimization of shim currents without ROI selection.(B) The top image shows an ROI (red) drawn manually for the brain region to be used for shim current optimization, while the bottom image shows the brain region automatically extracted also used for the same purpose.

Figure 2 .
Figure 2. Comparison of MR shimming approaches on the resulting field map.Linear, high order shimming (HOS), and automated HOS (autoHOS).AutoHOS and HOS improve global B0 homogeneity.(A-C) and (D-F) show representative B0 field maps at 3T and 7T, respectively.(G) and (I) display SD(B0) values across slices at 3T and 7T.(H) and (J) present boxplots of the global SD(B0) comparisons, where ** indicates p < 0.01 and *** means p < 0.001.(K,L) and (M,N) highlighted in red where HOS and autoHOS show significant improved shimming (p < 0.05) with results shown on the corresponding sagittal MRI images for 3T and 7T, respectively.

Figure 3 .
Figure 3. HOS and autoHOS reduce EPI image distortions at 3T. (A-D) and (G-J) display representative slices of the reference T2-weighted CUBE images (A,G), EPI images acquired with linear shimming (B,H), HOS (C,I), and autoHOS (D,J).The top row (A-D) and bottom row (G-J) showcase

Figure 2 .
Figure 2. Comparison of MR shimming approaches on the resulting field map.Linear, high order shimming (HOS), and automated HOS (autoHOS).AutoHOS and HOS improve global B 0 homogeneity.(A-C) and (D-F) show representative B 0 field maps at 3T and 7T, respectively.(G) and (I) display SD(B 0 ) values across slices at 3T and 7T.(H) and (J) present boxplots of the global SD(B 0 ) comparisons, where ** indicates p < 0.01 and *** means p < 0.001.(K,L) and (M,N) highlighted in red where HOS and autoHOS show significant improved shimming (p < 0.05) with results shown on the corresponding sagittal MRI images for 3T and 7T, respectively.

Figure 2 .
Figure 2. Comparison of MR shimming approaches on the resulting field map.Linear, high order shimming (HOS), and automated HOS (autoHOS).AutoHOS and HOS improve global B0 homogeneity.(A-C) and (D-F) show representative B0 field maps at 3T and 7T, respectively.(G) and (I) display SD(B0) values across slices at 3T and 7T.(H) and (J) present boxplots of the global SD(B0) comparisons, where ** indicates p < 0.01 and *** means p < 0.001.(K,L) and (M,N) highlighted in red where HOS and autoHOS show significant improved shimming (p < 0.05) with results shown on the corresponding sagittal MRI images for 3T and 7T, respectively.

Figure 3 .
Figure 3. HOS and autoHOS reduce EPI image distortions at 3T. (A-D) and (G-J) display representative slices of the reference T2-weighted CUBE images (A,G), EPI images acquired with linear shimming (B,H), HOS (C,I), and autoHOS (D,J).The top row (A-D) and bottom row (G-J) showcase

Figure 3 .
Figure 3. HOS and autoHOS reduce EPI image distortions at 3T. (A-D) and (G-J) display representative slices of the reference T2-weighted CUBE images (A,G), EPI images acquired with linear shimming (B,H), HOS (C,I), and autoHOS (D,J).The top row (A-D) and bottom row (G-J) showcase the EPI images and the corresponding outer edge of the brain in the EPI images, respectively.Panels (E) and (K) illustrate the PSNR and SSIM of the outer brain contour, respectively.Panels (F) and (L) depict comparisons of global PSNR and SSIM of outer brain contour, respectively.(** p < 0.01, *** p < 0.001).

Figure 4 .
Figure 4. HOS and autoHOS reduce EPI image distortions at 7T. (A-D) and (G-J) display representative slices of the reference T2-weighted CUBE images (A,G), EPI images acquired with linear shimming (B,H), HOS (C,I), and autoHOS (D,J).The top row (A-D) and bottom row (G-J) showcase the EPI images and the corresponding outer edge of the brain, respectively.Panels (E) and (K) illustrate the PSNR and SSIM of the outer brain contour, respectively.Panels (F) and (L) highlight slices in which both HOS and autoHOS demonstrate significant improvement, marked in red on sagittal view anatomical MRI, for PSNR and SSIM of the outer brain contour, respectively.

Figure 4 .
Figure 4. HOS and autoHOS reduce EPI image distortions at 7T. (A-D) and (G-J) display representative slices of the reference T2-weighted CUBE images (A,G), EPI images acquired with linear shimming (B,H), HOS (C,I), and autoHOS (D,J).The top row (A-D) and bottom row (G-J) showcase the EPI images and the corresponding outer edge of the brain, respectively.Panels (E) and (K) illustrate the PSNR and SSIM of the outer brain contour, respectively.Panels (F) and (L) highlight slices in which both HOS and autoHOS demonstrate significant improvement, marked in red on sagittal view anatomical MRI, for PSNR and SSIM of the outer brain contour, respectively.
: Selection of shim ROI.; Figure S2: HOS and autoHOS reduces EPI image distortions at 7T; Figure S3: Representative EPI images acquired with linear shimming, HOS, and autoHOS at 3T and 7T; Figure S4: Selected voxel of 31 P MRSI lineshape comparison between HOS and autoHOS at 7T. Author Contributions: Conceptualization, B.Y. and V.A.M.; methodology, J.X., B.Y., D.K. and V.A.M.; software, J.X. and B.Y.; investigation, J.X., B.Y. and V.A.M.; data curation, J.X. and V.A.M.; writingoriginal draft preparation, J.X. and V.A.M.; writing-review and editing, J.X., B.Y., D.K. and V.A.M.; visualization, J.X.; supervision, V.A.M.; project administration, V.A.M.; funding acquisition, V.A.M.All authors have read and agreed to the published version of the manuscript.Funding: This research used MR imaging equipment supported in part by NIH grants S10OD025025 and S10RR028821.VAM is supported in part by NIH grant P50HD103556.Institutional Review Board Statement: The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of the University of Iowa (IRB ID #: 201410788, 200810706) on 31 January 2022 and 2 February 2022, respectively.