Search

Purpose The aim of the study was to investigate whether incorrectly compensated eddy currents are the source of persistent X-nuclear spectroscopy and imaging artifacts, as well as methods to correct this. Methods Pulse-acquire spectra were collected for 1 H and X-nuclei (23 Na or 31 P) using the minimum TR permitted on a 3T clinical MRI system. Data were collected in 3 orientations (axial, sagittal, and coronal) with the spoiler gradient at the end of the TR applied along the slice direction for each. Modifications to system calibration files to tailor eddy current compensation for each X-nucleus were developed and applied, and data were compared with and without these corrections for: slice-selective MRS (for 23 Na and 31 P), 2D spiral trajectories (for 13 C), and 3D cones trajectories (for 23 Na). Results Line-shape distortions characteristic of eddy currents were demonstrated for X-nuclei, which were not seen for 1 H. The severity of these correlated with the amplitude of the eddy current frequency compensation term applied by the system along the axis of the applied spoiler gradient. A proposed correction to eddy current compensation, taking account of the gyromagnetic ratio, was shown to dramatically reduce these distortions. The same correction was also shown to improve data quality of non-Cartesian imaging (2D spiral and 3D cones trajectories). Conclusion A simple adaptation of the default compensation for eddy currents was shown to eliminate a range of artifacts detected on X-nuclear spectroscopy and imaging.

PURPOSE: To describe and implement a strategy for dynamic slice-by-slice and multiband B0 shimming using spherical harmonic shims in the human brain at 7T. THEORY: For thin axial slices, spherical harmonic shims can be divided into pairs of shims (z-degenerate and non-z-degenerate) that are spatially degenerate, such that only ½ of the shims (non-z-degenerate) are required for single slice optimizations. However, when combined, the pairs of shims can be used to simultaneously generate the same in-plane symmetries but with different amplitudes as a function of their z location. This enables multiband shimming equivalent to that achievable by single slice-by-slice optimization. METHODS: All data were acquired at 7T using a spherical harmonic shim insert enabling shimming up through 4th order with two additional 5th order shims (1st-4th+). Dynamic shim updating was achieved using a 10A shim power supply with 2 ms ramps and constrained optimizations to minimize eddy currents. RESULTS: In groups of eight subjects, we demonstrated that: 1) dynamic updating using 1st-4th+ order shims reduced the SD of the B0 field over the whole brain from 32.4 ± 2.6 and 24.9 ± 2 Hz with 1st-2nd and 1st-4th+ static global shimming to 15.1 ± 1.7 Hz; 2) near equivalent performance was achieved when dynamically updating only the non-z-degenerate shims (14.3 ± 1.5 Hz), or when a using multiband shim factor of 2, MBs = 2, and all shims (14.4 ± 2.0 Hz). CONCLUSION: High order spherical harmonics provide substantial improvements over static global shimming and enable dynamic multiband shimming with near equivalent performance to that of dynamic slice-by-slice shimming. This reduces distortion in echo planar imaging.

Purpose The impact of gradient imperfections on UTE images and UTE image-derived bone water quantification was investigated at 3 T field strength. Methods The effects of simple gradient time delays and eddy currents on UTE images, as well as the effects of gradient error corrections, were studied with simulation and phantom experiments. The k-space trajectory was mapped with a 2D sequence with phase encoding on both spatial axes by measuring the phase of the signal in small time increments during ramp-up of the read gradient. In vivo 3D UTE images were reconstructed with and without gradient error compensation to determine the bias in bone water quantification. Finally, imaging was performed on 2 equally configured Siemens TIM Trio systems (Siemens Medical Solutions, Erlangen, Germany) to investigate the impact of such gradient imperfections on inter-scanner measurement bias. Results Compared to values derived from UTE images with full gradient error compensation, total bone water was found to deviate substantially with no (up to 17%) or partial (delay-only) compensation (up to 10.8%). Bound water, obtained with inversion recovery-prepared UTE, was somewhat less susceptible to gradient errors (up to 2.2% for both correction strategies). Inter-scanner comparison indicated a statistically significant bias between measurements from the 2 MR systems for both total and bound water, which either vanished or was substantially reduced following gradient error correction. Conclusion Gradient imperfections impose spatially dependent artifacts on UTE images, which compromise not only bone water quantification accuracy but also inter-scanner measurement agreement if left uncompensated.

Purpose To propose an explicit Balanced steady-state free precession (bSSFP) signal model that predicts eddy current-induced steady-state disruptions and to provide a prospective, practical, and general eddy current compensation method. Theory and methods Gradient impulse response functions (GIRF) were used to simulate trajectory-specific eddy current-induced phase errors at the end of a repetition block. These phase errors were included in bloch simulations to establish a bSSFP signal model to predict steady-state disruptions and their corresponding image artifacts. The signal model was embedded in the MR system and used to compensate the phase errors by prospectively modifying the phase cycling scheme of the RF pulse. The signal model and eddy current compensation method were validated in phantom and in vivo experiments. In addition, the signal model was used to analyze pre-existing eddy current mitigation methods, such as 2D tiny golden angle radial and 3D paired phase encoded Cartesian acquisitions. Results The signal model predicted eddy current-induced image artifacts, with the zeroth-order GIRF being the primary factor to predict the steady-state disruption. Prospective RF phase cycling schemes were automatically computed online and considerably reduced eddy current-induced image artifacts. The signal model provides a direct relationship for the smoothness of k-space trajectories, which explains the effectiveness of phase encode pairing and tiny golden angle trajectory. Conclusions The proposed signal model can accurately predict eddy current-induced steady-state disruptions for bSSFP imaging. The signal model can be used to derive the eddy current-induced phase errors required for trajectory-specific RF phase cycling schemes, which considerably reduce eddy current-induced image artifacts.

