Advances in spiral fMRI: A high-resolution study with single-shot acquisition.

Spiral fMRI has been put forward as a viable alternative to rectilinear echo-planar imaging, in particular due to its enhanced average k-space speed and thus high acquisition efficiency. This renders spirals attractive for contemporary fMRI applications that require high spatiotemporal resolution, such as laminar or columnar fMRI. However, in practice, spiral fMRI is typically hampered by its reduced robustness and ensuing blurring artifacts, which arise from imperfections in both static and dynamic magnetic fields. Recently, these limitations have been overcome by the concerted application of an expanded signal model that accounts for such field imperfections, and its inversion by iterative image reconstruction. In the challenging ultra-high field environment of 7 Tesla, where field inhomogeneity effects are aggravated, both multi-shot and single-shot 2D spiral imaging at sub-millimeter resolution was demonstrated with high depiction quality and anatomical congruency. In this work, we further these advances towards a time series application of spiral readouts, namely, single-shot spiral BOLD fMRI at 0.8 mm in-plane resolution. We demonstrate that high-resolution spiral fMRI at 7 T is not only feasible, but delivers both excellent image quality, BOLD sensitivity, and spatial specificity of the activation maps, with little artifactual blurring. Furthermore, we show the versatility of the approach with a combined in/out spiral readout at a more typical resolution (1.5 mm), where the high acquisition efficiency allows to acquire two images per shot for improved sensitivity by echo combination.

Rapid anatomical brain imaging using spiral acquisition and an expanded signal model.

We report the deployment of spiral acquisition for high-resolution structural imaging at 7T. Long spiral readouts are rendered manageable by an expanded signal model including static off-resonance and B0 dynamics along with k-space trajectories and coil sensitivity maps. Image reconstruction is accomplished by inversion of the signal model using an extension of the iterative non-Cartesian SENSE algorithm. Spiral readouts up to 25 ms are shown to permit whole-brain 2D imaging at 0.5 mm in-plane resolution in less than a minute. A range of options is explored, including proton-density and T2* contrast, acceleration by parallel imaging, different readout orientations, and the extraction of phase images. Results are shown to exhibit competitive image quality along with high geometric consistency.

Single-shot spiral imaging at 7 T.

PURPOSE: The purpose of this work is to explore the feasibility and performance of single-shot spiral MRI at 7 T, using an expanded signal model for reconstruction. METHODS: Gradient-echo brain imaging is performed on a 7 T system using high-resolution single-shot spiral readouts and half-shot spirals that perform dual-image acquisition after a single excitation. Image reconstruction is based on an expanded signal model including the encoding effects of coil sensitivity, static off-resonance, and magnetic field dynamics. The latter are recorded concurrently with image acquisition, using NMR field probes. The resulting image resolution is assessed by point spread function analysis. RESULTS: Single-shot spiral imaging is achieved at a nominal resolution of 0.8 mm, using spiral-out readouts of 53-ms duration. High depiction fidelity is achieved without conspicuous blurring or distortion. Effective resolutions are assessed as 0.8, 0.94, and 0.98 mm in CSF, gray matter and white matter, respectively. High image quality is also achieved with half-shot acquisition yielding image pairs at 1.5-mm resolution. CONCLUSION: Use of an expanded signal model enables single-shot spiral imaging at 7 T with unprecedented image quality. Single-shot and half-shot spiral readouts deploy the sensitivity benefit of high field for rapid high-resolution imaging, particularly for functional MRI and arterial spin labeling.

Ultrafast Ligand Self-Exchanging Gadolinium Complexes in Ionic Liquids for NMR Field Probes.

Concurrent magnetic field monitoring in MRI with an array of NMR field probes allows for reducing image imperfections. High 19F concentrations together with short relaxation times are basic probe properties. We present the NMR properties of [Gd(NTf2)4]- dissolved in an ionic liquid consisting of an imidazolium cation and the [NTf2]- anion. These solutions achieve fluorine concentrations as high as 26 M and rapid relaxation times in the sub-ms time range. The second order self-exchange rate of coordinated [NTf2]- with [NTf2]- as determined with [Y(NTf2)4]- is 4.5 × 105·M-1·s-1. X-ray structure analyses confirm the coordination number of eight around Gd3+ and Y3+. These ionic-liquid solutions are thus excellent candidates for dynamically probing the NMR fields in MRI.

Gradient and shim pre-emphasis by inversion of a linear time-invariant system model.

