Improved motion correction capabilities for fast spin echo T1 FLAIR propeller using non‐Cartesian external calibration data driven parallel imaging

Patient motion is a common challenge in the clinical setting and fast spin echo longitudinal relaxation time fluid attenuating inversion recovery imaging method with motion correction would be highly desirable. The motion correction provided by transverse relaxation time‐ and diffusion‐weighted periodically rotated overlapping parallel lines with enhanced reconstruction methods has seen significant clinical adoption. However, periodically rotated overlapping parallel lines with enhanced reconstruction with fast spin echo longitudinal relaxation time fluid attenuating inversion recovery‐weighting has proved challenging since motion correction requires wide blades that are difficult to acquire while also maintaining short echo train lengths that are optimal for longitudinal relaxation time fluid attenuating inversion recovery‐weighting. Parallel imaging provides an opportunity to increase the effective blade width for a given echo train lengths. Coil‐by‐coil data‐driven autocalibrated parallel imaging methods provide greater robustness in the event of motion compared to techniques relying on accurate coil sensitivity maps. However, conventional internally calibrated data‐driven parallel imaging methods limit the effective acceleration possible for each blade. We present a method to share a single calibration dataset over all imaging blades on a slice by slice basis using the APPEAR non‐Cartesian parallel imaging method providing an effective blade width increase of 2.45×, enabling robust motion correction. Results comparing the proposed technique to conventional Cartesian and periodically rotated overlapping parallel lines with enhanced reconstruction methods demonstrate a significant improvement during subject motion and maintaining high image quality when no motion is present in normal and clinical volunteers. Magn Reson Med, 2012. © 2012 Wiley Periodicals, Inc.     

High‐resolution human diffusion tensor imaging using 2‐D navigated multishot SENSE EPI at 7 T

The combination of parallel imaging with partial Fourier acquisition has greatly improved the performance of diffusion‐weighted single‐shot EPI and is the preferred method for acquisitions at low to medium magnetic field strength such as 1.5 or 3 T. Increased off‐resonance effects and reduced transverse relaxation times at 7 T, however, generate more significant artifacts than at lower magnetic field strength and limit data acquisition. Additional acceleration of k‐space traversal using a multishot approach, which acquires a subset of k‐space data after each excitation, reduces these artifacts relative to conventional single‐shot acquisitions. However, corrections for motion‐induced phase errors are not straightforward in accelerated, diffusion‐weighted multishot EPI because of phase aliasing. In this study, we introduce a simple acquisition and corresponding reconstruction method for diffusion‐weighted multishot EPI with parallel imaging suitable for use at high field. The reconstruction uses a simple modification of the standard sensitivity‐encoding (SENSE) algorithm to account for shot‐to‐shot phase errors; the method is called image reconstruction using image‐space sampling function (IRIS). Using this approach, reconstruction from highly aliased in vivo image data using 2‐D navigator phase information is demonstrated for human diffusion‐weighted imaging studies at 7 T. The final reconstructed images show submillimeter in‐plane resolution with no ghosts and much reduced blurring and off‐resonance artifacts. Magn Reson Med, 2013. © 2012 Wiley Periodicals, Inc.     

Real‐time optical motion correction for diffusion tensor imaging

Head motion is a fundamental problem in brain MRI. The problem is further compounded in diffusion tensor imaging because of long acquisition times, and the sensitivity of the tensor computation to even small misregistration. To combat motion artifacts in diffusion tensor imaging, a novel real‐time prospective motion correction method was introduced using an in‐bore monovision system. The system consists of a camera mounted on the head coil and a self‐encoded checkerboard marker that is attached to the patient's forehead. Our experiments showed that optical prospective motion correction is more effective at removing motion artifacts compared to image‐based retrospective motion correction. Statistical analysis revealed a significant improvement in similarity between diffusion data acquired at different time points when prospective correction was used compared to retrospective correction (P < 0.001). Magn Reson Med, 2010. © 2011 Wiley‐Liss, Inc.     

