Gradient Switching Research Articles

Fast and silent MRI using nonlinear gradient fields at the ultrasonic gradient switching frequency of 20 kHz with a Point Spread Function framework reconstruction.

To demonstrate the feasibility of using a nonlinear gradient field for spatial encoding at the ultrasonic switching frequency of 20 kHz and present a framework to reconstruct data acquired in this way. Nonlinear encoding at 20 kHz was realized by using a single-axis silent gradient insert for imaging in the periphery, that, is the nonlinear region, of the gradient field. The gradient insert induces a rapidly oscillating gradient field in the phase-encode direction, which enables nonlinear encoding when combined with a Cartesian readout from the linear whole-body gradients. Data from a 2D gradient echo sequence were reconstructed using a point spread function (PSF) framework. Accelerated scans were also simulated via retrospective undersampling (R = 1 to R = 8) to determine the effectiveness of the PSF-framework for accelerated imaging. Using a nonlinear gradient field switched at 20 kHz and the PSF-framework resulted in images of comparable quality to images from conventional Cartesian linear encoding. At increased acceleration factors (R ≤ 8), the PSF-framework outperformed linear SENSE reconstructions by improved controlling of aliasing artifacts. Using the PSF-framework, images of comparable quality to conventional SENSE reconstructions are possible via combining traditional linear and ultrasonic oscillating nonlinear encoding fields. Using nonlinear gradient fields relaxes the demand for strictly linear gradient fields, enabling much higher slew rates with a reduced risk of peripheral nerve stimulation or cardiac stimulation, which could aid in extension to ultrasonic whole-body MRI. The lack of aliasing artifacts also highlights the potential of accelerated imaging using the PSF-framework.

Large-Scale Data Processing and Machine Learning Analysis Model Based on Distributed Algorithm

The model of large-scale data processing and ML(machine learning) analysis based on DA(distributed algorithm) is a powerful computing method, which aims at processing huge data sets and performing efficient ML analysis. In this paper, a cluster topology driver module based on gradient switching and aggregate communication is designed, and its core goal is to adapt the distributed system to various underlying network topologies. By designing decentralized gradient exchange algorithm and aggregate communication framework, the parallel transmission ability of multi-interface network can be fully exerted, thus improving the model synchronization efficiency of ML task. The experimental results show that the cluster topology driver module can provide better performance than the existing methods in terms of training convergence, cluster scalability and communication overhead. Large-scale data processing and ML analysis model based on DA is widely used in processing massive data and realizing complex analysis tasks.

Sophisticated Study of Time, Frequency and Statistical Analysis for Gradient-Switching-Induced Potentials during MRI.

Magnetic resonance imaging (MRI) is a standard procedure in medical imaging, on a par with echography and tomodensitometry. In contrast to radiological procedures, no harmful radiation is produced. The constant development of magnetic resonance imaging (MRI) techniques has enabled the production of higher resolution images. The switching of magnetic field gradients for MRI imaging generates induced voltages that strongly interfere with the electrophysiological signals (EPs) collected simultaneously. When the bandwidth of the collection amplifiers is higher than 150 Hz, these induced voltages are difficult to eliminate. Understanding the behavior of these artefacts contributes to the development of new digital processing tools for better quality EPs. In this paper, we present a study of induced voltages collected in vitro using a device (350 Hz bandwidth). The experiments were conducted on a 1.5T MRI machine with two MRI sequences (fast spin echo (FSE) and cine gradient echo (CINE)) and three slice orientations. The recorded induced voltages were then segmented into extract patterns called "artefact puffs". Two analysis series, "global" and "local", were then performed. The study found that the temporal and frequency characteristics were specific to the sequences and orientations of the slice and that, despite the pseudo-periodic character of the artefacts, the variabilities within the same recording were significant. These evolutions were confirmed by two stationarity tests: the Kwiatkowski-Phillips-Schmidt-Shin (KPSS) and the time-frequency approach. The induced potentials, all stationary at the global scale, are no longer stationary at the local scale, which is an important issue in the design of optimal filters adapted to reduce MRI artifacts contaminating a large bandwidth, which varies between 0 and 500 Hz.

Broadly applicable methods for the detection of artefacts in electroencephalography acquired simultaneously with hemodynamic recordings

