Referenceless interleaved echo-planar imaging

Unsupervised deep learning model for correcting Nyquist ghosts of single-shot spatiotemporal encoding

PURPOSE

To design an unsupervised deep learning (DL) model for correcting Nyquist ghosts of single-shot spatiotemporal encoding (SPEN) and evaluate the model for real MRI applications.

METHODS

The proposed method consists of three main components: (1) an unsupervised network that combines Residual Encoder and Restricted Subspace Mapping (RERSM-net) and is trained to generate a phase-difference map based on the even and odd SPEN images; (2) a spin physical forward model to obtain the corrected image with the learned phase difference map; and (3) cycle-consistency loss that is explored for training the RERSM-net.

RESULTS

The proposed RERSM-net could effectively generate smooth phase difference maps and correct Nyquist ghosts of single-shot SPEN. Both simulation and real in vivo MRI experiments demonstrated that our method outperforms the state-of-the-art SPEN Nyquist ghost correction method. Furthermore, the ablation experiments of generating phase-difference maps show the advantages of the proposed unsupervised model.

CONCLUSION

The proposed method can effectively correct Nyquist ghosts for the single-shot SPEN sequence.

Addressing concomitant gradient phase errors in time-interleaved chemical shift-encoded MRI fat fraction and R2 * mapping with a pass-specific phase fitting method

PURPOSE

Concomitant gradients induce phase errors that increase quadratically with distance from isocenter. This work proposes a complex-based fitting method that addresses concomitant gradient phase errors in chemical shift encoded (CSE) MRI estimation of proton density fat fraction (PDFF) and R₂ ^* through joint estimation of pass-specific phase terms. This method is applicable to time-interleaved multi-echo gradient-echo acquisitions (i.e., multi-pass acquisitions) and does not require prior knowledge of gradient waveforms typically needed to address concomitant gradient phase errors.

THEORY AND METHODS

A CSE-MRI spoiled gradient echo signal model, with pass-specific phase terms, is introduced for non-linear least squares estimation of PDFF and R₂ ^* in the presence of concomitant gradient phase errors. Cramér-Rao lower bound analysis was used to determine noise performance tradeoffs of the proposed fitting method, which was then validated in both phantom and in vivo experiments.

RESULTS

The proposed fitting method removed PDFF and R₂ ^* estimation errors up to 12% and 10 s^-1 , respectively, at ±12 cm off isocenter (S/I) in a water phantom. In healthy volunteers, PDFF and R₂ ^* bias was reduced by ~10% (12 cm off-isocenter) and ~30 s^-1 (16 cm off-isocenter), respectively. An evaluation in 29 clinical liver datasets demonstrated reduced PDFF bias and variability (8.4% improvement in the coefficient of variation), even with the imaging volume centered at isocenter.

CONCLUSION

Concomitant gradient induced phase errors in multi-pass CSE-MRI acquisitions can result in PDFF and R₂ ^* estimation biases away from isocenter. The proposed fitting method enables accurate PDFF and R₂ ^* quantification in the presence of concomitant gradient phase errors without knowledge of imaging gradient waveforms.

Accelerated k-space shift calibration for free-breathing stack-of-radial MRI quantification of liver fat and R 2 ∗

PURPOSE

To develop an accelerated k-space shift calibration method for free-breathing 3D stack-of-radial MRI quantification of liver proton-density fat fraction (PDFF) and R 2 ∗ .

METHODS

Accelerated k-space shift calibration was developed to partially skip acquisition of k-space shift data in the through-plane direction then interpolate in processing, as well as to reduce the in-plane averages. A multi-echo stack-of-radial sequence with the baseline calibration was evaluated on a phantom versus vendor-provided reference-standard PDFF and R 2 ∗ values at 1.5T, and in 13 healthy subjects and 5 clinical subjects at 3T with respect to reference-standard breath-hold Cartesian acquisitions. PDFF and R 2 ∗ maps were calculated with different calibration acceleration factors offline and compared to reference-standard values using Bland-Altman analysis. Bias and uncertainty were evaluated using normal distribution and Bayesian probability of difference (P < .05 considered significant).

RESULTS

Bland-Altman plots of phantom and in vivo data showed that substantial acceleration was highly feasible in both through-plane and in-plane directions. Compared to the baseline calibration without acceleration, Bayesian analysis revealed no significant differences on biases and uncertainties of PDFF and R 2 ∗ measurements with all acceleration methods in this study, except the method with through-plane acceleration equaling slices and averages equaling 20 for PDFF and R 2 ∗ (both P < .001) for the phantom. A six-fold reduction in equivalent calibration acquisition time (time saving ≥25 s and ≥80.7%) was achieved using recommended acceleration factors for the in vivo protocols in this study.

CONCLUSION

This proposed method may allow accelerated calibration for free-breathing stack-of-radial MRI PDFF and R 2 ∗ mapping.

Robust autocalibrated structured low-rank EPI ghost correction

PURPOSE

We propose and evaluate a new structured low-rank method for echo-planar imaging (EPI) ghost correction called Robust Autocalibrated LORAKS (RAC-LORAKS). The method can be used to suppress EPI ghosts arising from the differences between different readout gradient polarities and/or the differences between different shots. It does not require conventional EPI navigator signals, and is robust to imperfect autocalibration data.

METHODS

Autocalibrated LORAKS is a previous structured low-rank method for EPI ghost correction that uses GRAPPA-type autocalibration data to enable high-quality ghost correction. This method works well when the autocalibration data are pristine, but performance degrades substantially when the autocalibration information is imperfect. RAC-LORAKS generalizes Autocalibrated LORAKS in two ways. First, it does not completely trust the information from autocalibration data, and instead considers the autocalibration and EPI data simultaneously when estimating low-rank matrix structure. Second, it uses complementary information from the autocalibration data to improve EPI reconstruction in a multi-contrast joint reconstruction framework. RAC-LORAKS is evaluated using simulations and in vivo data, including comparisons to state-of-the-art methods.

RESULTS

RAC-LORAKS is demonstrated to have good ghost elimination performance compared to state-of-the-art methods in several complicated EPI acquisition scenarios (including gradient-echo brain imaging, diffusion-encoded brain imaging, and cardiac imaging).

CONCLUSIONS

RAC-LORAKS provides effective suppression of EPI ghosts and is robust to imperfect autocalibration data.

Prospective motion correction for diffusion weighted EPI of the brain using an optical markerless tracker

PURPOSE

To enable motion-robust diffusion weighted imaging of the brain using well-established imaging techniques.

METHODS

An optical markerless tracking system was used to estimate and correct for rigid body motion of the head in real time during scanning. The imaging coordinate system was updated before each excitation pulse in a single-shot EPI sequence accelerated by GRAPPA with motion-robust calibration. Full Fourier imaging was used to reduce effects of motion during diffusion encoding. Subjects were imaged while performing prescribed motion patterns, each repeated with prospective motion correction on and off.

