B0 shimming of the human heart at 7T

Implementing a two-stage, shim field-calibrated superconducting shimming method on a 7 T cryogen-free small animal MRI magnet

Ultrahigh field systems (≥ 7 T) can increase the signal-to-noise ratio of magnetic resonance imaging (MRI), improving imaging performance compared to systems with lower fields. However, these enhancements heavily rely on a high B0 magnetic field homogeneity level, which can be achieved through superconducting shimming. This paper presents a novel two-stage superconducting shimming method designed to achieve precise shimming for a 7 T MRI superconducting magnet. In the initial stage, detailed measurements and fittings were conducted to determine the current polarity and the axial or circumferential positions of the shim fields. Subsequently, an optimization strategy was implemented to determine the optimal shim currents with a flexible target field. The second stage involves an iterative process to fine-tune the current of a specific shim coil, identified as having the most significant impact on field homogeneity. The overall fitness of 99.5% underscores the precision in determining the current polarity and position of the shim fields. Significantly, the calibrated shim system substantially improves the peak-to-peak and Root Mean Square Error (RMSE) field homogeneities from 107.42 ppm and 37.00 ppm to 11.12 ppm and 3.26 ppm, respectively, representing improvements of 89.65% and 91.19%. Furthermore, the simulation results of the fine-tuning stage demonstrate additional enhancements in peak-to-peak field homogeneity, to 9.9 ppm by reducing the current of the Z2 shim coil by 51.3 mA. Additionally, the shimmed magnetic field exhibited high time stability, with a maximum variation of only 27 µT observed within 48 h. Thus, the proposed two-stage superconducting shimming framework effectively addresses the challenge of imperfect B0 magnetic fields, enhancing peak-to-peak and RMSE field homogeneity. The stepwise optimized approach also mitigates deviations caused by shim-to-shim coupling, demonstrating its efficacy in achieving precise shimming in ultrahigh-field MRI systems.

Rotary excitation of non-sinusoidal pulsed magnetic fields: Towards non-invasive direct detection of cardiac conduction

PURPOSE

There is a need for high resolution non-invasive imaging methods of physiologic magnetic fields. The purpose of this work is to develop a MRI detection approach for non-sinusoidal magnetic fields based on the rotary excitation (REX) mechanism which was previously successfully applied for the detection of oscillating magnetic fields in the sub-nT range.

METHODS

The new detection concept was examined by means of Bloch simulations, evaluating the interaction effect of spin-locked magnetization and low-frequency pulsed magnetic fields. The REX detection approach was validated under controlled conditions in phantom experiments at 3 T. Gaussian and sinc-shaped stimuli were investigated. In addition, the detection of artificial fields resembling a cardiac QRS complex, which is the most prominent peak visible on a magnetocardiogram, was tested.

RESULTS

Bloch simulations demonstrated that the REX method has a high sensitivity to pulsed fields in the resonance case, which is met when the spin-lock frequency coincides with a non-zero Fourier component of the stimulus field. In the experiments, we found that magnetic stimuli of different durations and waveforms can be distinguished by their characteristic REX response spectrum. The detected REX amplitude was proportional to the stimulus peak amplitude (R2 > 0.98) and the lowest field detection was 1 nT. Furthermore, the detection of QRS-like fields with varying QRS durations yielded significant results in a phantom setup (p < 0.001).

CONCLUSION

REX detection can be transferred to non-sinusoidal pulsed magnetic fields and could provide a non-invasive, quantitative tool for spatially resolved assessment of cardiac biomagnetism. Potential applications include the direct detection and characterization of cardiac conduction.

Quantification of the rotating frame relaxation time T2ρ: Comparison of balanced spin-lock and continuous-wave Malcolm-Levitt preparations

