Fast online-customized (FOCUS) parallel transmission pulses: A combination of universal pulses and individual optimization

Unsupervised deep learning with convolutional neural networks for static parallel transmit design: A retrospective study

PURPOSE

To mitigateB 1 + $$ {B}_1^{+} $$ inhomogeneity at 7T for multi-channel transmit arrays using unsupervised deep learning with convolutional neural networks (CNNs).

METHODS

Deep learning parallel transmit (pTx) pulse design has received attention, but such methods have relied on supervised training and did not use CNNs for multi-channelB 1 + $$ {B}_1^{+} $$ maps. In this work, we introduce an alternative approach that facilitates the use of CNNs with multi-channelB 1 + $$ {B}_1^{+} $$ maps while performing unsupervised training. The multi-channelB 1 + $$ {B}_1^{+} $$ maps are concatenated along the spatial dimension to enable shift-equivariant processing amenable to CNNs. Training is performed in an unsupervised manner using a physics-driven loss function that minimizes the discrepancy of the Bloch simulation with the target magnetization, which eliminates the calculation of reference transmit RF weights. The training database comprises 3824 2D sagittal, multi-channelB 1 + $$ {B}_1^{+} $$ maps of the healthy human brain from 143 subjects.B 1 + $$ {B}_1^{+} $$ data were acquired at 7T using an 8Tx/32Rx head coil. The proposed method is compared to the unregularized magnitude least-squares (MLS) solution for the target magnetization in static pTx design.

RESULTS

The proposed method outperformed the unregularized MLS solution for RMS error and coefficient-of-variation and had comparable energy consumption. Additionally, the proposed method did not show local phase singularities leading to distinct holes in the resulting magnetization unlike the unregularized MLS solution.

CONCLUSION

Proposed unsupervised deep learning with CNNs performs better than unregularized MLS in static pTx for speed and robustness.

Fast online spectral-spatial pulse design for subject-specific fat saturation in cervical spine and foot imaging at 1.5 T

OBJECTIVE

To compensate subject-specific field inhomogeneities and enhance fat pre-saturation with a fast online individual spectral-spatial (SPSP) single-channel pulse design.

METHODS

The RF shape is calculated online using subject-specific field maps and a predefined excitation k-space trajectory. Calculation acceleration options are explored to increase clinical viability. Four optimization configurations are compared to a standard Gaussian spectral selective pre-saturation pulse and to a Dixon acquisition using phantom and volunteer (N = 5) data at 1.5 T with a turbo spin echo (TSE) sequence. Measurements and simulations are conducted across various body parts and image orientations.

RESULTS

Phantom measurements demonstrate up to a 3.5-fold reduction in residual fat signal compared to Gaussian fat saturation. In vivo evaluations show improvements up to sixfold for dorsal subcutaneous fat in sagittal cervical spine acquisitions. The versatility of the tailored trajectory is confirmed through sagittal foot/ankle, coronal, and transversal cervical spine experiments. Additional measurements indicate that excitation field (B1) information can be disregarded at 1.5 T. Acceleration methods reduce computation time to a few seconds.

DISCUSSION

An individual pulse design that primarily compensates for main field (B0) inhomogeneities in fat pre-saturation is successfully implemented within an online "push-button" workflow. Both fat saturation homogeneity and the level of suppression are improved.

Dynamic parallel transmit diffusion MRI at 7T

Diffusion MRI (dMRI) is inherently limited by SNR. Scanning at 7 T increases intrinsic SNR but 7 T MRI scans suffer from regions of signal dropout, especially in the temporal lobes and cerebellum. We applied dynamic parallel transmit (pTx) to allow whole-brain 7 T dMRI and compared with circularly polarized (CP) pulses in 6 subjects. Subject-specific 2-spoke dynamic pTx pulses were designed offline for 8 slabs covering the brain. We used vendor-provided B0 and B1+ mapping. Spokes positions were set using the Fourier difference approach, and RF coefficients optimized with a Jacobi-matrix high-flip-angle optimizer. Diffusion data were analyzed with FSL. Comparing whole-brain averages for pTx against CP scans: mean flip angle error improved by 15% for excitation (2-spoke-VERSE 15.7° vs CP 18.4°, P = 0.012) and improved by 14% for refocusing (2-spoke-VERSE 39.7° vs CP 46.2°, P = 0.008). Computed spin-echo signal standard deviation improved by 14% (2-spoke-VERSE 0.185 vs 0.214 CP, P = 0.025). Temporal SNR increased by 5.4% (2-spoke-VERSE 8.47 vs CP 8.04, P = 0.004) especially in the inferior temporal lobes. Diffusion fitting uncertainty decreased by 6.2% for first fibers (2-spoke VERSE 0.0655 vs CP 0.0703, P < 0.001) and 1.3% for second fibers (2-spoke VERSE 0.139 vs CP 0.141, P = 0.01). In conclusion, dynamic parallel transmit improves the uniformity of 7 T diffusion-weighted imaging. In future, less restrictive SAR limits for parallel transmit scans are expected to allow further improvements.

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.