RESULTS

Prospective motion correction with dynamic ghost correction enabled high quality DWI in the presence of fast and continuous motion within a 10° range. Images acquired without motion were not degraded by the prospective correction. Calculated diffusion tensors tolerated the motion well, but ADC values were slightly increased.

CONCLUSIONS

Prospective correction by markerless optical tracking minimizes patient interaction and appears to be well suited for EPI-based DWI of patient groups unable to remain still including those who are not compliant with markers.

Dynamic distortion correction for functional MRI using FID navigators

PURPOSE

To develop a method for slice-wise dynamic distortion correction for EPI using rapid spatiotemporal B₀ field measurements from FID navigators (FIDnavs) and to evaluate the efficacy of this new approach relative to an established data-driven technique.

METHODS

A low-resolution reference image was used to create a forward model of FIDnav signal changes to enable estimation of spatiotemporal B₀ inhomogeneity variations up to second order from measured FIDnavs. Five volunteers were scanned at 3 T using a 64-channel coil with FID-navigated EPI. The accuracy of voxel shift measurements and geometric distortion correction was assessed for experimentally induced magnetic field perturbations. The temporal SNR was evaluated in EPI time-series acquired at rest and with a continuous nose-touching action, before and after image realignment.

RESULTS

Field inhomogeneity coefficients and voxel shift maps measured using FIDnavs were in excellent agreement with multi-echo EPI measurements. The FID-navigated distortion correction accurately corrected image geometry in the presence of induced magnetic field perturbations, outperforming the data-driven approach in regions with large field offsets. In functional MRI scans with nose touching, FIDnav-based correction yielded temporal SNR gains of 30% in gray matter. Following image realignment, which accounted for global image shifts, temporal SNR gains of 3% were achieved.

CONCLUSIONS

Our proposed application of FIDnavs enables slice-wise dynamic distortion correction with high temporal efficiency. We achieved improved signal stability by leveraging the encoding information from multichannel coils. This approach can be easily adapted to other EPI-based sequences to improve temporal SNR for a variety of clinical and research applications.

Effect of respiratory motion on free-breathing 3D stack-of-radial liver R 2 ∗ relaxometry and improved quantification accuracy using self-gating

PURPOSE

To develop an accurate free-breathing 3D liver R 2 ∗ mapping approach and to evaluate it in vivo.

METHODS

A free-breathing multi-echo stack-of-radial sequence was applied in 5 normal subjects and 6 patients at 3 Tesla. Respiratory motion compensation was implemented using the inherent self-gating signal. A breath-hold Cartesian acquisition was the reference standard. Proton density fat fraction and R 2 ∗ were measured and compared between radial and Cartesian methods using Bland-Altman plots. The normal subject results were fitted to a linear mixed model (P < .05 considered significant).

RESULTS

Free-breathing stack-of-radial without self-gating exhibited signal attenuation in echo images and artifactually elevated apparent R 2 ∗ values. In the Bland-Altman plots of normal subjects, compared to breath-hold Cartesian, free-breathing stack-of-radial acquisitions of 22, 30, 36, and 44 slices, had mean R 2 ∗ differences of 27.4, 19.4, 10.9, and 14.7 s^-1 with 800 radial views, and they had 18.4, 11.9, 9.7, and 27.7 s^-1 with 404 views, which were reduced to 0.4, 0.9, -0.2, and -0.7 s^-1 and to -1.7, -1.9, -2.1, and 0.5 s^-1 with self-gating, respectively. No substantial proton density fat fraction differences were found. The linear mixed model showed free-breathing radial R 2 ∗ results without self-gating were significantly biased by 17.2 s^-1 averagely (P = .002), which was eliminated with self-gating (P = .930). Proton density fat fraction results were not different (P > .234). For patients, Bland-Altman plots exhibited mean R 2 ∗ differences of 14.4 and 0.1 s^-1 for free-breathing stack-of-radial without self-gating and with self-gating, respectively, but no substantial proton density fat fraction differences.

CONCLUSION

The proposed self-gating method corrects the respiratory motion bias and enables accurate free-breathing stack-of-radial quantification of liver R 2 ∗ .

A comparison of static and dynamic ∆B0 mapping methods for correction of CEST MRI in the presence of temporal B0 field variations

PURPOSE

To assess the performance, in the presence of scanner instabilities, of three dynamic correction methods which integrate ∆B₀ mapping into the chemical exchange saturation transfer (CEST) measurement and three established static ∆B₀ -correction approaches.

METHODS

A homogeneous phantom and five healthy volunteers were scanned with a CEST sequence at 7 T. The in vivo measurements were performed twice: first with unaltered system frequency and again applying frequency shifts during the CEST acquisition. In all cases, retrospective voxel-wise ∆B₀ -correction was performed using one intrinsic and two extrinsic [prescans with dual-echo gradient-echo and water saturation shift referencing (WASSR)] static approaches. These were compared with two intrinsic [using phase data directly generated by single-echo or double-echo GRE (gradient-echo) CEST readout (CEST-GRE-2TE)] and one extrinsic [phase from interleaved dual-echo EPI (echo planar imaging) navigator (NAV-EPI-2TE)] dynamic ∆B₀ -correction approaches [allowing correction of each Z-spectral point before magnetization transfer ratio asymmetry (MTR_asym) analysis].

RESULTS

All three dynamic methods successfully mapped the induced drift. The intrinsic approaches were affected by the CEST labeling near water (∆ω < |0.3| ppm). The MTR_asym contrast was distorted by the frequency drift in the brain by up to 0.21%/Hz when static ∆B₀ -corrections were applied, whereas the dynamic ∆B₀ corrections reduced this to <0.01%/Hz without the need of external scans. The CEST-GRE-2TE and NAV-EPI-2TE resulted in highly consistent MTR_asym values with/without drift for all subjects.

CONCLUSION

Reliable correction of scanner instabilities is essential to establish clinical CEST MRI. The three dynamic approaches presented improved the ∆B₀ -correction performance significantly in the presence of frequency drift compared to established static methods. Among them, the self-corrected CEST-GRE-2TE was the most accurate and straightforward to implement.

A Wireless Radio Frequency Triggered Acquisition Device (WRAD) for Self-Synchronised Measurements of the Rate of Change of the MRI Gradient Vector Field for Motion Tracking

