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.     

Accelerated volumetric MRI with a SENSE/GRAPPA combination

PURPOSE: To combine the specific advantages of the generalized autocalibrating partially parallel acquisitions (GRAPPA) technique and sensitivity encoding (SENSE) with two-dimensional (2D) undersampling. MATERIALS AND METHODS: By splitting the 2D reconstruction process into multiple one-dimensional (1D) reconstructions, the normal 1D GRAPPA method can be used for image reconstruction. Due to this data-handling process, a GRAPPA reconstruction is performed along the phase-encoding (PE) direction and effectively a SENSE reconstruction is performed along the partition-encoding (PAE) direction. RESULTS: In vivo experiments demonstrate the successful implementation of the SENSE/GRAPPA combination. Experimental results with up to 9.6-fold acceleration using a prototype 32-channel receiver head coil array are presented. CONCLUSION: The proposed SENSE/GRAPPA combination for 3D imaging allows the GRAPPA method to be applied in combination with 2D undersampling. Because the SENSE/GRAPPA combination is not based on knowledge of spatial coil sensitivities, it should be the method of choice whenever it is difficult to extract the sensitivity information.     

Robust GRAPPA reconstruction and its evaluation with the perceptual difference model

PURPOSE: To develop and optimize a new modification of GRAPPA (generalized autocalibrating partially parallel acquisitions) MR reconstruction algorithm named "Robust GRAPPA." MATERIALS AND METHODS: In Robust GRAPPA, k-space data points were weighted before the reconstruction. Small or zero weights were assigned to "outliers" in k-space. We implemented a Slow Robust GRAPPA method, which iteratively reweighted the k-space data. It was compared to an ad hoc Fast Robust GRAPPA method, which eliminated (assigned zero weights to) a fixed percentage of k-space "outliers" following an initial estimation procedure. In comprehensive experiments the new algorithms were evaluated using the perceptual difference model (PDM), whereby image quality was quantitatively compared to the reference image. Independent variables included algorithm type, total reduction factor, outlier ratio, center filling options, and noise across multiple image datasets, providing 10,800 test images for evaluation. RESULTS: The Fast Robust GRAPPA method gave results very similar to Slow Robust GRAPPA, and showed significant improvements as compared to regular GRAPPA. Fast Robust GRAPPA added little computation time compared with regular GRAPPA. CONCLUSION: Robust GRAPPA was proposed and proved useful for improving the reconstructed image quality. PDM was helpful in designing and optimizing the MR reconstruction algorithms.     

Nonlinear GRAPPA: a kernel approach to parallel MRI reconstruction

GRAPPA linearly combines the undersampled k-space signals to estimate the missing k-space signals where the coefficients are obtained by fitting to some auto-calibration signals (ACS) sampled with Nyquist rate based on the shift-invariant property. At high acceleration factors, GRAPPA reconstruction can suffer from a high level of noise even with a large number of auto-calibration signals. In this work, we propose a nonlinear method to improve GRAPPA. The method is based on the so-called kernel method which is widely used in machine learning. Specifically, the undersampled k-space signals are mapped through a nonlinear transform to a high-dimensional feature space, and then linearly combined to reconstruct the missing k-space data. The linear combination coefficients are also obtained through fitting to the ACS data but in the new feature space. The procedure is equivalent to adding many virtual channels in reconstruction. A polynomial kernel with explicit mapping functions is investigated in this work. Experimental results using phantom and in vivo data demonstrate that the proposed nonlinear GRAPPA method can significantly improve the reconstruction quality over GRAPPA and its state-of-the-art derivatives.     

k-t GRAPPA: a k-space implementation for dynamic MRI with high reduction factor.