What can we gain from subpopulation universal pulses? A simulation-based study

PURPOSE

The aim of the study was to explore a novel methodology for designing universal pulses (UPs) that balances the benefits of a calibration-free approach with subject-specific online pulse design.

METHODS

The proposed method involves segmenting the population into subpopulations with variability in anatomical shapes and positions reduced to 75%, 50%, and 25% of their original values while keeping the mean values unchanged. An additional 25% extreme case with a large volume of interest and shifted position was included. For each group, a 5kT-points universal inversion pulse was designed and assessed by the normalized root mean square error (NRMSE) on the target longitudinal magnetization profile. The performance was compared to the conventional one-size-fits-all approach. A total of 132 electromagnetic simulations were executed to generate representative anatomies and specific absorption rate (SAR) distributions in a three-dimensional parameter space comprised of head breadth, head length, and Y-shift. The 99.9th percentile on the peak local SAR distribution was utilized to establish an intersubject variability safety margin.

RESULTS

UPs designed for subpopulations with decreased head shape and position variability reduced the anatomical safety margin by up to 20%. Furthermore, when a head was significantly different to the average case, the proposed approach improved the inversion homogeneity by up to 24%, compared to the conventional one-size-fits-all approach.

CONCLUSION

Subpopulation UPs present an opportunity to improve theB 1 + $$ {\mathrm{B}}_1^{+} $$ homogeneity and reduce anatomical SAR safety margins at 7T without additional acquisition time for calibration.

FastPtx: a versatile toolbox for rapid, joint design of pTx RF and gradient pulses using Pytorch's autodifferentiation

OBJECTIVE

With modern optimization methods, free optimization of parallel transmit pulses together with their gradient waveforms can be performed on-line within a short time. A toolbox which uses PyTorch's autodifferentiation for simultaneous optimization of RF and gradient waveforms is presented and its performance is evaluated.

METHODS

MR measurements were performed on a 9.4T MRI scanner using a 3D saturated single-shot turboFlash sequence for [Formula: see text] mapping. RF pulse simulation and optimization were done using a Python toolbox and a dedicated server. An RF- and Gradient pulse design toolbox was developed, including a cost function to balance different metrics and respect hardware and regulatory limits. Pulse performance was evaluated in GRE and MPRAGE imaging. Pulses for non-selective and for slab-selective excitation were designed.

RESULTS

Universal pulses for non-selective excitation reduced the flip angle error to an NRMSE of (12.3±1.7)% relative to the targeted flip angle in simulations, compared to (42.0±1.4)% in CP mode. The tailored pulses performed best, resulting in a narrow flip angle distribution with NRMSE of (8.2±1.0)%. The tailored pulses could be created in only 66 s, making it feasible to design them during an experiment. A 90° pulse was designed as preparation pulse for a satTFL sequence and achieved a NRMSE of 7.1%. We showed that both MPRAGE and GRE imaging benefited from the pTx pulses created with our toolbox.

CONCLUSION

The pTx pulse design toolbox can freely optimize gradient and pTx RF waveforms in a short time. This allows for tailoring high-quality pulses in just over a minute.

7 tricks for 7 T CEST: Improving the reproducibility of multipool evaluation provides insights into the effects of age and the early stages of Parkinson's disease

The objective of the current study was to optimize the postprocessing pipeline of 7 T chemical exchange saturation transfer (CEST) imaging for reproducibility and to prove this optimization for the detection of age differences and differences between patients with Parkinson's disease versus normal subjects. The following 7 T CEST MRI experiments were analyzed: repeated measurements of a healthy subject, subjects of two age cohorts (14 older, seven younger subjects), and measurements of 12 patients with Parkinson's disease. A slab-selective, B 1 + -homogeneous parallel transmit protocol was used. The postprocessing, consisting of motion correction, smoothing, B 0 -correction, normalization, denoising, B 1 + -correction and Lorentzian fitting, was optimized regarding the intrasubject and intersubject coefficient of variation (CoV) of the amplitudes of the amide pool and the aliphatic relayed nuclear Overhauser effect (rNOE) pool within the brain. Seven "tricks" for postprocessing accomplished an improvement of the mean voxel CoV of the amide pool and the aliphatic rNOE pool amplitudes of less than 5% and 3%, respectively. These postprocessing steps are: motion correction with interpolation of the motion of low-signal offsets (1) using the amide pool frequency offset image as reference (2), normalization of the Z-spectrum using the outermost saturated measurements (3), B 0 correction of the Z-spectrum with moderate spline smoothing (4), denoising using principal component analysis preserving the 11 highest intensity components (5), B 1 + correction using a linear fit (6) and Lorentzian fitting using the five-pool fit model (7). With the optimized postprocessing pipeline, a significant age effect in the amide pool can be detected. Additionally, for the first time, an aliphatic rNOE contrast between subjects with Parkinson's disease and age-matched healthy controls in the substantia nigra is detected. We propose an optimized postprocessing pipeline for CEST multipool evaluation. It is shown that by the use of these seven "tricks", the reproducibility and, thus, the statistical power of a CEST measurement, can be greatly improved and subtle changes can be detected.