In this paper, we present a device that is capable of wireless synchronization to the MRI pulse sequence time frame with sub-microsecond precision. This is achieved by detecting radio frequency pulses in the parent pulse sequence using a small resonant circuit. The device incorporates a 3-axis pickup coil, constructed using conventional printed circuit board (PCB) manufacturing techniques, to measure the rate of change of the gradient waveforms with respect to time. Using Maxwell's equations, assuming negligible rates of change of curl and divergence, a model of the expected gradient derivative (slew) vector field is presented. A 3-axis Hall effect magnetometer allows for the measurement of the direction of the static magnetic field in the device co-ordinate frame. By combining the magnetometer measurement with the pickup coil voltages and slew vector field model, the orientation and position can be determined to within a precision of 0.1 degrees and 0.1 mm, respectively, using a pulse series lasting 880 μs . The gradient pulses are designed to be sinusoidal, enabling the detection of a phase shift between the time frame of the pickup coil digitization circuit and the gradient amplifiers. The signal processing is performed by a low power micro-controller on the device and the results are transmitted out of the scanner bore using a low latency 2.4 GHz radio link. The device identified an unexpected 40 kHz oscillation relating to the pulse width modulation frequency of the gradient amplifiers that is predominantly in the direction of the static magnetic field. The proposed wireless radio frequency triggered acquisition device enables users to probe the scanner gradient slew vector field with minimal hardware set-up and shows promise for the future developments in the prospective motion correction.

Real-Time Filtering with Sparse Variations for Head Motion in Magnetic Resonance Imaging

Estimating a time-varying signal, such as head motion from magnetic resonance imaging data, becomes particularly challenging in the face of other temporal dynamics such as functional activation. This paper describes a new Kalman filter-like framework that includes a sparse residual term in the measurement model. This additional term allows the extended Kalman filter to generate real-time motion estimates suitable for prospective motion correction when such dynamics occur. An iterative augmented Lagrangian algorithm similar to the alterating direction method of multipliers implements the update step for this Kalman filter. This paper evaluates the accuracy and convergence rate of this iterative method for small and large motion in terms of its sensitivity to parameter selection. The included experiment on a simulated functional magnetic resonance imaging acquisition demonstrates that the resulting method improves the maximum Youden's J index of the time series analysis by 2-3% versus retrospective motion correction, while the sensitivity index increases from 4.3 to 5.4 when combining prospective and retrospective correction.

Navigator-Free EPI Ghost Correction With Structured Low-Rank Matrix Models: New Theory and Methods

Structured low-rank matrix models have previously been introduced to enable calibrationless MR image reconstruction from sub-Nyquist data, and such ideas have recently been extended to enable navigator-free echo-planar imaging (EPI) ghost correction. This paper presents a novel theoretical analysis which shows that, because of uniform subsampling, the structured low-rank matrix optimization problems for EPI data will always have either undesirable or non-unique solutions in the absence of additional constraints. This theory leads us to recommend and investigate problem formulations for navigator-free EPI that incorporate side information from either image-domain or k-space domain parallel imaging methods. The importance of using nonconvex low-rank matrix regularization is also identified. We demonstrate using phantom and in vivo data that the proposed methods are able to eliminate ghost artifacts for several navigator-free EPI acquisition schemes, obtaining better performance in comparison with the state-of-the-art methods across a range of different scenarios. Results are shown for both single-channel acquisition and highly accelerated multi-channel acquisition.

Free-breathing liver fat quantification using a multiecho 3D stack-of-radial technique

PURPOSE

The diagnostic gold standard for nonalcoholic fatty liver disease is an invasive biopsy. Noninvasive Cartesian MRI fat quantification remains limited to a breath-hold (BH). In this work, a novel free-breathing 3D stack-of-radial (FB radial) liver fat quantification technique is developed and evaluated in a preliminary study.

METHODS

Phantoms and healthy subjects (n = 11) were imaged at 3 Tesla. The proton-density fat fraction (PDFF) determined using FB radial (with and without scan acceleration) was compared to BH single-voxel MR spectroscopy (SVS) and BH 3D Cartesian MRI using linear regression (correlation coefficient ρ and concordance coefficient ρ_c ) and Bland-Altman analysis.

RESULTS

In phantoms, PDFF showed significant correlation (ρ > 0.998, ρ_c  > 0.995) and absolute mean differences < 2.2% between FB radial and BH SVS, as well as significant correlation (ρ > 0.999, ρ_c  > 0.998) and absolute mean differences < 0.6% between FB radial and BH Cartesian. In the liver and abdomen, PDFF showed significant correlation (ρ > 0.986, ρ_c  > 0.985) and absolute mean differences < 1% between FB radial and BH SVS, as well as significant correlation (ρ > 0.996, ρ_c  > 0.995) and absolute mean differences < 0.9% between FB radial and BH Cartesian.

CONCLUSION

Accurate 3D liver fat quantification can be performed in 1 to 2 min using a novel FB radial technique. Magn Reson Med 79:370-382, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Accelerated mapping of magnetic susceptibility using 3D planes-on-a-paddlewheel (POP) EPI at ultra-high field strength

With the advent of ultra-high field MRI scanners in clinical research, susceptibility based MRI has recently gained increasing interest because of its potential to assess subtle tissue changes underlying neurological pathologies/disorders. Conventional, but rather slow, three-dimensional (3D) spoiled gradient-echo (GRE) sequences are typically employed to assess the susceptibility of tissue. 3D echo-planar imaging (EPI) represents a fast alternative but generally comes with echo-time restrictions, geometrical distortions and signal dropouts that can become severe at ultra-high fields. In this work we assess quantitative susceptibility mapping (QSM) at 7 T using non-Cartesian 3D EPI with a planes-on-a-paddlewheel (POP) trajectory, which is created by rotating a standard EPI readout train around its own phase encoding axis. We show that the threefold accelerated non-Cartesian 3D POP EPI sequence enables very fast, whole brain susceptibility mapping at an isotropic resolution of 1 mm and that the high image quality has sufficient signal-to-noise ratio in the phase data for reliable QSM processing. The susceptibility maps obtained were comparable with regard to QSM values and geometric distortions to those calculated from a conventional 4 min 3D GRE scan using the same QSM processing pipeline. Copyright © 2016 John Wiley & Sons, Ltd.

New method to characterize and correct with sub-µs precision gradient delays in bipolar multispoke RF pulses

PURPOSE

Small gradient delays with respect to radiofrequency (RF) events can have disastrous effects on the performance of bipolar spokes RF pulses used in parallel transmission (pTx). In this work, we propose a new method to characterize and correct this delay with sub-µs precision.

METHODS

By determining experimentally the phase Δφ producing a 0 ° flip angle excitation in a α0°-α180°+Δφ bipolar two-spoke pulse configuration at multiple slice locations, we demonstrate the possibility of deducing the underlying gradient delay with precision. The technique also suggests prospectively compensating for the same delay by altering the phase of the second pulse. The approach was tested with a multislice gradient echo sequence on a phantom and on one healthy volunteer at 7 Tesla.

RESULTS

Application of the method returned an accuracy of approximately 50 ns on the gradient delay measurement, a performance shown in fact to be desirable for high-performance pTx 2D applications. Phase corrections of up to 180 ° on the second spoke RF pulse in the bipolar configuration allowed us to obtain similar performance as for unipolar designs, yet with significantly shorter excitations.