Electroencephalography (EEG) data, acquired simultaneously with magnetic resonance imaging (MRI), must be corrected for artefacts related to MR gradient switches (GS) and the cardioballistic (CB) effect. Canonical approaches require additional signal acquisition for artefact detection (e.g., MR volume onsets, ECG), without which the EEG data would be rendered uncleanable from these artefacts.•We present two broadly applicable methods for artefact detection based on peak detection combined with temporal constraints with respect to periodicity directly from the EEG data itself; no additional signals are required.We validated the performance of our methods versus the two canonical approaches for detection of GS/CB artefact, respectively, on 26 healthy human EEG-functional MRI resting-state datasets. Utilising various performance metrics, we found our methods to perform as well as – and sometimes better than - the canonical standard approaches. With as little as one EEG channel recording, our methods can be applied to detect GS/CB artefacts in EEG data acquired simultaneously with MRI in the absence of MR volume onsets and/or an ECG recording. The detected artefact onsets can then be fed into the standard artefact correction software.

A three‐dimensional positivity‐preserving and conservative multimoment finite‐volume transport model on a cubed‐sphere grid

AbstractA three‐dimensional positivity‐preserving and conservative transport model on a cubed‐sphere grid is presented using a multimoment finite‐volume (MMFV) method. The three‐point fourth‐order MMFV method with boundary gradient switching projection is utilized to ensure nonoscillatory numerical solutions in the horizontal, while the simple and computationally efficient semi‐Lagrangian piecewise rational method is used to allow a large time step in the vertical. By means of horizontal and vertical separation, a Strang‐type splitting technique is adopted to advance the numerical solutions in time. Non‐negative corrections are made to achieve positive‐definite preservation exactly. The resulting transport model is inherently conservative, high‐accuracy, and positive‐definite, and allows a large Courant number no matter how dense the vertical grid is. Various benchmark test cases, including three representative advection tests and an idealized terminator chemistry test, are conducted to validate the performance of the model presented. The numerical results indicate that the model is quite competitive in comparison with the existing models in the literature and looks very promising in real atmospheric simulations.

RF shielding designs for birdcage coils for preclinical MRI at 9.4 T

Birdcage coils are widely used in preclinical MRI as they perform well, allow for quadrature drive, and can provide a homogeneous transmit field. Unlike in larger bore scanners, an RF shield is essential to avoid strong cross-talk with gradient coils that are in close proximity. However, gradient switching induces eddy currents that heat the shield and coil and impair the temporal signal-to-noise ratio (tSNR). The motivation of this study is to investigate the performance of different designs of RF shields on a birdcage coil used for high resolution functional MRI of small primates at 9.4 T. We found the choice of materials for RF shields significantly affected ghosting and tSNR in fMRI scans. Both ultrathin foils and a slotted pattern reduce eddy currents and improve imaging quality. Our results also demonstrate that a 9-um-thick copper foil is sufficiently thin to reduce the eddy current effects for high-resolution fMRI scans and there is no need for high-cost 4-um-thick foil. For slotted shields, our results demonstrate that the number of slots should be carefully considered, and an excessive number of slots can lead to a lower SNR and tSNR. We believe the results from this study can be used as a reference to design future RF coil shields selection for preclinical scanners.

Tracking of rigid head motion during MRI using anEEG system.

PurposeTo demonstrate a novel method for tracking of head movements during MRI using electroencephalography (EEG) hardware for recording signals induced by native imaging gradients.Theory and MethodsGradient switching during simultaneous EEG–fMRI induces distortions in EEG signals, which depend on subject head position and orientation. When EEG electrodes are interconnected with high‐impedance carbon wire loops, the induced voltages are linear combinations of the temporal gradient waveform derivatives. We introduce head tracking based on these signals (CapTrack) involving 3 steps: (1) phantom scanning is used to characterize the target sequence and a fast calibration sequence; (2) a linear relation between changes of induced signals and head pose is established using the calibration sequence; and (3) induced signals recorded during target sequence scanning are used for tracking and retrospective correction of head movement without prolonging the scan time of the target sequence. Performance of CapTrack is compared directly to interleaved navigators.ResultsHead‐pose tracking at 27.5 Hz during echo planar imaging (EPI) was demonstrated with close resemblance to rigid body alignment (mean absolute difference: [0.14 0.38 0.15]‐mm translation, [0.30 0.27 0.22]‐degree rotation). Retrospective correction of 3D gradient‐echo imaging shows an increase of average edge strength of 12%/−0.39% for instructed/uninstructed motion with CapTrack pose estimates, with a tracking interval of 1561 ms and high similarity to interleaved navigator estimates (mean absolute difference: [0.13 0.33 0.12] mm, [0.28 0.15 0.22] degrees).ConclusionMotion can be estimated from recordings of gradient switching with little or no sequence modification, optionally in real time at low computational burden and synchronized to image acquisition, using EEG equipment already found at many research institutions.