A novel technique called "k-t GRAPPA" is introduced for the acceleration of dynamic magnetic resonance imaging. Dynamic magnetic resonance images have significant signal correlations in k-space and time dimension. Hence, it is feasible to acquire only a reduced amount of data and recover the missing portion afterward. Generalized autocalibrating partially parallel acquisitions (GRAPPA), as an important parallel imaging technique, linearly interpolates the missing data in k-space. In this work, it is shown that the idea of GRAPPA can also be applied in k-t space to take advantage of the correlations and interpolate the missing data in k-t space. For this method, no training data, filters, additional parameters, or sensitivity maps are necessary, and it is applicable for either single or multiple receiver coils. The signal correlation is locally derived from the acquired data. In this work, the k-t GRAPPA technique is compared with our implementation of GRAPPA, TGRAPPA, and sliding window reconstructions, as described in Methods. The experimental results manifest that k-t GRAPPA generates high spatial resolution reconstruction without significant loss of temporal resolution when the reduction factor is as high as 4. When the reduction factor becomes higher, there might be a noticeable loss of temporal resolution since k-t GRAPPA uses temporal interpolation. Images reconstructed using k-t GRAPPA have less residue/folding artifacts than those reconstructed by sliding window, much less noise than those reconstructed by GRAPPA, and wider temporal bandwidth than those reconstructed by GRAPPA with residual k-space. k-t GRAPPA is applicable to a wide range of dynamic imaging applications and is not limited to imaging parts with quasi-periodic motion. Since only local information is used for reconstruction, k-t GRAPPA is also preferred for applications requiring real time reconstruction, such as monitoring interventional MRI.     

Cross-validation-based kernel support selection for improved GRAPPA reconstruction.

The extended version of the generalized autocalibrating partially parallel acquisition (GRAPPA) technique incorporates multiple lines and multiple columns of measured k-space data to estimate missing data. For a given accelerated dataset, the selection of the measured data points for fitting a missing datum (i.e., the kernel support) that provides optimal reconstruction depends on coil array configuration, noise level in the acquired data, imaging configuration, and number and position of autocalibrating signal lines. In this work, cross-validation is used to select the kernel support that best balances the conflicting demands of fit accuracy and stability in GRAPPA reconstruction. The result is an optimized tradeoff between artifacts and noise. As demonstrated with experimental data, the method improves image reconstruction with GRAPPA. Because the method is simple and applied in postprocessing, it can be used with GRAPPA routinely.     

Self-calibration method for radial GRAPPA/k-t GRAPPA.

Generalized autocalibrating partially parallel acquisitions (GRAPPA), an important parallel imaging technique, can be easily applied to radial k-space data by segmenting the k-space. The previously reported radial GRAPPA method requires extra calibration data to determine the relative shift operators. In this work it is shown that pseudo-full k-space data can be generated from the partially acquired radial data by filtering in image space followed by inverse gridding. The relative shift operators can then be approximated from the pseudo-full k-space data. The self-calibration method using pseudo-full k-space data can be applied in both k and k-t space. This technique avoids the prescans and hence improves the applicability of radial GRAPPA to image static tissue, and makes k-t GRAPPA applicable to radial trajectory. Experiments show that radial GRAPPA calibrated with pseudo-full calibration data generates results similar to radial GRAPPA calibrated with the true full k-space data for that image. If motion occurs during acquisition, self-calibrated radial GRAPPA protects structural information better than externally calibrated GRAPPA. However, radial GRAPPA calibrated with pseudo-full calibration data suffers from residual streaking artifacts when the reduction factor is high. Radial k-t GRAPPA calibrated with pseudo-full calibration data generates reduced errors compared to the sliding-window method and temporal GRAPPA (TGRAPPA).     

Using GRAPPA to improve autocalibrated coil sensitivity estimation for the SENSE family of parallel imaging reconstruction algorithms

Two strategies are widely used in parallel MRI to reconstruct subsampled multicoil image data. SENSE and related methods employ explicit receiver coil spatial response estimates to reconstruct an image. In contrast, coil‐by‐coil methods such as GRAPPA leverage correlations among the acquired multicoil data to reconstruct missing k‐space lines. In self‐referenced scenarios, both methods employ Nyquist‐rate low‐frequency k‐space data to identify the reconstruction parameters. Because GRAPPA does not require explicit coil sensitivities estimates, it needs considerably fewer autocalibration signals than SENSE. However, SENSE methods allow greater opportunity to control reconstruction quality though regularization and thus may outperform GRAPPA in some imaging scenarios. Here, we employ GRAPPA to improve self‐referenced coil sensitivity estimation in SENSE and related methods using very few auto‐calibration signals. This enables one to leverage each methods' inherent strength and produce high quality self‐referenced SENSE reconstructions. Magn Reson Med 60:462–467, 2008. © 2008 Wiley‐Liss, Inc.     