CONCLUSIONS

A simple and accurate gradient-delay calibration method was proposed that offers the possibility of using bipolar multispoke pulses in multislice protocols. Magn Reson Med 78:2194-2202, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Referenceless one-dimensional Nyquist ghost correction in multicoil single-shot spatiotemporally encoded MRI

Single-shot spatiotemporally encoded (SPEN) MRI is a novel fast imaging method capable of retaining the time efficiency of single-shot echo planar imaging (EPI) but with distortion artifacts significantly reduced. Akin to EPI, the phase inconsistencies between mismatched even and odd echoes also result in the so-called Nyquist ghosts. However, the characteristic of the SPEN signals provides the possibility of obtaining ghost-free images directly from even and odd echoes respectively, without acquiring additional reference scans. In this paper, a theoretical analysis of the Nyquist ghosts manifested in single-shot SPEN MRI is presented, a one-dimensional correction scheme is put forward capable of maintaining definition of image features without blurring when the phase inconsistency along SPEN encoding direction is negligible, and a technique is introduced for convenient and robust correction of data from multi-channel receiver coils. The effectiveness of the proposed processing pipeline is validated by a series of experiments conducted on simulation data, in vivo rats and healthy human brains. The robustness of the method is further verified by implementing distortion correction on ghost corrected data.

K-space trajectory mapping and its application for ultrashort Echo time imaging

MR images are affected by system delays and gradient field imperfections which induce discrepancies between prescribed and actual k-space trajectories. This could be even more critical for non-Cartesian data acquisitions where even a small deviation from the assumed k-space trajectory results in severe image degradation and artifacts. Knowledge of the actual k-space trajectories is therefore crucial and can be incorporated in the reconstruction of high quality non-Cartesian images. A novel MR method for the calibration of actual gradient waveforms was developed using a combination of phase encoding increments and subsequent detection of the exact time point at which the corresponding trajectory is crossing the k-space origin. The measured sets of points were fitted to a parametrical model to calculate the complete actual acquisition trajectory. Measurements performed on phantoms and volunteers, positioned both in- and off-isocenter of the magnet, clearly demonstrate the improvement in reconstructed ultrashort echo time (UTE) images, when information from calibration of k-space sampling trajectories is employed in the MR image reconstruction procedure. The unique feature of the proposed method is its robustness and simple experimental setup, making it suitable for quick acquisition trajectory calibration procedures e.g. for non-Cartesian radial fast imaging.

Acquiring and processing ultrafast biomolecular 2D NMR experiments using a referenced-based correction

Thanks to their special spatiotemporal encoding/decoding scheme, ultrafast (UF) NMR sequences can deliver arbitrary 2D spectra following a single excitation. Regardless of their nature, these sequences have in common their tracing of a path in the [Formula: see text]-[Formula: see text] plane, that will deliver the spectrum being sought after a 1D Fourier transformation versus [Formula: see text]. This need to simultaneously digitize two domains, tends to impose bandwidth limitations along all spectral axes. Along the [Formula: see text]/[Formula: see text] dimension this problem is exacerbated by the fact that odd and even time points are not equispaced, and by additional artifacts such as time shifts between time points sampled while under the action of positive and negative decoding gradients. As a result, odd and even [Formula: see text] points are typically Fourier transformed separately, halving the potential spectral width along this dimension. While this halving of the [Formula: see text] span can be overcome by an interlaced Fourier transform, this post-processing is seldom used because of its sensitivity to hardware inaccuracies requiring even finer corrections of the even/odd [Formula: see text] data points. These corrections have so far been done manually, but are challenging to implement when dealing with low signal-to-noise ratio signals like those associated with biomolecular NMR experiments. This study introduces an algorithm for an automatic correction of all even/odd ultrafast NMR inconsistencies, based on the acquisition of a reference scan on the solvent. This algorithm was verified experimentally using an [Formula: see text]-[Formula: see text] UF-HSQC variant on ubiquitin at 600 MHz. Features of this method as well as of the interlaced Fourier transformation in general, are discussed.

A method for the dynamic correction of B0-related distortions in single-echo EPI at 7T

We propose a method to calculate field maps from the phase of each EPI in an fMRI time series. These field maps can be used to correct the corresponding magnitude images for distortion caused by inhomogeneity in the static magnetic field. In contrast to conventional static distortion correction, in which one ‘snapshot’ field map is applied to all subsequent fMRI time points, our method also captures dynamic changes to B₀ which arise due to motion and respiration. The approach is based on the assumption that the non-B₀-related contribution to the phase measured by each radio-frequency coil, which is dominated by the coil sensitivity, is stable over time and can therefore be removed to yield a field map from EPI.

Our solution addresses imaging with multi-channel coils at ultra-high field (7 T), where phase offsets vary rapidly in space, phase processing is non-trivial and distortions are comparatively large. We propose using dual-echo gradient echo reference scan for the phase offset calculation, which yields estimates with high signal-to-noise ratio. An extrapolation method is proposed which yields reliable estimates for phase offsets even where motion is large and a tailored phase unwrapping procedure for EPI is suggested which gives robust results in regions with disconnected tissue or strong signal decay.

Phase offsets are shown to be stable during long measurements (40 min) and for large head motions. The dynamic distortion correction proposed here is found to work accurately in the presence of large motion (up to 8.1°), whereas a conventional method based on single field map fails to correct or even introduces distortions (up to 11.2 mm). Finally, we show that dynamic unwarping increases the temporal stability of EPI in the presence of motion.

Our approach can be applied to any EPI measurements without the need for sequence modification.

Joint correction of Nyquist artifact and minuscule motion-induced aliasing artifact in interleaved diffusion weighted EPI data using a composite two-dimensional phase correction procedure

Diffusion-weighted imaging (DWI) obtained with interleaved echo-planar imaging (EPI) pulse sequence has great potential of characterizing brain tissue properties at high spatial-resolution. However, interleaved EPI based DWI data may be corrupted by various types of aliasing artifacts. First, inconsistencies in k-space data obtained with opposite readout gradient polarities result in Nyquist artifact, which is usually reduced with 1D phase correction in post-processing. When there exist eddy current cross terms (e.g., in oblique-plane EPI), 2D phase correction is needed to effectively reduce Nyquist artifact. Second, minuscule motion induced phase inconsistencies in interleaved DWI scans result in image-domain aliasing artifact, which can be removed with reconstruction procedures that take shot-to-shot phase variations into consideration. In existing interleaved DWI reconstruction procedures, Nyquist artifact and minuscule motion-induced aliasing artifact are typically removed subsequently in two stages. Although the two-stage phase correction generally performs well for non-oblique plane EPI data obtained from well-calibrated system, the residual artifacts may still be pronounced in oblique-plane EPI data or when there exist eddy current cross terms. To address this challenge, here we report a new composite 2D phase correction procedure, which effective removes Nyquist artifact and minuscule motion induced aliasing artifact jointly in a single step. Our experimental results demonstrate that the new 2D phase correction method can much more effectively reduce artifacts in interleaved EPI based DWI data as compared with the existing two-stage artifact correction procedures. The new method robustly enables high-resolution DWI, and should prove highly valuable for clinical uses and research studies of DWI.

