New Theory and Faster Computations for Subspace-Based Sensitivity Map Estimation in Multichannel MRI.

Sensitivity map estimation is important in many multichannel MRI applications. Subspace-based sensitivity map estimation methods like ESPIRiT are popular and perform well, though can be computationally expensive and their theoretical principles can be nontrivial to understand. In the first part of this work, we present a novel theoretical derivation of subspace-based sensitivity map estimation based on a linear-predictability/structured low-rank modeling perspective. This results in an estimation approach that is equivalent to ESPIRiT, but with distinct theory that may be more intuitive for some readers. In the second part of this work, we propose and evaluate a set of computational acceleration approaches (collectively known as PISCO) that can enable substantial improvements in computation time (up to  âˆ¼ 100× in the examples we show) and memory for subspace-based sensitivity map estimation.

Single breath-hold CINE imaging with combined simultaneous multislice and region-optimized virtual coils.

PURPOSE: To investigate the feasibility of combining simultaneous multislice (SMS) and region-optimized virtual coils (ROVir) for single breath-hold CINE imaging. METHOD: ROVir is a recent virtual coil approach that allows reduced-field of view (FOV) imaging by localizing the signal from a region-of-interest (ROI) and/or suppressing the signal from unwanted spatial regions. In this work, ROVir is used for reduced-FOV SMS bSSFP CINE imaging, which enables whole heart CINE with a single breath-hold acquisition. RESULTS: Reduced-FOV CINE with either SMS-only or ROVir-only resulted in significant aliasing, with severely reduced image quality when compared to the full FOV reference CINE, while the visual appearance of aliasing was substantially reduced with the proposed SMS+ROVir. The end diastolic volume, end systolic volume, and ejection fraction obtained using the proposed approach were similar to the clinical reference (correlations of 0.92, 0.94, and 0.88, respectively with p < 0 . 05 $$ p<0.05 $$ in each case, and biases of 0.1, 1.6 mL, and - 0 . 6 % $$ -0.6\% $$ , respectively). No statistically significant differences for these parameters were found with a Wilcoxon rank test (p = 0.96, 0.20, and 0.40, respectively). CONCLUSION: We demonstrated that reduced-FOV CINE imaging with SMS+ROVir enables single breath-hold whole-heart imaging without compromising visual image quality or quantitative cardiac function parameters.

High-fidelity, high-spatial-resolution diffusion magnetic resonance imaging of ex vivo whole human brain at ultra-high gradient strength with structured low-rank echo-planar imaging ghost correction.

Diffusion magnetic resonance imaging (dMRI) of whole ex vivo human brain specimens enables three-dimensional (3D) mapping of structural connectivity at the mesoscopic scale, providing detailed evaluation of fiber architecture and tissue microstructure at a spatial resolution that is difficult to access in vivo. To account for the short T2 and low diffusivity of fixed tissue, ex vivo dMRI is often acquired using strong diffusion-sensitizing gradients and multishot/segmented 3D echo-planar imaging (EPI) sequences to achieve high spatial resolution. However, the combination of strong diffusion-sensitizing gradients and multishot/segmented EPI readout can result in pronounced ghosting artifacts incurred by nonlinear spatiotemporal variations in the magnetic field produced by eddy currents. Such ghosting artifacts cannot be corrected with conventional correction solutions and pose a significant roadblock to leveraging human MRI scanners with ultrahigh gradients for ex vivo whole-brain dMRI. Here, we show that ghosting-correction approaches that correct for either polarity-related ghosting or shot-to-shot variations in a separate manner are suboptimal for 3D multishot diffusion-weighted EPI experiments in fixed human brain specimens using strong diffusion-sensitizing gradients on the 3-T Connectom MRI scanner, resulting in orientationally biased dMRI estimates. We apply a recently developed advanced k-space reconstruction method based on structured low-rank matrix (SLM) modeling that handles both polarity-related ghosting and shot-to-shot variation simultaneously, to mitigate artifacts in high-angular resolution multishot dMRI data acquired in several fixed human brain specimens at 0.7-0.8-mm isotropic spatial resolution using b-values up to 10,000 s/mm2 and gradient strengths up to 280 mT/m. We demonstrate the improved mapping of diffusion tensor imaging and fiber orientation distribution functions in key neuroanatomical areas distributed across the whole brain using SLM-based EPI ghost correction compared with alternative techniques.

On the shape of convolution kernels in MRI reconstruction: Rectangles versus ellipsoids.

PURPOSE: Many MRI reconstruction methods (including GRAPPA, SPIRiT, ESPIRiT, LORAKS, and convolutional neural network [CNN] methods) involve shift-invariant convolution models. Rectangular convolution kernel shapes are often chosen by default, although ellipsoidal kernel shapes have potentially appealing theoretical characteristics. In this work, we systematically investigate the differences between different kernel shape choices in several contexts. THEORY: It is well-understood that a rectangular region of k-space is associated with anisotropic spatial resolution, while ellipsoidal regions can be associated with more isotropic resolution. Further, for a fixed spatial resolution, ellipsoidal kernels are associated with substantially fewer parameters than rectangular kernels. These characteristics suggest that ellipsoidal kernels may have certain advantages over rectangular kernels. METHODS: We used real retrospectively undersampled k-space data to empirically study the characteristics of rectangular and ellipsoidal kernels in the context of seven methods (GRAPPA, SPIRiT, ESPIRiT, SAKE, LORAKS, AC-LORAKS, and CNN-based reconstructions). RESULTS: Empirical results suggest that both kernel shapes can produce reconstructed images with similar error metrics, although the ellipsoidal shape can often achieve this with reduced computation time and memory usage and/or fewer model parameters. CONCLUSION: Ellipsoidal kernel shapes may offer advantages over rectangular kernel shapes in various MRI applications.

