Dynamic Field Mapping Research Articles

FPGA-Based Magnetic Field Camera for Dynamic Magnetic Field Mapping

FPGA-Based Magnetic Field Camera for Dynamic Magnetic Field Mapping

Dynamic field mapping and distortion correction using single-shot blip-rewound EPI and joint multi-echo reconstruction.

To develop a method for dynamic mapping and distortion correction. A blip-rewound EPI trajectory was developed to acquire multiple 2D EPI images in a single readout with an interleaved order, which allows a short TE difference. A joint multi-echo reconstruction was utilized to exploit the shared information between EPI images. The reconstructed images from each readout are combined to produce a final magnitude image. A map is calculated from the phase of these images for distortion correction. The efficacy of the proposed method is assessed with phantom and in vivo experiments. The performance of the proposed method in the presence of subject motion is also investigated. Compared to conventional multi-echo EPI, the proposed method allows dynamic mapping at matched resolution with a much shorter TR. Phantom and in vivo results show that the proposed method can provide a comparable magnitude image as conventional single-shot EPI. The maps calculated from the proposed method are consistent with conventional multi-echo EPI in the phantom experiment. For in vivo experiments, the proposed method provides a more accurate estimation of than conventional multi-echo EPI, which is prone to phase wrapping problems due to the long TE difference. In-vivo scan with subject motion shows the proposed dynamic field mapping method can improve the temporal stability of EPI time series compared to gradient echo (GRE) based static field mapping. The proposed method allows accurate dynamic mapping for robust distortion correction without compromising spatial or temporal resolution.

Predicting dynamic, motion-related changes in B0 field inthe brain at a 7T MRI using a subject-specific fine-trainedU-net.

Subject movement during the MR examination is inevitable and causes not only image artifacts but also deteriorates the homogeneity of the main magnetic field (B0 ), which is a prerequisite for high quality data. Thus, characterization of changes to B0 , for example induced by patient movement, is important for MR applications that are prone to B0 inhomogeneities. We propose a deep learning based method to predict such changes within the brain from the change of the head position to facilitate retrospective or even real-time correction. A 3D U-net was trained on in vivo gradient-echo brain 7T MRI data. The input consisted of B0 maps and anatomical images at an initial position, and anatomical images at a different head position (obtained by applying a rigid-body transformation on the initial anatomical image). The output consisted of B0 maps at the new head positions. We further fine-trained the network weights to each subject by measuring a limited number of head positions of the given subject, and trained the U-net with these data. Our approach was compared to established dynamic B0 field mapping via interleaved navigators, which suffer from limited spatial resolution and the need for undesirable sequence modifications. Qualitative and quantitative comparison showed similar performance between an interleaved navigator-equivalent method and proposed method. It is feasible to predict B0 maps from rigid subject movement and, when combined with external tracking hardware, this information could be used to improve the quality of MR acquisitions without the use of navigators.

Patient specific distortion detection and mitigation in MR images used for stereotactic radiosurgery

Objective. In MRI-based radiation therapy planning, mitigating patient-specific distortion with standard high bandwidth scans can result in unnecessary sacrifices of signal to noise ratio. This study investigates a technique for distortion detection and mitigation on a patient specific basis. Approach. Fast B0 mapping was performed using a previously developed technique for high-resolution, large dynamic range field mapping without the need for phase unwrapping algorithms. A phantom study was performed to validate the method. Distortion mitigation was validated by reducing geometric distortion with increased acquisition bandwidth and confirmed by both the B0 mapping technique and manual measurements. Images and contours from 25 brain stereotactic radiosurgery patients and 95 targets were analyzed to estimate the range of geometric distortions expected in the brain and to estimate bandwidth required to keep all treatment targets within the ±0.5 mm iso-distortion contour. Main Results. The phantom study showed, at 3 T, the technique can measure distortions with a mean absolute error of 0.12 mm (0.18 ppm), and a maximum error of 0.37 mm (0.6 ppm). For image acquisition at 3 T and 1.0 mm resolution, mean absolute distortion under 0.5 mm in patients required bandwidths from 109 to 200 Hz px−1 for patients with the least and most distortion, respectively. Maximum absolute distortion under 0.5 mm required bandwidths from 120 to 390 Hz px−1. Significance. The method for B0 mapping was shown to be valid and may be applied to assess distortion clinically. Future work will adapt the readout bandwidth to prospectively mitigate distortion with the goal to improve radiosurgery treatment outcomes by reducing healthy tissue exposure.