PURPOSE: The goal of this contribution is to enhance the fidelity and switching speed of gradient and shim fields by advancing pre-emphasis toward broadband and full cross-term correction. THEORY AND METHODS: The proposed approach is based on viewing gradient and shim chains as linear, time-invariant (LTI) systems. Pre-emphasis is accomplished by inversion of a broadband digital system model. In the multiple-channel case, it amounts to a matrix of broadband filters that perform concerted self- and cross-term correction. This approach is demonstrated with gradients and shims up to the third order in a 7 Tesla whole-body MR system. RESULTS: Pre-emphasis by LTI model inversion is first verified by studying settling speeds and response behavior without and with the correction. It is then demonstrated for rapid shim updating, achieving substantially enhanced fidelity of field dynamics and shim settling within approximately 1 ms. In single-shot echo-planar imaging (EPI) acquisitions in vivo, this benefit is shown to translate into enhanced geometric fidelity. CONCLUSIONS: The fidelity of gradient and shim dynamics can be greatly enhanced by pre-emphasis based on inverting a general LTI system model. Permitting shim settling on the millisecond scale, broadband multiple-channel pre-emphasis promises to render higher-order shimming fully versatile at the level of MRI sequence design. Magn Reson Med 78:1607-1622, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Single-shot spiral imaging enabled by an expanded encoding model: Demonstration in diffusion MRI.

PURPOSE: The purpose of this work was to improve the quality of single-shot spiral MRI and demonstrate its application for diffusion-weighted imaging. METHODS: Image formation is based on an expanded encoding model that accounts for dynamic magnetic fields up to third order in space, nonuniform static B0 , and coil sensitivity encoding. The encoding model is determined by B0 mapping, sensitivity mapping, and concurrent field monitoring. Reconstruction is performed by iterative inversion of the expanded signal equations. Diffusion-tensor imaging with single-shot spiral readouts is performed in a phantom and in vivo, using a clinical 3T instrument. Image quality is assessed in terms of artefact levels, image congruence, and the influence of the different encoding factors. RESULTS: Using the full encoding model, diffusion-weighted single-shot spiral imaging of high quality is accomplished both in vitro and in vivo. Accounting for actual field dynamics, including higher orders, is found to be critical to suppress blurring, aliasing, and distortion. Enhanced image congruence permitted data fusion and diffusion tensor analysis without coregistration. CONCLUSION: Use of an expanded signal model largely overcomes the traditional vulnerability of spiral imaging with long readouts. It renders single-shot spirals competitive with echo-planar readouts and thus deploys shorter echo times and superior readout efficiency for diffusion imaging and further prospective applications. Magn Reson Med 77:83-91, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Image reconstruction using a gradient impulse response model for trajectory prediction.

PURPOSE: Gradient imperfections remain a challenge in MRI, especially for sequences relying on long imaging readouts. This work aims to explore image reconstruction based on k-space trajectories predicted by an impulse response model of the gradient system. THEORY AND METHODS: Gradient characterization was performed twice with 3 years interval on a commercial 3 Tesla (T) system. The measured gradient impulse response functions were used to predict actual k-space trajectories for single-shot echo-planar imaging (EPI), spiral and variable-speed EPI sequences. Image reconstruction based on the predicted trajectories was performed for phantom and in vivo data. Resulting images were compared with reconstructions based on concurrent field monitoring, separate trajectory measurements, and nominal trajectories. RESULTS: Image reconstruction using model-based trajectories yielded high-quality images, comparable to using separate trajectory measurements. Compared with using nominal trajectories, it strongly reduced ghosting, blurring, and geometric distortion. Equivalent image quality was obtained with the recent characterization and that performed 3 years prior. CONCLUSION: Model-based trajectory prediction enables high-quality image reconstruction for technically challenging sequences such as single-shot EPI and spiral imaging. It thus holds great promise for fast structural imaging and advanced neuroimaging techniques, including functional MRI, diffusion tensor imaging, and arterial spin labeling. The method can be based on a one-time system characterization as demonstrated by successful use of 3-year-old calibration data. Magn Reson Med 76:45-58, 2016. © 2015 Wiley Periodicals, Inc.

A field camera for MR sequence monitoring and system analysis.