High‐order conservative transport on Yin‐Yang grids using the multi‐moment constrained finite volume method

AbstractA numerical constraint was developed to improve the global and local conservation for the Yin‐Yang grid system, which has been known as one of the quasi‐uniform grids on a sphere. Two‐dimensional cubic mass distribution within an individual mesh was assumed, to describe the subgrid‐scale structure of local properties and to ensure high‐order‐accuracy mass flux specification for the Yin‐Yang boundary. A three‐point Multi‐moment Constrained finite Volume scheme, in cooperation with a fourth‐order Runge–Kutta scheme, was selected for numerical transport with the help of Boundary Gradient Switching for oscillation suppression. The new scheme was tested with a couple of idealized numerical experiments in advection and shallow‐water models on the Yin‐Yang grid to verify its performance. Numerical results confirmed the exact mass conservation in spherical advection problems, and the numerical convergence rate reached fourth order in both advection and shallow‐water models. Computational stability, shape‐preserving and numerical oscillation‐free properties were also revealed in the nonlinear testing problems.

Acoustic Control through Gradient Coil Design Using a Finite-Difference-Based Method for MRI

An acoustic control scheme is proposed in this paper through the process of gradient coil design for magnetic resonance imaging (MRI). With a finite-difference-based method, the stream function and coil displacement caused by fast gradient switching can be unified by a simplified momentum equation, which can be incorporated into the conventional gradient coil design. A three-dimensional transverse gradient coil with an edge-connected cylindrical structure is used as a design example to verify the proposed design method. In addition, an acoustic model is established to simulate the sound pressure level (SPL).In the model, two hemispherical air volumes are added flush with the ends of the cylindrical main magnet to mimic the free propagation of sound waves on the boundaries. The simulation results show that by optimizing coil displacement, the overall SPL can be attenuated by 4 dB over the frequency range from 0 to 3000 Hz with the displacement reduced by about 50%, at the cost of a figure of merit (FOM) loss by about 8%. Therefore, the proposed acoustic control scheme can be used as a complement to conventional acoustic control methods for further noise reduction.

Sparse and continuous sampling approaches to fMRI of overt vocalization tasks

Sparse temporal sampling has become the dominant paradigm for functional magnetic resonance imaging studies of auditory stimuli or verbal responses, as it allows the presentation or production of stimuli during the relatively quiet periods when there is no gradient switching and ensures that task-related movements are not occurring during scan acquisitions. To date, however, there has been no direct comparison between sparse and continuous acquisition protocols for overt auditory-verbal studies (i.e., speaking or singing). The aim of this study was to determine whether sparse temporal sampling would reduce movement artefacts and show better network activation for overt singing compared to continuous imaging. Fifteen healthy adults performed the same overt singing task under both sparse and continuous scanning conditions. We noted significant variations in signal intensity between adjacent slices in our sparse acquisition, with (odd-numbered) slices acquired in the second half of each volume acquisition being of lower intensity and showing less reliable task-related activation, and thus requiring the removal of these slices prior to preprocessing. Edge artefacts, presumably due to movement, were observed in both acquisition types at a subthreshold level, although ventricular space artefacts were more apparent in the continuous data. However, statistical comparison revealed no significant differences in functional activation nor in motion correction parameters. Our results show that sparse imaging has the potential to introduce significant image artefacts affecting downstream analyses. While sparse sampling provides benefits that may be essential for certain studies (e.g., periods free from scanner noise), the technical requirements of such sequences should not be overlooked and inspection of raw data is essential. Our data also show that continuous imaging can be used for overt response auditory-verbal studies and may be of wider utility than previously appreciated.

Quasi-Stable Salt Gradient and Resistive Switching in Solid-State Nanopores.