Dynamic off-resonance correction for spiral real-time MRI of speech.

To improve the depiction and tracking of vocal tract articulators in spiral real-time MRI (RT-MRI) of speech production by estimating and correcting for dynamic changes in off-resonance. The proposed method computes a dynamic field map from the phase of single-TE dynamic images after a coil phase compensation where complex coil sensitivity maps are estimated from the single-TE dynamic scan itself. This method is tested using simulations and in vivo data. The depiction of air-tissue boundaries is evaluated quantitatively using a sharpness metric and visual inspection. Simulations demonstrate that the proposed method provides robust off-resonance correction for spiral readout durations up to 5 ms at 1.5T. In -vivo experiments during human speech production demonstrate that image sharpness is improved in a majority of data sets at air-tissue boundaries including the upper lip, hard palate, soft palate, and tongue boundaries, whereas the lower lip shows little improvement in the edge sharpness after correction. Dynamic off-resonance correction is feasible from single-TE spiral RT-MRI data, and provides a practical performance improvement in articulator sharpness when applied to speech production imaging.

Fast rotary nonlinear spatial acquisition (FRONSAC) imaging.

Nonlinear spatial encoding magnetic fields (SEMs) have been studied to reconstruct images from a minimum number of echoes. Previous work has also explored single shot trajectories in nonlinear SEMs. However, the search continues for optimal schemes that apply nonlinear SEMs to improve spatial encoding efficiency and image quality. We enhance the encoding efficiency of standard linear gradient trajectories by adding a rapidly rotating nonlinear SEM of moderate amplitude, the so called FRONSAC (Fast ROtary Nonlinear Spatial ACquisition) imaging. This additional gradient greatly improves the image quality of highly undersampled single-shot trajectories, including EPI, Spiral, and Rosette trajectories. Our simulations, including noise and dephasing effects, test the effect of adding FRONSAC gradients, demonstrating the applicability of this approach. Performance is explained by demonstrating the additional k-space sampling the nonlinear gradient provides. Studies of the optimal amplitude and frequency of the additional FRONSAC field are presented, and the role of enhanced sampling during the readout demonstrated. Dynamic field mapping in a second-order gradient system shows the proposed gradient waveforms are feasible. Images resulting from highly undersampled existing k-space trajectories, such as EPI, Spiral, and Rosette, are greatly enhanced simply by adding a rotating nonlinear SEM field.

Note: Dynamic strain field mapping with synchrotron X-ray digital image correlation.

We present a dynamic strain field mapping method based on synchrotron X-ray digital image correlation (XDIC). Synchrotron X-ray sources are advantageous for imaging with exceptional spatial and temporal resolutions, and X-ray speckles can be produced either from surface roughness or internal inhomogeneities. Combining speckled X-ray imaging with DIC allows one to map strain fields with high resolutions. Based on experiments on void growth in Al and deformation of a granular material during Kolsky bar/gas gun loading at the Advanced Photon Source beamline 32ID, we demonstrate the feasibility of dynamic XDIC. XDIC is particularly useful for dynamic, in-volume, measurements on opaque materials under high strain-rate, large, deformation.

High‐resolution, large dynamic range field map estimation