PURPOSE: MR image formation and interpretation relies on highly accurate dynamic magnetic fields of high fidelity. A range of mechanisms still limit magnetic field fidelity, including magnet drifts, eddy currents, and finite linearity and stability of power amplifiers used to drive gradient and shim coils. Addressing remaining errors by means of hardware, sequence, or signal processing optimizations, calls for immediate observation by magnetic field monitoring. The present work presents a stand-alone monitoring system delivering insight into such field imperfections for MR sequence and system analysis. METHODS: A flexible NMR field probe-based stand-alone monitoring system, built on a software-defined-radio approach, is introduced and used to sense field dynamics up to third-order in space in a selection of situations with different time scales. RESULTS: Highly sensitive trajectories are measured and successfully used for image reconstruction. Further field perturbations due to mechanical oscillations and thermal field drifts following demanding gradient use and external interferences are studied. CONCLUSION: A flexible and versatile monitoring system is presented, delivering camera-like access to otherwise hardly accessible field dynamics with nanotesla resolution. Its stand-alone nature enables field analysis even during unknown MR system states.

Concurrent recording of RF pulses and gradient fields - comprehensive field monitoring for MRI.

Reconstruction of MRI data is based on exact knowledge of all magnetic field dynamics, since the interplay of RF and gradient pulses generates the signal, defines the contrast and forms the basis of resolution in spatial and spectral dimensions. Deviations caused by various sources, such as system imperfections, delays, eddy currents, drifts or externally induced fields, can therefore critically limit the accuracy of MRI examinations. This is true especially at ultra-high fields, because many error terms scale with the main field strength, and higher available SNR renders even smaller errors relevant. Higher baseline field also often requires higher acquisition bandwidths and faster signal encoding, increasing hardware demands and the severity of many types of hardware imperfection. To address field imperfections comprehensively, in this work we propose to expand the concept of magnetic field monitoring to also encompass the recording of RF fields. In this way, all dynamic magnetic fields relevant for spin evolution are covered, including low- to audio-frequency magnetic fields as produced by main magnets, gradients and shim systems, as well as RF pulses generated with single- and multiple-channel transmission systems. The proposed approach permits field measurements concurrently with actual MRI procedures on a strict common time base. The combined measurement is achieved with an array of miniaturized field probes that measure low- to audio-frequency fields via (19) F NMR and simultaneously pick up RF pulses in the MRI system's (1) H transmit band. Field recordings can form the basis of system calibration, retrospective correction of imaging data or closed-loop feedback correction, all of which hold potential to render MRI more robust and relax hardware requirements. The proposed approach is demonstrated for a range of imaging methods performed on a 7 T human MRI system, including accelerated multiple-channel RF pulses. Copyright © 2015 John Wiley & Sons, Ltd.

Real-time motion correction using gradient tones and head-mounted NMR field probes.

PURPOSE: Sinusoidal gradient oscillations in the kilohertz range are proposed for position tracking of NMR probes and prospective motion correction for arbitrary imaging sequences without any alteration of sequence timing. The method is combined with concurrent field monitoring to robustly perform image reconstruction in the presence of potential dynamic field deviations. METHODS: Benchmarking experiments were done to assess the accuracy and precision of the method and to compare it with theoretical predictions based on the field probe's time-dependent signal-to-noise ratio. An array of four field probes was used to perform real-time prospective motion correction in vivo. Images were reconstructed based on both predetermined and concurrently measured k-space trajectories. RESULTS: For observation windows of 4.8 ms, the precision of probe position determination was found to be 35 to 62 µm, and the maximal measurement error was 595 µm root-mean-square on a single axis. Sequence update per repetition time on this basis yielded images free of conspicuous artifacts despite substantial head motion. Predetermined and concurrently observed k-space trajectories yielded equivalent image quality. CONCLUSION: NMR field probes in conjunction with gradient tones permit the tracking and prospective correction of rigid-body motion. Relying on gradient oscillations in the kilohertz range, the method allows for concurrent motion detection and image encoding.

Single-shot imaging with higher-dimensional encoding using magnetic field monitoring and concomitant field correction.

PURPOSE: PatLoc (Parallel Imaging Technique using Localized Gradients) accelerates imaging and introduces a resolution variation across the field-of-view. Higher-dimensional encoding employs more spatial encoding magnetic fields (SEMs) than the corresponding image dimensionality requires, e.g. by applying two quadratic and two linear spatial encoding magnetic fields to reconstruct a 2D image. Images acquired with higher-dimensional single-shot trajectories can exhibit strong artifacts and geometric distortions. In this work, the source of these artifacts is analyzed and a reliable correction strategy is derived. METHODS: A dynamic field camera was built for encoding field calibration. Concomitant fields of linear and nonlinear spatial encoding magnetic fields were analyzed. A combined basis consisting of spherical harmonics and concomitant terms was proposed and used for encoding field calibration and image reconstruction. RESULTS: A good agreement between the analytical solution for the concomitant fields and the magnetic field simulations of the custom-built PatLoc SEM coil was observed. Substantial image quality improvements were obtained using a dynamic field camera for encoding field calibration combined with the proposed combined basis. CONCLUSION: The importance of trajectory calibration for single-shot higher-dimensional encoding is demonstrated using the combined basis including spherical harmonics and concomitant terms, which treats the concomitant fields as an integral part of the encoding.