Fast four-dimensional coronary MR angiography with k-t GRAPPA

PURPOSE: To investigate the effectiveness of k-t GRAPPA for accelerating four-dimensional (4D) coronary MRA in comparison with GRAPPA and the feasibility of combining variable density undersampling with conventional k-t GRAPPA (k-t(2) GRAPPA) to alleviate the overhead of acquiring autocalibration signals. MATERIALS AND METHODS: The right coronary artery of nine healthy volunteers was scanned at 1.5 Tesla. The 4D k-space datasets were fully acquired and subsequently undersampled to simulate partially parallel acquisitions, namely, GRAPPA, k-t GRAPPA, and k-t(2) GRAPPA. Comparisons were made between the images reconstructed from full k-space datasets and those reconstructed from undersampled k-space datasets. RESULTS: k-t GRAPPA significantly reduced artifacts compared with GRAPPA and high acceleration factors were achieved with only minimal sacrifices in vessel depiction. k-t(2) GRAPPA could further increase imaging speed without significant losses in image quality. CONCLUSION: By exploiting high-degree spatiotemporal correlations during the rest period of a cardiac cycle, k-t GRAPPA and k-t(2) GRAPPA can greatly increase data acquisition efficiency and, therefore, are promising solutions for fast 4D coronary MRA.     

Parallel magnetic resonance imaging using the GRAPPA operator formalism.

In this article it is shown that GRAPPA reconstruction can be reformulated as a matrix operator, similar to ladder or propagator operators used in quantum mechanics, that shifts data in k-space. Using this formalism, it is shown that there exists an infinitesimal GRAPPA operator that shifts data in k-space by arbitrarily small amounts. Other desired k-space shifts can then be accomplished through repeated applications of this infinitesimal GRAPPA operator. Implications of these ideas are described.     

2D partially parallel imaging with k-space surrounding neighbors-based data reconstruction.

Partially parallel imaging (PPI) achieves imaging acceleration by replacing partial phase encoding (PE) with the spatially localized sensitivity encoding of a receiver surface coil array. Further accelerations can be achieved through 2D PPI along two PE directions in 3D MRI. This paper is to explore the k-space-based PPI acquisition and reconstruction strategies for 3D MRI. A surrounding neighbors-based autocalibrating PPI (SNAPPI) was first presented by generalizing the 2D multicolumn multiline interpolation method. Several 2D PPI reconstruction methods were then provided by applying SNAPPI to recover the partially skipped k-space data along two PE directions separately or nonseparately, in k-space or in the hybrid k and image space. An optimal 2D PPI sampling-based reconstruction approach was also presented for applying PPI along certain spatial direction along which the array coil has not sufficient sensitivity variation for a valid PPI reconstruction. Both simulated and in vivo 2D PPI data were used to evaluate the proposed methods.     

Non-Cartesian data reconstruction using GRAPPA operator gridding (GROG).

A novel approach that uses the concepts of parallel imaging to grid data sampled along a non-Cartesian trajectory using GRAPPA operator gridding (GROG) is described. GROG shifts any acquired data point to its nearest Cartesian location, thereby converting non-Cartesian to Cartesian data. Unlike other parallel imaging methods, GROG synthesizes the net weight for a shift in any direction from a single basis set of weights along the logical k-space directions. Given the vastly reduced size of the basis set, GROG calibration and reconstruction requires fewer operations and less calibration data than other parallel imaging methods for gridding. Instead of calculating and applying a density compensation function (DCF), GROG requires only local averaging, as the reconstructed points fall upon the Cartesian grid. Simulations are performed to demonstrate that the root mean square error (RMSE) values of images gridded with GROG are similar to those for images gridded using the gold-standard convolution gridding. Finally, GROG is compared to the convolution gridding technique using data sampled along radial, spiral, rosette, and BLADE (a.k.a. periodically rotated overlapping parallel lines with enhanced reconstruction [PROPELLER]) trajectories.     

Self-calibrating GRAPPA operator gridding for radial and spiral trajectories.