Partial fourier partially parallel imaging.

The techniques of partial Fourier (PF) and partially parallel imaging have been combined using a constrained reconstruction technique. The benefits compared with the individual techniques are reduced imaging time and/or an increase in signal-to-noise ratio. Low-resolution phase maps and coil sensitivities may be obtained using autocalibration or from a prescan followed by additional processing. Minor phase artifacts that are introduced by relying on conjugate symmetry can be reduced using a novel regularization scheme to vary the degree to which PF is used in the reconstruction. A nonrectilinear reconstruction algorithm is presented and the potential for motion artifact reduction is investigated using robust reconstruction.     

3D radial sampling and 3D affine transform‐based respiratory motion correction technique for free‐breathing whole‐heart coronary MRA with 100% imaging efficiency

The navigator gating and slice tracking approach currently used for respiratory motion compensation during free‐breathing coronary magnetic resonance angiography (MRA) has low imaging efficiency (typically 30–50%), resulting in long imaging times. In this work, a novel respiratory motion correction technique with 100% scan efficiency was developed for free‐breathing whole‐heart coronary MRA. The navigator signal was used as a reference respiratory signal to segment the data into six bins. 3D projection reconstruction k‐space sampling was used for data acquisition and enabled reconstruction of low resolution images within each respiratory bin. The motion between bins was estimated by image registration with a 3D affine transform. The data from the different respiratory bins was retrospectively combined after motion correction to produce the final image. The proposed method was compared with a traditional navigator gating approach in nine healthy subjects. The proposed technique acquired whole‐heart coronary MRA with 1.0 mm³ isotropic spatial resolution in a scan time of 6.8 ± 0.9 min, compared with 16.2 ± 2.8 min for the navigator gating approach. The image quality scores, and length, diameter and sharpness of the right coronary artery (RCA), left anterior descending coronary artery (LAD), and left circumflex coronary artery (LCX) were similar for both approaches (P > 0.05 for all), but the proposed technique reduced scan time by a factor of 2.5. Magn Reson Med, 2011. © 2011 Wiley‐Liss, Inc.     

Image reconstruction: An overview for clinicians

Image reconstruction plays a critical role in the clinical use of magnetic resonance imaging (MRI). The MRI raw data is not acquired in image space and the role of the image reconstruction process is to transform the acquired raw data into images that can be interpreted clinically. This process involves multiple signal processing steps that each have an impact on the image quality. This review explains the basic terminology used for describing and quantifying image quality in terms of signal‐to‐noise ratio and point spread function. In this context, several commonly used image reconstruction components are discussed. The image reconstruction components covered include noise prewhitening for phased array data acquisition, interpolation needed to reconstruct square pixels, raw data filtering for reducing Gibbs ringing artifacts, Fourier transforms connecting the raw data with image space, and phased array coil combination. The treatment of phased array coils includes a general explanation of parallel imaging as a coil combination technique. The review is aimed at readers with no signal processing experience and should enable them to understand what role basic image reconstruction steps play in the formation of clinical images and how the resulting image quality is described . J. Magn. Reson. Imaging 2014 . © 2014 Wiley Periodicals, Inc. J. Magn. Reson. Imaging 2015;41:573–585. © 2014 Wiley Periodicals, Inc.     

Independent slab‐phase modulation combined with parallel imaging in bilateral breast MRI