Real-time feedback for spatiotemporal field stabilization in MR systems.

PURPOSE: MR imaging and spectroscopy require a highly stable, uniform background field. The field stability is typically limited by hardware imperfections, external perturbations, or field fluctuations of physiological origin. The purpose of the present work is to address these issues by introducing spatiotemporal field stabilization based on real-time sensing and feedback control. METHODS: An array of NMR field probes is used to sense the field evolution in a whole-body MR system concurrently with regular system operation. The field observations serve as inputs to a proportional-integral controller that governs correction currents in gradient and higher-order shim coils such as to keep the field stable in a volume of interest. RESULTS: The feedback system was successfully set up, currently reaching a minimum latency of 20 ms. Its utility is first demonstrated by countering thermal field drift during an EPI protocol. It is then used to address respiratory field fluctuations in a T2 *-weighted brain exam, resulting in substantially improved image quality. CONCLUSION: Feedback field control is an effective means of eliminating dynamic field distortions in MR systems. Third-order spatial control at an update time of 100 ms has proven sufficient to largely eliminate thermal and breathing effects in brain imaging at 7 Tesla.

Retrospective correction of physiological field fluctuations in high-field brain MRI using concurrent field monitoring.

PURPOSE: Magnetic field fluctuations caused by subject motion, such as breathing or limb motion, can degrade image quality in brain MRI, especially at high field strengths. The purpose of this study was to investigate the feasibility of retrospectively correcting for such physiological field perturbations based on concurrent field monitoring. THEORY AND METHODS: High-resolution T2*-weighted gradient-echo images of the brain were acquired at 7T with subjects performing different breathing and hand movement patterns. Field monitoring with a set of (19) F NMR probes distributed around the head was performed in two variants: concurrently with imaging or as a single field measurement per readout. The measured field fluctuations were then accounted for in the image reconstruction. RESULTS: Significant field fluctuations due to motion were observed in all subjects, resulting in severe artifacts in uncorrected images. The artifacts were largely removed by reconstruction based on field monitoring. Accounting for field perturbations up to the 1st spatial order was generally sufficient to recover good image quality. CONCLUSIONS: It has been demonstrated that artifacts due to physiologically induced dynamic field perturbations can be greatly reduced by retrospective image correction based on field monitoring. The necessity to perform such correction is greatest at high fields and for field-sensitive techniques such as T2*-weighted imaging.

Diffusion MRI with concurrent magnetic field monitoring.

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.

Matched-filter acquisition for BOLD fMRI.

We introduce matched-filter fMRI, which improves BOLD (blood oxygen level dependent) sensitivity by variable-density image acquisition tailored to subsequent image smoothing. Image smoothing is an established post-processing technique used in the vast majority of fMRI studies. Here we show that the signal-to-noise ratio of the resulting smoothed data can be substantially increased by acquisition weighting with a weighting function that matches the k-space filter imposed by the smoothing operation. We derive the theoretical SNR advantage of this strategy and propose a practical implementation of 2D echo-planar acquisition matched to common Gaussian smoothing. To reliably perform the involved variable-speed trajectories, concurrent magnetic field monitoring with NMR probes is used. Using this technique, phantom and in vivo measurements confirm reliable SNR improvement in the order of 30% in a "resting-state" condition and prove robust in different regimes of physiological noise. Furthermore, a preliminary task-based visual fMRI experiment equally suggests a consistent BOLD sensitivity increase in terms of statistical sensitivity (average t-value increase of about 35%). In summary, our study suggests that matched-filter acquisition is an effective means of improving BOLD SNR in studies that rely on image smoothing at the post-processing level.

Field camera measurements of gradient and shim impulse responses using frequency sweeps.

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.

Feedback field control improves linewidths in in vivo magnetic resonance spectroscopy.