Self-calibrating GRAPPA operator gridding (GROG) is a method by which non-Cartesian MRI data can be gridded using spatial information from a multichannel coil array without the need for an additional calibration dataset. Using self-calibrating GROG, the non-Cartesian datapoints are shifted to nearby k-space locations using parallel imaging weight sets determined from the datapoints themselves. GROG employs the GRAPPA Operator, a special formulation of the general reconstruction method GRAPPA, to perform these shifts. Although GROG can be used to grid undersampled datasets, it is important to note that this method uses parallel imaging only for gridding, and not to reconstruct artifact-free images from undersampled data. The innovation introduced here, namely, self-calibrating GROG, allows the shift operators to be calculated directly out of the non-Cartesian data themselves. This eliminates the need for an additional calibration dataset, which reduces the imaging time and also makes the GROG reconstruction more robust by removing possible inconsistencies between the calibration and non-Cartesian datasets. Simulated and in vivo examples of radial and spiral datasets gridded using self-calibrating GROG are compared to images gridded using the standard method of convolution gridding.     

Artifact reduction in moving-table acquisitions using parallel imaging and multiple averages.

Two-dimensional (2D) axial continuously-moving-table imaging has to deal with artifacts due to gradient nonlinearity and breathing motion, and has to provide the highest scan efficiency. Parallel imaging techniques (e.g., generalized autocalibrating partially parallel acquisition GRAPPA)) are used to reduce such artifacts and avoid ghosting artifacts. The latter occur in T(2)-weighted multi-spin-echo (SE) acquisitions that omit an additional excitation prior to imaging scans for presaturation purposes. Multiple images are reconstructed from subdivisions of a fully sampled k-space data set, each of which is acquired in a single SE train. These images are then averaged. GRAPPA coil weights are estimated without additional measurements. Compared to conventional image reconstruction, inconsistencies between different subsets of k-space induce less artifacts when each k-space part is reconstructed separately and the multiple images are averaged afterwards. These inconsistencies may lead to inaccurate GRAPPA coil weights using the proposed intrinsic GRAPPA calibration. It is shown that aliasing artifacts in single images are canceled out after averaging. Phantom and in vivo studies demonstrate the benefit of the proposed reconstruction scheme for free-breathing axial continuously-moving-table imaging using fast multi-SE sequences.     

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.     

Iterative GRAPPA (iGRAPPA) for improved parallel imaging reconstruction.

In this work an iterative reconstruction method based on generalized autocalibrating partially parallel acquisitions (GRAPPA) reconstruction is introduced. In the new method the reconstructed lines are used to reestimate and refine the weights from all the acquired data by applying the GRAPPA procedure iteratively with regularization. Both phantom and in vivo MRI experiments demonstrated that, compared to GRAPPA, the iterative approach reduces parallel imaging artifacts and permits high-quality image reconstruction with a relatively small number of calibration lines and slight changes of GRAPPA weights.     

Highly k‐t‐space–accelerated phase‐contrast MRI

The purpose of this study was to combine a recently introduced spatiotemporal parallel imaging technique, PEAK‐GRAPPA (parallel MRI with extended and averaged generalized autocalibrating partially parallel acquisition), with two‐dimensional (2D) cine phase‐contrast velocity mapping. Phase‐contrast MRI was applied to measure the blood flow in the thoracic aorta and the myocardial motion of the left ventricle. To evaluate the performance of different reconstruction methods, fully acquired k‐space data sets were used to compare conventional parallel imaging using GRAPPA with reduction factors of R = 2–6 and PEAK‐GRAPPA as well as sliding window reconstruction with reduction factors R = 2–12 (net acceleration factors up to 5.2). PEAK‐GRAPPA reconstruction resulted in improved image quality with considerably reduced artifacts, which was also supported by error analysis. To analyze potential blurring or low‐pass filtering effects of spatiotemporal PEAK‐GRAPPA, the velocity time courses of aortic flow and myocardial tissue motion were evaluated and compared with conventional image reconstructions. Quantitative comparisons of blood flow velocities and pixel‐wise correlation analysis of velocities highlight the potential of PEAK‐GRAPPA for highly accelerated dynamic phase‐contrast velocity mapping. Magn Reson Med 60:1169–1177, 2008. © 2008 Wiley‐Liss, Inc.     

