Diffusion-weighted interleaved echo-planar imaging with a pair of orthogonal navigator echoes

This work describes a diffusion-weighted (DW) interleaved echo-planar imaging (IEPI) method for use on either conventional whole-body scanners or scanners equipped with highspeed gradient and receiver hardware. In combination with cardiac gating and motion correction with a pair of orthogonal navigator echoes, the presented method is time-efficient, compensates for patient motions, and is less sensitive to image distortions than single-shot methods. The motion-correction scheme consists of correction for constant and linear phase terms found from the orthogonal navigator echoes. The correction for the linear phase term in the phase-encode direction includes gridding the data to the Cartesian grid. The DW IEPI was used to image a phantom rotating about the slice-select direction, and motion correction was performed to eliminate ghost artifacts arising from motion in either the readout- or phase-encoding directions. DW IEPI with motion correction was performed on a normal volunteer and on a patient with a 26-day-old region of ischemia over much of the right hemisphere.     

Reduced field-of-view MRI using outer volume suppression for spinal cord diffusion imaging.

A spin-echo single-shot echo-planar imaging (SS-EPI) technique with a reduced field of view (FOV) in the phase-encoding direction is presented that simultaneously reduces susceptibility effects and motion artifacts in diffusion-weighted (DW) imaging (DWI) of the spinal cord at a high field strength (3T). To minimize aliasing, an outer volume suppression (OVS) sequence was implemented. Effective fat suppression was achieved with the use of a slice-selection gradient-reversal technique. The OVS was optimized by numerical simulations with respect to T(1) relaxation times and B(1) variations. The optimized sequence was evaluated in vitro and in vivo. In simulations the optimized OVS showed suppression to <0.25% and approximately 3% in an optimal and worst-case scenario, respectively. In vitro measurements showed a mean residual signal of <0.95% +/- 0.42 for all suppressed areas. In vivo acquisition with 0.9 x 1.05 mm(2) in-plane resolution resulted in artifact-free images. The short imaging time of this technique makes it promising for clinical studies.     

A Velocity k-Space Analysis of Flow Effects in Echo-Planar and Spiral Imaging

A velocity k-space formalism facilitates the analysis of flow effects for imaging sequences involving time-varying gradients such as echo-planar and spiral. For each sequence, the velocity k-space trajectory can be represented by k_v (k)_r; that is, its velocity-frequency (k_r) position as a function of spatial-frequency (k_r) position. In an echo-planar sequence, k_r is discontinuous and asymmetric. However, in a spiral sequence, k_r is smoothly varying, circularly symmetric, and small near the k_r origin. To compare the effects of these trajectory differences, simulated images were generated by computing the k-space values for an in-plane vessel with parabolic flow. Whereas the resulting echo-planar images demonstrate distortions and ghosting that depend on the vessel orientation, the spiral images exhibit minimal artifacts.     

Analysis of flow effects in echo-planar imaging

The effects on the phase of spins moving during echo-planar imaging (EPI) acquisition were studied. Standard single-shot and interleaved multishot blipped EPI acquisitions were considered, assuming either high gradient strength and slew rates or standard gradient strength and slew rates. A spiral k-space trajectory was also considered. Flow components in the section-select and phase- and frequency-encoding directions were analyzed separately. While the effect of flow in the section-select direction is identical to that in a standard two-dimensional Fourier transform (2DFT) acquisition, flow in the phase- or frequency-encoding directions can have substantial effects on the image, different from that in 2DFT imaging. The magnitude of these effects, which include displacement, distortion, and/or ghosting of vascular structures, is analyzed and predicted for a given velocity and direction of flow, the specific acquisition sequence, and the strength and slew rate of the gradients. For example, 50-cm/sec flow along the phase-encoding direction can cause a blurring of 1.25 cm full width at half maximum for blipped EPI with high-strength gradients, assuming a 40-cm field of view and 64 × 64 matrix.     

Abbreviated moment-compensated phase encoding