Independent slab‐phase modulation allows three‐dimensional imaging of multiple volumes without encoding the space between volumes, thus reducing scan time. Parallel imaging further accelerates data acquisition by exploiting coil sensitivity differences between volumes. This work compared bilateral breast image quality from self‐calibrated parallel imaging reconstruction methods such as modified sensitivity encoding, generalized autocalibrating partially parallel acquisitions and autocalibrated reconstruction for Cartesian sampling (ARC) for data with and without slab‐phase modulation. A study showed an improvement of image quality by incorporating slab‐phase modulation. Geometry factors measured from phantom images were more homogenous and lower on average when slab‐phase modulation was used for both mSENSE and GRAPPA reconstructions. The resulting improved signal‐to‐noise ratio (SNR) was validated for in vivo images as well using ARC instead of GRAPPA, illustrating average SNR efficiency increases in mSENSE by 5% and ARC by 8% based on region of interest analysis. Furthermore, aliasing artifacts from mSENSE reconstruction were reduced when slab‐phase modulation was used. Overall, slab‐phase modulation with parallel imaging improved image quality and efficiency for 3D bilateral breast imaging. Magn Reson Med, 2009. © 2009 Wiley‐Liss, Inc.     

Diffusion tensor imaging (DTI) with retrospective motion correction for large-scale pediatric imaging

Purpose:

To develop and implement a clinical DTI technique suitable for the pediatric setting that retrospectively corrects for large motion without the need for rescanning and/or reacquisition strategies, and to deliver high‐quality DTI images (both in the presence and absence of large motion) using procedures that reduce image noise and artifacts.

Materials and Methods:

We implemented an in‐house built generalized autocalibrating partially parallel acquisitions (GRAPPA)‐accelerated diffusion tensor (DT) echo‐planar imaging (EPI) sequence at 1.5T and 3T on 1600 patients between 1 month and 18 years old. To reconstruct the data, we developed a fully automated tailored reconstruction software that selects the best GRAPPA and ghost calibration weights; does 3D rigid‐body realignment with importance weighting; and employs phase correction and complex averaging to lower Rician noise and reduce phase artifacts. For select cases we investigated the use of an additional volume rejection criterion and b‐matrix correction for large motion.

Results:

The DTI image reconstruction procedures developed here were extremely robust in correcting for motion, failing on only three subjects, while providing the radiologists high‐quality data for routine evaluation.

Conclusion:

This work suggests that, apart from the rare instance of continuous motion throughout the scan, high‐quality DTI brain data can be acquired using our proposed integrated sequence and reconstruction that uses a retrospective approach to motion correction. In addition, we demonstrate a substantial improvement in overall image quality by combining phase correction with complex averaging, which reduces the Rician noise that biases noisy data. J. Magn. Reson. Imaging 2012;36:961–971. © 2012 Wiley Periodicals, Inc.     

Noise and artifact comparison for Fourier and polynomial phase correction used with Fourier reconstruction of asymmetric data sets.

Two processes of phase correction, Fourier phase mapping and second-order polynomial phase fitting, are compared in combination with four different schemes for Fourier reconstruction of asymmetric data, using one-dimensional simulations and two-dimensional human head data. Polynomial phase correction provides systematically less image noise and is much less affected by localized phase differences caused by object edges and motion.     

Combination of compressed sensing and parallel imaging for highly accelerated first‐pass cardiac perfusion MRI

First‐pass cardiac perfusion MRI is a natural candidate for compressed sensing acceleration since its representation in the combined temporal Fourier and spatial domain is sparse and the required incoherence can be effectively accomplished by k‐t random undersampling. However, the required number of samples in practice (three to five times the number of sparse coefficients) limits the acceleration for compressed sensing alone. Parallel imaging may also be used to accelerate cardiac perfusion MRI, with acceleration factors ultimately limited by noise amplification. In this work, compressed sensing and parallel imaging are combined by merging the k‐t SPARSE technique with sensitivity encoding (SENSE) reconstruction to substantially increase the acceleration rate for perfusion imaging. We also present a new theoretical framework for understanding the combination of k‐t SPARSE with SENSE based on distributed compressed sensing theory. This framework, which identifies parallel imaging as a distributed multisensor implementation of compressed sensing, enables an estimate of feasible acceleration for the combined approach. We demonstrate feasibility of 8‐fold acceleration in vivo with whole‐heart coverage and high spatial and temporal resolution using standard coil arrays. The method is relatively insensitive to respiratory motion artifacts and presents similar temporal fidelity and image quality when compared to Generalized autocalibrating partially parallel acquisitions (GRAPPA) with 2‐fold acceleration. Magn Reson Med, 2010. © 2010 Wiley‐Liss, Inc.     

