Direct parallel image reconstructions for spiral trajectories using GRAPPA.

The use of spiral trajectories is an efficient way to cover a desired k-space partition in magnetic resonance imaging (MRI). Compared to conventional Cartesian k-space sampling, it allows faster acquisitions and results in a slight reduction of the high gradient demand in fast dynamic scans, such as in functional MRI (fMRI). However, spiral images are more susceptible to off-resonance effects that cause blurring artifacts and distortions of the point-spread function (PSF), and thereby degrade the image quality. Since off-resonance effects scale with the readout duration, the respective artifacts can be reduced by shortening the readout trajectory. Multishot experiments represent one approach to reduce these artifacts in spiral imaging, but result in longer scan times and potentially increased flow and motion artifacts. Parallel imaging methods are another promising approach to improve image quality through an increase in the acquisition speed. However, non-Cartesian parallel image reconstructions are known to be computationally time-consuming, which is prohibitive for clinical applications. In this study a new and fast approach for parallel image reconstructions for spiral imaging based on the generalized autocalibrating partially parallel acquisitions (GRAPPA) methodology is presented. With this approach the computational burden is reduced such that it becomes comparable to that needed in accelerated Cartesian procedures. The respective spiral images with two- to eightfold acceleration clearly benefit from the advantages of parallel imaging, such as enabling parallel MRI single-shot spiral imaging with the off-resonance behavior of multishot acquisitions.     

Fat/water separation using a concentric rings trajectory

The concentric rings two‐dimensional (2D) k‐space trajectory enables flexible trade‐offs between image contrast, signal‐to‐noise ratio (SNR), spatial resolution, and scan time. However, to realize these benefits for in vivo imaging applications, a robust method is desired to deal with fat signal in the acquired data. Multipoint Dixon techniques have been shown to achieve uniform fat suppression with high SNR‐efficiency for Cartesian imaging, but application of these methods for non‐Cartesian imaging is complicated by the fact that fat off‐resonance creates significant blurring artifacts in the reconstruction. In this work, two fat–water separation algorithms are developed for the concentric rings. A retracing design is used to sample rings near the center of k‐space through multiple revolutions to characterize the fat–water phase evolution difference at multiple time points. This acquisition design is first used for multipoint Dixon reconstruction, and then extended to a spectroscopic approach to account for the trajectory's full evolution through 3D k‐t space. As the trajectory is resolved in time, off‐resonance effects cause shifts in frequency instead of spatial blurring in 2D k‐space. The spectral information can be used to assess field variation and perform robust fat–water separation. In vivo experimental results demonstrate the effectiveness of both algorithms. Magn Reson Med, 2009. © 2008 Wiley‐Liss, Inc.     

3D magnetization‐prepared imaging using a stack‐of‐rings trajectory

Efficient acquisition strategies for magnetization‐prepared imaging based on the three‐dimensional (3D) stack‐of‐rings k‐space trajectory are presented in this work. The 3D stack‐of‐rings can be acquired with centric ordering in all three dimensions for greater efficiency in capturing the desired contrast. In addition, the 3D stack‐of‐rings naturally supports spherical coverage in k‐space for shorter scan times while achieving isotropic spatial resolution. While non‐Cartesian trajectories generally suffer from greater sensitivity to system imperfections, the 3D stack‐of‐rings can enhance magnetization‐prepared imaging with a high degree of robustness to timing delays and off‐resonance effects. As demonstrated with phantom scans, timing errors and gradient delays only cause a bulk rotation of the 3D stack‐of‐rings reconstruction. Furthermore, each ring can be acquired with a time‐efficient retracing design to resolve field inhomogeneities and enable fat/water separation. To demonstrate its effectiveness, the 3D stack‐of‐rings are considered for the case of inversion‐recovery‐prepared structural brain imaging. Experimental results show that the 3D stack‐of‐rings can achieve higher signal‐to‐noise ratio and higher contrast‐to‐noise ratio within a shorter scan time when compared to the standard inversion‐recovery‐prepared sequence based on 3D Cartesian encoding. The design principles used for this specific case of inversion‐recovery‐prepared brain imaging can be applied to other magnetization‐prepared imaging applications. Magn Reson Med 63:1210–1218, 2010. © 2010 Wiley‐Liss, Inc.     

Novel interleaved spiral imaging motion correction technique using orbital navigators.