For the quantification of rotating frame relaxation times, the T2ρ relaxation pathway plays an essential role. Nevertheless, T2ρ imaging has been studied only to a small extent compared with T1ρ, and preparation techniques for T2ρ have so far been adapted from T1ρ methods. In this work, two different preparation concepts are compared specifically for the use of T2ρ mapping. The first approach involves transferring the balanced spin-locking (B-SL) concept of T1ρ imaging. The second and newly proposed approach is a continuous-wave Malcolm-Levitt (CW-MLEV) pulse train with zero echo times and was motivated from T2 preparation strategies. The modules are tested in Bloch simulations for their intrinsic sensitivity to field inhomogeneities and validated in phantom experiments. In addition, myocardial T2ρ mapping was performed in mice as an exemplary application. Our results demonstrate that the CW-MLEV approach provides superior robustness and thus suggest that established methods of T1ρ imaging are not best suited for T2ρ experiments. In the presence of field inhomogeneities, the simulations indicated an increased banding compensation by a factor of 4.1 compared with B-SL. Quantification of left ventricular T2ρ time in mice yielded more consistent results, and values in the range of 59.2-61.1 ms (R2 = 0.986-0.992) were observed at 7 T.

Tailored and universal parallel transmit broadband pulses for homogeneous 3D excitation of the human heart at 7T

PURPOSE

To research and evaluate the performance of broadband tailored kT-point pulses (TP) and universal pulses (UP) for homogeneous excitation of the human heart at 7T.

METHODS

Relative 3DB 1 + $$ {\mathrm{B}}_1^{+} $$ -maps of the thorax were acquired from 29 healthy volunteers. TP and UP were designed using the small-tip-angle approximation for a different composition of up to seven resonance frequencies. TP were computed for each of the 29B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps, and UPs were calculated using 22B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps and tested in seven testcases. The performance of the pulses was analyzed using the coefficient of variation (CV) in the 3D heart volumes. The 3D gradient-echo (GRE) scans were acquired for the seven testcases to qualitatively validate theB 1 + $$ {\mathrm{B}}_1^{+} $$ -predictions.

RESULTS

Single- and double-frequency optimized pulses achieved homogeneity in flip angle (FA) for the frequencies they were optimized for, while the broadband pulses achieved uniformity in FA across a 1300 Hz frequency range.

CONCLUSION

Broadband TP and UP can be used for homogeneous excitation of the heart volume across a 1300 Hz frequency range, including the water and the main six fat peaks, or with longer pulse durations and higher FAs for a smaller transmit bandwidth. Moreover, despite large inter-volunteer variations, broadband UP can be used for calibration-free 3D heart FA homogenization in time-critical situations.

Precision imaging of cardiac function and scar size in acute and chronic porcine myocardial infarction using ultrahigh-field MRI

BACKGROUND

7 T cardiac magnetic resonance imaging (MRI) studies may enable higher precision in clinical metrics like cardiac function, ventricular mass, and more. Higher precision may allow early detection of functional impairment and early evaluation of treatment responses in clinical practice and pre-clinical studies.

METHODS

Seven female German Landrace pigs were scanned prior to and at three time points (3-4 days, 7-10 days, and ~60 days) post myocardial infarction using a whole body 7 T system and three radiofrequency (RF) coils developed and built in-house to accompany animal growth.

RESULTS

The combination of dedicated RF hardware and 7 T MRI enables a longitudinal study in a pig model of acute and chronic infarction, providing consistent blood tissue contrast and high signal-to-noise ratio (SNR) in measurements of cardiac function, as well as low coefficients of variation (CoV) for ejection fraction (CoVintra-observer: 2%, CoVinter-observer: 3.8%) and infarct size (CoVintra-observer: 8.4%, CoVinter-observer: 3.8%), despite drastic animal growth.

CONCLUSIONS

Best results are achieved via manual segmentation. We define state-of-the-art procedures for large animal studies at 7 T.

A modular motion compensation pipeline for prospective respiratory motion correction of multi-nuclear MR spectroscopy