Robust autocalibrated structured low-rank EPI ghost correction.

PURPOSE: We propose and evaluate a new structured low-rank method for echo-planar imaging (EPI) ghost correction called Robust Autocalibrated LORAKS (RAC-LORAKS). The method can be used to suppress EPI ghosts arising from the differences between different readout gradient polarities and/or the differences between different shots. It does not require conventional EPI navigator signals, and is robust to imperfect autocalibration data. METHODS: Autocalibrated LORAKS is a previous structured low-rank method for EPI ghost correction that uses GRAPPA-type autocalibration data to enable high-quality ghost correction. This method works well when the autocalibration data are pristine, but performance degrades substantially when the autocalibration information is imperfect. RAC-LORAKS generalizes Autocalibrated LORAKS in two ways. First, it does not completely trust the information from autocalibration data, and instead considers the autocalibration and EPI data simultaneously when estimating low-rank matrix structure. Second, it uses complementary information from the autocalibration data to improve EPI reconstruction in a multi-contrast joint reconstruction framework. RAC-LORAKS is evaluated using simulations and in vivo data, including comparisons to state-of-the-art methods. RESULTS: RAC-LORAKS is demonstrated to have good ghost elimination performance compared to state-of-the-art methods in several complicated EPI acquisition scenarios (including gradient-echo brain imaging, diffusion-encoded brain imaging, and cardiac imaging). CONCLUSIONS: RAC-LORAKS provides effective suppression of EPI ghosts and is robust to imperfect autocalibration data.

Autoregression and Structured Low-Rank Modeling of Sinogram Neighborhoods.

Sinograms are commonly used to represent the raw data from tomographic imaging experiments. Although it is already well-known that sinograms posess some amount of redundancy, in this work, we present novel theory suggesting that sinograms will often possess substantial additional redundancies that have not been explicitly exploited by previous methods. Specifically, we derive that sinograms will often satisfy multiple simple data-dependent autoregression relationships. This kind of autoregressive structure enables missing/degraded sinogram samples to be linearly predicted using a simple shift-invariant linear combination of neighboring samples. Our theory also further implies that if sinogram samples are assembled into a structured Hankel/Toeplitz matrix, then the matrix will be expected to have low-rank characteristics. As a result, sinogram restoration problems can be formulated as structured low-rank matrix recovery problems. Illustrations of this approach are provided using several different (real and simulated) X-ray imaging datasets, including comparisons against a state-of-the-art deep learning approach. Results suggest that structured low-rank matrix methods for sinogram recovery can have comparable performance to state-of-the-art approaches. Although our evaluation focuses on competitive comparisons against other approaches, we believe that autoregressive constraints are actually complementary to existing approaches with strong potential synergies.

Navigator-Free EPI Ghost Correction With Structured Low-Rank Matrix Models: New Theory and Methods.

Structured low-rank matrix models have previously been introduced to enable calibrationless MR image reconstruction from sub-Nyquist data, and such ideas have recently been extended to enable navigator-free echo-planar imaging (EPI) ghost correction. This paper presents a novel theoretical analysis which shows that, because of uniform subsampling, the structured low-rank matrix optimization problems for EPI data will always have either undesirable or non-unique solutions in the absence of additional constraints. This theory leads us to recommend and investigate problem formulations for navigator-free EPI that incorporate side information from either image-domain or k-space domain parallel imaging methods. The importance of using nonconvex low-rank matrix regularization is also identified. We demonstrate using phantom and in vivo data that the proposed methods are able to eliminate ghost artifacts for several navigator-free EPI acquisition schemes, obtaining better performance in comparison with the state-of-the-art methods across a range of different scenarios. Results are shown for both single-channel acquisition and highly accelerated multi-channel acquisition.

ROBUST AUTOCALIBRATED LORAKS FOR EPI GHOST CORRECTION.

Nyquist ghosts are a longstanding problem in a variety of fast MRI experiments that use echo-planar imaging (EPI). Recently, several structured low-rank matrix modeling approaches have been proposed that achieve state-of-the-art ghost-elimination, although the performance of these approaches is still inadequate in several important scenarios. We present a new structured low-rank matrix recovery ghost correction method that we call Robust Autocalibrated LORAKS (RAC-LORAKS). RAC-LORAKS incorporates constraints from autocalibration data to avoid ill-posedness, but allows adaptation of these constraints to gain robustness against possible autocalibration imperfections. RAC-LORAKS is tested in two challenging scenarios: highly-undersampled multi-channel EPI of the brain, and cardiac EPI with a double-oblique slice orientation. Results show that RAC-LORAKS can provide substantial improvements over existing ghost correction methods, and potentially enables new imaging applications that were previously confounded by ghost artifacts.