Thermal variation in gradient response: measurement and modeling.

PURPOSE: Many aspects and imperfections of gradient dynamics in MRI have been successfully captured by linear time-invariant (LTI) models. Changes in gradient behavior due to heating, however, violate time invariance. The goal of this work is to study such changes at the level of transfer functions and model them by thermal extension of the LTI framework. METHODS: To study the impact of gradient heating on transfer functions, a clinical MR system was heated using a range of high-amplitude DC and AC waveforms, each followed by measuring transfer functions in rapid succession while the system cooled down. Simultaneously, gradient temperature was monitored with an array of temperature sensors positioned according to initial infrared recordings of the gradient tube. The relation between temperatures and transfer functions is cast into local and global linear models. The models are analysed in terms of self-consistency, conditioning, and prediction performance. RESULTS: Pronounced thermal effects are observed in the time resolved transfer functions, largely attributable to in-coil eddy currents and mechanical resonances. Thermal modeling is found to capture these effects well. The keys to good model performance are well-placed temperature sensors and suitable training data. CONCLUSION: Heating changes gradient response, violating time invariance. The utility of LTI modeling can nevertheless be recovered by a linear thermal extension, relying on temperature sensing and adequate one-time training.

Motion detection with NMR markers using real-time field tracking in the laboratory frame.

PURPOSE: To enhance the utility of motion detection with nuclear magnetic resonance (NMR) markers by removing the need for sequence-dependent calibration. METHODS: Two sets of NMR markers are used for simultaneous observation of magnetic field dynamics during imaging procedures. A set of stationary markers at known positions in the laboratory frame serves to determine the field evolution in that frame. Concurrent recording from a set of head-mounted markers then permits calculating their lab-frame positions and derived rigid-body motion parameters. The precision and accuracy of this approach are evaluated relative to current calibration-based solutions. Use for prospective motion correction is then demonstrated in high-resolution imaging of long scan duration. RESULTS: Motion detection with real-time field tracking overcomes the need for explicit calibration without compromising precision, which is assessed at 10 to 30 µm. Relative to full conventional calibration, it is found to offer superior robustness against thermal drift. Relative to more economical modes of calibration, it achieves substantially higher accuracy. Prospective motion correction based on real-time field tracking resulted in consistently high image quality even when head motion exceeded the image resolution by one order of magnitude. CONCLUSION: Real-time field tracking enables motion detection with NMR markers without calibration overhead and thus overcomes a key obstacle toward routine use. In addition, it renders this mode of motion tracking more robust against system imperfections.

High-resolution short-T2 MRI using a high-performance gradient.

PURPOSE: To achieve high resolution in imaging of short-T2 materials and tissues by using a high-performance human-sized gradient insert with strength up to 200 mT/m and 100% duty cycle. METHODS: Dedicated short-T2 methodology and hardware are used, such as the pointwise encoding time reduction with radial acquisition (PETRA) technique with modulated excitation pulses, optimized radio-frequency hardware, and a high-performance gradient insert. A theoretical analysis of actual spatial resolution and SNR is provided to support the choice of scan parameters and interpretation of the results. Imaging is performed in resolution phantoms, animal specimen, and human volunteers at both conventional and maximum available gradient strengths and compared using image subtraction. RESULTS: Calculations suggest that increasing gradient strength beyond conventional values considerably improves both actual resolution and SNR efficiency in short-T2 imaging. Resolution improvements are confirmed in all investigated samples, in particular 2 mm slots were resolved in a hard-plastic plate with T2 ≈ 10 µs and in vivo musculoskeletal images were acquired at isotropic 200 µm resolution. Expected improvements in signal yield are realized in fine structures benefitting from high resolution but to less extent in regions of low contrast where decay-related blurring leads to signal overlap between neighboring locations. CONCLUSION: Strong gradients with high duty cycle enable short-T2 imaging at unprecedentedly high resolution, holding the potential for improving MRI of, eg, bone, tendon, lung, or teeth. Moreover, it allows direct access of tissues with T2 of tens of microseconds such as myelin or collagen.

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.

Filling the dead-time gap in zero echo time MRI: Principles compared.

PURPOSE: MRI of tissues with short coherence lifetimes T2 or T2* can be performed efficiently using zero echo time (ZTE) techniques such as algebraic ZTE, pointwise encoding time reduction with radial acquisition (PETRA), and water- and fat-suppressed proton projection MRI (WASPI). They share the principal challenge of recovering data in central k-space missed due to an initial radiofrequency dead time. The purpose of this study was to compare the three techniques directly, with a particular focus on their behavior in the presence of ultra-short-lived spins. METHODS: The most direct comparison was enabled by aligning acquisition and reconstruction strategies of the three techniques. Image quality and short- T2* performance were investigated using point spread functions, 3D simulations, and imaging of phantom and bone samples with short (<1 ms) and ultra-short (<100 µs) T2*. RESULTS: Algebraic ZTE offers favorable properties but is limited to k-space gaps up to approximately three Nyquist dwells. At larger gaps, PETRA enables robust imaging with little compromise in image quality, whereas WASPI may be prone to artifacts from ultra-short T2* species. CONCLUSION: For small k-space gaps (<4 dwells) and T2* much larger than the dead time, all techniques enable artifact-free short- T2* MRI. However, if these requirements are not fulfilled careful consideration is needed and PETRA will generally achieve better image quality. Magn Reson Med 79:2036-2045, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Prospective motion correction with NMR markers using only native sequence elements.