Magnetic resonance (MR) acquisitions of the torso are frequently affected by respiratory motion with detrimental effects on signal quality. The motion of organs inside the body is typically decoupled from surface motion and is best captured using rapid MR imaging (MRI). We propose a pipeline for prospective motion correction of the target organ using MR image navigators providing absolute motion estimates in millimeters. Our method is designed to feature multi-nuclear interleaving for non-proton MR acquisitions and to tolerate local transmit coils with inhomogeneous field and sensitivity distributions. OpenCV object tracking was introduced for rapid estimation of in-plane displacements in 2D MR images. A full three-dimensional translation vector was derived by combining displacements from slices of multiple and arbitrary orientations. The pipeline was implemented on 3 T and 7 T MR scanners and tested in phantoms and volunteers. Fast motion handling was achieved with low-resolution 2D MR image navigators and direct implementation of OpenCV into the MR scanner's reconstruction pipeline. Motion-phantom measurements demonstrate high tracking precision and accuracy with minor processing latency. The feasibility of the pipeline for reliable in-vivo motion extraction was shown on heart and kidney data. Organ motion was manually assessed by independent operators to quantify tracking performance. Object tracking performed convincingly on 7774 navigator images from phantom scans and different organs in volunteers. In particular the kernelized correlation filter (KCF) achieved similar accuracy (74%) as scored from inter-operator comparison (82%) while processing at a rate of over 100 frames per second. We conclude that fast 2D MR navigator images and computer vision object tracking can be used for accurate and rapid prospective motion correction. This and the modular structure of the pipeline allows for the proposed method to be used in imaging of moving organs and in challenging applications like cardiac magnetic resonance spectroscopy (MRS) or magnetic resonance imaging (MRI) guided radiotherapy.

Reliable Off-Resonance Correction in High-Field Cardiac MRI Using Autonomous Cardiac B0 Segmentation with Dual-Modality Deep Neural Networks

B0 field inhomogeneity is a long-lasting issue for Cardiac MRI (CMR) in high-field (3T and above) scanners. The inhomogeneous B0 fields can lead to corrupted image quality, prolonged scan time, and false diagnosis. B0 shimming is the most straightforward way to improve the B0 homogeneity. However, today's standard cardiac shimming protocol requires manual selection of a shim volume, which often falsely includes regions with large B0 deviation (e.g., liver, fat, and chest wall). The flawed shim field compromises the reliability of high-field CMR protocols, which significantly reduces the scan efficiency and hinders its wider clinical adoption. This study aims to develop a dual-channel deep learning model that can reliably contour the cardiac region for B0 shim without human interaction and under variable imaging protocols. By utilizing both the magnitude and phase information, the model achieved a high segmentation accuracy in the B0 field maps compared to the conventional single-channel methods (Dice score: 2D-mag = 0.866, 3D-mag = 0.907, and 3D-mag-phase = 0.938, all p < 0.05). Furthermore, it shows better generalizability against the common variations in MRI imaging parameters and enables significantly improved B0 shim compared to the standard method (SD(B0Shim): Proposed = 15 ± 11% vs. Standard = 6 ± 12%, p < 0.05). The proposed autonomous model can boost the reliability of cardiac shimming at 3T and serve as the foundation for more reliable and efficient high-field CMR imaging in clinical routines.

High resolution simulation and measurement demonstrate oscillatory spatiotemporal B0 fluctuations across the human cardiac cycle

PURPOSE

Functional cardiac MRI scans employing balanced steady-state free precession sequences suffer from dark band artifacts in the myocardium due to B0 inhomogeneity. We recently introduced a novel method for the theoretical derivation of B0 distributions in the human heart. This study aims to simulate the B0 distributions in the heart across the cardiac cycle using structural MR images and validate the simulations via in vivo measured cardiac phase-specific B0 maps on the same subjects at 3T.

METHODS

Cardiac phase-specific B0 field maps were acquired from eight healthy subjects at 3T. B0 conditions were simulated based on tissue masks created from the cardiac-phase specific structural images from the in vivo B0 map scan and anatomical images from a thoracic MRI scan, adopting our recently published approach. The simulations and in vivo measurements were compared by calculating the spatial correlation of their B0 distributions and temporal correlation of the derived spherical harmonic coefficients throughout the cardiac cycle.

RESULTS