Understanding and control of ion transport in a fluidic channel is of crucial importance for iontronics. The present study reports on quasi-stable ionic current characteristics in a SiNx nanopore under a salinity gradient. An intriguing interplay between electro-osmotic flow and local ion density distributions in a solid-state pore is found to induce highly asymmetric ion transport to negative differential resistance behavior under a 100-fold difference in the cross-membrane salt concentrations. Meanwhile, a subtle change in the salinity gradient profile led to observations of resistive switching. This peculiar characteristic was suggested to stem from quasi-stable local ion density around the channel that can be switched between two distinct states via the electro-osmotic flow under voltage control. The present findings may be useful for neuromorphic devices based on micro- and nanofluidic channels.

Temporal Signal-to-Noise Changes in Combined Multislice- and In-Plane-Accelerated Echo-Planar Imaging with a 20- and 64-Channel Coil

Echo-planar imaging (EPI) is the most common method of functional MRI for acquiring the blood oxygenation level-dependent (BOLD) contrast, allowing the acquisition of an entire brain volume within seconds. However, because imaging protocols are limited by hardware (e.g., fast gradient switching), researchers must compromise between spatial resolution, temporal resolution, or whole-brain coverage. Earlier attempts to circumvent this problem included developing protocols in which slices of a volume were acquired faster (i.e., in-plane acceleration (S)) or simultaneously (i.e., multislice acceleration (M)). However, applying acceleration methods can lead to a reduction in the temporal signal-to-noise ratio (tSNR): a critical measure of signal stability over time. Using a 20- and 64-channel receiver coil, we show that enabling S-acceleration consistently yielded a substantial decrease in tSNR, regardless of the receiver coil, whereas M-acceleration yielded less pronounced tSNR decrease. Moreover, tSNR losses tended to occur in temporal, insular, and medial brain regions and were more noticeable with the 20-channel coil, while with the 64-channel coil, the tSNR in lateral frontoparietal regions remained relatively stable up to six-fold M-acceleration producing comparable tSNR to that of no acceleration. Such methodological explorations can guide researchers and clinicians in optimizing imaging protocols depending on the brain regions under investigation.

Reduction of vibration‐induced signal loss by matching mechanical vibrational states: Application in high b‐value diffusion‐weighted MRS

Diffusion encoding gradients are known to yield vibrations of the typical clinical MR scanner hardware with a frequency of 20 to 30 Hz, which may lead to signal loss in diffusion-weighted MR measurements. This work proposes to mitigate vibration-induced signal loss by introducing a vibration-matching gradient (VMG) to match vibrational states during the 2 diffusion gradient pulses. A theoretical description of displacements induced by gradient switching was introduced and modeled by a 2-mass-spring-damper system. An additional preceding VMG mimicking timing and properties of the diffusion encoding gradients was added to a high b-value diffusion-weighted MR spectroscopy sequence. Laser interferometry was employed to measure 3D displacements of a phantom surface. Lipid ADC was assessed in water-fat phantoms and in vivo in the tibial bone marrow of 3 volunteers. The modeling and the laser interferometer measurements revealed that the displacement curves are more similar during the 2 diffusion gradients with the VMG compared to the standard sequence, resulting in less signal loss of the diffusion-weighted signal. Phantom results showed lipid ADC overestimation up to 119% with the standard sequence and an error of 5.5% with the VMG. An 18% to 35% lower coefficient of variation was obtained for in vivo lipid ADC measurement when employing the VMG. The application of the VMG reduces the signal loss introduced by hardware vibrations in a high b-value diffusion-weighted MRS sequence in phantoms and in vivo. Reference measurements based on laser interferometry and mechanical modelling confirmed the findings.

Artifacts in Simultaneous hdEEG/fMRI Imaging: A Nonlinear Dimensionality Reduction Approach.

Simultaneous recordings of electroencephalogram (EEG) and functional magnetic resonance imaging (fMRI) are at the forefront of technologies of interest to physicians and scientists because they combine the benefits of both modalities—better time resolution (hdEEG) and space resolution (fMRI). However, EEG measurements in the scanner contain an electromagnetic field that is induced in leads as a result of gradient switching slight head movements and vibrations, and it is corrupted by changes in the measured potential because of the Hall phenomenon. The aim of this study is to design and test a methodology for inspecting hidden EEG structures with respect to artifacts. We propose a top-down strategy to obtain additional information that is not visible in a single recording. The time-domain independent component analysis algorithm was employed to obtain independent components and spatial weights. A nonlinear dimension reduction technique t-distributed stochastic neighbor embedding was used to create low-dimensional space, which was then partitioned using the density-based spatial clustering of applications with noise (DBSCAN). The relationships between the found data structure and the used criteria were investigated. As a result, we were able to extract information from the data structure regarding electrooculographic, electrocardiographic, electromyographic and gradient artifacts. This new methodology could facilitate the identification of artifacts and their residues from simultaneous EEG in fMRI.