Although spiral imaging seldom produces apparent artifacts related to flow, it remains sensitive to rapid object motion. In this article, a new correction method is presented for rapid rigid body motion in interleaved spiral imaging. With this technique, an identical circular navigator k-space trajectory is linked to each spiral trajectory. Data inconsistency due to both rotation and translation among spiral interleaves can be corrected by evaluating the magnitudes and phases of the data contained in the navigator "ring." Further, it is difficult to create a frequency field map for off-resonance correction when an object moves during a scan, because there is motion-dependent misregistration between the two images acquired with different TEs. However, this difficulty can be overcome by combining the motion-correction method with a recently proposed technique (off-resonance correction using variable-density spirals (ORC-VDS)), thereby enabling both motion compensation and off-resonance correction with no additional scanning.     

Fast three-dimensional k-space trajectory design using missile guidance ideas.

Three-dimensional (3D) k-space trajectories are needed to acquire volumetric images in MRI. While scan time is determined by the trajectory efficiency, image quality and distortions depend on the shape of the trajectories. There are several 3D trajectory strategies for sampling the k-space using rectilinear or curve schemes. Since there is no evidence about their optimality in terms of image quality and acquisition time, a new design method based on missile guidance ideas is explored. Since air-to-air missile guidance shares similar goals and constraints with the problem of k-space trajectory design, a control approach for missiles is used to design a 3D trajectory. The k-space is divided into small cubes, and each one is treated as a target to be sampled. The main goal is to cover the entire space as quickly and efficiently as possible, with good performance under different conditions. This novel design method is compared to other trajectories using simulated and real data. As an example, a trajectory that requires 0.11 times the number of shots needed by the cylindrical 3DFT acquisition was designed. This trajectory requires more shots (1.66 times) than the stack of spirals, but behaves better under nonideal conditions, such as off-resonance and motion.     

Inversion recovery radial MRI with interleaved projection sets.

The radial trajectory has found applications in cardiac imaging because of its resilience to undersampling and motion artifacts. Recent work has shown that interleaved and weighted radial imaging can produce images with multiple contrasts from a single data set. This feature was investigated for inversion recovery imaging of scar using a radial technique. The 2D radial imaging method was modified to acquire quadruply interleaved projection sets within each acquisition window of the cardiac cycle. These data were reconstructed using k-space weightings that used a smaller segment of the acquisition window for the central k-space data, the determinant of image contrast. This method generates four images with different T1 weightings. The novel approach was compared with noninterleaved radial imaging, interleaved radial without weightings, and Cartesian imaging in simulations, phantoms, and seven subjects with clinical myocardial infarction. The results show that during a typical acquisition window after an inversion pulse, magnetization changes rapidly. The interleaved acquisition provided better image quality than the noninterleaved radial acquisition. Interleaving with weighting provided better quality when the inversion time (TI) was shorter than optimal; otherwise, interleaving without weighting was superior. These methods enable a radial trajectory to be employed in conjunction with preparation pulses for viability imaging.     

Design and analysis of a practical 3D cones trajectory.

The 3D Cones k-space trajectory has many desirable properties for rapid and ultra-short echo time magnetic resonance imaging. An algorithm is presented that generates the 3D Cones gradient waveforms given a desired field of view and resolution. The algorithm enables a favorable trade-off between increases in readout time and decreases in the total number of required readouts. The resulting trajectory is very signal-to-noise ratio (SNR) efficient and has excellent aliasing properties. A rapid high-resolution ultra-short echo time imaging sequence is used to compare the 3D Cones trajectory to 3D projection reconstruction (3DPR) sampling schemes. For equivalent scan times, the 3D Cones trajectory has better SNR performance and fewer aliasing artifacts as compared to the 3DPR trajectory.     

Spectroscopic imaging using concentrically circular echo-planar trajectories in vivo

An alternative to the standard echo-planar spectroscopic imaging technique is presented, spectroscopic imaging using concentrically circular echo-planar trajectories (SI-CONCEPT). In contrast to the conventional chemical shift imaging data, the sampled data from each set of concentric rings were regridded into Cartesian space. Usage of concentric k-space trajectories has the advantage of requiring significantly reduced slew rates than echo-planar spectroscopic imaging, allowing for the collection of higher spectral bandwidths and opening the door for high-bandwidth echo-planar styled spectroscopic imaging at higher magnetic fields. Before two-dimensional spatial and one-dimensional spectral encoding, the volume of interest was localized using the standard point-resolved spectroscopy sequence. The feasibility of using concentric k-space trajectories is demonstrated, and the spatial profiles and representative spectra are compared with the standard echo-planar spectroscopic imaging technique in a gray matter phantom containing metabolites at physiological concentrations and healthy human brain in vivo. The symmetric nature of the concentric circles also reduces the number of required excitations for a given resolution by a factor of two. Possible artifacts and limitations are discussed.     