The effects of RF coils and SAR supervision strategies for clinically applicable nonselective parallel-transmit pulses at 7 T

PURPOSE

To investigate the effects of using different parallel-transmit (pTx) head coils and specific absorption rate (SAR) supervision strategies on pTx pulse design for ultrahigh-field MRI using a 3D-MPRAGE sequence.

METHODS

The PTx universal pulses (UPs) and fast online-customized (FOCUS) pulses were designed with pre-acquired data sets (B₀ , B₁ ⁺ maps, specific absorption rate [SAR] supervision data) from two different 8 transmit/32 receive head coils on two 7T whole-body MR systems. For one coil, the SAR supervision model consisted of per-channel RF power limits. In the other coil, SAR estimations were done with both per-channel RF power limits as well as virtual observation points (VOPs) derived from electromagnetic field (EMF) simulations using three virtual human body models at three different positions. All pulses were made for nonselective excitation and inversion and evaluated on 132 B₀ , B₁ ⁺ , and SAR supervision datasets obtained with one coil and 12 from the other. At both sites, 3 subjects were examined using MPRAGE sequences that used UP/FOCUS pulses generated for both coils.

RESULTS

For some subjects, the UPs underperformed when simulated on a different coil from which they were derived, whereas FOCUS pulses still showed acceptable performance in that case. FOCUS inversion pulses outperformed adiabatic pulses when scaled to the same local SAR level. For the self-built coil, the use of VOPs showed reliable overestimation compared with the ground-truth EMF simulations, predicting about 52% lower local SAR for inversion pulses compared with per-channel power limits.

CONCLUSION

FOCUS inversion pulses offer a low-SAR alternative to adiabatic pulses and benefit from using EMF-based VOPs for SAR estimation.

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.

Rapid estimation of 2D relative B 1 + -maps from localizers in the human heart at 7T using deep learning

PURPOSE

Subject-tailored parallel transmission pulses for ultra-high fields body applications are typically calculated based on subject-specific B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps of all transmit channels, which require lengthy adjustment times. This study investigates the feasibility of using deep learning to estimate complex, channel-wise, relative 2D B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps from a single gradient echo localizer to overcome long calibration times.

METHODS

126 channel-wise, complex, relative 2D B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps of the human heart from 44 subjects were acquired at 7T using a Cartesian, cardiac gradient-echo sequence obtained under breath-hold to create a library for network training and cross-validation. The deep learning predicted maps were qualitatively compared to the ground truth. Phase-only B 1 + $$ {\mathrm{B}}_1^{+} $$ -shimming was subsequently performed on the estimated B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps for a region of interest covering the heart. The proposed network was applied at 7T to 3 unseen test subjects.

RESULTS

The deep learning-based B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps, derived in approximately 0.2 seconds, match the ground truth for the magnitude and phase. The static, phase-only pulse design performs best when maximizing the mean transmission efficiency. In-vivo application of the proposed network to unseen subjects demonstrates the feasibility of this approach: the network yields predicted B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps comparable to the acquired ground truth and anatomical scans reflect the resulting B 1 + $$ {\mathrm{B}}_1^{+} $$ -pattern using the deep learning-based maps.

CONCLUSION

The feasibility of estimating 2D relative B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps from initial localizer scans of the human heart at 7T using deep learning is successfully demonstrated. Because the technique requires only sub-seconds to derive channel-wise B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps, it offers high potential for advancing clinical body imaging at ultra-high fields.

Parkinson's disease or multiple system atrophy: potential separation by quantitative susceptibility mapping

Background

Due to the absence of robust biomarkers, and the low sensitivity and specificity of routine imaging techniques, the differential diagnosis between Parkinson's disease (PD) and multiple system atrophy (MSA) is challenging. High-field magnetic resonance imaging (MRI) opened up new possibilities regarding the analysis of pathological alterations associated with neurodegenerative processes. Recently, we have shown that quantitative susceptibility mapping (QSM) enables visualization and quantification of two major histopathologic hallmarks observed in MSA: reduced myelin density and iron accumulation in the basal ganglia of a transgenic murine model of MSA. It is therefore emerging as a promising imaging modality on the differential diagnosis of Parkinsonian syndromes.

Objectives

To assess QSM on high-field MRI for the differential diagnosis of PD and MSA.

Methods

We assessed 23 patients (nine PDs and 14 MSAs) and nine controls using QSM on 3T and 7T MRI scanners at two academic centers.

Results