Simultaneous EEG-fMRI: Evaluating the Effect of the EEG Cap-Cabling Configuration on the Gradient Artifact.

Electroencephalography (EEG) data recorded during simultaneous EEG-fMRI experiments are contaminated by large gradient artifacts (GA). The amplitude of the GA depends on the area of the wire loops formed by the EEG leads, as well as on the rate of switching of the magnetic field gradients, which are essential for MR imaging. Average artifact subtraction (AAS), the most commonly used method for GA correction, relies on the EEG amplifier having a large enough dynamic range to characterize the artifact voltages. Low-pass filtering (250 Hz cut-off) is generally used to attenuate the high-frequency voltage fluctuations of the GA, but even with this precaution channel saturation can occur, particularly during acquisition of high spatial resolution MRI data. Previous work has shown that the ribbon cable, used to connect the EEG cap and amplifier, makes a significant contribution to the GA, since the cable geometry produces large effective wire-loop areas. However, by appropriately connecting the wires of the ribbon cable to the EEG cap it should be possible to minimize the overall range and root mean square (RMS) amplitude of the GA by producing partial cancelation of the cap and cable contributions. Here by modifying the connections of the EEG cap to a 1 m ribbon cable we were able to reduce the range of the GA for a high-resolution coronal echo planar Imaging (EPI) acquisition by a factor of ∼ 1.6 and by a factor of ∼ 1.15 for a standard axial EPI acquisition. These changes could potentially be translated into a reduction in the required dynamic range, an increase in the EEG bandwidth or an increase in the achievable image resolution without saturation, all of which could be beneficially exploited in EEG-fMRI studies. The re-wiring could also prevent the system from saturating when small subject movements occur using the standard recording bandwidth.

A cone-shaped gradient coil design for high-resolution MRI head imaging

Insertable head gradient coils offer significant advantages such as high gradient strength and fast gradient switching speed owing to shorter distances to the target region of interest than whole-body cylindrical coils. To produce superior gradient performance, the local head coil is typically designed with an asymmetric configuration to accommodate both the shoulders and head of a patient, leading to tough dimensional constraints and practical limits to the coil implementation. In this paper, we propose a new cone-shaped model to improve the performance of the asymmetric head coils and to mitigate patient claustrophobia. The primary coils are designed with a larger diameter at the patient end for access and a smaller diameter at the service end to bring wires closer to the human head, while the secondary coils are arranged on a cylindrical former to improve coil efficiency. Two cases are studied in this paper. Case I: inner bore size at the patient end (diameter 42 cm) is fixed as the design reference. In this case, inner diameters at any other position vary with the conical tilting angles. Compared with a set of conical gradient coils designed with tilting angles ranging from 0 to 14°, it is found that the optimal coil performance is achieved at the tilting angle of 14°. The key performance parameters have been improved by 100%–200% for the transverse coils, and about 50% for the longitudinal coils compared with the cylindrical counterpart with the reference bore size (that is, the same diameter of 42 cm). The conical coils also produce less heat in the gradient structure and lower acoustic noise in the field of view. Case II: inner bore size at the iso-centre (diameter 34 cm) is set as the design reference. It is also found that, compared with 34 cm diameter cylindrical coils, the conical transverse coil performance has been improved at an angle of 14°. The key coil performance increases by 20%–50% for transverse coil but decreases by 20%–40% for the longitudinal coil. However, compared with the tight cylindrical structure (e.g. 34 cm diameter), the tilting angle will provide patient-friendly space for imaging and handling, which can be critical for fMRI and other brain studies.

Studying the effects of metallic components of PET-insert on PET and MRI performance due to gradient switching