Real-time cardiac cine imaging with SPIDER: steady-state projection imaging with dynamic echo-train readout.

Steady-state projection imaging with dynamic echo-train readout (SPIDER) is a multiecho radial k-space trajectory TrueFISP sequence developed for real-time cine imaging of the heart. This new pulse sequence combines the superior SNR and blood-to-myocardium contrast of TrueFISP with the increased scan time efficiency of EPI and undersampled projection reconstruction. SPIDER sequence RF repetition time (TR) was minimized by limiting the echo-train to a length of three while acquiring the first and third echoes asymmetrically. A temporal resolution of 45 ms was achieved with TR/TE1/TE2/TE3 of 3.24/0.6/1.6/2.6 ms and a factor of 2 view sharing scheme. Phantom experiments showed little difference between the weighting of the signals acquired at each of the echo times but did show considerable off-resonance modulation between them. In vivo experiments demonstrated the feasibility of using the SPIDER sequence for real-time imaging in the cardiac short axis orientation.     

Centric ordering is superior to gradient moment nulling for motion artifact reduction in EPI

Echo-planar imaging (EPI) is sensitive to motion despite its rapid data acquisition rate. Compared with traditional imaging techniques, it is more sensitive to motion or flow in the phase-encode direction, which can cause image artifacts such as ghosting, misregistration, and loss of spatial resolution. Consequently, EPI of dynamic structures (eg, the cardiovascular system) could benefit from methods that eliminate these artifacts. In this paper, two methods of artifact reduction for motion in the phase-encode direction are evaluated. First, the k-space trajectory is evaluated by comparing centric with top-down ordered sequences. Next, velocity gradient moment nulling (GMN) of the phase-encode direction is evaluated for each trajectory. Computer simulations and experiments in flow phantoms and rabbits in vivo show that uncompensated centric ordering produces the highest image quality. This is probably due to a shorter readout duration, which reduces T2* relaxation losses and off-resonance effects, and to the linear geometry of phantoms and vessels, which can obscure centric blurring artifacts.     

Fast magnetic resonance coronary angiography with a three-dimensional stack of spirals trajectory.

In this work, three-dimensional (3D) spiral imaging has been utilized for magnetic resonance coronary angiography. Spiral-based 3D techniques can dramatically reduce imaging time requirements compared with 3D Fourier Transform imaging. The method developed here utilized a "stack of spirals" trajectory, to traverse 3D k-space rapidly. Both thick-slab volumes encompassing the entire coronary tree with isotropic resolution and thin-slab volumes targeted to a particular vessel of interest were acquired. Respiratory compensation was achieved using the diminishing variance algorithm. T2-prepared contrast was also applied in some cases to improve contrast between vessel and myocardium, while off-resonance blurring was minimized by applying a linear correction to the acquired data. Images from healthy volunteers were displayed using a curved reformatting technique to view long segments of vessel in a single projection. The results demonstrate that this 3D spiral technique is capable of producing high-quality coronary magnetic resonance angiograms.     

Generalized k-space decomposition with chemical shift correction for non-Cartesian water-fat imaging.

Chemical-shift artifacts associated with non-Cartesian imaging are more complex to model and less clinically acceptable than the bulk fat shift that occurs with conventional spin-warp Cartesian imaging. A novel k-space based iterative decomposition of water and fat with echo asymmetry and least-squares estimation (IDEAL) approach is introduced that decomposes multiple species while simultaneously correcting distortion of off-resonant species. The new signal model accounts for the additional phase accumulated by off-resonant spins at each point in the k-space acquisition trajectory. This phase can then be corrected by adjusting the decomposition matrix for each k-space point during the final IDEAL processing step with little increase in reconstruction time. The technique is demonstrated with water-fat decomposition using projection reconstruction (PR)/radial, spiral, and Cartesian spin-warp imaging of phantoms and human subjects, in each case achieving substantial correction of chemical-shift artifacts. Simulations of the point-spread-function (PSF) for off-resonant spins are examined to show the nature of the chemical-shift distortion for each acquisition. Also introduced is an approach to improve the signal model for species which have multiple resonant peaks. Many chemical species, including fat, have multiple resonant peaks, although such species are often approximated as a single peak. The improved multipeak decomposition is demonstrated with water-fat imaging, showing a substantial improvement in water-fat separation.     