To achieve correct spatial location of blood vessels, first order gradient moment nulling applied to the phase encoding axes can be used. However, gradient moment nulling prolongs echo time (TE), which may degrade the flow image in regions of complex flow. The fact that abbreviated moment compensated phase-encoding (AMCPE) can be used to apply partial flow compensation to the phase-encoding axes to prevent spatial misregistration of vessels without requiring the use of long echo times or using arbitrary chosen TE is demonstrated. AMCPE defines two cutoff lines in k-space. The flow-induced phase is completely compensated for values between the cutoff lines and partially compensated beyond the cutoff lines. The AMCPE technique has been tested on both a flow phantom and a human volunteer. The AMCPE images from both the in vivo and the in vitro study demonstrate correctly imaged flow. Computer simulations have been performed to analyze the penalty caused by the incomplete flow compensation. The result shows that the ripple artifacts due to the incomplete flow compensation are unobservable when 60%–70% of k-space is completely flow compensated.     

Ghost phase cancellation with phase-encoding gradient modulation

Motion artifacts are a dominant cause of magnetic resonance image quality degradation. Periodic or nearly periodic motion results in image replicates of the moving structures in spin-warp Fourier imaging. The replicates, or ghosts, propagate in the image in the phase encoding, or y, direction. These ghosted images can be considered to consist of the time-averaged spin density I₀and a ghost mask g. A set of J ghosted images I_j may be acquired in which the ghost mask is intentionally phase shifted by varying amounts relative to I₀ with interleaved acquisitions that have shifted phaseencoding orders or by acquiring multiple images during a single readout period in the presence of an oscillating phase-encoding gradient. The resulting complex images I_j have the same time-averaged spin density I₀ but have ghost contributions g_j that, on a pixel-by-pixel basis, trace part of a circle around I₀. The source images I_j can then be used to estimate I₀. Simulations and experiments with the phase-encoding gradient modulation method show good general ghost suppression for a variety of quasi-periodic motion sources including both respiratory-type artifacts and flow artifacts. The primary limitation of the method is the need for rapid gradient switching.     

Effect of and correction for in-plane myocardial motion on estimates of coronary-volume flow rates