The spatial comparison of B0 maps between the simulation and in vivo measurement indicates an overall average correlation coefficient of 0.91 across the cardiac cycle in all subjects. Both groups show consistent high-level B0 patterns. Temporal variations of B0 conditions exhibit sinusoidal characteristics and are strongly correlated between simulation and in vivo.

CONCLUSION

Theoretical simulations employing regional anatomical features were validated by direct in vivo B0 mapping in the same subjects. The spatial B0 condition throughout the cardiac cycle exhibits oscillatory characteristics due to structural distortions of cardiac motion.

Ultra-high field cardiac MRI in large animals and humans for translational cardiovascular research

A key step in translational cardiovascular research is the use of large animal models to better understand normal and abnormal physiology, to test drugs or interventions, or to perform studies which would be considered unethical in human subjects. Ultrahigh field magnetic resonance imaging (UHF-MRI) at 7 T field strength is becoming increasingly available for imaging of the heart and, when compared to clinically established field strengths, promises better image quality and image information content, more precise functional analysis, potentially new image contrasts, and as all in-vivo imaging techniques, a reduction of the number of animals per study because of the possibility to scan every animal repeatedly. We present here a solution to the dual use problem of whole-body UHF-MRI systems, which are typically installed in clinical environments, to both UHF-MRI in large animals and humans. Moreover, we provide evidence that in such a research infrastructure UHF-MRI, and ideally combined with a standard small-bore UHF-MRI system, can contribute to a variety of spatial scales in translational cardiovascular research: from cardiac organoids, Zebra fish and rodent hearts to large animal models such as pigs and humans. We present pilot data from serial CINE, late gadolinium enhancement, and susceptibility weighted UHF-MRI in a myocardial infarction model over eight weeks. In 14 pigs which were delivered from a breeding facility in a national SARS-CoV-2 hotspot, we found no infection in the incoming pigs. Human scanning using CINE and phase contrast flow measurements provided good image quality of the left and right ventricle. Agreement of functional analysis between CINE and phase contrast MRI was excellent. MRI in arrested hearts or excised vascular tissue for MRI-based histologic imaging, structural imaging of myofiber and vascular smooth muscle cell architecture using high-resolution diffusion tensor imaging, and UHF-MRI for monitoring free radicals as a surrogate for MRI of reactive oxygen species in studies of oxidative stress are demonstrated. We conclude that UHF-MRI has the potential to become an important precision imaging modality in translational cardiovascular research.

Germany's journey toward 14 Tesla human magnetic resonance

Multiple sites within Germany operate human MRI systems with magnetic fields either at 7 Tesla or 9.4 Tesla. In 2013, these sites formed a network to facilitate and harmonize the research being conducted at the different sites and make this technology available to a larger community of researchers and clinicians not only within Germany, but also worldwide. The German Ultrahigh Field Imaging (GUFI) network has defined a strategic goal to establish a 14 Tesla whole-body human MRI system as a national research resource in Germany as the next progression in magnetic field strength. This paper summarizes the history of this initiative, the current status, the motivation for pursuing MR imaging and spectroscopy at such a high magnetic field strength, and the technical and funding challenges involved. It focuses on the scientific and science policy process from the perspective in Germany, and is not intended to be a comprehensive systematic review of the benefits and technical challenges of higher field strengths.

High-resolution simulation of B0 field conditions in the human heart from segmented computed tomography images

B₀ inhomogeneity leads to imaging artifacts in cardiac magnetic resonance imaging (MRI), in particular dark band artifacts with steady-state free precession pulse sequences. The limited spatial resolution of MR-derived in vivo B₀ maps and the lack of population data prevent systematic analysis of the problem at hand and the development of optimized B₀ shim strategies. We used readily available clinical computed tomography (CT) images to simulate the B₀ conditions in the human heart at high spatial resolution. Calculated B₀ fields showed consistency with MRI-based B₀ measurements. The B₀ maps for both the simulations and in vivo measurements showed local field inhomogeneities in the vicinity of lung tips with dominant Z3 spherical harmonic terms in the field distribution. The presented simulation approach allows for the derivation of B₀ field conditions at high spatial resolution from CT images and enables the development of subject- and population-specific B₀ shim strategies for the human heart.