Two-dimensional-NGC-SENSE-GRAPPA for fast, ghosting-robust reconstruction of in-plane and slice-accelerated blipped-CAIPI echo planar imaging

Purpose

Ghosting‐robust reconstruction of blipped‐CAIPI echo planar imaging simultaneous multislice data with low computational load.

Methods

To date, Slice‐GRAPPA, with “odd–even” kernels that improve ghosting performance, has been the framework of choice for such reconstructions due to its predecessor SENSE‐GRAPPA being deemed unsuitable for blipped‐CAIPI data. Modifications to SENSE‐GRAPPA are used to restore CAIPI compatibility and to make it robust against ghosting. Two implementations are tested, one where slices and in‐plane unaliasing are dealt in the same serial manner as in Slice‐GRAPPA [referred to as one‐dimensional (1D)‐NGC‐SENSE‐GRAPPA, where NGC stands for Nyquist Ghost Corrected] and one where both are unaliased in a single step (2D‐NGC‐SENSE‐GRAPPA), which is analytically and experimentally shown to be computationally cheaper.

Results

The 1D‐NGC‐SENSE‐GRAPPA and odd‐even Slice‐GRAPPA perform identically, whereas 2D‐NGC‐SENSE‐GRAPPA shows reduced error propagation, less residual ghosting when reliable reference data were available. When the latter was not the case, error propagation was increased.

Conclusion

Unlike Slice‐GRAPPA, SENSE‐GRAPPA operates fully within the GRAPPA framework, for which improved reconstructions (e.g., iterative, nonlinear) have been developed over the past decade. It could, therefore, bring benefit to the reconstruction of SMS data as an attractive alternative to Slice‐GRAPPA. Magn Reson Med 77:998–1009, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Foundations of MRI phase imaging and processing for Quantitative Susceptibility Mapping (QSM)

Quantitative Susceptibility Mapping (QSM) is a novel MRI based technique that relies on estimates of the magnetic field distribution in the tissue under examination. Several sophisticated data processing steps are required to extract the magnetic field distribution from raw MRI phase measurements. The objective of this review article is to provide a general overview and to discuss several underlying assumptions and limitations of the pre-processing steps that need to be applied to MRI phase data before the final field-to-source inversion can be performed. Beginning with the fundamental relation between MRI signal and tissue magnetic susceptibility this review covers the reconstruction of magnetic field maps from multi-channel phase images, background field correction, and provides an overview of state of the art QSM solution strategies.

Correcting dynamic distortions in 7T echo planar imaging using a jittered echo time sequence

Purpose

To develop a distortion correction method for echo planar imaging (EPI) that is able to measure dynamic changes in B₀.

Theory and Methods

The approach we propose is based on single‐echo EPI with a jittering of the echo time between two values for alternate time points. Field maps are calculated between phase images from adjacent volumes and are used to remove distortion from corresponding magnitude images. The performance of our approach was optimized using an analytical model and by comparison with field maps from dual‐echo EPI. The method was tested in functional MRI experiments at 7T with motor tasks and compared with the conventional static approach.

Results

Unwarping using our method was accurate even for head rotations up to 8.2°, where the static approach introduced errors up to 8.2 mm. Jittering the echo time between 19 and 25 ms had no measurable effect on blood oxygenation level–dependent (BOLD) sensitivity. Our approach reduced the distortions in activated regions to <1 mm and repositioned active voxels correctly.

Conclusion

This method yields accurate distortion correction in the presence of motion. No reduction in BOLD sensitivity was observed. As such, it is suitable for application in a wide range of functional MRI experiments. Magn Reson Med 76:1388–1399, 2016. © 2015 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

New spatiotemporal approaches for fully refocused, multislice ultrafast 2D MRI

PURPOSE

Single-scan multislice acquisition schemes play key roles in magnetic resonance imaging. Central among these "ultrafast" experiments stands echo-planar imaging, a technique that although of optimal sampling is challenged by T2* artifacts. Recent studies described alternatives based on spatiotemporal encoding (SPEN), which are particularly robust if implemented in a "full-refocusing" mode. This work extends this modality from the single-slice acquisitions in which it has hitherto been implemented, by introducing a variety of multislice schemes scanning 3D volumes.

METHODS

Multislice SPEN employing either inversion or stimulated echo pulses and timed to fulfill the demands of full refocusing, are demonstrated. The performance of the ensuing methods was examined in "Hybrid" modalities encoding data in k- and direct-space, in low-specific absorption rate stimulated-echo approaches, and in direct-space SPEN approaches.

RESULTS

When applied in phantoms and in in vivo systems, the ensuing single-shot sequences evidenced similar robustness, sensitivity, and resolution qualities as previously discussed 2D single-slice schemes, while enabling a rapid sampling of the third dimension via multislicing.

CONCLUSION

The unique benefits deriving from fully refocused, multislice, single-scan SPEN sequences were corroborated by phantom tests, as well as by in vivo scans at 3 and 7 T. Low specific absorption rate multislice SPEN variants compatible with human studies were demonstrated.

Referenceless reconstruction of spatiotemporally encoded imaging data: principles and applications to real-time MRI

PURPOSE

Ultrafast sequences based on "Hybrid" spatiotemporal encoding (SPEN) replace echo-planar imaging's phase encoding "blips," while retaining a k-space readout acquisition. Hardware imperfections during acquisition may lead to ghosts and striped artifacts along the SPEN dimension; akin to echo-planar imaging's Nyquist ghosts, but weaker. A referenceless method to eliminate these artifacts in Hybrid SPEN is demonstrated.

THEORY AND METHODS

Owing to its encoding in direct space, rather than reciprocal space, undersampling in SPEN does not generate an echo-planar-imaging-like aliasing, but instead lowers the spatial resolution. Hybrid SPEN data can be split into two undersampled signals: a reference one comprised of the odd-echos, and an even-echo set that has to be "corrected" for consistency with the former. A simple way of implementing such a correction that enables a joint high-resolution reconstruction is proposed.

RESULTS

The referenceless algorithm is demonstrated with various examples, including oblique scans, large in vivo datasets from real-time dynamic contrast-enhanced perfusion experiments, and human brain imaging.

CONCLUSIONS

The referenceless correction enables robust single-scan imaging under changing conditions-such as patient motion and changes in shimming over time-without the need of ancillary navigators. This opens new options for real-time MRI and interactive scanning.

Turboprop+: enhanced Turboprop diffusion-weighted imaging with a new phase correction