The sensitivities of phase-difference (PD) and complex-difference (CD) processing strategies to in-plane motion were examined theoretically and experimentally. Errors in velocity and volume flow rate (VFR) estimates were attributed to (a) motion between different velocity encodings and, in the case of segmented k-space acquisition strategies, (b) motion over the segment duration. PD estimates were found to be insensitive to in-plane motion between velocity encodings, whereas CD VFR estimates were found to be sensitive to this motion. PD estimates, however, were affected by partial volume effects. A corrected CD (CD') scheme was developed that minimizes both partial-volume and in-plane motion effects. Segmented k-space acquisitions with sequential offset and sequential interleaved offset (or centric) phase-encoding schemes were studied. Images obtained using these techniques were found to include both blurring and replication artifacts. The amount of artifact generally increased with the number of views per segment (vps) and the in-plane velocity. PD, CD, and CD' VFR estimates were found to be degraded by these artifacts. The sequential offset phase-encoding scheme generally had acceptable VFR errors (at 4 vps, a CD' VFR error of 7.0%) when averaged over the physiologic range of myocardial motion (>12 cm second^?1); however, larger errors were observed outside this range. VFR estimates obtained using the sequential interleaved phase-encoding scheme at 4 vps were unacceptable. More accurate VFR measurements were obtained using a revised segmented PC strategy, which reversed the order in which the velocity and phase encodings were interleaved. The weighted average CD' VFR error obtained using the revised strategy was 24.5% (for 4 vps). Using displacement information obtained from the two velocity-encoded images, an estimate of the in-plane velocity was obtained and used to correct the acquired data. This decreased the VFR error (weighted average CD' error at 4 vps decreased from 24.5% to ?6.3%); however, the implemented correction algorithm could potentially introduced other artifacts in the images.     

Elimination of oblique flow artifacts in magnetic resonance imaging.

We present an analysis of how flow oblique to the frequency-encoding direction generates displacement artifacts in MR imaging and show that for flow which has constant velocity between the start of the phase encoding and the center of the echo it is possible to eliminate these artifacts by gradient moment nulling in the phase-encoding direction. However, unlike the standard moment nulling calculations for flow compensating the frequency-encode and slice-selection gradients, the phase-encoding first moment must be nulled specifically with respect to the echo center. Limitations of this method imposed by finite gradient strengths are analyzed. In 3D volume acquisitions with two axes phase encoded it is possible to correct for oblique flow in all directions, and this is demonstrated in images of a human volunteer. Correction for oblique flow displacement artifacts may be particularly useful in quantitative flow and angiographic applications.     

GRASE imaging at 3 Tesla with template interactive phase–encoding

A new method for ordering the phase-encoding gradient is proposed, and an application for short effective TE gradient-and spin-echo (GRASE) imaging is demonstrated. The proposed method calculates the phase-encoding order from the signal decay of a template scan (hence “template interactive phase-encoding” or TIPE). Computer simulations are used to compare the point spread functions of different phase-encoding orders giving short effective echo times (kb centric GRASE, centric GRASE, centric TIPE). The conventional centric phase-encoding order is also considered for GRASE. The conventional centric method is sensitive to both amplitude and phase modulation of the signal in K-space. The centric TIPE method gives the least amplitude modulation artifacts but is vulnerable to phase artifacts. The TIPE experiment was implemented on a 3 Tesla system. To the best of our knowledge, we present the first in vivo GRASE images at this field strength.     

Reducing oblique flow effects in interleaved EPI with a centric reordering technique.

Segmented interleaved echo planar imaging offers a fast and efficient approach to magnetic resonance angiography. Unfortunately, this technique is particularly sensitive to oblique flow in the imaging plane. In this work, a mathematical analysis of oblique flow effects for several types of k-space coverage is presented. The conventional linear acquisition scheme, an alternating centric and a nonalternating centric encoding scheme are compared with respect to their flow properties. It is shown both by simulations and imaging experiments that artifacts from oblique in-plane flow are effectively reduced by both centric reordered phase-encoding schemes. The nonalternating centric acquisition scheme is preferred to the alternating centric scheme due to the smoother phase error transition in k-space in the presence of obliquely-angled flow. Magn Reson Med 45:623-629, 2001.     

Time of flight quantification of coronary flow with echo-planar MRI

Detection and quantification of flow of the left anterior descending (LAD) coronary artery in healthy volunteers are demonstrated using echo-planar imaging (EPI). A time-of-flight TOF) model was used to derive coronary flow velocities from wash-in curves, free of cardiac wall motion contamination. Short-axis cardiac studies were performed using a gated, gradient echo EPI technique to limit the effect of cardiac wall motion on coronary vessel imaging. A series of 10 to 20 single or multislice images were acquired within a single breath-hold. Real-time cine series showed the LAD coronary artery with a detectability of 91% (n = 23) and revealed beat-to-beat variability in vessel position of a magnitude equal to or greater than its diameter. Flow velocity was measured in the proximal portion of the artery at rest and during exercise. The data demonstrated the known phasic pattern of LAD flow:-V_systole ≤ 5 cm/s and peak V_diastole = 14 ± 3 cm/s (n = 11, V = mean laminar flow velocity). During isometric exercise, a LAD flow velocity increase (52 ± 24%) was detected in eight of nine subjects. The capacity of the EPI TOF method to detect flow velocity changes should prove clinically useful for future assessment of coronary flow reserve.     

Theoretical aspects of motion sensitivity and compensation in echo-planar imaging.

Magnetic resonance (MR) imaging can be performed on or below the time scale of most anatomic motion via echo-planar imaging (EPI) techniques and their derivatives. The goal is to image rapidly and reduce artifacts that typically result from view-to-view changes in the spatial distribution of spins due to motion. However, the required time-dependent magnetic field gradient waveforms remain sensitive to the dephasing effects of motion. Sources of motion artifact are simulated for spins moving along the imaging axes and are shown to be an important source of reduced image quality in EPI. A novel method of EPI is proposed that (a) refocuses single or multiple derivatives of motion at all echoes and (b) prevents accumulation of velocity (or higher derivative)--induced dephasing along the phase-encoding axis by moment nulling all phase-encoding-step waveforms about a single instant of time. Theoretical EPI sequences with considerable reductions in ghosts, blurring, and signal loss due to motion sensitivity are produced and compared with other EPI methods. Their time efficiency is presented as a function of available (relative) gradient strength for a variety of sequence waveforms.     

Isotropic diffusion-weighted and spiral-navigated interleaved EPI for routine imaging of acute stroke

An interleaved echo-planar imaging (EPI) technique is presented for the rapid acquisition of isotropic diffusion-weighted images of stroke patients. Sixteen isotropic diffusion-weighted images at three b values are acquired in less than 3 min. A spiral navigator echo is used to measure the constant and linear phase shifts across the head in both the x and y directions which result from motion during the isotropic diffusion- sensitizing gradients. The measured k-space errors are corrected during a gridding reconstruction. The gridding kernel has a constant width in k_x and a variable width in k_y which eliminates variable data-density ghosts. The resulting isotropic diffusion-weighted images have excellent lesion-to-normal brain contrast, very good spatial resolution, and little sensitivity to susceptibility effects in the base of the brain. Examples of diffusion-weighted images and ADC maps from several stroke patients are shown.     

Phase-Encode reordering to minimize errors caused by motion

A new method for suppressing the effects of motion in MR images by reordering the acquisition of k space has been developed. Existing reordering methods suffer from image blurring. The method presented here applies specifically to translation along the phase-encoding direction, in which case it reduces both ghosting and blurring. The method is conceptually similar to a linear phase shift of k space, except that it has the advantage of not corrupting stationary structures inadvertently. The method is intended for anatomic sites in which substantial translational motion occurs along the phase-encoding direction, such as the cranial-caudal motion of the liver and kidneys. The reordering method is motivated from an analysis of factors affecting the severity of motion artifacts. The theory behind the reordering method is described and validated experimentally by imaging a moving phantom.     

Consistent fat suppression with compensated spectral-spatial pulses

Reliable fat suppression is especially important with fast imaging techniques such as echo-planar (EPI), spiral, and fast spin-echo (FSE) T₂-weighted imaging. Spectral-spatial excitation has a number of advantages over spectrally selective presaturation techniques, including better resilience to B₀ and B₁, inhomogeneity. In this paper, a FSE sequence using a spectral-spatial excitation pulse for superior fat suppression is presented. Previous problems maintaining the CPMG condition are solved using simple methods to accurately program radio-frequency (RF) phase. Next an analysis shows how B₀ eddy currents can reduce fat suppression effectiveness for spectral-spatial pulses designed for conventional gradient systems. Three methods to compensate for the degradation are provided. Both the causes of the degradation and the compensation techniques apply equally to gradient-recalled applications using these pulses. These problems do not apply to pulses designed for high-speed gradient systems. The spectral-spatial FSE sequence delivers clinically lower fat signal with better uniformity than spectrally selective pre-saturation techniques.     

Usefulness of free-breathing readout-segmented echo-planar imaging (RESOLVE) for detection of malignant liver tumors: Comparison with single-shot echo-planar imaging (SS-EPI)

Purpose

We aimed to clarify the usefulness of free-breathing readout-segmented echo-planar imaging (RESOLVE), which is multi-shot echo-planar imaging based on a 2D-navigator-based reacquisition technique, for detecting malignant liver tumor.

Materials and methods

In 77 patients with malignant liver tumors, free-breathing RESOLVE and respiratory-triggered single-shot echo-planar imaging (SS-EPI) at 3-T MR unit were performed. We set a scan time up to approximately 5 min (300 s) before examination, measured actual scan time and assessed (1) susceptibility and (2) motion artifacts in the right and left liver lobes (3, no artifact; 1, marked), and (3) detectability of malignant liver tumors (3, good; 1, poor) using a 3-point scale.

Results

The median actual scan time of RESOLVE/SS-EPI was 365/423 s. The median scores of each factor in RESOLVE/SS-EPI were as following in this order: (1) 3/2 (right lobe); 3/3 (left lobe), (2) 2/3 (right lobe); 1/2 (left lobe), and (3) 3/3, respectively. Significant differences were noted between RESOLVE and SS-EPI in all evaluated factors (P < 0.05) except for susceptibility of left lobe and detectability of the lesions.

Conclusion

Despite the effect of motion artifacts, RESOLVE provides a comparable detectability of the lesion and the advantage of reducing scanning time compared with SS-EPI.     

Reduction of MR imaging time by the hybrid fast-scan technique

The time taken to collect high-resolution and high signal-to-noise (S/N) data in magnetic resonance (MR) imaging may be the limiting factor in patient throughput and in reducing patient motion. A hybrid fast-scan technique combining static and oscillatory phase-encoding gradients from two-dimensional Fourier transform (2DFT) and echo-planar imaging can reduce the time needed to collect data at the expense of loss in S/N. The flexibility of this technique is that any amplitude or frequency of oscillation of the phase-encoding gradient can be used. The technique was used for different frequencies and amplitudes, and images are presented that were acquired in one-half and one-quarter the time required with standard 2DFT techniques. The images illustrate that the hybrid and 2DFT techniques produce comparable resolution and contrast under identical conditions.     

Measurement of the point spread function in MRI using constant time imaging

The point spread function is a fundamental property of magnetic resonance imaging methods that affects image quality and spatial resolution. The point spread function is difficult to measure precisely in magnetic resonance even with the use of carefully designed phantoms, and it is difficult to calculate this function for complex sequences such as echo-planar imaging. This report describes a method that measures the point spread function with high spatial resolution at each pixel in samples of uniform intensity distribution. This method uses additional phase encoding gradients before the echo-planar acquisition that are constant in length but vary in amplitude. The additional gradients are applied to image the contents within each individual voxel. This method has been used to measure the point spread function for echo-planar imaging to demonstrate the effects of limited k-space sampling and transverse relaxation, as well as the effects of object motion. By considering the displacement of the point spread function, local distortions due to susceptibility and chemical shift effects have been quantified and corrected. The method allows rapid assessment of the point spread function in echo-planar imaging, in vivo, and may also be applied to other rapid imaging sequences that can be modified to include these additional phase encoding gradients.     

Improved clinical echo-planar MRI using spatial-spectral excitation

Echo-planar imaging (EPI) is markedly susceptible to B₀ field distortions and to frequency differences caused by chemical shift, because the phase of the signals is accumulating during the acquisition train. Thus, only water proton signals are usually recorded after frequency-selective suppression of the fat portions of magnetization. Otherwise, a shifted signal frequency from fat results in ghosting artifacts. In this article, a technique is presented working with spatial-spectral excitation for highly selective water or fat EPI. The proposed method allows recording in multislice operation on EPI scanners without irregular gradient or RF pulse shapes. Examples of gradient-echo and spin-echo EPI using spatial-spectral excitation by series of two to eight single slice-selective RF pulses are demonstrated. The method is not sensitive to misadjustments or inhomogeneities of the B₁ field, but sufficient homogeneity of the static magnetic field B₀ is required. Especially the quality of diffusion-weighted echo-planar images can be markedly improved by the new technique compared to conventional EPI, because artifacts from undesired chemical shift components are completely avoided.     

Real time blood flow imaging by spiral scan phase velocity mapping

The work describes the development of a novel sequence that uses rapid spiral k-space sampling, combined with phase velocity mapping, for real time flow velocity imaging. The performance of the technique is assessed on phantoms for both through-plane and in-plane flows. The flow measurements compared well with those measured using a bucket and stopwatch. One advantage of the technique is that flow related signal loss is minimal due to the early acquisition of the center of k-space data. Flow artifacts were observed for in-plane flow and these were understood with the aid of computer simulations. In vivo studies involved cine velocity mapping in normal volunteers; aortic blood flow waveforms acquired by spiral scanning in two cardiac cycles compared well with data from a conventional gradient-echo sequence. Potential applications of the method are demonstrated by studying the response of aortic flow to physical exercise and the real time monitoring of aortic flow during a valsalver maneuver.