PurposeWe present a theory and a corresponding method to compute high‐resolution field maps over a large dynamic range.Theory and MethodsWe derive a closed‐form expression for the error in the field map value when computed from two echoes. We formulate an optimization problem to choose three echo times which result in a pair of maximally distinct error distributions. We use standard field mapping sequences at the prescribed echo times. We then design a corresponding estimation algorithm which takes advantage of the optimized echo times to disambiguate the field offset value.ResultsWe validate our method using high‐resolution images of a phantom at 7T. The resulting field maps demonstrate robust mapping over both a large dynamic range, and in low SNR regions. We also present high‐resolution offset maps in vivo using both, GRE and multiecho gradient echo sequences. Even though the proposed echo time spacings are larger than the well known phase aliasing cutoff, the resulting field maps exhibit a large dynamic range without the use of phase unwrapping or spatial regularization techniques.ConclusionWe demonstrate a novel three‐echo field map estimation method which overcomes the traditional noise‐dynamic range trade‐off. Magn Reson Med 71:105–117, 2014. © 2013 Wiley Periodicals, Inc.

Quantification and correction of respiration induced dynamic field map changes in fMRI using 3D single shot techniques

Respiration induced dynamic field map changes in the brain are quantified and the influence on the magnitude signal (physiological noise) is investigated. Dynamic off-resonance correction allows to reduce the signal fluctuations overlaying the blood oxygenation level dependent signal in T2*-weighted functional imaging. A single-shot whole brain imaging technique with 100 ms temporal resolution was used to measure dynamic off-resonance maps that were calculated from the incremental changes of the image phase. These off-resonance maps are then used to dynamically update the off-resonance corrected reconstruction. A global resonance offset and a pronounced gradient in head-foot direction were identified as the main components of the change during a respiration cycle. On average, correction for these fluctuations decreases the magnitude fluctuations by around 30%. Single shot 3D imaging allows for a robust quantification of dynamic off-resonance changes in the brain. Correction for these fluctuations removes the physiological noise component associated with dynamic point spread function changes.

High‐resolution, large dynamic range field map estimation

We present a theory and a corresponding method to compute high-resolution field maps over a large dynamic range. We derive a closed-form expression for the error in the field map value when computed from two echoes. We formulate an optimization problem to choose three echo times which result in a pair of maximally distinct error distributions. We use standard field mapping sequences at the prescribed echo times. We then design a corresponding estimation algorithm which takes advantage of the optimized echo times to disambiguate the field offset value. We validate our method using high-resolution images of a phantom at 7T. The resulting field maps demonstrate robust mapping over both a large dynamic range, and in low SNR regions. We also present high-resolution offset maps in vivo using both, GRE and multiecho gradient echo sequences. Even though the proposed echo time spacings are larger than the well known phase aliasing cutoff, the resulting field maps exhibit a large dynamic range without the use of phase unwrapping or spatial regularization techniques. We demonstrate a novel three-echo field map estimation method which overcomes the traditional noise-dynamic range trade-off.

Joint Estimation and Correction of Geometric Distortions for EPI Functional MRI Using Harmonic Retrieval

Magnetic resonance imaging (MRI) uses applied spatial variations in the magnetic field to encode spatial position. Therefore, nonuniformities in the main magnetic field can cause image distortions. In order to correct the image distortions, it is desirable to simultaneously acquire data with a field map in registration. We propose a joint estimation (JE) framework with a fast, noniterative approach using harmonic retrieval (HR) methods, combined with a multi-echo echo-planar imaging (EPI) acquisition. The connection with HR establishes an elegant framework to solve the JE problem through a sequence of 1-D HR problems in which efficient solutions are available. We also derive the condition on the smoothness of the field map in order for HR techniques to recover the image with high signal-to-noise ratio. Compared to other dynamic field mapping methods, this method is not constrained by the absolute level of the field inhomogeneity over the slice, but relies on a generous pixel-to-pixel smoothness. Moreover, this method can recover image, field map, and T2* map simultaneously.

A new EPI‐based dynamic field mapping method: Application to retrospective geometrical distortion corrections