Faster periodically rotated overlapping parallel lines with enhanced reconstruction (PROPELLER) diffusion-weighted imaging acquisitions, such as Turboprop and X-prop, remain subject to phase errors inherent to a gradient echo readout, which ultimately limits the applied turbo factor (number of gradient echoes between each pair of radiofrequency refocusing pulses) and, thus, scan time reductions. This study introduces a new phase correction to Turboprop, called Turboprop+. This technique employs calibration blades, which generate 2-D phase error maps and are rotated in accordance with the data blades, to correct phase errors arising from off-resonance and system imperfections. The results demonstrate that with a small increase in scan time for collecting calibration blades, Turboprop+ had a superior immunity to the off-resonance-related artifacts when compared to standard Turboprop and recently proposed X-prop with the high turbo factor (turbo factor = 7). Thus, low specific absorption rate and short scan time can be achieved in Turboprop+ using a high turbo factor, whereas off-resonance related artifacts are minimized.

Ultra-high resolution diffusion tensor imaging of the microscopic pathways of the medial temporal lobe

Diseases involving the medial temporal lobes (MTL) such as Alzheimer's disease and mesial temporal sclerosis pose an ongoing diagnostic challenge because of the difficulty in identifying conclusive imaging features, particularly in pre-clinical states. Abnormal neuronal connectivity may be present in the circuitry of the MTL, but current techniques cannot reliably detect those abnormalities. Diffusion tensor imaging (DTI) has shown promise in defining putative abnormalities in connectivity, but DTI studies of the MTL performed to date have shown neither dramatic nor consistent differences across patient populations. Conventional DTI methodology provides an inadequate depiction of the complex microanatomy present in the medial temporal lobe because of a typically employed low isotropic resolution of 2.0-2.5 mm, a low signal-to-noise ratio (SNR), and echo-planar imaging (EPI) geometric distortions that are exacerbated by the inhomogeneous magnetic environment at the skull base. In this study, we pushed the resolving power of DTI to near-mm isotropic voxel size to achieve a detailed depiction of mesial temporal microstructure at 3 T. High image fidelity and SNR at this resolution are achieved through several mechanisms: (1) acquiring multiple repetitions of the minimum field of view required for hippocampal coverage to boost SNR; (2) utilizing a single-refocused diffusion preparation to enhance SNR further; (3) performing a phase correction to reduce Rician noise; (4) minimizing distortion and maintaining left-right distortion symmetry with axial-plane parallel imaging; and (5) retaining anatomical and quantitative accuracy through the use of motion correction coupled with a higher-order eddy-current correction scheme. We combined this high-resolution methodology with a detailed segmentation of the MTL to identify tracks in all subjects that may represent the major pathways of the MTL, including the perforant pathway. Tractography performed on a subset of the data identified similar tracks, although they were lesser in number. This detailed analysis of MTL substructure may have applications to clinical populations.

Two-dimensional phase cycled reconstruction for inherent correction of echo-planar imaging Nyquist artifacts

The inconsistency of k-space trajectories results in Nyquist artifacts in echo-planar imaging (EPI). Traditional techniques often only correct for phase errors along the frequency-encoding direction (one-dimensional correction), which may leave significant residual artifacts, particularly for oblique-plane EPI or in the presence of cross-term eddy currents. As compared with one-dimensional correction, two-dimensional (2D) phase correction can be much more effective in suppressing Nyquist artifacts. However, most existing 2D correction methods require reference scans and may not be generally applicable to different imaging protocols. Furthermore, EPI reconstruction with these 2D phase correction methods is susceptible to error amplification due to subject motion. To address these limitations, we report an inherent and general 2D phase correction technique for EPI Nyquist removal. First, a series of images are generated from the original dataset, by cycling through different possible values of phase errors using a 2D reconstruction framework. Second, the image with the lowest artifact level is identified from images generated in the first step using criteria based on background energy in sorted and sigmoid-weighted signals. In this report, we demonstrate the effectiveness of our new method in removing Nyquist ghosts in single-shot, segmented and parallel EPI without acquiring additional reference scans and the subsequent error amplifications.

Single breathhold cardiac CINE imaging with multi-echo three-dimensional hybrid radial SSFP acquisition

PURPOSE

To achieve single breathhold whole heart cardiac CINE imaging with improved spatial resolution and temporal resolution by using a multi-echo three-dimensional (3D) hybrid radial SSFP acquisition.

MATERIALS AND METHODS

Multi-echo 3D hybrid radial SSFP acquisitions were used to acquire cardiac CINE imaging within a single breathhold. An optimized interleaving scheme was developed for view ordering throughout the cardiac cycle.

RESULTS

Whole heart short axis views were acquired with a spatial resolution of 1.3 x 1.3 x 8.0 mm(3) and temporal resolution of 45 ms, within a single 17 s breathhold. The technique was validated on eight healthy volunteers by measuring the left ventricular volume throughout the cardiac cycle and comparing with the conventional 2D multiple breathhold technique. The left ventricle functional measurement bias of our proposed 3D technique from the conventional 2D technique: end diastolic volume -3.3 mL +/- 13.7 mL, end systolic volume 1.4 mL +/- 6.1 mL, and ejection fraction -1.7% +/- 4.3%, with high correlations 0.94, 0.97, and 0.91, accordingly.

CONCLUSION

A multi-echo 3D hybrid radial SSFP acquisition was developed to allow for a whole heart cardiac CINE exam in a single breathhold. Cardiac function measurements in volunteers compared favorably with the standard multiple breathhold exams.

Robust 2D phase correction for echo planar imaging under a tight field-of-view

Nyquist ghost artifacts are a serious issue in echo planar imaging. These artifacts primarily originate from phase difference between even and odd echo images and can be removed or reduced using phase correction methods. The commonly used 1D phase correction can only correct phase difference along readout axis. 2D correction is, therefore, necessary when phase difference presents along both readout and phase encoding axes. However, existing 2D methods have several unaddressed issues that affect their practicality. These issues include uncharacterized noise behavior, image artifact due to unoptimized phase estimation, Gibbs ringing artifact when directly applying to partial k(y) data, and most seriously a new image artifact under tight field-of-view (i.e., field-of-view slightly smaller than object size). All these issues are addressed in this article. Specifically, theoretical analysis of noise amplification and effect of phase estimation error is provided, and tradeoff between noise and ghost is studied. A new 2D phase correction method with improved polynomial fitting, joint homodyne processing and phase correction, compatibility with tight field-of-view is then proposed. Various results show that the proposed method can robustly generate images free of Nyquist ghosts and other image artifacts even in oblique scans or when cross-term eddy current terms are significant.

Phase and amplitude correction for multi-echo water-fat separation with bipolar acquisitions

PURPOSE

To address phase and amplitude errors for multi-point water-fat separation with "bipolar" acquisitions, which efficiently collect all echoes with alternating read-out gradient polarities in one repetition.