Purpose To perform direct, selective MRI of short-T2 tissues using zero echo time (ZTE) imaging with weighted echo subtraction (WSUB). Methods Radial imaging was performed at 7T, acquiring both ZTE and gradient echo (GRE) signals created by bipolar gradients. Long-T2 suppression was achieved by weighted subtraction of ZTE and GRE images. Special attention was given to imperfections of gradient dynamics, to which radial GRE imaging is particularly susceptible. To compensate for gradient errors, matching of gradient history was combined with data correction based on trajectory measurement. The proposed approach was first validated in phantom experiments and then demonstrated in musculoskeletal (MSK) imaging. Results Trajectory analysis and phantom imaging demonstrated that gradient imperfections were successfully addressed. Gradient history matching enabled consistency between antiparallel projections as required for deriving zeroth-order eddy current dynamics. Trajectory measurement provided individual echo times per projection that showed considerable variation between gradient directions. In in vivo imaging of knee, ankle, and tibia, the proposed approach enabled high-resolution 3D depiction of bone, tendons, and ligaments. Distinct contrast of these structures indicates excellent selectivity of long-T2 suppression. Clarity of depiction also confirmed sufficient SNR of short-T2 tissues, achieved by high baseline sensitivity at 7T combined with high SNR efficiency of ZTE acquisition. Conclusion Weighted subtraction of ZTE and GRE data reconciles robust long-T2 suppression with fastest k-space coverage and high SNR efficiency. This approach enables high-resolution imaging with excellent selectivity to short-T2 tissues, which are of major interest in MSK and neuroimaging applications.

Purpose To overcome existing challenges in dynamic B0 shimming by implementing a shim optimization algorithm which limits shim current amplitudes and their temporal variation through the application of constraints and regularization terms. Theory and methods Spherical harmonic dynamic B0 shimming is complicated by eddy currents, ill-posed optimizations, and the need for strong power supplies. Based on the fact that eddy current amplitudes are proportional to the magnitude of the shim current changes, and assuming a smoothness of the B0 inhomogeneity variation in the slice direction, a novel algorithm was implemented to reduce eddy current generation by limiting interslice shim current changes. Shim degeneracy issues and resulting high current amplitudes are additionally addressed by penalizing high solution norms. Applicability of the proposed algorithm was validated in simulations and in phantom and in vivo measurements. Results High-order dynamic shimming simulations and measurements have shown that absolute shim current amplitudes and their temporal variation can be substantially reduced with negligible loss in achievable B0 homogeneity. Whereas conventional dynamic shim updating optimizations improve the B0 homogeneity, on average, by a factor of 2.1 over second-order static solutions, our proposed routine reached a factor of 2.0, while simultaneously providing a 14-fold reduction of the average maximum shim current changes. Conclusions The proposed algorithm substantially reduces the shim amplitudes and their temporal variation, while only marginally affecting the achievable B0 homogeneity. As a result, it has the potential to mitigate the remaining challenges in dynamic B0 shimming and help in making its application more readily available.

PURPOSE Velocity-selective saturation (VSS) pulse trains provide a viable alternative to the spatially selective methods for measuring cerebral blood volume (CBV) by reducing the sensitivity to arterial transit time. This study is to compare the Fourier-transform-based velocity-selective saturation (FT-VSS) pulse trains with the conventional flow-dephasing VSS techniques for CBV quantification. METHODS The proposed FT-VSS label and control modules were compared with VSS pulse trains utilizing double refocused hyperbolic tangent (DRHT) and 8-segment B1-insensitive rotation (BIR-8). This was done using both numerical simulations and phantom studies to evaluate their sensitivities to gradient imperfections such as eddy currents. DRHT, BIR-8, and FT-VSS prepared CBV mapping was further compared for velocity-encoding gradients along 3 orthogonal directions in healthy subjects at 3T. RESULTS The phantom studies exhibited more consistent immunity to gradient imperfections for the utilized FT-VSS pulse trains. Compared to DRHT and BIR-8, FT-VSS delivered more robust CBV results across the 3 VS encoding directions with significantly reduced artifacts along the superior-inferior direction and improved temporal signal-to-noise ratio (SNR) values. Average CBV values obtained from FT-VSS based sequences were 5.3 mL/100 g for gray matter and 2.3 mL/100 g for white matter, comparable to literature expectations. CONCLUSION Absolute CBV quantification utilizing advanced FT-VSS pulse trains had several advantages over the existing approaches using flow-dephasing VSS modules. A greater immunity to gradient imperfections and the concurrent tissue background suppression of FT-VSS pulse trains enabled more robust CBV measurements and higher SNR than the conventional VSS pulse trains.

Purpose B0 eddy currents are a subtle but important source of artifacts in spiral MRI. This study illustrates the importance of addressing these artifacts and presents a system response-based eddy current correction strategy using B0 eddy current phase measurements on a phantom. Methods B0 and linear eddy current system response measurements were estimated from phantom-based measurement and used to predict residual eddy current effects in spiral acquisitions. The measurements were evaluated across multiple systems and gradient sets. The corresponding eddy current corrections were studied in both axial spiral-in/out TSE and sagittal spiral-out MPRAGE volunteer data. Results Correction of B0 eddy currents using the proposed method mitigated blurriness in the axial spiral-in/out images and artifacts in the sagittal spiral-out images. The system response measurement was found to yield repeatable results over time with some variation in the B0 eddy current responses measured between different systems. Conclusions The proposed eddy current correction framework was effective in mitigating the effects of residual B0 and linear eddy currents. Any spiral acquisition should take residual eddy currents into account. This is particularly important in spiral-in/out acquisitions.

PURPOSE Sensitivity to imperfections of image-encoding gradient fields may significantly impair widespread use of radial MR data acquisition. Such imperfections can cause individual echo shifts for each spoke acquired in the k-space and may produce severe image artifacts. Therefore, fast and robust methods to quantify and correct for those echo shifts are required. THEORY AND METHODS Echo shifts can be induced by inhomogeneities of the static magnetic field (δnB ) and by imbalances of the imaging gradients (δnG ) mainly caused by eddy currents. However, mismatch between data acquisition and gradient switching may additionally play a role. From the position of the echo maxima of 2 opposing spokes, δnG and δnB can be determined and corrected by adapting the read-dephasing gradient accordantly. This approach was implemented on MR-systems of different field strengths, gradient systems, and vendors, and the dependencies of echo shift and acquisition parameters were analyzed. Data sets of phantoms and of mice under in vivo conditions were obtained using RF-spoiled 2D radial-FLASH. RESULTS The presented method allowed for echo-shift detection and correction of