Simultaneous phase correction and SENSE reconstruction for navigated multi-shot DWI with non-cartesian k-space sampling.

Phase-navigated multi-shot acquisition and parallel imaging are two techniques that have been applied to diffusion-weighted imaging (DWI) to diminish distortions and to enhance spatial resolution. Specifically, sensitivity encoding (SENSE) has been combined with single-shot echo planar imaging (EPI). Thus far, it has been difficult to apply parallel imaging methods, like SENSE, to multi-shot DWI because motion-induced phase error varies from shot to shot and interferes with sensitivity encoding. Although direct phase subtraction methods have been introduced to correct this phase error, they generally are not suitable for SENSE reconstruction, and they cannot remove all the motion artifacts even if the phase error is fully known. Here, an effective algorithm is proposed to correct the motion-induced phase error using an iterative reconstruction. In this proposed conjugate-gradient (CG) algorithm, the phase error is treated as an image encoding function. Given the complex perturbation terms, diffusion-weighted images can be reconstructed using an augmented sensitivity map. The mathematical formulation and image reconstruction procedures of this algorithm are similar to the SENSE reconstruction. By defining a dynamic composite sensitivity, the CG phase correction method can be conveniently incorporated with SENSE reconstruction for the application of multi-shot SENSE DWI. Effective phase correction and multi-shot SENSE DWI (R = 1 to 3) are demonstrated on both simulated and in vivo data acquired with PROPELLER and SNAILS.     

Motion correction using coil arrays (MOCCA) for free‐breathing cardiac cine MRI

In this study, we present a motion correction technique using coil arrays (MOCCA) and evaluate its application in free‐breathing respiratory self‐gated cine MRI. Motion correction technique using coil arrays takes advantages of the fact that motion‐induced changes in k‐space signal are modulated by individual coil sensitivity profiles. In the proposed implementation of motion correction technique using coil arrays self‐gating for free‐breathing cine MRI, the k‐space center line is acquired at the beginning of each k‐space segment for each cardiac cycle with 4 repetitions. For each k‐space segment, the k‐space center line acquired immediately before was used to select one of the 4 acquired repetitions to be included in the final self‐gated cine image by calculating the cross correlation between the k‐space center line with a reference line. The proposed method was tested on a cohort of healthy adult subjects for subjective image quality and objective blood‐myocardium border sharpness. The method was also tested on a cohort of patients to compare the left and right ventricular volumes and ejection fraction measurements with that of standard breath‐hold cine MRI. Our data indicate that the proposed motion correction technique using coil arrays method provides significantly improved image quality and sharpness compared with free‐breathing cine without respiratory self‐gating and provides similar volume measurements compared with breath‐hold cine MRI. Magn Reson Med, 2011. © 2011 Wiley‐Liss, Inc.     

Volumetric navigators for real-time motion correction in diffusion tensor imaging

A method for motion correction in multicoil imaging applications, involving both data collection and reconstruction, is presented. A bit‐reversed radial acquisition scheme, in conjunction with a rapid self‐calibrated parallel imaging method, Generalized auto‐calibrating partial parallel acquisition (GRAPPA) operator for wider radial bands (GROWL), is used to achieve motion correction at a high temporal resolution. View‐by‐view in‐plane motion correction is achieved in 2D imaging, while 3D motion correction is achieved for every two consecutive slice‐encoding planes in 3D imaging. In the proposed technique, GROWL contributes in two aspects: First, a central k‐space circle/cylinder used as the motion‐free reference is generated from a small number of radial lines/planes; Second, undersampled k‐space regions resulting from rotation and inconsistent (e.g. intraview and nonrigid body) motion can be filled in. When compared with navigator‐based motion correction methods, the proposed method does not prolong scan time and can be applied to short‐TR sequences. Magn Reson Med, 2011. © 2011 Wiley‐Liss, Inc.     