PURPOSE: Magnetic resonance spectroscopy (MRS) experiments rely on a homogeneous and stable magnetic field within the sample. Field homogeneity is typically optimized by static B0 shimming while reproducible effects from dynamic field variation are commonly diminished by means of gradient system calibration as well as calibration based on non-water suppressed reference data. However, residual encoding deficiencies from incomplete calibration and nonreproducible field perturbations deteriorate the quality of the obtained data. To overcome this problem, we propose to adapt higher-order feedback field control based on NMR field probes for its application in MRS. METHODS: To allow for field measurements simultaneously with the spectroscopy readout, radiofrequency-shielded field probes were employed. The setup was evaluated in vitro and tested in vivo for single-voxel MRS at 7T to correct for field perturbations that occur due to subject breathing and limb motion. RESULTS: The in vitro experiments showed an effective field control during the MRS sequence. The resulting spectroscopy data were free of spurious signal and the achieved field stabilization improved the spectral resolution in vitro and in vivo. CONCLUSION: High-field MRS is limited by nonreproducible field perturbations for which spatiotemporal field feedback provides a solution without compromising sequence timing and efficiency.

Analysis of temperature dependence of background phase errors in phase-contrast cardiovascular magnetic resonance.

BACKGROUND: The accuracy of phase-contrast cardiovascular magnetic resonance (PC-CMR) can be compromised by background phase errors. It is the objective of the present work to provide an analysis of the temperature dependence of background phase errors in PC-CMR by means of gradient mount temperature sensing and magnetic field monitoring. METHODS: Background phase errors were measured for various temperatures of the gradient mount using magnetic field monitoring and validated in a static phantom. The effect of thermal changes during k-space acquisition was simulated and confirmed with measurements in a stationary phantom. RESULTS: The temperature of the gradient mount was found to increase by 20-30 K during PC-CMR measurements of 6-12 min duration. Associated changes in background phase errors of up to 11% or 0.35 radian were measured at 10 cm from the magnet's iso-center as a result of first order offsets. Zeroth order phase errors exhibited little thermal dependence. CONCLUSIONS: It is concluded that changes in gradient mount temperature significantly modify background phase errors during PC-CMR with high gradient duty cycle. Since temperature increases significantly during the first minutes of scanning the results presented are also of relevance for single-slice or multi-slice PC-CMR scans. The findings prompt for further studies to investigate advanced correction methods taking into account gradient temperature and/or the use of concurrent field-monitoring to map gradient-induced fields throughout the scan.

Monitoring and compensating phase imperfections in cine balanced steady-state free precession.

PURPOSE: To analyze and correct for eddy current-induced phase imperfections in cardiac cine balanced steady-state free precession (bSSFP) imaging. METHODS: Eddy current-induced phase offsets were measured for different phase-encoding schemes using a higher order dynamic field camera. Based on these measurements, offset phases were corrected for in postprocessing and by run-time phase compensation applying radiofrequency phase increments and additional compensatory gradient areas. The findings were validated using numerical simulations, phantom experiments, and in vivo cardiac scans. RESULTS: Depending on the phase-encoding scheme, significant eddy current-induced phase offsets were detected. Time-varying phase offsets were observed at subsequent excitations leading to steady-state distortions and hence to profile-dependent amplitude modulations in k-space. Taking into account measured k-space trajectories algebraic image reconstruction allowed compensating imperfect spatial encoding. Correction of amplitude modulations was successfully accomplished by run-time phase compensation. CONCLUSION: Using magnetic field monitoring, artifacts in cine balanced steady-state free precession caused by uncompensated eddy current fields can be significantly reduced.

Gradient system characterization by impulse response measurements with a dynamic field camera.

This work demonstrates a fast, sensitive method of characterizing the dynamic performance of MR gradient systems. The accuracy of gradient time-courses is often compromised by field imperfections of various causes, including eddy currents and mechanical oscillations. Characterizing these perturbations is instrumental for corrections by pre-emphasis or post hoc signal processing. Herein, a gradient chain is treated as a linear time-invariant system, whose impulse response function is determined by measuring field responses to known gradient inputs. Triangular inputs are used to probe the system and response measurements are performed with a dynamic field camera consisting of NMR probes. In experiments on a whole-body MR system, it is shown that the proposed method yields impulse response functions of high temporal and spectral resolution. Besides basic properties such as bandwidth and delay, it also captures subtle features such as mechanically induced field oscillations. For validation, measured response functions were used to predict gradient field evolutions, which was achieved with an error below 0.2%. The field camera used records responses of various spatial orders simultaneously, rendering the method suitable also for studying cross-responses and dynamic shim systems. It thus holds promise for a range of applications, including pre-emphasis optimization, quality assurance, and image reconstruction.