PURPOSE: To develop a method of tracking active NMR markers that requires no alterations of common imaging sequences and can be used for prospective motion correction (PMC) in brain MRI. METHODS: Localization of NMR markers is achieved by acquiring short signal snippets in rapid succession and evaluating them jointly. To spatially encode the markers, snippets are timed such that signal phase is accrued during sequence intervals with suitably diverse gradient actuation. For motion tracking and PMC in brain imaging, the markers are mounted on a lightweight headset. PMC is then demonstrated with high-resolution T2 *- and T1 -weighted imaging sequences in the presence of instructed as well as residual unintentional head motion. RESULTS: With both unaltered sequences, motion tracking was achieved with precisions on the order of 10 µm and 0.01° and temporal resolution of 48 and 39 ms, respectively. On this basis, PMC improved image quality significantly throughout. CONCLUSION: The proposed approach permits high-precision motion tracking and PMC with standard imaging sequences. It does so without altering sequence design and thus overcomes a key hindrance to routine motion tracking with NMR markers. Magn Reson Med 79:2046-2057, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Physiology recording with magnetic field probes for fMRI denoising.

Physiological noise originating in cardiovascular and respiratory processes is a substantial confound in BOLD fMRI. When unaccounted for it reduces the temporal SNR and causes error in inferred brain activity and connectivity. Physiology correction typically relies on auxiliary measurements with peripheral devices such as ECG, pulse oximeters, and breathing belts. These require direct skin contact or at least a tight fit, impairing subject comfort and adding to the setup time. In this work, we explore a touch-free alternative for physiology recording, using magnetic detection with NMR field probes. Placed close to the chest such probes offer high sensitivity to cardiovascular and respiratory dynamics without mechanical contact. This is demonstrated by physiology regression in a typical fMRI scenario at 7T, including validation against standard devices. The study confirms essentially equivalent performance of noise models based on conventional recordings and on field probes. It is shown that the field probes may be positioned in the subject's back such that they could be readily integrated in the patient table.

Analysis and correction of field fluctuations in fMRI data using field monitoring.

This work investigates the role of magnetic field fluctuations as a confound in fMRI. In standard fMRI experiments with single-shot EPI acquisition at 3 Tesla the uniform and gradient components of the magnetic field were recorded with NMR field sensors. By principal component analysis it is found that differences of field evolution between the EPI readouts are explainable by few components relating to slow and within-shot field dynamics of hardware and physiological origin. The impact of fluctuating field components is studied by selective data correction and assessment of its influence on image fluctuation and SFNR. Physiological field fluctuations, attributed to breathing, were found to be small relative to those of hardware origin. The dominant confounds were hardware-related and attributable to magnet drift and thermal changes. In raw image time series, field fluctuation caused significant SFNR loss, reflected by a 67% gain upon correction. Large part of this correction can be accomplished by traditional image realignment, which addresses slow and spatially uniform field changes. With realignment, explicit field correction increased the SFNR on the order of 6%. In conclusion, field fluctuations are a relevant confound in fMRI and can be addressed effectively by retrospective data correction. Based on the physics involved it is anticipated that the advantage of full field correction increases with field strength, with non-Cartesian readouts, and upon phase-sensitive BOLD analysis.

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.

Utility of real-time field control in T2 *-Weighted head MRI at 7T.

PURPOSE: Real-time field control can serve to reduce respiratory field perturbations during T2 * imaging at high fields. This work investigates the effectiveness of this approach in relation to key variables such as patient physique, breathing patterns, slice location, and the choice of sequence. METHODS: To cover variation in physical constitution and breathing behavior, volunteers with a wide range of body-mass-indices were asked to breathe either normally or deeply during T2 *-weighted image acquisition at 7T. Ensuing field fluctuation was countered by real-time field control or merely recorded in reference experiments. The impact of the control system on image quality was assessed by classifying and grading artifacts related to field fluctuation. RESULTS: The amplitude of respiratory field changes and related artifacts were generally stronger for subjects with higher body-mass-index and for lower slices. Field control was found effective at mitigating all five types of artifacts that were studied. Overall image quality was systematically improved. Residual artifacts in low slices are attributed to insufficient spatial order of the control system. CONCLUSION: Real-time field control was found to be a robust means of countering respiratory field perturbations in variable conditions encountered in high-field brain imaging. Reducing net fluctuation, it generally expands the feasibility of high-field T2 * imaging toward challenging patients and brain regions. Magn Reson Med 76:430-439, 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.

Monitoring, analysis, and correction of magnetic field fluctuations in echo planar imaging time series.

PURPOSE: To assess the utility of concurrent magnetic field monitoring for observing and correcting for variations in k-space trajectories and global background fields that occur in single-shot echo planar imaging (EPI) time series as typically used in functional MRI (fMRI). METHODS: Field monitoring was performed using an array of NMR field probes operated concurrently with series of single-shot EPI acquisitions from a static phantom. The observed fluctuations in field evolution were analyzed in terms of their temporal and spatial behavior at the field level as well as at the level of reconstructed image series. The potential to correct for such fluctuations was assessed by accounting for them upon image reconstruction. An indication of the number and relative magnitude of underlying effects was obtained via principal component analysis. RESULTS: Trajectory and global field variations were found to induce substantial image fluctuations. Global field fluctuations induced standard deviations in image intensity up to 31%. Fluctuations in the trajectory induced ghosting artifacts with standard deviations up to 2%. Concurrent magnetic field monitoring reduced the fluctuations in the EPI time series to a maximum of 1.2%. CONCLUSION: Concurrent magnetic field monitoring holds the potential to improve the net sensitivity of fMRI by reducing signal fluctuations unrelated to brain activity.

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.

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.