MATERIALS AND METHODS

With the bipolar acquisitions, eddy currents and other system nonidealities can induce inconsistent phase errors between echoes, disrupting water-fat separation. Previous studies have addressed phase correction in the read-out direction. However, the bipolar acquisitions may be subject to spatially high order phase errors as well as an amplitude modulation in the read-out direction. A method to correct for the 2D phase and amplitude errors is introduced. Low resolution reference data with reversed gradient polarities are collected. From the pair of low-resolution data collected with opposite gradient polarities, the two-dimensional phase errors are estimated and corrected. The pair of data are then combined for water-fat separation.

RESULTS

We demonstrate that the proposed method can effectively remove the high order errors with phantom and in vivo experiments, including obliquely oriented scans.

CONCLUSION

For bipolar multi-echo acquisitions, uniform water-fat separation can be achieved by removing high order phase errors with the proposed method.

Robust EPI Nyquist ghost elimination via spatial and temporal encoding

Nyquist ghosts are an inherent artifact in echo planar imaging acquisitions. An approach to robustly eliminate Nyquist ghosts is presented that integrates two previous Nyquist ghost correction techniques: temporal domain encoding (phase labeling for additional coordinate encoding: PLACE and spatial domain encoding (phased array ghost elimination: PAGE). Temporal encoding modulates the echo planar imaging acquisition trajectory from frame to frame, enabling one to interleave data to remove inconsistencies that occur between sampling on positive and negative gradient readouts. With PLACE, one can coherently combine the interleaved data to cancel residual Nyquist ghosts. If the level of ghosting varies significantly from image to image, however, the signal cancellation that occurs with PLACE can adversely affect SNR-sensitive applications such as perfusion imaging with arterial spin labeling. This work proposes integrating PLACE into a PAGE-based reconstruction process to yield significantly better Nyquist ghost correction that is more robust than PLACE or PAGE alone. The robustness of this method is demonstrated in the presence of magnetic field drift with an in-vivo arterial spin labeling perfusion experiment.

Characterizing and correcting gradient errors in non-cartesian imaging: Are gradient errors linear time-invariant (LTI)?

Non-Cartesian and rapid imaging sequences are more sensitive to scanner imperfections such as gradient delays and eddy currents. These imperfections vary between scanners and over time and can be a significant impediment to successful implementation and eventual adoption of non-Cartesian techniques by scanner manufacturers. Differences between the k-space trajectory desired and the trajectory actually acquired lead to misregistration and reduction in image quality. While early calibration methods required considerable scan time, more recent methods can work more quickly by making certain approximations. We examine a rapid gradient calibration procedure applied to multiecho three-dimensional projection reconstruction (3DPR) acquisitions in which the calibration runs as part of every scan. After measuring the trajectories traversed for excitations on each of the orthogonal gradient axes, trajectories for the oblique projections actually acquired during the scan are synthesized as linear combinations of these measurements. The ability to do rapid calibration depends on the assumption that gradient errors are linear and time-invariant (LTI). This work examines the validity of these assumptions and shows that the assumption of linearity is reasonable, but that gradient errors can vary over short time periods (due to changes in gradient coil temperature) and thus it is important to use calibration data matched to the scan data.

Turboprop IDEAL: a motion-resistant fat-water separation technique

Suppression of the fat signal in MRI is very important for many clinical applications. Multi-point water-fat separation methods, such as IDEAL (Iterative Decomposition of water and fat with Echo Asymmetry and Least-squares estimation), can robustly separate water and fat signal, but inevitably increase scan time, making separated images more easily affected by patient motions. PROPELLER (Periodically Rotated Overlapping ParallEL Lines with Enhanced Reconstruction) and Turboprop techniques offer an effective approach to correct for motion artifacts. By combining these techniques together, we demonstrate that the new TP-IDEAL method can provide reliable water-fat separation with robust motion correction. The Turboprop sequence was modified to acquire source images, and motion correction algorithms were adjusted to assure the registration between different echo images. Theoretical calculations were performed to predict the optimal shift and spacing of the gradient echoes. Phantom images were acquired, and results were compared with regular FSE-IDEAL. Both T1- and T2-weighted images of the human brain were used to demonstrate the effectiveness of motion correction. TP-IDEAL images were also acquired for pelvis, knee, and foot, showing great potential of this technique for general clinical applications.

Least-squares chemical shift separation for (13)C metabolic imaging

PURPOSE

To describe a new least-squares chemical shift (LSCSI) method for separation of chemical species with widely spaced peaks in a sparse spectrum. The ability to account for species with multiple peaks is addressed.

MATERIALS AND METHODS

This method is applied to imaging of (13)C-labeled pyruvate and its metabolites alanine, pyruvate, and lactate. The method relies on a priori knowledge of the resonant frequencies of the different chemical species, as well as the relative signal from the two pyruvate peaks, one of which lies near the alanine peak. With this information a least-squares method was utilized for separation of signal from the three metabolites, facilitating tremendous reductions in the amount of data required to decompose the different chemical species. Optimization of echo spacing for maximum noise performance of the signal separation is also described.

RESULTS

Imaging an enriched (13)C phantom at 3.0T, the LSCSI method demonstrates excellent metabolite separation, very similar to echo planar spectroscopic imaging (EPSI), while only using 1/16th as much data.

CONCLUSION

This approach may be advantageous for in vivo hyperpolarized (13)C metabolic applications for reduced scan time compared with EPSI.

Fast decomposition of water and lipid using a GRASE technique with the IDEAL algorithm

Three-point Dixon techniques achieve good lipid-water separation by estimating the phase due to field inhomogeneities. Recently it was demonstrated that the combination of an iterative algorithm (iterative decomposition of water and fat with echo asymmetry and least-squares estimation (IDEAL)) with a fast spin-echo (FSE) three-point Dixon method yielded robust lipid-water decomposition. As an alternative to FSE, the gradient- and spin-echo (GRASE) technique has been developed for efficient data collection. In this work we present a method for lipid-water separation by combining IDEAL with the GRASE technique. An approach to correct for errors in the lipid-water decomposition caused by phase distortions due to the switching of the readout gradient polarities inherent to GRASE is presented. The IDEAL-GRASE technique is demonstrated in phantoms and in vivo for various applications, including pelvic, musculoskeletal, and (breath-hold) cardiac imaging.

Robust prescan calibration for multiple spin-echo sequences: application to FSE and b-SSFP

The collection of fast imaging techniques that use multiple spin-echo (MSE) sequences relies on a precise phase relationship between spin echoes and stimulated echoes that form along the radiofrequency refocusing pulse train. Failure to achieve this condition produces dark banding artifacts that result from destructive interference between signal coherence pathways. Satisfying this condition on the microsecond timescale required is technically challenging for conditions involving strong diffusion-weighted gradients, for arbitrary orientation acquisitions and at large field strengths with high-resolution acquisitions. Two clinically significant MSE sequences, fast spin echo (FSE) and balanced steady-state free precession (b-SSFP), are investigated in this work using a 4-T whole-body scanner. We developed a readout-projection-based prescan technique that ensures coherent signal formation by utilizing banding artifacts to automatically adjust gradient balance. Subsequent image acquisition using the results of this prescan permits the formation of coherent-echo images, which are robust under challenging imaging conditions. The robustness of this approach is demonstrated for FSE and b-SSFP images obtained from the knees of human volunteers. We believe that the use of this prescan calibration technique for the alignment of signal pools in MSE sequences is critical at high fields and will facilitate the implementation of high-quality clinically significant sequences such as FSE and b-SSFP.

Turboprop: Improved PROPELLER imaging

A variant of periodically rotated overlapping parallel lines with enhanced reconstruction (PROPELLER) MRI, called turboprop, is introduced. This method employs an oscillating readout gradient during each spin echo of the echo train to collect more lines of data per echo train, which reduces the minimum scan time, motion-related artifact, and specific absorption rate (SAR) while increasing sampling efficiency. It can be applied to conventional fast spin-echo (FSE) imaging; however, this article emphasizes its application in diffusion-weighted imaging (DWI). The method is described and compared with conventional PROPELLER imaging, and clinical images collected with this PROPELLER variant are shown.

Propeller EPI in the other direction

A new propeller EPI pulse sequence with reduced sensitivity to field inhomogeneities is proposed. Image artifacts such as blurring due to Nyquist ghosting and susceptibility gradients are investigated and compared with those obtained in previous propeller EPI studies. The proposed propeller EPI sequence uses a readout that is played out along the short axis of the propeller blade, orthogonal to the readout used in previous propeller methods. In contrast to long-axis readout propeller EPI, this causes the echo spacing between two consecutive phase-encoding (PE) lines to decrease, which in turn increases the k-space velocity in this direction and hence the pseudo-bandwidth. Long- and short-axis propeller EPI, and standard single-shot EPI sequences were compared on phantoms and a healthy volunteer. Diffusion-weighted imaging (DWI) was also performed on the volunteer. Short-axis propeller EPI produced considerably fewer image artifacts compared to the other two sequences. Further, the oblique blades for the long-axis propeller EPI were also prone to one order of magnitude higher residual ghosting than the proposed short-axis propeller EPI.

Calibration of gradient propagation delays for accurate two-dimensional radiofrequency pulses

Hardware-related delays between the requested and actual start times of the gradient waveforms on each physical axis are of particular importance for multidimensional selective excitation in which the synchronization of gradient and radiofrequency (RF) waveforms is critical. A method is proposed for the accurate calibration of gradient propagation delays to optimize the spatial accuracy of 2D RF pulses, although the results may also be used to reduce artifacts in other MR techniques. The sensitivity of 2D RF pulses to uncorrected time shifts between the gradient and RF waveforms was exploited to calibrate accurately the propagation delays on each physical gradient axis. This was achieved using a technique that relates the effect of gradient delays in the component waveforms of a constant-angular rate spiral k-space trajectory 2D RF pulse to the spatial location of the subsequent excitation profile. Comparison was also made with a procedure based on a previously described k-space plotting method, showing broad agreement, but with some discrepancies that illustrate the value of a self-referenced correction method for multidimensional RF pulses.

Removal of EPI Nyquist ghost artifacts with two-dimensional phase correction

Odd-even echo inconsistencies result in Nyquist ghost artifacts in the reconstructed EPI images. The ghost artifacts reduce the image signal-to-noise ratio and make it difficult to correctly interpret the EPI data. In this article a new 2D phase mapping protocol and a postprocessing algorithm are presented for an effective Nyquist ghost artifacts removal. After an appropriate k-space data regrouping, a 2D map accurately encoding low- and high-order phase errors is derived from two phase-encoded reference scans, which were originally proposed by Hu and Le (Magn Reson Med 36:166-171;1996) for their 1D nonlinear correction method. The measured phase map can be used in the postprocessing algorithm developed to remove ghost artifacts in subsequent EPI experiments. Experimental results from phantom, animal, and human studies suggest that the new technique is more effective than previously reported methods and has a better tolerance to signal intensity differences between reference and actual EPI scans. The proposed method may potentially be applied to repeated EPI measurements without subject movements, such as functional MRI and diffusion coefficient mapping.

Centering the projection reconstruction trajectory: reducing gradient delay errors

The projection reconstruction (PR) trajectory was investigated for the effect of gradient timing delays between the actual and requested start time of each physical gradient. Radial trajectories constructed with delayed gradients miss the center of k-space in an angularly dependent manner, causing effective echo times to vary with projection angle. The gradient timing delays were measured in phantoms, revealing delays on the x, y, and z gradients which differed by as much as 5 micro sec. Using this one-time calibration measurement, the trajectories were corrected for gradient delays by addition of compensatory gradient areas to the prephasers of the logical x and y readout gradients. Effective projection-to-projection echo time variability was reduced to less than 1 micro sec for all imaging orientations. Using corrected trajectories, artifacts were reduced in phantom images and in volunteer studies. This correction should potentiate greater clinical use of the PR trajectory.

Elimination of Nyquist ghosting caused by read-out to phase-encode gradient cross-terms in EPI

Echo-planar imaging (EPI) commonly suffers from ghosting artifacts caused by zero- or first-order phase differences between the odd and even echoes that constitute an EPI dataset. Small-bore imaging systems with shielded gradients may suffer significantly from cross-term eddy currents due to the high degree of manufacturing precision required in such systems compared to larger whole-body coils. A Nyquist ghost caused by cross-term eddy current contributions from the read-out to the phase-encode axis was identified in a small-bore system and characterized using a modified EPI experiment and a computer simulation. The artifact was corrected for using both a postprocessing approach and compensation blips along the phase-encode axis. Correction using compensation blips proved to be a more effective strategy to reduce this artifact than the postprocessing method used.

Ultrafast pulse sequence techniques for cardiac magnetic resonance imaging

Cardiac magnetic resonance imaging is a rapidly emerging field that has seen tremendous advances in the past decade. Central to the development of effective imaging strategies has been the advent of high-performance gradient hardware and the exploitation of their speed characteristics through specialized pulse sequences well suited for cardiac imaging. These advances have facilitated unprecedented acquisition times that now approach echocardiographic frame rates, while maintaining excellent image quality. This article provides a detailed overview of advanced pulse sequence technology and approaches currently taken to maximize speed performance and image quality. In particular, segmented K-space techniques that include single-echo and multiecho spoiled gradient-echo imaging as well as steady-state free precession imaging are discussed. Finally, spiral and fast spin-echo techniques are explored. Examples of common applications of these pulse sequences are presented.