3D diffusion tensor MRI with isotropic resolution using a steady‐state radial acquisition

Purpose

To obtain diffusion tensor images (DTI) over a large image volume rapidly with 3D isotropic spatial resolution, minimal spatial distortions, and reduced motion artifacts, a diffusion‐weighted steady‐state 3D projection (SS 3DPR) pulse sequence was developed.

Materials and Methods

A diffusion gradient was inserted in a SS 3DPR pulse sequence. The acquisition was synchronized to the cardiac cycle, linear phase errors were corrected along the readout direction, and each projection was weighted by measures of consistency with other data. A new iterative parallel imaging reconstruction method was also implemented for removing off‐resonance and undersampling artifacts simultaneously.

Results

The contrast and appearance of both the fractional anisotropy and eigenvector color maps were substantially improved after all correction techniques were applied. True 3D DTI datasets were obtained in vivo over the whole brain (240 mm field of view in all directions) with 1.87 mm isotropic spatial resolution, six diffusion encoding directions in under 19 minutes.

Conclusion

A true 3D DTI pulse sequence with high isotropic spatial resolution was developed for whole brain imaging in under 20 minutes. To minimize the effects of brain motion, a cardiac synchronized, multiecho, DW‐SSFP pulse sequence was implemented. Motion artifacts were further reduced by a combination of linear phase correction, corrupt projection detection and rejection, sampling density reweighting, and parallel imaging reconstruction. The combination of these methods greatly improved the quality of 3D DTI in the brain. J. Magn. Reson. Imaging 2009;29:1175–1184. © 2009 Wiley‐Liss, Inc.     

Partially parallel imaging with localized sensitivities (PILS).

In this study a novel partially parallel acquisition method is presented, which can be used to accelerate image acquisition using an RF coil array for spatial encoding. In this technique, Parallel Imaging with Localized Sensitivities (PILS), it is assumed that the individual coils in the array have localized sensitivity patterns, in that their sensitivity is restricted to a finite region of space. Within the PILS model, a detailed, highly accurate RF field map is not needed prior to reconstruction. In PILS, each coil in the array is fully characterized by only two parameters: the center of coil's sensitive region in the FOV and the width of the sensitive region around this center. In this study, it is demonstrated that the incorporation of these coil parameters into a localized Fourier transform allows reconstruction of full FOV images in each of the component coils from data sets acquired with a reduced number of phase encoding steps compared to conventional imaging techniques. After the introduction of the PILS technique, primary focus is given to issues related to the practical implementation of PILS, including coil parameter determination and the SNR and artifact power in the resulting images. Finally, in vivo PILS images are shown which demonstrate the utility of the technique.     

Robust GRAPPA‐accelerated diffusion‐weighted readout‐segmented (RS)‐EPI

Readout segmentation (RS‐EPI) has been suggested as a promising variant to echo‐planar imaging (EPI) for high‐resolution imaging, particularly when combined with parallel imaging. This work details some of the technical aspects of diffusion‐weighted (DW)‐RS‐EPI, outlining a set of reconstruction methods and imaging parameters that can both minimize the scan time and afford high‐resolution diffusion imaging with reduced distortions. These methods include an efficient generalized autocalibrating partially parallel acquisition (GRAPPA) calibration for DW‐RS‐EPI data without scan time penalty, together with a variant for the phase correction of partial Fourier RS‐EPI data. In addition, the role of pulsatile and rigid‐body brain motion in DW‐RS‐EPI was assessed. Corrupt DW‐RS‐EPI data arising from pulsatile nonlinear brain motion had a prevalence of ～7% and were robustly identified via k‐space entropy metrics. For DW‐RS‐EPI data corrupted by rigid‐body motion, we showed that no blind overlap was required. The robustness of RS‐EPI toward phase errors and motion, together with its minimized distortions compared with EPI, enables the acquisition of exquisite 3 T DW images with matrix sizes close to 512². Magn Reson Med, 2009. © 2009 Wiley‐Liss, Inc.     