A novel antisymmetric 16-element transceiver dipole antenna array for parallel transmit cardiac MRI in pigs at 7 T

To improve parallel transmit (pTx) and receive performance for cardiac MRI (cMRI) in pigs at 7 T, a dedicated transmit/receive (Tx/Rx), 16-element antisymmetric dipole antenna array, which combines L-shaped and straight dipoles, was designed, implemented, and evaluated in both cadavers and animals in vivo. Electromagnetic-field simulations were performed with the new 16-element dipole antenna array loaded with a pig thorax-shaped phantom and compared with an eight-element array of straight dipoles. The new dipole array was interfaced to a 7 T scanner in pTx mode (8Tx/16Rx). Imaging performance of the novel array was validated through MRI measurements in a pig phantom, an 85 kg pig cadaver, and two pigs in vivo (74 and 81 kg). Due to the improved decoupling between interleaved L-shaped and straight dipole elements, the 16-element dipole array fits within the same outer dimensions as an eight-element array of straight dipoles. This provides improvement of both transmit and receive characteristics and additional degrees of freedom for B1+ shimming. The antisymmetric dipole array demonstrated efficient suppression of destructive interferences in the B1+ field, with up to 25% improvement in the B1+ homogeneity achieved using static pTx-RFPA B1+ shimming in comparison with the hardware-adjusted state, which was optimized for single transmit. High-resolution (0.5 × 0.5 × 4 mm³ ) anatomical images of the heart after cardiac arrest proved good transmit and receive characteristics of the novel array design. Parallel imaging with an acceleration factor up to R = 6 was possible while maintaining a mean g factor of 1.55 within the pig heart. CINE images acquired in vivo in two pigs demonstrated SNR and parallel imaging capabilities similar to those of a reference 8Tx/16Rx dedicated loop array for cMRI in pigs.

In Vivo Magnetic Resonance Spectroscopy Methods for Investigating Cardiac Metabolism

Magnetic resonance spectroscopy (MRS) is a non-invasive and non-ionizing technique, enabling in vivo investigation of cardiac metabolism in normal and diseased hearts. In vivo measurement tools are critical for studying mechanisms that regulate cardiac energy metabolism in disease developments and to assist in early response assessments to novel therapies. For cardiac MRS, proton (1H), phosphorus (31P), and hyperpolarized 13-carbon (13C) provide valuable metabolic information for diagnosis and treatment assessment purposes. Currently, low sensitivity and some technical limitations limit the utility of MRS. An essential step in translating MRS for clinical use involves further technological improvements, particularly in coil design, improving the signal-to-noise ratios, field homogeneity, and optimizing radiofrequency sequences. This review addresses the recent advances in metabolic imaging by MRS from primarily the literature published since 2015.

Motion-compensated fat-water imaging for 3D cardiac MRI at ultra-high fields

PURPOSE

Respiratory motion-compensated (MC) 3D cardiac fat-water imaging at 7T.

METHODS

Free-breathing bipolar 3D triple-echo gradient-recalled-echo (GRE) data with radial phase-encoding (RPE) trajectory were acquired in 11 healthy volunteers (7M\4F, 21-35 years, mean: 30 years) with a wide range of body mass index (BMI; 19.9-34.0 kg/m² ) and volunteer tailored B 1 + shimming. The bipolar-corrected triple-echo GRE-RPE data were binned into different respiratory phases (self-navigation) and were used for the estimation of non-rigid motion vector fields (MF) and respiratory resolved (RR) maps of the main magnetic field deviations (ΔB₀ ). RR ΔB₀ maps and MC ΔB₀ maps were compared to a reference respiratory phase to assess respiration-induced changes. Subsequently, cardiac binned fat-water images were obtained using a model-based, respiratory motion-corrected image reconstruction.

RESULTS