Inserting positron emission tomography (PET) detection modules inside an MRI bore imposes extra challenges owing to the behavior of metallic materials in a strong magnetic field. The metallic parts even when placed outside an MRI field of view may not only disturb MRI performance, but could also increase temperature and vibrations, leading to premature failure of PET electronics. To investigate the compatibility of detection modules inside 3 T, 7 T and 9.4 T MRI bore, a theoretical study of the metal induced artifacts originating from component materials of electronic circuit is presented. The LabPET II detection module and a modified version of it in which the connector was replaced by ball grid array (BGA) were studied. In addition, the effect of eddy current and the associated heat loss on the PET detection module have been examined using COMSOL Multiphysics® simulations for 10 kHz and 100 kHz gradient switching. Results show that displacement artifacts resulting from the presence of small amounts of ferromagnetic metal and the heating effects of metal due to gradient switching can be compensated by using the slightly modified LabPET II detection module. Thus, the LabPET II system would be MR-compatible with some minor adjustments to operate effectively inside an MRI bore without interfering with its performance.

MRI-Interactions with magnetically active and electrically conductive material

Patients undergoing MRI examinations are exposed to astrong static magnetic field and powerful electromagnetic alternating fields. Undesired or even dangerous effects could be caused if implants or objects with magnetic or electrically conductive elements are accidentally brought into the examination area. Relevant interactions in MRI between magnetic/electric fields and body tissue as well as foreign materials are systematically presented, based on proven physical principles. Natural components of the human body are mainly diamagnetic leading to only hardly perceptible magnetic forces in MRI. In contrast, ferromagnetic items as iron show translational forces of more than hundred times their weight force when brought to the entry of the bore. Lengthy ferromagnetic items are additionally subjected to torque. Materials with high electrical conductivity as metals and carbon fibre-reinforced plastic are also safety relevant. Especially long conductive structures as often present in implants are prone to induced strong electrical currents and high voltages at their end portions. Maximum voltages occurring at the implants and current density in adjacent tissue which might cause significant heating are hardly predictable for individual cases. Implants providing extended conductive loops for ring currents often show strong vibrations due to gradient switching. Counter forces must be considered when tilting conductive plates or ring structures inside the magnetic field area.

Imperfect magnetic field gradients in radial k-space encoding-Quantification, correction, and parameter dependency.

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. 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. The presented method allowed for echo-shift detection and correction of<1 data point, significantly improving the image quality in vitro and in vivo. Moreover, the method separated the effect of imbalanced gradients from those of magnetic inhomogeneities. The observed echo shifts were MR system-specifically dependent on acquisition parameters such as gradient strengths and dwell time. By acquiring 12 spokes for a certain set of acquisition parameters, the proposed method successfully corrects echo shift-related image artifacts independently of the MR system.

Technical Note: Interleaved bipolar acquisition and low-rank reconstruction for water-fat separation in MRI.

To acquire interleaved bipolar data and reconstruct the full data using low-rank property for water-fat separation. Bipolar acquisition suffers from issues related to gradient switching, the opposite gradient polarities, and other system imperfections, which prevent accurate water-fat separation. In this study, an interleaved bipolar acquisition scheme and a low-rank reconstruction method were proposed to reduce issues from the bipolar gradients while achieving a short imaging time. The proposed interleaved bipolar acquisition scheme collects echo-time signals from both gradient polarities; however, the sequence increases the imaging time. To reduce the imaging time, the signals were subsampled at every dimension of k-space. The low-rank property of the bipolar acquisition was defined and exploited to estimate the full data from the acquired subsampled data. To eliminate the bipolar issues, in the proposed method, the water-fat separation was performed separately for each gradient polarity, and the results for the positive and negative gradient polarities were combined after the water-fat separation. A phantom study and invivo experiments were conducted on a 3T Siemens Verio system. The results for the proposed method were compared with the results of the fully sampled interleaved bipolar acquisition and the Soliman's method, which was the previous water-fat separation approach for reducing the issues of bipolar gradients and accelerating the interleaved bipolar acquisition. The proposed method provided accurate water and fat images without the issues of bipolar gradients and demonstrated a better performance compared with the results of the previous methods. The water-fat separation using the bipolar acquisition has several benefits including a short echo-spacing time. However, it suffers from bipolar-gradient issues such as strong gradient switching, system imperfection, and eddy current effects. This study demonstrated that accurate water-fat separated images can be obtained using the proposed interleaved bipolar acquisition and low-rank reconstruction by using the benefits of the bipolar acquisition while reducing the bipolar-gradient issues with a short imaging time.