We observed increased susceptibility in MSA at 3T in prototypical subcortical and brainstem regions. Susceptibility measures of putamen, pallidum, and substantia nigra reached excellent diagnostic accuracy to separate both synucleinopathies. Increase toward 100% sensitivity and specificity was achieved using 7T MRI in a subset of patients. Magnetic susceptibility correlated with age in all groups, but not with disease duration in MSA. Sensitivity and specificity were particularly high for possible MSA, and reached 100% in the putamen.

Conclusion

Putaminal susceptibility measures, in particular on ultra-high-field MRI, may distinguish MSA patients from both, PD and controls, allowing an early and sensitive diagnosis of MSA.

Ultra-high field MRI: parallel-transmit arrays and RF pulse design

This paper reviews the field of multiple or parallel radiofrequency (RF) transmission for magnetic resonance imaging (MRI). Currently the use of ultra-high field (UHF) MRI at 7 tesla and above is gaining popularity, yet faces challenges with non-uniformity of the RF field and higher RF power deposition. Since its introduction in the early 2000s, parallel transmission (pTx) has been recognized as a powerful tool for accelerating spatially selective RF pulses and combating the challenges associated with RF inhomogeneity at UHF. We provide a survey of the types of dedicated RF coils used commonly for pTx and the important modeling of the coil behavior by electromagnetic (EM) field simulations. We also discuss the additional safety considerations involved with pTx such as the specific absorption rate (SAR) and how to manage them. We then describe the application of pTx with RF pulse design, including a practical guide to popular methods. Finally, we conclude with a description of the current and future prospects for pTx, particularly its potential for routine clinical use.

Clipped DeepControl: Deep neural network two-dimensional pulse design with an amplitude constraint layer

Advanced radio-frequency pulse design used in magnetic resonance imaging has recently been demonstrated with deep learning of (convolutional) neural networks and reinforcement learning. For two-dimensionally selective radio-frequency pulses, the (convolutional) neural network pulse prediction time (a few milliseconds) was in comparison more than three orders of magnitude faster than the conventional optimal control computation. The network pulses were from the supervised training capable of compensating scan-subject dependent inhomogeneities of B₀ and B₁⁺ fields. Unfortunately, the network presented with a small percentage of pulse amplitude overshoots in the test subset, despite the optimal control pulses used in training were fully constrained. Here, we have extended the convolutional neural network with a custom-made clipping layer that completely eliminates the risk of pulse amplitude overshoots, while preserving the ability to compensate for the inhomogeneous field conditions.

Parallel transmission 2D RARE imaging at 7T with transmit field inhomogeneity mitigation and local SAR control

PURPOSE

We develop and test a parallel transmit (pTx) pulse design framework to mitigate transmit field inhomogeneity with control of local specific absorption rate (SAR) in 2D rapid acquisition with relaxation enhancement (RARE) imaging at 7T.

METHODS

We design large flip angle RF pulses with explicit local SAR constraints by numerical simulation of the Bloch equations. Parallel computation and analytical expressions for the Jacobian and the Hessian matrices are employed to reduce pulse design time. The refocusing-excitation "spokes" pulse pairs are designed to satisfy the Carr-Purcell-Meiboom-Gill (CPMG) condition using a combined magnitude least squares-least squares approach.

RESULTS

In a simulated dataset, the proposed approach reduced peak local SAR by up to 56% for the same level of refocusing uniformity error and reduced refocusing uniformity error by up to 59% (from 32% to 7%) for the same level of peak local SAR compared to the circularly polarized birdcage mode of the pTx array. Using explicit local SAR constraints also reduced peak local SAR by up to 46% compared to an RF peak power constrained design. The excitation and refocusing uniformity error were reduced from 20%-33% to 4%-6% in single slice phantom experiments. Phantom experiments demonstrated good agreement between the simulated excitation and refocusing uniformity profiles and experimental image shading.

CONCLUSION

PTx-designed excitation and refocusing CPMG pulse pairs can mitigate transmit field inhomogeneity in the 2D RARE sequence. Moreover, local SAR can be decreased significantly using pTx, potentially leading to better slice coverage, enabling larger flip angles or faster imaging.

Clinically compatible subject-specific dynamic parallel transmit pulse design for homogeneous fat saturation and water-excitation at 141657T: Proof-of-concept for 14165CEST MRI of the brain

PURPOSE

To evaluate the benefits and challenges of dynamic parallel transmit (pTx) pulses for fat saturation (FS) and water-excitation (WE), in the context of CEST MRI.

METHODS

"Universal" k_T -points (for FS) and spiral non-selective (for WE) trajectories were optimized offline for flip angle (FA) homogeneity. Routines to optimize the pulse shape online, based on the subject's fields maps, were implemented (target FA of 110°/0° for FS, 0°/5° for WE at fat/water frequencies). The pulses were inserted in a CEST sequence with a pTx readout. The different fat suppression schemes and their effects on CEST contrasts were compared in 12 volunteers at 7T.

RESULTS