The 3D cardiac fat-water imaging at 7T was successfully demonstrated. Local respiration-induced frequency shifts in MC ΔB₀ maps are small compared to the chemical shifts used in the multi-peak model. Compared to the reference exhale ΔB₀ map these changes are in the order of 10 Hz on average. Cardiac binned MC fat-water reconstruction reduced respiration induced blurring in the fat-water images, and flow artifacts are reduced in the end-diastolic fat-water separated images.

CONCLUSION

This work demonstrates the feasibility of 3D fat-water imaging at UHF for the entire human heart despite spatial and temporal B 1 + and B₀ variations, as well as respiratory and cardiac motion.

On the way to routine cardiac MRI at 7 Tesla - a pilot study on consecutive 84 examinations

INTRODUCTION

Cardiac magnetic resonance (CMR) at ultrahigh field (UHF) offers the potential of high resolution and fast image acquisition. Both technical and physiological challenges associated with CMR at 7T require specific hardware and pulse sequences. This study aimed to assess the current status and existing, publicly available technology regarding the potential of a clinical application of 7T CMR.

METHODS

Using a 7T MRI scanner and a commercially available radiofrequency coil, a total of 84 CMR examinations on 72 healthy volunteers (32 males, age 19-70 years, weight 50-103 kg) were obtained. Both electrocardiographic and acoustic triggering were employed. The data were analyzed regarding the diagnostic image quality and the influence of patient and hardware dependent factors. 50 complete short axis stacks and 35 four chamber CINE views were used for left ventricular (LV) and right ventricular (RV), mono-planar LV function, and RV fractional area change (FAC). Twenty-seven data sets included aortic flow measurements that were used to calculate stroke volumes. Subjective acceptance was obtained from all volunteers with a standardized questionnaire.

RESULTS

Functional analysis showed good functions of LV (mean EF 56%), RV (mean EF 59%) and RV FAC (mean FAC 52%). Flow measurements showed congruent results with both ECG and ACT triggering. No significant influence of experimental parameters on the image quality of the LV was detected. Small fractions of 5.4% of LV and 2.5% of RV segments showed a non-diagnostic image quality. The nominal flip angle significantly influenced the RV image quality.

CONCLUSION

The results demonstrate that already now a commercially available 7T MRI system, without major methods developments, allows for a solid morphological and functional analysis similar to the clinically established CMR routine approach. This opens the door towards combing routine CMR in patients with development of advanced 7T technology.

Sensitivity analysis for interpretation of machine learning based segmentation models in cardiac MRI

BACKGROUND

Image segmentation is a common task in medical imaging e.g., for volumetry analysis in cardiac MRI. Artificial neural networks are used to automate this task with performance similar to manual operators. However, this performance is only achieved in the narrow tasks networks are trained on. Performance drops dramatically when data characteristics differ from the training set properties. Moreover, neural networks are commonly considered black boxes, because it is hard to understand how they make decisions and why they fail. Therefore, it is also hard to predict whether they will generalize and work well with new data. Here we present a generic method for segmentation model interpretation. Sensitivity analysis is an approach where model input is modified in a controlled manner and the effect of these modifications on the model output is evaluated. This method yields insights into the sensitivity of the model to these alterations and therefore to the importance of certain features on segmentation performance.

RESULTS

We present an open-source Python library (misas), that facilitates the use of sensitivity analysis with arbitrary data and models. We show that this method is a suitable approach to answer practical questions regarding use and functionality of segmentation models. We demonstrate this in two case studies on cardiac magnetic resonance imaging. The first case study explores the suitability of a published network for use on a public dataset the network has not been trained on. The second case study demonstrates how sensitivity analysis can be used to evaluate the robustness of a newly trained model.

CONCLUSIONS

Sensitivity analysis is a useful tool for deep learning developers as well as users such as clinicians. It extends their toolbox, enabling and improving interpretability of segmentation models. Enhancing our understanding of neural networks through sensitivity analysis also assists in decision making. Although demonstrated only on cardiac magnetic resonance images this approach and software are much more broadly applicable.