Reconstruction of MRI data encoded with arbitrarily shaped,curvilinear, nonbijective magnetic fields

A basic framework for image reconstruction from spatial encoding by curvilinear, nonbijective magnetic encoding fields in combination with multiple receivers is presented. The theory was developed in the context of the recently introduced parallel imaging technique using localized gradients (PatLoc) approach. In this new imaging modality, the linear gradient fields are generalized to arbitrarily shaped, nonbijective spatial encoding magnetic fields, which lead to ambiguous encoding. Ambiguities are resolved by adaptation of concepts developed for parallel imaging. Based on theoretical considerations, a practical algorithm for Cartesian trajectories is derived in the case that the conventional gradient coils are replaced by coils for PatLoc. The reconstruction method extends Cartesian sensitivity encoding (SENSE) reconstruction with an additional voxelwise intensity‐correction step. Spatially varying resolution, signal‐to‐noise ratio, and truncation artifacts are described and analyzed. Theoretical considerations are validated by two‐dimensional simulations based on multipolar encoding fields and they are confirmed by applying the reconstruction algorithm to initial experimental data. Magn Reson Med, 2010. © 2010 Wiley‐Liss, Inc.     

Sensitivity profiles from an array of coils for encoding and reconstruction in parallel (SPACE RIP).

A new parallel imaging technique was implemented which can result in reduced image acquisition times in MRI. MR data is acquired in parallel using an array of receiver coils and then reconstructed simultaneously with multiple processors. The method requires the initial estimation of the 2D sensitivity profile of each coil used in the receiver array. These sensitivity profiles are then used to partially encode the images of interest. A fraction of the total number of k-space lines is consequently acquired and used in a parallel reconstruction scheme, allowing for a substantial reduction in scanning and display times. This technique is in the family of parallel acquisition schemes such as simultaneous acquisition of spatial harmonics (SMASH) and sensitivity encoding (SENSE). It extends the use of the SMASH method to allow the placement of the receiver coil array around the object of interest, enabling imaging of any plane within the volume of interest. In addition, this technique permits the arbitrary choice of the set of k-space lines used in the reconstruction and lends itself to parallel reconstruction, hence allowing for real-time rendering. Simulated results with a 16-fold increase in temporal resolution are shown, as are experimental results with a 4-fold increase in temporal resolution. Magn Reson Med 44:301-308, 2000.     

Motion correction using an enhanced floating navigator and GRAPPA operations

A method for motion correction in multicoil imaging applications, involving both data collection and reconstruction, is presented. The floating navigator method, which acquires a readout line off center in the phase‐encoding direction, is expanded to detect translation/rotation and inconsistent motion. This is done by comparing floating navigator data with a reference k‐space region surrounding the floating navigator line, using a correlation measure. The technique of generalized autocalibrating partially parallel acquisition is further developed to correct for a fully sampled, motion‐corrupted dataset. The flexibility of generalized autocalibrating partially parallel acquisition kernels is exploited by extrapolating readout lines to fill in missing “pie slices” of k‐space caused by rotational motion and regenerating full k‐space data from multiple interleaved datasets, facilitating subsequent rigid‐body motion correction or proper weighting of inconsistent data (e.g., with through‐plane and nonrigid motion). Phantom and in vivo imaging experiments with turbo spin‐echo sequence demonstrate the correction of severe motion artifacts. Magn Reson Med, 2010. © 2009 Wiley‐Liss, Inc.     

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.