With a 25%-shorter pulse duration, pTx FS largely improved the FA homogeneity (root-mean-square-error (RMSE) = 12.3° vs. 53.4° with circularly-polarized mode, at the fat frequency). However, the spectral selectivity was degraded mainly in the cerebellum and close to the sinuses (RMSE = 5.8° vs. 0.2° at the water frequency). Similarly, pTx WE showed a trade-off between FA homogeneity and spectral selectivity compared to pTx non-selective pulses (RMSE = 0.9° and 1.1° at the fat and water frequencies, vs. 4.6° and 0.5°). In the brain, CEST metrics were reduced by up to 31.9% at -3.3 ppm with pTx FS, suggesting a mitigated lipid-induced bias.

CONCLUSION

This clinically compatible implementation of dynamic pTx pulses improved the fat suppression homogeneity at 7T taking into account the subject-specific B₀ heterogeneities online. This study highlights the lipid-induced biases on the CEST z-spectrum. The results are promising for body applications where B₀ heterogeneities and fat are more substantial.

Optimized ultrahigh field parallel transmission workflow using rapid presaturated TurboFLASH transmit field mapping with a three-dimensional centric single-shot readout

PURPOSE

To evaluate the usage of three-dimensional (3D) presaturated TurboFLASH (satTFL) for B 1 + $$ {\mathrm{B}}_1^{+} $$ and B 0 $$ {\mathrm{B}}_0 $$ mapping on single channel and parallel transmission (pTx) systems.

METHODS

B 1 + $$ {\mathrm{B}}_1^{+} $$ maps recorded with 3D satTFL were compared to maps from three other 3D B 1 + $$ {\mathrm{B}}_1^{+} $$ mapping sequences in an agar phantom. Furthermore, individual-channel B 1 + $$ {\mathrm{B}}_1^{+} $$ maps of 18 human subjects were recorded with 3D satTFL using B 1 + $$ {\mathrm{B}}_1^{+} $$ interferometry. A neural network was trained for masking of the maps.

RESULTS

Out of the sequences compared satTFL was the only one with a mapping range exceeding well over 90°. In regions with lower flip angles there was high correspondence between satTFL and AFI. DREAM and double angle method also showed high qualitative similarity, however the magnitude differed from the other two measurements. The individual-channel B 1 + $$ {\mathrm{B}}_1^{+} $$ maps were successfully used for pTx pulse calculation in a separate study.

CONCLUSION

3D satTFL can record high-quality B 1 + $$ {\mathrm{B}}_1^{+} $$ maps with a high dynamic range in a short time. Correspondence with AFI maps is high, while measurement duration is reduced drastically.

Interleaved binomial kT-Points for water-selective imaging at 7T

PURPOSE

We present a time-efficient water-selective, parallel transmit RF excitation pulse design for ultra-high field applications.

METHODS

The proposed pulse design method achieves flip angle homogenization at ultra-high fields by employing spatially nonselective k T $$ {\mathrm{k}}_T $$ -points pulses. In order to introduce water-selection, the concept of binomial pulses is applied. Due to the composite nature of k T $$ {\mathrm{k}}_T $$ -points, the pulse can be split into multiple binomial subpulse blocks shorter than half the precession period of fat, that are played out successively. Additional fat precession turns, that would otherwise impair the spectral response, can thus be avoided. Bloch simulations of the proposed interleaved binomial k T $$ {\mathrm{k}}_T $$ -points pulses were carried out and compared in terms of duration, homogeneity, fat suppression and pulse energy. For validation, in vivo MP-RAGE and 3D-EPI data were acquired.

RESULTS

Simulation results show that interleaved binomial k T $$ {\mathrm{k}}_T $$ -points pulses achieve shorter total pulse durations, improved flip angle homogeneity and more robust fat suppression compared to available methods. Interleaved binomial k T $$ {\mathrm{k}}_T $$ -points can be customized by changing the number of k T $$ {\mathrm{k}}_T $$ -points, the subpulse duration and the order of the binomial pulse. Using shorter subpulses, the number of k T $$ {\mathrm{k}}_T $$ -points can be increased and hence better homogeneity is achieved, while still maintaining short total pulse durations. Flip angle homogenization and fat suppression of interleaved binomial k T $$ {\mathrm{k}}_T $$ -points pulses is demonstrated in vivo at 7T, confirming Bloch simulation results.

CONCLUSION

In this work, we present a time efficient and robust parallel transmission technique for nonselective water excitation with simultaneous flip angle homogenization at ultra-high field.

OTUP workflow: target specific optimization of the transmit k-space trajectory for flexible universal parallel transmit RF pulse design

PURPOSE

To optimize transmit k-space trajectories for a wide range of excitation targets and to design "universal pTx RF pulses" based on these trajectories.

METHODS

Transmit k-space trajectories (stack of spirals and SPINS) were optimized to best match different excitation targets using the parameters of the analytical equations of spirals and SPINS. The performances of RF pulses designed based on optimized and non-optimized trajectories were compared. The optimized trajectories were utilized for universal pulse design. The universal pulse performances were compared with subject specific tailored pulse performances. The OTUP workflow (optimization of transmit k-space trajectories and universal pulse calculation) was tested on three test target excitation patterns. For one target (local excitation of a central area in the human brain) the pulses were tested in vivo at 9.4 T.