Parallel MRI with extended and averaged GRAPPA kernels (PEAK-GRAPPA): optimized spatiotemporal dynamic imaging

Purpose

To evaluate an optimized k‐t‐space related reconstruction method for dynamic magnetic resonance imaging (MRI), a method called PEAK‐GRAPPA (Parallel MRI with Extended and Averaged GRAPPA Kernels) is presented which is based on an extended spatiotemporal GRAPPA kernel in combination with temporal averaging of coil weights.

Materials and Methods

The PEAK‐GRAPPA kernel consists of a uniform geometry with several spatial and temporal source points from acquired k‐space lines and several target points from missing k‐space lines. In order to improve the quality of coil weight estimation sets of coil weights are averaged over the temporal dimension.

Results

The kernel geometry leads to strongly decreased reconstruction times compared to the recently introduced k‐t‐GRAPPA using different kernel geometries with only one target point per kernel to fit. Improved results were obtained in terms of the root mean square error and the signal‐to‐noise ratio as demonstrated by in vivo cardiac imaging.

Conclusion

Using a uniform kernel geometry for weight estimation with the properties of uncorrelated noise of different acquired timeframes, optimized results were achieved in terms of error level, signal‐to‐noise ratio, and reconstruction time. J. Magn. Reson. Imaging 2008;28:1226–1232. © 2008 Wiley‐Liss, Inc.     

Fast method for 1D non-cartesian parallel imaging using GRAPPA.

MRI with non-Cartesian sampling schemes can offer inherent advantages. Radial acquisitions are known to be very robust, even in the case of vast undersampling. This is also true for 1D non-Cartesian MRI, in which the center of k-space is oversampled or at least sampled at the Nyquist rate. There are two main reasons for the more relaxed foldover artifact behavior: First, due to the oversampling of the center, high-energy foldover artifacts originating from the center of k-space are avoided. Second, due to the non-equidistant sampling of k-space, the corresponding field of view (FOV) is no longer well defined. As a result, foldover artifacts are blurred over a broad range and appear less severe. The more relaxed foldover artifact behavior and the densely sampled central k-space make trajectories of this type an ideal complement to autocalibrated parallel MRI (pMRI) techniques, such as generalized autocalibrating partially parallel acquisitions (GRAPPA). Although pMRI can benefit from non-Cartesian trajectories, this combination has not yet entered routine clinical use. One of the main reasons for this is the need for long reconstruction times due to the complex calculations necessary for non-Cartesian pMRI. In this work it is shown that one can significantly reduce the complexity of the calculations by exploiting a few specific properties of k-space-based pMRI.     

Discrepancy-based adaptive regularization for GRAPPA reconstruction

PURPOSE: To develop a novel regularization method for GRAPPA by which the regularization parameters can be optimally and adaptively chosen. MATERIALS AND METHODS: In the fit procedures in GRAPPA, the discrepancy principle, which chooses the regularization parameter based on a priori information about the noise level in the autocalibrating signals (ACS), is used with the truncated singular value decomposition (TSVD) regularization and the Tikhonov regularization, and its performance is compared with the singular value (SV) threshold method and the L-curve method, respectively by axial and sagittal head imaging experiments. RESULTS: In both axial and sagittal reconstructions, normal GRAPPA reconstruction results exhibit a relatively high level of noise. With discrepancy-based choices of parameters, regularization can improve the signal-to-noise ratio (SNR) with only a very modest increase in aliasing artifacts. The L-curve method in all of the reconstructions leads to overregularization, which causes severe residual aliasing artifacts. The 10% SV threshold method yields good overall image quality in the axial case, but in the sagittal case it also leads to an obvious increase in aliasing artifacts. CONCLUSION: Neither a fixed SV threshold nor the L-curve are robust means of choosing the appropriate parameters in GRAPPA reconstruction. However, with the discrepancy-based parameter-choice strategy, adaptively regularized GRAPPA can be used to automatically choose nearly optimal parameters for reconstruction and achieve an excellent compromise between SNR and artifacts.