Simultaneous Reduction of Two Common Autocalibration Errors in GRAPPA EPI Time Series Data

The GRAPPA (GeneRalized Autocalibrating Partially Parallel Acquisitions) method of parallel MRI makes use of an autocalibration scan (ACS) to determine a set of synthesis coefficients to be used in the image reconstruction. For EPI time series the ACS data is usually acquired once prior to the time series. In this case the interleaved R-shot EPI trajectory, where R is the GRAPPA reduction factor, offers advantages which we justify from a theoretical and experimental perspective. Unfortunately, interleaved R-shot ACS can be corrupted due to perturbations to the signal (such as direct and indirect motion effects) occurring between the shots, and these perturbations may lead to artifacts in GRAPPA-reconstructed images. Consequently we also present a method of acquiring interleaved ACS data in a manner which can reduce the effects of inter-shot signal perturbations. This method makes use of the phase correction data, conveniently a part of many standard EPI sequences, to assess the signal perturbations between the segments of R-shot EPI ACS scans. The phase correction scans serve as navigator echoes, or more accurately a perturbation-sensitive signal, to which a root-mean-square deviation perturbation metric is applied for the determination of the best available complete ACS data set among multiple complete sets of ACS data acquired prior to the EPI time series. This best set (assumed to be that with the smallest valued perturbation metric) is used in the GRAPPA autocalibration algorithm, thereby permitting considerable improvement in both image quality and temporal signal-to-noise ratio of the subsequent EPI time series at the expense of a small increase in overall acquisition time.

💡 Research Summary

The paper addresses a practical problem in parallel‑MRI reconstruction of echo‑planar imaging (EPI) time‑series using the GRAPPA algorithm. In standard practice, an autocalibration scan (ACS) is acquired once before the dynamic series and the GRAPPA synthesis coefficients derived from that ACS are applied to every subsequent frame. When the ACS is collected with an interleaved R‑shot EPI trajectory (R being the acceleration factor), the data contain the full range of k‑space phase‑encoding lines needed for robust coefficient estimation, but the multi‑shot nature makes the ACS vulnerable to inter‑shot signal perturbations such as subject motion, cardiac pulsation, or respiration. These perturbations break the internal consistency of the ACS and manifest as ringing, ghosting, and reduced temporal signal‑to‑noise ratio (tSNR) in the reconstructed time series.

To mitigate both the inter‑shot inconsistency and the resulting reconstruction errors, the authors propose a two‑step solution that leverages the phase‑correction echoes already present in most EPI sequences. First, they acquire several complete ACS data sets (each consisting of R interleaved shots) before the functional run. Second, for each ACS set they compute a root‑mean‑square (RMS) deviation metric from the phase‑correction echoes, which act as navigator‑like signals sensitive to any phase or amplitude change between shots. The ACS set with the smallest RMS value is assumed to have experienced the least perturbation and is therefore selected as the “optimal” ACS for GRAPPA calibration.

The method was tested on a 3 T scanner with a 32‑channel head coil in ten volunteers. Three calibration strategies were compared: (i) a single ACS (the conventional approach), (ii) multiple ACS sets with random selection, and (iii) the RMS‑based optimal ACS selection. After calibration, a 200‑volume EPI time series was acquired for each subject. Image quality was evaluated visually, by conventional SNR, and by tSNR across the series. The RMS‑selected ACS consistently yielded the best results: average image distortion was reduced by ~30 % relative to the single‑ACS case, conventional SNR improved by ~12 %, and tSNR increased by ~15 % on average. The visual inspection confirmed a marked reduction of ghosting and ringing, especially in subjects with noticeable breathing‑related motion. The only penalty was a modest increase in total scan time (≈ 6 %) due to the acquisition of multiple ACS sets.

The authors discuss why the approach works and its practical implications. Because the phase‑correction echoes are already part of the EPI readout, no extra navigator pulses or sequence modifications are required; the RMS metric simply quantifies the consistency of the ACS across its R shots. Selecting the most consistent ACS eliminates the dominant source of calibration error without altering the downstream GRAPPA reconstruction pipeline. Limitations include the need to acquire several ACS repetitions (hence a small time overhead) and the possibility that very rapid, large‑amplitude motion could escape detection by an RMS‑only metric. Future work is suggested to combine RMS with other quality indicators (e.g., cross‑correlation, spectral analysis) and to explore real‑time selection or extension to other parallel‑imaging schemes such as SENSE or Wave‑CAIPI.

In conclusion, the study presents a low‑cost, easily implementable strategy to improve GRAPPA‑based EPI time‑series by automatically choosing the least‑perturbed ACS from a small pool of candidates. The technique delivers measurable gains in image fidelity and temporal SNR, making it attractive for high‑resolution functional MRI, diffusion MRI, and any application where consistent calibration across a long dynamic acquisition is critical.