Importance of k-space trajectory in echo-planar myocardial tagging at rest and during dobutamine stress.

Hybrid fast gradient echo/echo-planar imaging (FGRE-EPI) can be used to increase temporal resolution, enhance tag contrast, and/or decrease scan time for breathhold myocardial tagging. However, off-resonance effects and motion can lead to local phase discontinuities in FGRE-EPI raw data when a conventional interleaved bottom-up k-space trajectory is used. These discontinuities can be particularly problematic for myocardial tagging, where the image energy is not only concentrated near the k-space origin, but is also concentrated in multiple spectral peaks centered throughout k-space. In this study, tag distortion artifacts in FGRE-EPI tagging due to off-resonance and velocity-induced phase discontinuities were characterized at rest and dobutamine stress, and the flyback and gradient moment smoothing (GMS) methods were shown to reduce these artifacts. For the specific parameters used in this study, flyback and GMS resulted in improved image quality at rest and stress, increased myocardium-tag contrast-to-noise ratio (11.4 +/- 2.1 vs. 10.0 +/- 2.9, P < 0.01 at rest; 11.1 +/- 1.8 vs. 8.1 +/- 2.4, P < 0.01 at stress), and reduced full width at half maximum of the tag profile (3.6 vs. 3.8 pixels at rest; 4.0 vs. 5.1 pixels at stress) compared to the conventional trajectory. A limitation of the improved trajectory is a parameter-dependent decrease in data acquisition efficiency. For the specific imaging protocol used, the repetition time of the improved trajectory increased by 36% compared to the conventional trajectory.     

Spiral scanning with anisotropic field of view

Spiral scanning has been used to achieve much shorter scan times than conventional techniques for a wide range of applications. The major drawback with spiral scans is blurring from off-resonant spins, which is proportional to the readout time. Blurring limits maximal spatial resolution and minimal scan time potentially achievable with spiral scanning. Anisotropic field of view is used in conventional scanning to improve image quality by matching /c-space trajectory to object characteristics. Anisotropic field of view improves spatial resolution in spiral scanning without increasing scan time or blurring. The resolution improvement results from increased maximal k-space radius allowed by the lower field of view. A field of view reduction by a factor of 2 in one direction provides up to 60% resolution improvement in that direction. Reduced SNR also results from non-uniform /c-space sampling.     

Fast 3D imaging using variable-density spiral trajectories with applications to limb perfusion.

Variable-density k-space sampling using a stack-of-spirals trajectory is proposed for ultra fast 3D imaging. Since most of the energy of an image is concentrated near the k-space origin, a variable-density k-space sampling method can be used to reduce the sampling density in the outer portion of k-space. This significantly reduces scan time while introducing only minor aliasing artifacts from the low-energy, high-spatial-frequency components. A stack-of-spirals trajectory allows control over the density variations in both the k(x)-k(y) plane and the k(z) direction while fast k-space coverage is provided by spiral trajectories in the k(x)-k(y) plane. A variable-density stack-of-spirals trajectory consists of variable-density spirals in each k(x)-k(y) plane that are located in varying density in the k(z) direction. Phantom experiments demonstrate that reasonable image quality is preserved with approximately half the scan time. This technique was then applied to first-pass perfusion imaging of the lower extremities which demands very rapid volume coverage. Using a variable-density stack-of-spirals trajectory, 3D images were acquired at a temporal resolution of 2.8 sec over a large volume with a 2.5 x 2.5 x 8 mm(3) spatial resolution. These images were used to resolve the time-course of muscle intensity following contrast injection.     

Direct comparison of 3D spiral vs. Cartesian gradient-echo coronary magnetic resonance angiography.