PURPOSE Uneven flow distribution in patients with Fontan circulation is suspected to lead to complications. 4D flow MRI offers evaluation using time-resolved pathlines; however, the potential error is not well understood. The aim of this study was to systematically assess variability in flow distribution caused by well-known sources of error. METHODS 4D flow MRI was acquired in 14 patients with Fontan circulation. Flow distribution was quantified by the % of caval venous flow pathlines reaching the left and right pulmonary arteries. Impact of data acquisition and data processing uncertainties were investigated by (1) probabilistic 4D blood flow tracking at varying noise levels, (2) down-sampling to mimic acquisition at different spatial resolutions, (3) pathline calculation with and without eddy current correction, and (4) varied segmentation of the Fontan geometry to mimic analysis errors. RESULTS Averaged among the cohort, uncertainties accounted for flow distribution errors from noise ≤3.2%, low spatial resolution ≤2.3% to 3.8%, eddy currents ≤6.4%, and inaccurate segmentation ≤3.9% to 9.1% (dilation and erosion, respectively). In a worst-case scenario (maximum additive errors for all 4 sources), flow distribution errors were as high as 22.5%. CONCLUSION Inaccuracies related to postprocessing (segmentation, eddy currents) resulted in the largest potential error (≤15.5% combined) whereas errors related to data acquisition (noise, low spatial resolution) had a lower impact (≤5.5%-7.0% combined). Whereas it is unlikely that these errors will be additive or affect the identification of severe asymmetry, these results illustrate the importance of eddy current correction and accurate segmentation to minimize Fontan flow distribution errors.

PURPOSE To combine the metabolite-cycling technique with diffusion-weighted H-MR spectroscopy and to use the inherent water reference for compensation of motion-related signal loss for improved estimation of metabolite apparent diffusion coefficients (ADCs). METHODS Diffusion-weighted spectra of water and metabolites were acquired simultaneously using metabolite-cycling at 3 T. The water information was used for signal correction of phase, frequency, and eddy currents, as well as for compensation of motion-induced signal loss. ADCs were estimated by 2D simultaneous fitting. The quality of ADC restoration was investigated in vitro. Subsequently, the new approach was applied in 13 subjects for enhanced metabolite ADC estimation in gray matter. RESULTS Metabolite-cycled diffusion H-MRS is suitable to measure metabolite and water ADCs simultaneously. The water reference facilitates signal amplitude restoration, compensating for motion-related artefacts. 2D fitting stabilizes the fitting procedure and allows the estimation of ADCs even for low signal-to-noise metabolites. Use of the motion-compensation scheme leads to estimation of smaller ADCs for virtually all metabolites (44% smaller ADC on average), to a reduction of fitting uncertainties for metabolite ADCs in individual subjects and reduced variance over the cohort (45% smaller SD on average). CONCLUSION Using the simultaneously acquired water signal as internal reference allows not only for compensation of phase and frequency fluctuations but also for signal amplitude restoration, and thus improved metabolite ADC estimation. Combination with 2D simultaneous fitting promises access to the diffusion properties even for low signal-to-noise metabolites. The combination of both techniques increases the specificity and sensitivity of estimated metabolite ADC values in the cohort.