RESULTS

The workflow produced appropriate transmit k-space trajectories for each test target. Utilization of an optimized trajectory was crucial for the pulse performance. Using unsuited trajectories diminished the performance. It was possible to create target specific universal pulses. However, not every test target is equally well suited for universal pulse design. There was no significant difference in the in vivo performance between subject specific tailored pulses and a universal pulse at 9.4 T.

CONCLUSIONS

The proposed workflow further exploited and improved the universal pulse concept by combining it with gradient trajectory optimization for stack of spirals and SPINS. It emphasized the importance of a well suited trajectory for pTx RF pulse design. Universal and tailored pulses performed with a sufficient degree of similarity in simulations and a high degree of similarity in vivo. The implemented OTUP workflow and the B₀ /B₁⁺ map data from 18 subjects measured at 9.4 T are available as open source (https://github.com/ole1965/workflow_OTUP.git).

The Impact of Respiratory Motion on Electromagnetic Fields and Specific Absorption Rate in Cardiac Imaging at 7T

PURPOSE

To present electromagnetic simulation setups for detailed analyses of respiration's impact on B 1 + $$ {B}_1^{+} $$ and E-fields, local specific absorption rate (SAR) and associated safety-limits for 7T cardiac imaging.

METHODS

Finite-difference time-domain electromagnetic field simulations were performed at five respiratory states using a breathing body model and a 16-element 7T body transceiver RF-coil array. B 1 + $$ {B}_1^{+} $$ and SAR are analyzed for fixed and moving coil configurations. SAR variations are investigated using phase/amplitude shimming considering (i) a local SAR-controlled mode (here SAR calculations consider RF amplitudes and phases) and (ii) a channel-wise power-controlled mode (SAR boundary calculation is independent of the channels' phases, only dependent on the channels' maximum amplitude).

RESULTS

Respiration-induced variations of both B 1 + $$ {B}_1^{+} $$ amplitude and phase are observed. The flip angle homogeneity depends on the respiratory state used for B 1 + $$ {B}_1^{+} $$ shimming; best results were achieved for shimming on inhale and exhale simultaneously ( | Δ C V | < 35 % $$ \mid \Delta CV\mid <35\% $$ ). The results reflect that respiration impacts position and amplitude of the local SAR maximum. With the local-SAR-control mode, a safety factor of up to 1.4 is needed to accommodate for respiratory variations while the power control mode appears respiration-robust when the coil moves with respiration (SAR peak decrease: 9% exhale→inhale). Instead, a spatially fixed coil setup yields higher SAR variations with respiration.

CONCLUSION

Respiratory motion does not only affect the B 1 + $$ {B}_1^{+} $$ distribution and hence the image contrast, but also location and magnitude of the peak spatial SAR. Therefore, respiration effects may need to be included in safety analyses of RF coils applied to the human thorax.

B1 field map synthesis with generative deep learning used in the design of parallel-transmit RF pulses for ultra-high field MRI

Spoke trajectory parallel transmit (pTX) excitation in ultra-high field MRI enables B₁⁺ inhomogeneities arising from the shortened RF wavelength in biological tissue to be mitigated. To this end, current RF excitation pulse design algorithms either employ the acquisition of field maps with subsequent non-linear optimization or a universal approach applying robust pre-computed pulses. We suggest and evaluate an intermediate method that uses a subset of acquired field maps combined with generative machine learning models to reduce the pulse calibration time while offering more tailored excitation than robust pulses (RP). The possibility of employing image-to-image translation and semantic image synthesis machine learning models based on generative adversarial networks (GANs) to deduce the missing field maps is examined. Additionally, an RF pulse design that employs a predictive machine learning model to find solutions for the non-linear (two-spokes) pulse design problem is investigated. As a proof of concept, we present simulation results obtained with the suggested machine learning approaches that were trained on a limited data-set, acquired in vivo. The achieved excitation homogeneity based on a subset of half of the B₁⁺ maps acquired in the calibration scans and half of the B₁⁺ maps synthesized with GANs is comparable with state of the art pulse design methods when using the full set of calibration data while halving the total calibration time. By employing RP dictionaries or machine-learning RF pulse predictions, the total calibration time can be reduced significantly as these methods take only seconds or milliseconds per slice, respectively.

Standardized universal pulse: A fast RF calibration approach to improve flip angle accuracy in parallel transmission

PURPOSE

In parallel transmission (pTX), subject-tailored RF pulses allow achieving excellent flip angle (FA) accuracy but often require computationally extensive online optimizations, precise characterization of the static field ( Δ B 0 ), and the transmit RF field ( B 1 + ) distributions. This costs time and requires expertise from the MR user. Universal pulses (UPs) have been proposed to reduce this burden, yet, with a penalty in FA accuracy. This study introduces the concept of standardized universal pulses (SUPs), where pulses are designed offline and adjusted to the subject through a fast online calibration scan.

METHODS

A SUP is designed offline using a so-called standardized database, wherein each B 1 + map has been normalized to a reference transmit RF field distribution. When scanning a new subject, a 3-slice B 1 + acquisition (scan time < 10  s) is performed and used to adjust the SUP to the subject through a linear transform. SUP performance was assessed at 7T with simulations by computing the FA-normalized root mean square error (FA-NRMSE) and the FA pattern stability as measured by the average and coefficient of variation of the FA across 15 control subjects, along with in vivo experiments using an MP2RAGE sequence implementing the SUP variant for the FLASH readout.

RESULTS

Adjusted SUP improved the FA-NRMSE (8.8 % for UP vs. 7.1 % for adjusted SUP). Experimentally in vivo, this translated in an improved signal homogeneity and more accurate T 1 quantification using MP2RAGE.

CONCLUSION

The proposed SUP approach improves excitation accuracy (FA-NRMSE) while preserving the same offline pulse design principle as offered by UPs.

Rigid motion-resolved B 1 + prediction using deep learning for real-time parallel-transmission pulse design

Purpose

Tailored parallel-transmit (pTx) pulses produce uniform excitation profiles at 7 T, but are sensitive to head motion. A potential solution is real-time pulse redesign. A deep learning framework is proposed to estimate pTx B₁⁺ distributions following within-slice motion, which can then be used for tailored pTx pulse redesign.

Methods

Using simulated data, conditional generative adversarial networks were trained to predict B₁⁺ distributions in the head following a displacement. Predictions were made for two virtual body models that were not included in training. Predicted maps were compared with groundtruth (simulated, following motion) B₁⁺ maps. Tailored pTx pulses were designed using B₁⁺ maps at the original position (simulated, no motion) and evaluated using simulated B₁⁺ maps at displaced position (ground-truth maps) to quantify motion-related excitation error. A second pulse was designed using predicted maps (also evaluated on ground-truth maps) to investigate improvement offered by the proposed method.

Results

Predicted B₁⁺ maps corresponded well with ground-truth maps. Error in predicted maps was lower than motion-related error in 99% and 67% of magnitude and phase evaluations, respectively. Worst-case flip-angle normalized RMS error due to motion (76% of target flip angle) was reduced by 59% when pulses were redesigned using predicted maps.

Conclusion

We propose a framework for predicting B₁⁺ maps online with deep neural networks. Predicted maps can then be used for real-time tailored pulse redesign, helping to overcome head motion–related error in pTx.

7 Tesla and Beyond: Advanced Methods and Clinical Applications in Magnetic Resonance Imaging

ABSTRACT

Ultrahigh magnetic fields offer significantly higher signal-to-noise ratio, and several magnetic resonance applications additionally benefit from a higher contrast-to-noise ratio, with static magnetic field strengths of B0 ≥ 7 T currently being referred to as ultrahigh fields (UHFs). The advantages of UHF can be used to resolve structures more precisely or to visualize physiological/pathophysiological effects that would be difficult or even impossible to detect at lower field strengths. However, with these advantages also come challenges, such as inhomogeneities applying standard radiofrequency excitation techniques, higher energy deposition in the human body, and enhanced B0 field inhomogeneities. The advantages but also the challenges of UHF as well as promising advanced methodological developments and clinical applications that particularly benefit from UHF are discussed in this review article.

Hybrid B 1 + -shimming and gradient adaptions for improved pseudo-continuous arterial spin labeling at 7 Tesla

PURPOSE

To improve pseudo-continuous arterial spin labeling (pcASL) at 7T by exploiting a hybrid homogeneity- and efficiency-optimized B 1 + -shim with adapted gradient strength as well as background suppression.

METHODS

The following three experiments were performed at 7T, each employing five volunteers: (1) A hybrid (ie, homogeneity-efficiency optimized) B 1 + -shim was introduced and evaluated for variable-rate selective excitation pcASL labeling. Therefore, B 1 + -maps in the V3 segment and time-of-flight images were acquired to identify the feeding arteries. For validation, a gradient-echo sequence was applied in circular polarized (CP) mode and with the hybrid B 1 + -shim. Additionally, the gray matter (temporal) signal-to-noise ratio (tSNR) in pcASL perfusion images was evaluated. (2) Bloch simulations for the pcASL labeling were conducted and validated experimentally, with a focus on the slice-selective gradients. (3) Background suppression was added to the B 1 + -shimmed, gradient-adapted 7T sequence and this was then compared to a matched sequence at 3T.

RESULTS

The B 1 + -shim improved the signal within the labeling plane (23.6%) and the SNR/tSNR increased (+11%/+11%) compared to its value in CP mode; however, the increase was not significant. In accordance with the simulations, the adapted gradients increased the tSNR (35%) and SNR (45%) significantly. Background suppression further improved the perfusion images at 7T, and this protocol performed as well as a resolution-matched protocol at 3T.

CONCLUSION

The combination of the proposed hybrid B 1 + -phase-shim with the adapted slice-selective gradients and background suppression shows great potential for improved pcASL labeling under suboptimal B 1 + conditions at 7T.

Calibration-free pTx of the human heart at 7T via 3D universal pulses

PURPOSE

MRI at ultra-high fields in the human body is highly challenging and requires lengthy calibration times to compensate for spatially heterogeneous B 1 + profiles. This study investigates the feasibility of using pre-computed universal pulses for calibration-free homogeneous 3D flip angle distribution in the human heart at 7T.

METHODS

Twenty-two channel-wise 3D B 1 + data sets were acquired under free-breathing in 19 subjects to generate a library for an offline universal pulse (UP) design (group 1: 12 males [M] and 7 females [F], 21-66 years, 19.8-28.3 kg/m² ). Three of these subjects (2M/1F, 21-33 years, 20.8-23.6 kg/m² ) were re-scanned on different days. A 4kT-points UP optimized for the 22 channel-wise 3D B 1 + data sets in group 1 (UP22-4kT) is proposed and applied at 7T in 9 new and unseen subjects (group 2: 4M/5F, 25-56 years, 19.5-35.3 kg/m² ). Multiple tailored and universal static and dynamic parallel-transmit (pTx) pulses were designed and evaluated for different permutations of the B 1 + data sets in group 1 and 2.

RESULTS

The proposed UP22-4kT provides low B 1 + variation in all subjects, seen and unseen, without severe signal drops. Experimental data at 7T acquired with UP22-4kT shows comparable image quality as data acquired with tailored-4kT pulses and demonstrates successful calibration-free pTx of the human heart.

CONCLUSION

UP22-4kT allows for calibration-free homogeneous flip angle distributions across the human heart at 7T. Large inter-subject variations because of sex, age, and body mass index are well tolerated. The proposed universal pulse removes the need for lengthy (10-15 min) calibration scans and therefore has the potential to bring body imaging at 7T closer to the clinical application.

Brain tissues have single-voxel signatures in multi-spectral MRI

Since the seminal works by Brodmann and contemporaries, it is well-known that different brain regions exhibit unique cytoarchitectonic and myeloarchitectonic features. Transferring the approach of classifying brain tissues - and other tissues - based on their intrinsic features to the realm of magnetic resonance (MR) is a longstanding endeavor. In the 1990s, atlas-based segmentation replaced earlier multi-spectral classification approaches because of the large overlap between the class distributions. Here, we explored the feasibility of performing global brain classification based on intrinsic MR features, and used several technological advances: ultra-high field MRI, q-space trajectory diffusion imaging revealing voxel-intrinsic diffusion properties, chemical exchange saturation transfer and semi-solid magnetization transfer imaging as a marker of myelination and neurochemistry, and current neural network architectures to analyze the data. In particular, we used the raw image data as well to increase the number of input features. We found that a global brain classification of roughly 97 brain regions was feasible with gross classification accuracy of 60%; and that mapping from voxel-intrinsic MR data to the brain region to which the data belongs is possible. This indicates the presence of unique MR signals of different brain regions, similar to their cytoarchitectonic and myeloarchitectonic fingerprints.

Whole-brain quantitative CEST MRI at 7T using parallel transmission methods and B 1 + correction

PURPOSE

To enable whole-brain quantitative CEST MRI at ultra-high magnetic field strengths (B₀ ≥ 7T) within short acquisition times.

METHODS

Multiple interleaved mode saturation (MIMOSA) was combined with fast online-customized (FOCUS) parallel transmission (pTx) excitation pulses and B 1 + correction to achieve homogenous whole-brain coverage. Examinations of 13 volunteers were performed on a 7T MRI system with 3 different types of pulse sequences: (1) saturation in circular polarized (CP) mode and CP mode readout, (2) MIMOSA and CP readout, and (3) MIMOSA and FOCUS readout. For comparison, the inverse magnetic transfer ratio metric for relayed nuclear Overhauser effect and amide proton transfer were calculated. To investigate the number of required acquisitions for a good B 1 + correction, 4 volunteers were measured with 6 different B₁ amplitudes. Finally, time point repeatability was investigated for 6 volunteers.

RESULTS

MIMOSA FOCUS sequence using B 1 + correction, with both single and multiple points, reduced inhomogeneity of the CEST contrasts around the occipital lobe and cerebellum. Results indicate that the most stable inter-subject coefficient of variation was achieved using the MIMOSA FOCUS sequence. Time point repeatability of MIMOSA FOCUS with single-point B 1 + correction showed a maximum coefficient of variation below 8% for 3 measurements in a single volunteer.

CONCLUSION

A combination of MIMOSA FOCUS with a single-point B 1 + correction can be used to achieve quantitative CEST measurements at ultra-high magnetic field strengths. Compared to previous B 1 + correction methods, acquisition time can be reduced as additional scans required for B 1 + correction can be omitted.