While 3D thin-slab coronary magnetic resonance angiography (MRA) has traditionally been performed using a Cartesian acquisition scheme, spiral k-space data acquisition offers several potential advantages. However, these strategies have not been directly compared in the same subjects using similar methodologies. Thus, in the present study a comparison was made between 3D coronary MRA using Cartesian segmented k-space gradient-echo and spiral k-space data acquisition schemes. In both approaches the same spatial resolution was used and data were acquired during free breathing using navigator gating and prospective slice tracking. Magnetization preparation (T(2) preparation and fat suppression) was applied to increase the contrast. For spiral imaging two different examinations were performed, using one or two spiral interleaves, during each R-R interval. Spiral acquisitions were found to be superior to the Cartesian scheme with respect to the signal-to-noise ratio (SNR) and contrast-to-noise-ratio (CNR) (both P < 0.001) and image quality. The single spiral per R-R interval acquisition had the same total scan duration as the Cartesian acquisition, but the single spiral had the best image quality and a 2.6-fold increase in SNR. The double-interleaf spiral approach showed a 50% reduction in scanning time, a 1.8-fold increase in SNR, and similar image quality when compared to the standard Cartesian approach. Spiral 3D coronary MRA appears to be preferable to the Cartesian scheme. The increase in SNR may be "traded" for either shorter scanning times using multiple consecutive spiral interleaves, or for enhanced spatial resolution.     

Reconstruction of undersampled non-Cartesian data sets using pseudo-Cartesian GRAPPA in conjunction with GROG.

Most k-space-based parallel imaging reconstruction techniques, such as Generalized Autocalibrating Partially Parallel Acquisitions (GRAPPA), necessitate the acquisition of regularly sampled Cartesian k-space data to reconstruct a nonaliased image efficiently. However, non-Cartesian sampling schemes offer some inherent advantages to the user due to their better coverage of the center of k-space and faster acquisition times. On the other hand, these sampling schemes have the disadvantage that the points acquired generally do not lie on a grid and have complex k-space sampling patterns. Thus, the extension of Cartesian GRAPPA to non-Cartesian sequences is nontrivial. This study introduces a simple, novel method for performing Cartesian GRAPPA reconstructions on undersampled non-Cartesian k-space data gridded using GROG (GRAPPA Operator Gridding) to arrive at a nonaliased image. Because the undersampled non-Cartesian data cannot be reconstructed using a single GRAPPA kernel, several Cartesian patterns are selected for the reconstruction. This flexibility in terms of both the appearance and number of patterns allows this pseudo-Cartesian GRAPPA to be used with undersampled data sets acquired with any non-Cartesian trajectory. The successful implementation of the reconstruction algorithm using several different trajectories, including radial, rosette, spiral, one-dimensional non-Cartesian, and zig-zag trajectories, is demonstrated.     

Fast Spiral two-point Dixon technique using block regional off-resonance correction.

The Spiral two-point Dixon (Spiral 2PD) technique has recently been proposed as a method for unambiguous water-fat decomposition in spiral imaging. It also corrects for off-resonance blurring artifacts using only two data sets. In the Spiral 2PD technique, several predetermined off-resonance frequencies are tested to both separate water and fat signals and deblur the decomposed images. Unfortunately, the algorithm is computationally quite intensive since the range of tested frequencies must be set sufficiently large to span the full range of anticipated B(0) variation over the scanned objects. The block regional off-resonance correction (BRORC) algorithm corrects for off-resonance blurring artifacts block by block through the reconstructed image and usually provides several times higher computational efficiency than the conventional frequency-segmented off-resonance correction algorithm. This work shows that both water-fat decomposition and blurring artifact correction can be performed block by block using two spiral images with different TEs and that this new technique (BRORC-Spiral2PD technique) significantly improves the computational efficiency of other Spiral 2PD algorithms, opening new opportunities for spiral imaging.     

Double average parallel steady-state free precession imaging: optimized eddy current and transient oscillation compensation.

Balanced steady-state free precession (SSFP) imaging is sensitive to off-resonance effects, which can lead to considerable artifacts during a transient phase following magnetization preparation or steady-state interruption. In addition, nonlinear k-space encoding is required if contrast-relevant k-space regions need to be acquired at specific delays following magnetization preparation or for transient artifact reduction in cardiac-gated k-space segmented CINE imaging. Such trajectories are problematic for balanced SSFP imaging due to nonconstant eddy current effects and resulting disruption of the steady state.In this work, a novel acquisition strategy for balanced SSFP imaging is presented that utilizes scan time reduction by parallel imaging for optimized "double average" eddy current compensation and artifact reduction during the transient phase following steady-state storage and magnetization preparation. Double average parallel SSFP imaging was applied to k-space segmented CINE SSFP tagging as well as nongated centrically encoded SSFP imaging. Phantom and human studies exhibit substantial reduction in steady-state storage and eddy current artifacts while maintaining spatial resolution, signal-to-noise ratio, and similar total scan time of a standard SSFP acquisition. The proposed technique can easily be extended to other acquisition schemes that would benefit from nonlinear reordering schemes and/or rely on interruption of the balanced SSFP steady state.