Purpose To investigate the effect of the asymmetric gradient concomitant fields (CF) with zeroth and first-order spatial dependence on fast/turbo spin-echo acquisitions, and to demonstrate the effectiveness of their real-time compensation. Methods After briefly reviewing the CF produced by asymmetric gradients, the effects of the additional zeroth and first-order CFs on these systems are investigated using extended-phase graph simulations. Phantom and in vivo experiments are performed to corroborate the simulation. Experiments are performed before and after the real-time compensations using frequency tracking and gradient pre-emphasis to demonstrate their effectiveness in correcting the additional CFs. The interaction between the CFs and prescan-based correction to compensate for eddy currents is also investigated. Results It is demonstrated that, unlike the second-order CFs on conventional gradients, the additional zeroth/first-order CFs on asymmetric gradients cause substantial signal loss and dark banding in fast spin-echo acquisitions within a typical brain-scan field of view. They can confound the prescan correction for eddy currents and degrade image quality. Performing real-time compensation successfully eliminates the artifacts. Conclusions We demonstrate that the zeroth/first-order CFs specific to asymmetric gradients can cause substantial artifacts, including signal loss and dark bands for brain imaging. These effects can be corrected using real-time compensation. Magn Reson Med, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Purpose Magnetic resonance–compatible medical devices operate within the MR environment while benefitting from the superior anatomic information of MRI. Avoiding electromagnetic interference between such instrumentation and the MR system is crucial. In this work, various shielding configurations for positron emission tomography (PET) detectors were studied and analyzed regarding radiofrequency (RF) shielding effectiveness and gradient-induced eddy current performances. However, the results of this work apply to shielding considerations for any MR-compatible devices. Methods Six shielding enclosure configurations with various thicknesses, patterns, and materials were designed: solid and segmented copper, phosphor bronze mesh (PBM), and carbon fiber composite (CFC). A series of tests was performed on RF shielding effectiveness and the gradient-induced eddy current. Results For the shielding effectiveness, the solid copper with various thickness and PBM configurations yield significantly better shielding effectiveness (>15 dB) compared with CFC and segmented configurations. For the gradient-induced eddy current performance, the solid copper shielding configurations with different thicknesses showed significantly worse results, up to a factor of 3.89 dB, compared with the segmented copper, PBM, and the CFC configurations. Conclusions We evaluated the RF shielding effectiveness and the gradient-induced eddy current artifacts of several shielding designs, and only the PBM showed positive outcomes for both aspects. Magn Reson Med 79:1745–1752, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Purpose To design and evaluate eddy current–nulled convex optimized diffusion encoding (EN-CODE) gradient waveforms for efficient diffusion tensor imaging (DTI) that is free of eddy current–induced image distortions. Methods The EN-CODE framework was used to generate diffusion-encoding waveforms that are eddy current–compensated. The EN-CODE DTI waveform was compared with the existing eddy current–nulled twice refocused spin echo (TRSE) sequence as well as monopolar (MONO) and non–eddy current–compensated CODE in terms of echo time (TE) and image distortions. Comparisons were made in simulations, phantom experiments, and neuro imaging in 10 healthy volunteers. Results The EN-CODE sequence achieved eddy current compensation with a significantly shorter TE than TRSE (78 versus 96 ms) and a slightly shorter TE than MONO (78 versus 80 ms). Intravoxel signal variance was lower in phantoms with EN-CODE than with MONO (13.6 ± 11.6 versus 37.4 ± 25.8) and not different from TRSE (15.1 ± 11.6), indicating good robustness to eddy current–induced image distortions. Mean fractional anisotropy values in brain edges were also significantly lower with EN-CODE than with MONO (0.16 ± 0.01 versus 0.24 ± 0.02, P

Purpose Design, implement, integrate, and characterize a customized coil system that allows for generating spatial encoding magnetic fields (SEMs) in a highly-flexible fashion. Methods A gradient coil with a high number of individual elements was designed. Dimensions of the coil were chosen to mimic a whole-body gradient system, scaled down to a head insert. Mechanical shape and wire layout of each element were optimized to increase the local gradient strength while minimizing eddy current effects and simultaneously considering manufacturing constraints. Results Resulting wire layout and mechanical design is presented. A prototype matrix gradient coil with 12 × 7 = 84 elements consisting of two element types was realized and characterized. Measured eddy currents are

Purpose To develop single and double acquisition methods to compensate for artifacts from eddy currents and transient oscillations in balanced steady-state free precession (bSSFP) with centric phase-encoding (PE) order for magnetization-prepared bSSFP imaging. Theory and Methods A single and four different double acquisition methods were developed and evaluated with Bloch equation simulations, phantom/in vivo experiments, and quantitative analyses. For the single acquisition method, multiple PE groups, each of which was composed of N linearly changing PE lines, were ordered in a pseudocentric manner for optimal contrast and minimal signal fluctuations. Double acquisition methods used complex averaging of two images that had opposite artifact patterns from different acquisition orders or from different numbers of dummy scans. Results Simulation results showed high sensitivity of eddy-current and transient-oscillation artifacts to off-resonance frequency and PE schemes. The artifacts were reduced with the PE-grouping with N values from 3 to 8, similar to or better than the conventional pairing scheme of N = 2. The proposed double acquisition methods removed the remaining artifacts significantly. The proposed methods conserved detailed structures in magnetization transfer imaging well, compared with the conventional methods. Conclusion The proposed single and double acquisition methods can be useful for artifact-free magnetization-prepared bSSFP imaging with desired contrast and minimized dummy scans. Magn Reson Med 78:254–263, 2017. © 2016 International Society for Magnetic Resonance in Medicine

Purpose To develop a gradient pre-emphasis scheme that prospectively counteracts the effects of the first-order concomitant fields for any arbitrary gradient waveform played on asymmetric gradient systems, and to demonstrate the effectiveness of this approach using a real-time implementation on a compact gradient system. Methods After reviewing the first-order concomitant fields that are present on asymmetric gradients, we developed a generalized gradient pre-emphasis model assuming arbitrary gradient waveforms to counteract their effects. A numerically straightforward, easily implemented approximate solution to this pre-emphasis problem was derived that was compatible with the current hardware infrastructure of conventional MRI scanners for eddy current compensation. The proposed method was implemented on the gradient driver subsystem, and its real-time use was tested using a series of phantom and in vivo data acquired from two-dimensional Cartesian phase-difference, echo-planar imaging, and spiral acquisitions. Results The phantom and in vivo results demonstrated that unless accounted for, first-order concomitant fields introduce considerable phase estimation error into the measured data and result in images with spatially dependent blurring/distortion. The resulting artifacts were effectively prevented using the proposed gradient pre-emphasis. Conclusion We have developed an efficient and effective gradient pre-emphasis framework to counteract the effects of first-order concomitant fields of asymmetric gradient systems. Magn Reson Med 77:2250–2262, 2017. © 2016 International Society for Magnetic Resonance in Medicine

Purpose To develop a new Nyquist ghost and geometric distortion correction method in echo planar imaging (EPI) using parallel imaging. Methods Two frames of EPI data are acquired with normal and phase-labeled sequence. The phase label is applied by modifying the PE prephase gradient to shift the central echo by one echo spacing. GRAPPA weights are trained from both frames and used to reconstruct images from positive or negative echoes in each frame to remove Nyquist ghost. Geometric distortion is then corrected by the B0 field map generated from the phase difference between positive and negative echo images. Phantom and in vivo experiments at 7 Tesla (T) and 3T were performed to evaluate the proposed method. Results Nyquist ghost was greatly reduced in all images even under oblique imaging and poor eddy current conditions, yielding significant improvements over the existing reference scan and image entropy minimization based methods. Image geometries were fully restored after distortion correction. Phantom results indicated that the signal-to-noise ratio efficiency was largely preserved while fMRI results showed no apparent degradation of temporal resolution. Conclusion The proposed method provides robust correction of both Nyquist ghost and geometric distortion in EPI, and it is particularly suitable for dynamic EPI applications. Magn Reson Med 77:1749-1761, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Purpose To shorten 4D flow acquisitions by shortening TRs with fast RF pulses and gradient waveforms. Real-time convex optimization is used to generate these gradients waveforms on the scanner. Theory and methods RF and slab-select waveforms were shortened with a minimum phase SLR excitation and the time-optimal variable-rate selective excitation method. Real-time convex optimization was used to shorten bipolar and spoiler gradients by finding the shortest gradient waveforms that satisfied constraints on scan parameters, gradient hardware, M0 , M1 , and peripheral nerve stimulation. Waveforms were calculated and TE and/or TR values were compared for a range of scan parameters and compared to a conventional 4D flow sequence. The method was tested in flow phantoms, and in the aorta and neurovasculature of volunteers (N = 10). Additionally, eddy current error was measured in a large phantom. Results TEs and TRs were shortened by 21-32% and 20-34%, respectively, compared to the conventional sequence over a range of scan parameters. Bland-Altman analysis of 2 flow phantom configurations showed flow rate bias of 0.3 mL/s and limits of agreement (LOA) of [-6.9, 7.5] mL/s for a cardiac phantom and a bias of -0.1 mL/s with LOA = [-0.4, 0.2] mL/s for a neuro phantom. Similar agreement was also seen for flow measurements in volunteers (bias = -1.0 and -0.1 mL/s, LOA = [-34.9, 33.0] and [-0.7, 0.6] mL/s). Measured eddy currents were 39% larger with the CVX + mpVERSE method. Conclusion The real-time optimized 4D flow gradients and fast slab-selection excitation methods produced up to 34% faster TRs with excellent flow measurement agreement compared to a conventional 4D flow sequence.

Purpose To develop a Fourier-transform based velocity-selective inversion (FT-VSI) pulse train for velocity-selective arterial spin labeling (VSASL). Methods This new pulse contains paired and phase cycled refocusing pulses. Its sensitivities to B0/B1 inhomogeneity and gradient imperfections such as eddy currents were evaluated through simulation and phantom studies. Cerebral blood flow (CBF) quantification using FT-VSI prepared VSASL was compared with conventional VSASL and pseudocontinuous ASL (PCASL) at 3 Tesla. Results Simulation and phantom results of the proposed FT-VSI pulse train demonstrated excellent robustness to B0/B1 field inhomogeneity and eddy currents. The estimated CBF of gray matter and white matter for the FT-VSI prepared VSASL, averaged among eight healthy volunteers, were 49.5 ± 7.5 mL/100 g/min and 14.8 ± 2.4 mL/100 g/min, respectively. Excellent correlation and agreement between the FT-VSI method and conventional VSASL and PCASL were found. The averaged signal-to-noise ratio (SNR) value in gray matter of the FT-VSI method was 39% higher than VSASL using conventional double refocused hyperbolic tangent pulses and 9% lower than PCASL. Conclusion A novel FT-VSI pulse train was demonstrated to be a suitable labeling module for VSASL with robustness of velocity-selective profile to B0/B1 field inhomogeneity and gradient imperfections. Compared with conventional VSASL, FT-VSI prepared VSASL produced consistent CBF maps with higher SNR values. Magn Reson Med 76:1136-1148, 2016. © 2015 Wiley Periodicals, Inc.

PURPOSE: Knowledge on low frequency (LF) tissue conductivity is relevant for various biomedical purposes. To obtain this information, LF phase maps arising from time-varying imaging gradients have been demonstrated to create a LF conductivity contrast. Essential in this methodology is the subtraction of phase images acquired with opposite gradient polarities to separate LF and RF phase effects. Here we demonstrate how sensitive these subtractions are with respect to geometrical distortions. THEORY AND METHODS: The effect of geometrical distortions on LF phase maps is mathematically defined. After quantifying typical geometrical distortions, their effects on LF phase maps are evaluated using conductive phantoms. For validation, electromagnetic simulations of LF phase maps were performed. RESULTS: Even sub-voxel distortions of 10% of the voxel size, measured for a typical LF MR sequence, cause leakage of RF phase into LF phase of several milli-radians, leading to a misleading pattern of LF phase maps. This leakage is mathematically confirmed, while simulations indicate that the expected LF phase should be in order of micro-radians. CONCLUSION: The conductivity scaling of LF phase maps is attributable to the RF phase leakage, thus dependent on the RF conductivity. In fact, simulations show that the LF phase is not measurable. Magn Reson Med 76:905-912, 2016. © 2015 Wiley Periodicals, Inc.

Purpose Tailored excitation (TEx) based on interspersing multiple radio frequency pulses with linear gradient and higher-order shim pulses can be used to obtain uniform flip angle in the presence of large radio frequency transmission (B ) inhomogeneity. Here, an implementation of dynamic, multislice tailored excitation using the recently developed multi-coil nonlinear shim hardware (MC-DTEx) is reported. Methods MC-DTEx was developed and tested both in a phantom and in vivo at 7 T, and its efficacy was quantitatively assessed. Predicted outcomes of MC-DTEx and DTEx based on spherical harmonic shims (SH-DTEx) were also compared. Results For a planned 30 flip angle, in a phantom, the standard deviation in excitation improved from 28% (regular excitation) to 12% with MC-DTEx. The SD in in vivo excitation improved from 22 to 12%. The improvements achieved with experimental MC-DTEx closely matched the theoretical predictions. Simulations further showed that MC-DTEx outperforms SH-DTEx for both scenarios. Conclusion Successful implementation of multislice MC-DTEx is presented and is shown to be capable of homogenizing excitation over more than twofold B variations. Its benefits over SH-DTEx are also demonstrated. A distinct advantage of MC hardware over SH shim hardware is the absence of significant eddy current effects, which allows for a straightforward, multislice implementation of MC-DTEx. Magn Reson Med, 2015. © 2015 Wiley Periodicals, Inc.

Purpose Diffusion MRI is compromised by unknown field perturbation during image encoding. The purpose of this study was to address this problem using the recently described approach of concurrent magnetic field monitoring. Methods Magnetic field dynamics were monitored during the echo planar imaging readout of a common diffusion-weighted MRI sequence using an integrated magnetic field camera setup. The image encoding including encoding changes over the duration of entire scans were quantified and analyzed. Field perturbations were corrected by accounting for them in generalized image reconstruction. The impact on image quality along with geometrical congruence among different diffusion-weighted images was assessed both qualitatively and quantitatively. Results The most significant field perturbations were found to be related to higher-order eddy currents from diffusion-weighting gradients and B0 field drift as well as gradual changes of short-term eddy current behavior and mechanical oscillations during the scan. All artifacts relating to dynamic field perturbations were eliminated by incorporating the measured encoding in image reconstruction. Conclusion Concurrent field monitoring combined with generalized reconstruction enhances depiction fidelity in diffusion imaging. In addition to artifact reduction, it improves geometric congruence and thus facilitates image combination for quantitative diffusion analysis. Magn Reson Med 74:925–933, 2015. © 2015 Wiley Periodicals, Inc.

Purpose The combination of fully balanced SSFP sequences with iterative golden angle radial sparse parallel (iGRASP) MRI leads to strong image artifacts due to eddy currents caused by the large angular increment of the golden angle ordering. The purpose of this work is to enable the combination of iterative golden angle radial sparse parallel MRI with balanced SSFP using the recently presented tiny golden angles. Methods The tiny golden angle trajectories are analyzed for their incoherence properties in relation to sparse imaging using the time-resolved point-spread functions. Tiny golden angle radial sparse parallel (tyGRASP) MRI is introduced and evaluated with applications in cardiac imaging and dynamic imaging of the temporomandibular joint. The results are analyzed in detail for 3 T and verified for 1.5 T. Results The incoherence properties of the tiny golden angle trajectory are comparable to the incoherence properties of the golden angle trajectory and are well suited for sparse MRI reconstruction. The proposed tiny golden angle radial sparse parallel MRI method strongly reduces eddy current related artifacts for both applications. Conclusion This work enables sparse, golden-ratio-based imaging with balanced SSFP sequences. Magn Reson Med, 2015. © 2015 Wiley Periodicals, Inc.

Purpose To accurately correct diffusion-encoding direction-dependent eddy-current-induced geometric distortions in diffusion-weighted echo-planar imaging (DW-EPI) and to minimize the calibration time at 7 Tesla (T). Methods A point spread function (PSF) mapping based eddy-current calibration method is newly presented to determine eddy-current-induced geometric distortions even including nonlinear eddy-current effects within the readout acquisition window. To evaluate the temporal stability of eddy-current maps, calibration was performed four times within 3 months. Furthermore, spatial variations of measured eddy-current maps versus their linear superposition were investigated to enable correction in DW-EPIs with arbitrary diffusion directions without direct calibration. For comparison, an image-based eddy-current correction method was additionally applied. Finally, this method was combined with a PSF-based susceptibility-induced distortion correction approach proposed previously to correct both susceptibility and eddy-current-induced distortions in DW-EPIs. Results Very fast eddy-current calibration in a three-dimensional volume is possible with the proposed method. The measured eddy-current maps are very stable over time and very similar maps can be obtained by linear superposition of principal-axes eddy-current maps. High resolution in vivo brain results demonstrate that the proposed method allows more efficient eddy-current correction than the image-based method. Conclusion The combination of both PSF-based approaches allows distortion-free images, which permit reliable analysis in diffusion tensor imaging applications at 7T. Magn Reson Med, 2015. © 2015 Wiley Periodicals, Inc.

Purpose Lipid signals measured by 1H MR spectroscopy cannot be adequately quantified by common fitting routines like VARPRO or AMARES, if lipid spectra are distorted by irregular spatial and temporal inhomogeneities of the static magnetic field during readout. A fully automatic reference deconvolution algorithm is presented that eliminates these distortions before application of fitting routines. Methods The measured signal of the dominant methyl resonance is isolated with aid of a spectral estimator (estimation of parameters via rotational invariance techniques) and used as reference signal for estimation of distortions. A Wiener filter is applied to deconvolve those distortions in the lipid spectrum. Performance of the algorithm is assessed for different bandwidths and shapes of distortions, using artificially distorted as well as measured data. Results Application of the fully automatic reference deconvolution algorithm on simulated spectra yields a distinct increase in quantification accuracy. Deconvolved in vivo spectra of subcutaneous fat indicate reduced spectral overlap after application of the proposed strategy. Conclusion The proposed method is helpful for in vivo magnetic resonance spectroscopy of adipose tissue to correct for effects of field inhomogeneities within the voxel and for inevitable eddy current effects. Quantification accuracy is improved by eliminating distortions before application of fitting routines. Magn Reson Med 75:1391–1401, 2016. © 2015 Wiley Periodicals, Inc.

Purpose A preliminary assessment of the MRI-compatibility of metallic object possibly embedded within the patient is required before conducting the MRI examination. The Magnetic Iron Detector (MID) is a highly sensitive susceptometer that uses a weak magnetic field to measure iron overload in the liver. MID might be used to perform a screening procedure for MRI by determining the ferromagnetic/conductive properties of embedded metallic objects. Methods The study was composed by: (i) definition of MID sensitivity threshold; (ii) application of MID in a procedure to characterize the ferromagnetic/conductive properties of metallic foreign objects in 958 patients scheduled for MID examination. Results The detection threshold for ferromagnetic objects was found to be the equivalent of a piece of wire of length 2 mm and gauge 0.8 mm2 and, representing purely conductive objects, an aluminum sheet of area 2 × 2 cm2. Of 958 patients, 165 had foreign bodies of unknown nature. MID was able to detect those with ferromagnetic and/or conducting properties based on fluctuations in the magnetic and eddy current signals versus control. Conclusion The high sensitivity of MID makes it suitable for assessing the ferromagnetic/conductive properties of metallic foreign objects embedded within the body of patients scheduled for MRI. Magn Reson Med 73:2030–2037, 2015. © 2014 Wiley Periodicals, Inc.

Purpose We investigated the temporal variation of the induced magnetic field due to the transverse and the longitudinal gradient coils in tungsten collimators arranged in hexagonal and pentagonal geometries with and without gaps between the collimators. Methods We modeled x-, y-, and z-gradient coils and different arrangements of single-photon emission computed tomography (SPECT) collimators using FEKO, a three-dimensional electromagnetic simulation tool. A time analysis approach was used to generate the pulsed magnetic field gradient. The approach was validated with measurements using a 7T MRI scanner. Results Simulations showed an induced magnetic field representing 4.66% and 0.87% of the applied gradient field (gradient strength = 500 mT/m) for longitudinal and transverse gradient coils, respectively. These values can be reduced by 75% by adding gaps between the collimators for the pentagonal arrangement, bringing the maximum induced magnetic field to less than 2% of the applied gradient for all of the gradient coils. Conclusion Characterization of the maximum induced magnetic field shows that by adding gaps between the collimators for an integrated SPECT/MRI system, eddy currents can be corrected by the MRI system to avoid artifact. The numerical model was validated and was proposed as a tool for studying the effect of a SPECT collimator within the MRI gradient coils. Magn Reson Med, 2014. © 2014 Wiley Periodicals, Inc.

Purpose To optimize the design of radiofrequency (RF) shielding of transmit coils at 7T and reduce eddy currents generated on the RF shielding when imaging with rapid gradient waveforms. Methods One set of a four-element, 2 × 2 Tic-Tac-Toe head coil structure was selected and constructed to study eddy currents on the RF coil shielding. The generated eddy currents were quantitatively studied in the time and frequency domains. The RF characteristics were studied using the finite difference time domain method. Five different kinds of RF shielding were tested on a 7T MRI scanner with phantoms and in vivo human subjects. Results The eddy current simulation method was verified by the measurement results. Eddy currents induced by solid/intact and simple-structured slotted RF shielding significantly distorted the gradient fields. Echo-planar images, maps, and S matrix measurements verified that the proposed slot pattern suppressed the eddy currents while maintaining the RF characteristics of the transmit coil. Conclusion The presented dual-optimization method could be used to design RF shielding and reduce the gradient field-induced eddy currents while maintaining the RF characteristics of the transmit coil. Magn Reson Med 74:1461–1469, 2015. © 2014 Wiley Periodicals, Inc.

Purpose Velocity-selective arterial spin labeling (VSASL) tags arterial blood on a velocity-selective (VS) basis and eliminates the tagging/imaging gap and associated transit delay sensitivity observed in other ASL tagging methods. However, the flow-weighting gradient pulses in VS tag preparation can generate eddy currents (ECs), which may erroneously tag the static tissue and create artificial perfusion signal, compromising the accuracy of perfusion quantification. Methods A novel VS preparation design is presented using an eight-segment B1 insensitive rotation with symmetric radio frequency and gradient layouts (sym-BIR-8), combined with delays after gradient pulses to optimally reduce ECs of a wide range of time constants while maintaining B0 and B1 insensitivity. Bloch simulation, phantom, and in vivo experiments were carried out to determine robustness of the new and existing pulse designs to ECs, B0, and B1 inhomogeneity. Results VSASL with reduced EC sensitivity across a wide range of EC time constants was achieved with the proposed sym-BIR-8 design, and the accuracy of cerebral blood flow measurement was improved. Conclusion The sym-BIR-8 design performed the most robustly among the existing VS tagging designs, and should benefit studies using VS preparation with improved accuracy and reliability. Magn Reson Med 73:1085–1094, 2015. © 2014 Wiley Periodicals, Inc.

Purpose Phase-contrast MRI (PC-MRI) is a versatile tool allowing evaluation of in vivo motion, but is sensitive to eddy current induced phase offsets, causing errors in the measured velocities. In high-resolution PC-MRI, these offsets can be sufficiently large to cause wrapping in the baseline phase, rendering conventional eddy current compensation (ECC) inadequate. The purpose of this study was to develop an improved ECC technique (unwrapping ECC) able to handle baseline phase discontinuities. Theory and Methods Baseline phase discontinuities are unwrapped by minimizing the spatiotemporal standard deviation of the static-tissue phase. Computer simulations were used for demonstrating the theoretical foundation of the proposed technique. The presence of baseline wrapping was confirmed in high-resolution myocardial PC-MRI of a normal rat heart at 9.4 Tesla (T), and the performance of unwrapping ECC was compared with conventional ECC. Results Areas of phase wrapping in static regions were clearly evident in high-resolution PC-MRI. The proposed technique successfully eliminated discontinuities in the baseline, and resulted in significantly better ECC than the conventional approach. Conclusion We report the occurrence of baseline phase wrapping in PC-MRI, and provide an improved ECC technique capable of handling its presence. Unwrapping ECC offers improved correction of eddy current induced baseline shifts in high-resolution PC-MRI. Magn Reson Med 72:1096–1102, 2014. © 2013 Wiley Periodicals, Inc.

Purpose Delta relaxation enhanced magnetic resonance (dreMR) imaging requires an auxiliary B0 electromagnet capable of shifting the main magnetic field within a clinical 1.5 Tesla (T) MR system. In this work, the main causes of interaction between an actively shielded, insertable resistive B0 electromagnet and a 1.5T superconducting system are systematically identified and mitigated. Methods The effects of nonideal fabrication of the field-shifting magnet are taken into consideration through careful measurement during winding and improved accuracy in the design of the associated active shield. The shielding performance of the resultant electromagnet is compared against a previously built system in which the shield design was based on an ideal primary coil model. Hardware and software approaches implemented to eliminate residual image artifacts are presented in detail. Results The eddy currents produced by the newly constructed dreMR system are shown to have a significantly smaller “long-time-constant” component, consistent with the hypothesis that less energy is deposited into the cryostat of the MR system. Conclusion With active compensation, the dreMR imaging system is capable of 0.22T field shifts within a clinical 1.5T MRI with no significant residual eddy-current fields. Magn Reson Med 72:1182–1190, 2014. © 2013 Wiley Periodicals, Inc.

Purpose Applications of dynamic shimming require high field fidelity, and characterizing the shim field dynamics is therefore necessary. Modeling the system as linear and time-invariant, the purpose of this work was to measure the impulse response function with optimal sensitivity. Theory and Methods Frequency-swept pulses as inputs are analyzed theoretically, showing that the sweep speed is a key factor for the measurement sensitivity. By adjusting the sweep speed it is possible to achieve any prescribed noise profile in the measured system response. Impulse response functions were obtained for the third-order shim system of a 7 Tesla whole-body MR scanner. Measurements of the shim fields were done with a dynamic field camera, yielding also cross-term responses. Results The measured shim impulse response functions revealed system characteristics such as response bandwidth, eddy currents and specific resonances, possibly of mechanical origin. Field predictions based on the shim characterization were shown to agree well with directly measured fields, also in the cross-terms. Conclusion Frequency sweeps provide a flexible tool for shim or gradient system characterization. This may prove useful for applications involving dynamic shimming by yielding accurate estimates of the shim fields and a basis for setting shim pre-emphasis. Magn Reson Med 72:570–583, 2014. © 2013 Wiley Periodicals, Inc.

Purpose The measurement of velocities based on phase contrast MRI can be subject to different phase offset errors which can affect the accuracy of velocity data. The purpose of this study was to determine the impact of these inaccuracies and to evaluate different correction strategies on three-dimensional visualization. Methods Phase contrast MRI was performed on a 3 T system (Siemens Trio) for in vitro (curved/straight tube models; venc: 0.3 m/s) and in vivo (aorta/intracranial vasculature; venc: 1.5/0.4 m/s) data. For comparison of the impact of different magnetic field gradient designs, in vitro data was additionally acquired on a wide bore 1.5 T system (Siemens Espree). Different correction methods were applied to correct for eddy currents, Maxwell terms, and gradient field inhomogeneities. Results The application of phase offset correction methods lead to an improvement of three-dimensional particle trace visualization and count. The most pronounced differences were found for in vivo/in vitro data (68%/82% more particle traces) acquired with a low venc (0.3 m/s/0.4 m/s, respectively). In vivo data acquired with high venc (1.5 m/s) showed noticeable but only minor improvement. Conclusion This study suggests that the correction of phase offset errors can be important for a more reliable visualization of particle traces but is strongly dependent on the velocity sensitivity, object geometry, and gradient coil design. Magn Reson Med 72:33–40, 2014. © 2013 Wiley Periodicals, Inc.

Purpose This article aims to present a fast, efficient and accurate multi-layer integral method (MIM) for the evaluation of complex spatiotemporal eddy currents in nonmagnetic and thin volumes of irregular geometries induced by arbitrary arrangements of gradient coils. Methods The volume of interest is divided into a number of layers, wherein the thickness of each layer is assumed to be smaller than the skin depth and where one of the linear dimensions is much smaller than the remaining two dimensions. The diffusion equation of the current density is solved both in time-harmonic and transient domain. Results The experimentally measured magnetic fields produced by the coil and the induced eddy currents as well as the corresponding time-decay constants were in close agreement with the results produced by the MIM. Relevant parameters such as power loss and force induced by the eddy currents in a split cryostat were simulated using the MIM. Conclusion The proposed method is capable of accurately simulating the current diffusion process inside thin volumes, such as the magnet cryostat. The method permits the priori-calculation of optimal pre-emphasis parameters. The MIM enables unified designs of gradient coil - magnet structures for an optimal mitigation of deleterious eddy current effects.

Purpose Time-varying system delays and eddy currents can substantially reduce the image quality of spiral images. A new method is proposed to estimate variable system delays for spiral-based trajectories. Methods This approach requires a minor modification of a conventional stack-of-spirals sequence and analyzes data collected on three overlapping orthogonal cylinders. Results Initial results are presented for acquired and synthesized phantom data and in vivo data. Conclusion The proposed method includes gradient coupling effects, requires no phantom measurements or specialized hardware, is robust to off-resonance effects, and estimates independent continuous delays for each gradient channel over the data-sampling period. Magn Reson Med 71:773–782, 2014. © 2013 Wiley Periodicals, Inc.

Purpose To correct gradient-induced phase errors in radial MRI. Methods Gradient-induced eddy currents affect the MRI data acquisition by gradient delays and phase errors that may lead to severe image artifacts for non-Cartesian imaging scenarios such as radial trajectories. While gradient delays are dealt with by respective shifts of the acquisition window during radial image acquisition, this work introduces a simple method for quantifying and correcting phase errors from the actual data prior to image reconstruction. For a given gradient system, the approach yields a specific phase error per gradient that can be used for correcting the raw data. Results Phantom studies at 9.4 T demonstrated marked improvements in radial image quality. It could be shown that the phase correction is not compromised by data undersampling. Moreover, the selective correction of gradient-induced phase errors retained the phase information caused by different concentrations of a paramagnetic contrast agent. Conclusion The proposed method does not require additional reference measurements and separately corrects for phase errors induced by eddy currents, while retaining the residual phase of the object which may carry physiologic information. Magn Reson Med 71:308–312, 2014. © 2013 Wiley Periodicals, Inc.

Purpose The twice-refocused spin echo sequence is widely used in diffusion imaging due to its excellent performance in reducing eddy currents. The three radio frequency pulses give rise to eight separate signal pathways. Because there is no general solution for the size and arrangement for crusher gradients, with constant size and orientation, that is effective for all arbitrary diffusion-sensitizing b-values and directions, this article introduces and validates a solution whereby the crusher and diffusion-encoding gradients are always kept orthogonal, thus ensuring their independence. Methods The cancellation of the crusher and diffusion gradients was demonstrated. Subsequently, crusher gradients were implemented in such a way that they were always orthogonal to the diffusion gradient. Phantom and in-vivo experiments were performed to ascertain that orthogonally implemented crusher gradients alleviate the problem without lowering image quality. Results In all experiments, when the crusher gradients' action was cancelled by the diffusion-encoding gradients artifactual signal modulation was observed. When orthogonal gradients were implemented the artifacts were eliminated without detrimental effects on image quality. Conclusions Orthogonal crushers are easy to implement and can be used for any variant of diffusion imaging sequences (e.g., diffusion tensor imaging, fiber diameter mapping) where the twice-refocused scheme is used. Magn Reson Med 71:506–515, 2014. © 2013 Wiley Periodicals, Inc.