To retrospectively correct for geometrical distortions, a new dynamic field mapping method suitable for dynamic single-shot gradient-echo type echo-planar imaging (GRE-EPI) is proposed. The method requires a single volume additional acquisition and allows the extraction of a field map from each phase volume, assuming invariance across time of the echo time-independent phase component. Performances of the method are assessed using three sets of experiments: the first tests the prerequisite and the modeling; the second tests the method with time-dependent geometrical distortions; and the third presents a comparison with two other methods. Our results legitimize the modeling procedure and demonstrate that the dynamic method is less sensitive to noise than the other methods. A theoretical explanation for this is proposed in the discussion section. Given the minor increase in the acquisition time, this method is well suited for functional magnetic resonance imaging; prospective direction.

Dynamic Field Mapping for Obtaining the Current Distribution in High-Temperature Superconducting Tapes

The magnetic field profile on the surface of Bi-2223/Ag tapes has been dynamically measured across the width of the samples. The experimental technique uses a Hall-probe array with 7 sensors connected to a multiple channel lock-in amplifier especially programmed for fast and synchronous data acquisition measurements. The speed of the system is high enough to measure real-time profiles with 7 probes and 50 Hz sine current through the sample. A numerical method to estimate the current distribution inside the tapes using the measured field profile data is proposed. The inverse problem has been solved using certain assumptions on the current distribution in the superconductor. Validation of the results has been done by comparison with finite element method simulations.

Dynamic field map estimation using a spiral-in/spiral-out acquisition.

The long readout times of single-shot acquisitions and the high field strengths desired for functional MRI (fMRI) using blood oxygenation level-dependent (BOLD) contrast make functional scans sensitive to magnetic field inhomogeneity. If it is not corrected during image reconstruction, field inhomogeneity can cause geometric distortions in the images when Cartesian k-space trajectories are used or blurring with spiral acquisitions. Many traditional methods to correct for field inhomogeneity distortions rely on a static field map measured with the use of images that are themselves distorted. In this work, we employ a regularized least-squares approach to jointly estimate both the undistorted image and field map at each acquisition using a spiral-in/spiral-out pulse sequence. Simulation and phantom studies show that this method is accurate and stable over a time series. Human functional studies show that the jointly estimated field map may be more accurate than standard field map estimates in the presence of respiration-induced phase oscillations, leading to better detection of functional activation. The proposed method measures a dynamic field map that accurately tracks magnetic field drift and respiration-induced phase oscillations during the course of a functional study.

Multispectral image analysis method to determine dynamic fluid saturation distribution in two-dimensional three-fluid phase flow laboratory experiments

The need for measuring dynamic fluid saturation distribution in multi-dimensional three-fluid phase flow experiments is hampered by lack of appropriate techniques to monitor full field transient flow phenomena. There is no conventional technique able to measure dynamic three-fluid phase saturation at several array points of the flow field at the same time. A multispectral image analysis technique was developed to determine dynamic NAPL, water and air saturation distribution in two-dimensional three-fluid phase laboratory experiments. Using a digital near-infrared camera, images of sand samples with various degrees of NAPL, water and air saturation were taken, under constant lighting conditions and within three narrow spectral bands of the visible and near-infrared spectrum. It was shown that the optical density defined for the reflected luminous intensity was a linear function of the NAPL and the water saturation for each spectral band and for any two and three-fluid phase systems. This allowed the definition of dimensionless lump reflection coefficients for the NAPL and the water phase within each spectral band. Consequently, at any given time, two images taken within two different spectral bands provided two linear equations which could be solved for the water and the NAPL saturation. The method was applied to two-dimensional three-phase flow experiments, which were conducted to investigate the migration and the distribution of LNAPL in the vadose zone. The method was used to obtain continuous, quantitative and dynamic full field mapping of the NAPL saturation as well as the variation of the water and the air saturation during NAPL flow. The method provides a non-destructive and non-intrusive tool for studying multiphase flow for which rapid changes in fluid saturation in the entire flow domain is difficult to measure using conventional techniques.