The "hidden noise" problem in MR image reconstruction.

PURPOSE: The performance of modern image reconstruction methods is commonly judged using quantitative error metrics like root mean squared-error and the structural similarity index, which are calculated by comparing reconstructed images against fully sampled reference data. In practice, the reference data will contain noise and is not a true gold standard. In this work, we demonstrate that the "hidden noise" present in reference data can substantially confound standard approaches for ranking different image reconstruction results. METHODS: Using both experimental and simulated k-space data and several different image reconstruction techniques, we examined whether there was correlation between performance metrics obtained with typical noisy reference data versus those obtained with higher-quality reference data. RESULTS: For conventional performance metrics, the reconstructions that matched best with the higher-quality reference data were substantially different from the reconstructions that matched best with typical noisy reference data. This leads to suboptimal reconstruction results if the performance with respect to noisy reference data is used to select which reconstruction methods/parameters to employ. These issues were reduced when employing alternative error metrics that better account for noise. CONCLUSION: Reference data containing hidden noise can substantially mislead the ranking of image reconstruction methods when using conventional error metrics, but this issue can be mitigated with alternative error metrics.

Diffusion tensor brain imaging at 0.55T: A feasibility study.

PURPOSE: To investigate the feasibility of diffusion tensor brain imaging at 0.55T with comparisons against 3T. METHODS: Diffusion tensor imaging data with 2 mm isotropic resolution was acquired on a cohort of five healthy subjects using both 0.55T and 3T scanners. The signal-to-noise ratio (SNR) of the 0.55T data was improved using a previous SNR-enhancing joint reconstruction method that jointly reconstructs the entire set of diffusion weighted images from k-space using shared-edge constraints. Quantitative diffusion tensor parameters were estimated and compared across field strengths. We also performed a test-retest assessment of repeatability at each field strength. RESULTS: After applying SNR-enhancing joint reconstruction, the diffusion tensor parameters obtained from 0.55T data were strongly correlated ( R 2 ≥ 0 . 70 $$ {R}^2\ge 0.70 $$ ) with those obtained from 3T data. Test-retest analysis showed that SNR-enhancing reconstruction improved the repeatability of the 0.55T diffusion tensor parameters. CONCLUSION: High-resolution in vivo diffusion MRI of the human brain is feasible at 0.55T when appropriate noise-mitigation strategies are applied.

Spherical echo-planar time-resolved imaging (sEPTI) for rapid 3D quantitative T 2 * and susceptibility imaging.

PURPOSE: To develop a 3D spherical EPTI (sEPTI) acquisition and a comprehensive reconstruction pipeline for rapid high-quality whole-brain submillimeter T 2 * $$ {\mathrm{T}}_2^{\ast } $$ and QSM quantification. METHODS: For the sEPTI acquisition, spherical k-space coverage is utilized with variable echo-spacing and maximum kx ramp-sampling to improve efficiency and signal incoherency compared to existing EPTI approaches. For reconstruction, an iterative rank-shrinking B0 estimation and odd-even high-order phase correction algorithms were incorporated into the reconstruction to better mitigate artifacts from field imperfections. A physics-informed unrolled network was utilized to boost the SNR, where 1-mm and 0.75-mm isotropic whole-brain imaging were performed in 45 and 90 s at 3 T, respectively. These protocols were validated through simulations, phantom, and in vivo experiments. Ten healthy subjects were recruited to provide sufficient data for the unrolled network. The entire pipeline was validated on additional five healthy subjects where different EPTI sampling approaches were compared. Two additional pediatric patients with epilepsy were recruited to demonstrate the generalizability of the unrolled reconstruction. RESULTS: sEPTI achieved 1.4 × $$ \times $$ faster imaging with improved image quality and quantitative map precision compared to existing EPTI approaches. The B0 update and the phase correction provide improved reconstruction performance with lower artifacts. The unrolled network boosted the SNR, achieving high-quality T 2 * $$ {\mathrm{T}}_2^{\ast } $$ and QSM quantification with single average data. High-quality reconstruction was also obtained in the pediatric patients using this network. CONCLUSION: sEPTI achieved whole-brain distortion-free multi-echo imaging and T 2 * $$ {\mathrm{T}}_2^{\ast } $$ and QSM quantification at 0.75 mm in 90 s which has the potential to be useful for wide clinical applications.

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 Ambiguity in Linear Inverse Problems: Entrywise Bounds on Nearly Data-Consistent Solutions and Entrywise Condition Numbers.

Ill-posed linear inverse problems appear frequently in various signal processing applications. It can be very useful to have theoretical characterizations that quantify the level of ill-posedness for a given inverse problem and the degree of ambiguity that may exist about its solution. Traditional measures of ill-posedness, such as the condition number of a matrix, provide characterizations that are global in nature. While such characterizations can be powerful, they can also fail to provide full insight into situations where certain entries of the solution vector are more or less ambiguous than others. In this work, we derive novel theoretical lower- and upper-bounds that apply to individual entries of the solution vector, and are valid for all potential solution vectors that are nearly data-consistent. These bounds are agnostic to the noise statistics and the specific method used to solve the inverse problem, and are also shown to be tight. In addition, our results also lead us to introduce an entrywise version of the traditional condition number, which provides a substantially more nuanced characterization of scenarios where certain elements of the solution vector are less sensitive to perturbations than others. Our results are illustrated in an application to magnetic resonance imaging reconstruction, and we include discussions of practical computation methods for large-scale inverse problems, connections between our new theory and the traditional Cramér-Rao bound under statistical modeling assumptions, and potential extensions to cases involving constraints beyond just data-consistency.

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.

Local perturbation responses and checkerboard tests: Characterization tools for nonlinear MRI methods.

PURPOSE: Modern methods for MR image reconstruction, denoising, and parameter mapping are becoming increasingly nonlinear, black-box, and at risk of "hallucination." These trends mean that traditional tools for judging confidence in an image (visual quality assessment, point-spread functions (PSFs), g-factor maps, etc.) are less helpful than before. This paper describes and evaluates an approach that can help with assessing confidence in images produced by arbitrary nonlinear methods. THEORY AND METHODS: We propose to characterize nonlinear methods by examining the images they produce before and after applying controlled perturbations to the measured data. This results in functions known as local perturbation responses (LPRs) that can provide useful insight into sensitivity, spatial resolution, and aliasing characteristics. LPRs can be viewed as generalizations of classical PSFs, and are are very flexible-they can be applied to arbitary nonlinear methods and arbitrary datasets across a range of different reconstruction, denoising, and parameter mapping applications. Importantly, LPRs do not require a ground truth image. RESULTS: Impulse-based and checkerboard-pattern LPRs are demonstrated in image reconstruction and denoising scenarios. We observe that these LPRs provide insights into spatial resolution, signal leakage, and aliasing that are not available with other methods. We also observe that popular reference-based image quality metrics (eg, mean-squared error and structural similarity) do not always correlate with good LPR characteristics. CONCLUSIONS: LPRs are a useful tool that can be used to characterize and assess confidence in nonlinear MR methods, and provide insights that are distinct from and complementary to existing quality assessments.

Region-optimized virtual (ROVir) coils: Localization and/or suppression of spatial regions using sensor-domain beamforming.

PURPOSE: In many MRI scenarios, magnetization is often excited from spatial regions that are not of immediate interest. Excitation of uninteresting magnetization can complicate the design of efficient imaging methods, leading to either artifacts or acquisitions that are longer than necessary. While there are many hardware- and sequence-based approaches for suppressing unwanted magnetization, this paper approaches this longstanding problem from a different and complementary angle, using beamforming to suppress signals from unwanted regions without modifying the acquisition hardware or pulse sequence. Unlike existing beamforming approaches, we use a spatially invariant sensor-domain approach that can be applied directly to raw data to facilitate image reconstruction. THEORY AND METHODS: We use beamforming to linearly mix a set of original coils into a set of region-optimized virtual (ROVir) coils. ROVir coils optimize a signal-to-interference ratio metric, are easily calculated using simple generalized eigenvalue decomposition methods, and provide coil compression. RESULTS: ROVir coils were compared against existing coil-compression methods, and were demonstrated to have substantially better signal suppression capabilities. In addition, examples were provided in a variety of different application contexts (including brain MRI, vocal tract MRI, and cardiac MRI; accelerated Cartesian and non-Cartesian imaging; and outer volume suppression) that demonstrate the strong potential of this kind of approach. CONCLUSION: The beamforming-based ROVir framework is simple to implement, has promising capabilities to suppress unwanted MRI signal, and is compatible with and complementary to existing signal suppression methods. We believe that this general approach could prove useful across a wide range of different 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.

Combined diffusion-relaxometry microstructure imaging: Current status and future prospects.

Microstructure imaging seeks to noninvasively measure and map microscopic tissue features by pairing mathematical modeling with tailored MRI protocols. This article reviews an emerging paradigm that has the potential to provide a more detailed assessment of tissue microstructure-combined diffusion-relaxometry imaging. Combined diffusion-relaxometry acquisitions vary multiple MR contrast encodings-such as b-value, gradient direction, inversion time, and echo time-in a multidimensional acquisition space. When paired with suitable analysis techniques, this enables quantification of correlations and coupling between multiple MR parameters-such as diffusivity, T1 , T2 , and T2∗ . This opens the possibility of disentangling multiple tissue compartments (within voxels) that are indistinguishable with single-contrast scans, enabling a new generation of microstructural maps with improved biological sensitivity and specificity.

Corticospinal Tract Microstructure Predicts Distal Arm Motor Improvements in Chronic Stroke.

BACKGROUND AND PURPOSE: The corticospinal tract (CST) is a crucial brain pathway for distal arm and hand motor control. We aimed to determine whether a diffusion tensor imaging (DTI)-derived CST metric predicts distal upper extremity (UE) motor improvements in chronic stroke survivors. METHODS: We analyzed clinical and neuroimaging data from a randomized controlled rehabilitation trial. Participants completed clinical assessments and neuroimaging at baseline and clinical assessments 4 months later, postintervention. Using univariate linear regression analysis, we determined the linear relationship between the DTI-derived CST fractional anisotropy asymmetry (FAasym) and the percentage of baseline change in log-transformed average Wolf Motor Function Test time for distal items (ΔlnWMFT-distal_%). The least absolute shrinkage and selection operator (LASSO) linear regressions with cross-validation and bootstrapping were used to determine the relative weighting of CST FAasym, other brain metrics, clinical outcomes, and demographics on distal motor improvement. Logistic regression analyses were performed to test whether the CST FAasym can predict clinically significant UE motor improvement. RESULTS: lnWMFT-distal significantly improved at the group level. Baseline CST FAasym explained 26% of the variance in ΔlnWMFT-distal_%. A multivariate LASSO model including baseline CST FAasym, age, and UE Fugl-Meyer explained 39% of the variance in ΔlnWMFT-distal_%. Further, CST FAasym explained more variance in ΔlnWMFT-distal_% than the other significant predictors in the LASSO model. DISCUSSION AND CONCLUSIONS: CST microstructure is a significant predictor of improvement in distal UE motor function in the context of an UE rehabilitation trial in chronic stroke survivors with mild-to-moderate motor impairment.Video Abstract available for more insight from the authors (see the Video, Supplemental Digital Content 1, available at: http://links.lww.com/JNPT/A350).

Impact of (k,t) sampling on DCE MRI tracer kinetic parameter estimation in digital reference objects.

PURPOSE: To evaluate the impact of (k,t) data sampling on the variance of tracer-kinetic parameter (TK) estimation in high-resolution whole-brain dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) using digital reference objects. We study this in the context of TK model constraints, and in the absence of other constraints. METHODS: Three anatomically and physiologically realistic brain-tumor digital reference objects were generated. Data sampling strategies included uniform and variable density; zone-based, lattice, pseudo-random, and pseudo-radial; with 50-time frames and 4-fold to 25-fold undersampling. In all cases, we assume a fully sampled first time frame, and prior knowledge of the arterial input function. TK parameters were estimated by indirect estimation (i.e., image-time-series reconstruction followed by model fitting), and direct estimation from the under-sampled data. We evaluated methods based on the Cramér-Rao bound and Monte-Carlo simulations, over the range of signal-to-noise ratio (SNR) seen in clinical brain DCE-MRI. RESULTS: Lattice-based sampling provided the lowest SDs, followed by pseudo-random, pseudo-radial, and zone-based. This ranking was consistent for the Patlak and extended Tofts model. Pseudo-random sampling resulted in 19% higher averaged SD compared to lattice-based sampling. Zone-based sampling resulted in substantially higher SD at undersampling factors above 10. CRB analysis showed only a small difference between uniform and variable density for both lattice-based and pseudo-random sampling up to undersampling factors of 25. CONCLUSION: Lattice sampling provided the lowest SDs, although the differences between sampling schemes were not substantial at low undersampling factors. The differences between lattice-based and pseudo-random sampling strategies with both uniform and variable density were within the range of error induced by other sources, at up to 25-fold undersampling.

Fast submillimeter diffusion MRI using gSlider-SMS and SNR-enhancing joint reconstruction.

PURPOSE: We evaluate a new approach for achieving diffusion MRI data with high spatial resolution, large volume coverage, and fast acquisition speed. THEORY AND METHODS: A recent method called gSlider-SMS enables whole-brain submillimeter diffusion MRI with high signal-to-noise ratio (SNR) efficiency. However, despite the efficient acquisition, the resulting images can still suffer from low SNR due to the small size of the imaging voxels. This work proposes to mitigate the SNR problem by combining gSlider-SMS with a regularized SNR-enhancing reconstruction approach. RESULTS: Illustrative results show that, from gSlider-SMS data acquired over a span of only 25 minutes on a 3T scanner, the proposed method is able to produce 71 MRI images (64 diffusion encoding orientations with b = 1500 s/ mm2 , and 7 images without diffusion weighting) of the entire in vivo human brain with nominal 0.66 mm spatial resolution. Using data acquired from 75 minutes of acquisition as a gold standard reference, we demonstrate that the proposed SNR-enhancement procedure leads to substantial improvements in estimated diffusion parameters compared to conventional gSlider reconstruction. Results also demonstrate that the proposed method has advantages relative to denoising methods based on low-rank matrix modeling. A theoretical analysis of the trade-off between spatial resolution and SNR suggests that the proposed approach has high efficiency. CONCLUSIONS: The combination of gSlider-SMS with advanced regularized reconstruction enables high-resolution quantitative diffusion MRI from a relatively fast acquisition.

Multidimensional correlation spectroscopic imaging of exponential decays: From theoretical principles to in vivo human applications.

Multiexponential modeling of relaxation or diffusion MR signal decays is a popular approach for estimating and spatially mapping different microstructural tissue compartments. While this approach can be quite powerful, it is also limited by the fact that one-dimensional multiexponential modeling is an ill-posed inverse problem with substantial ambiguities. In this article, we present an overview of a recent multidimensional correlation spectroscopic imaging approach to this problem. This approach helps to alleviate ill-posedness by making advantageous use of multidimensional contrast encoding (e.g., 2D diffusion-relaxation encoding or 2D relaxation-relaxation encoding) combined with a regularized spatial-spectral estimation procedure. Theoretical calculations, simulations, and experimental results are used to illustrate the benefits of this approach relative to classical methods. In addition, we demonstrate an initial proof-of-principle application of this kind of approach to in vivo human MRI experiments.

Wave-LORAKS: Combining wave encoding with structured low-rank matrix modeling for more highly accelerated 3D imaging.

PURPOSE: Wave-CAIPI is a novel acquisition approach that enables highly accelerated 3D imaging. This paper investigates the combination of Wave-CAIPI with LORAKS-based reconstruction (Wave-LORAKS) to enable even further acceleration. METHODS: LORAKS is a constrained image reconstruction framework that can impose spatial support, smooth phase, sparsity, and/or parallel imaging constraints. LORAKS requires minimal prior information, and instead uses the low-rank subspace structure of the raw data to automatically learn which constraints to impose and how to impose them. Previous LORAKS implementations addressed 2D image reconstruction problems. In this work, several recent advances in structured low-rank matrix recovery were combined to enable large-scale 3D Wave-LORAKS reconstruction with improved quality and computational efficiency. Wave-LORAKS was investigated by retrospective subsampling of two fully sampled Wave-encoded 3D MPRAGE datasets, and comparisons were made against existing Wave reconstruction approaches. RESULTS: Results show that Wave-LORAKS can yield higher reconstruction quality with 16×-accelerated data than is obtained by traditional Wave-CAIPI with 9×-accerated data. CONCLUSIONS: There are strong synergies between Wave encoding and LORAKS, which enables Wave-LORAKS to achieve higher acceleration and more flexible sampling compared to Wave-CAIPI.

Childhood Music Training Induces Change in Micro and Macroscopic Brain Structure: Results from a Longitudinal Study.

Several studies comparing adult musicians and nonmusicians have shown that music training is associated with structural brain differences. It is not been established, however, whether such differences result from pre-existing biological traits, lengthy musical training, or an interaction of the two factors, or if comparable changes can be found in children undergoing music training. As part of an ongoing longitudinal study, we investigated the effects of music training on the developmental trajectory of children's brain structure, over two years, beginning at age 6. We compared these children with children of the same socio-economic background but either involved in sports training or not involved in any systematic after school training. We established at the onset that there were no pre-existing structural differences among the groups. Two years later we observed that children in the music group showed (1) a different rate of cortical thickness maturation between the right and left posterior superior temporal gyrus, and (2) higher fractional anisotropy in the corpus callosum, specifically in the crossing pathways connecting superior frontal, sensory, and motor segments. We conclude that music training induces macro and microstructural brain changes in school-age children, and that those changes are not attributable to pre-existing biological traits.

Improving parallel imaging by jointly reconstructing multi-contrast data.

PURPOSE: To develop parallel imaging techniques that simultaneously exploit coil sensitivity encoding, image phase prior information, similarities across multiple images, and complementary k-space sampling for highly accelerated data acquisition. METHODS: We introduce joint virtual coil (JVC)-generalized autocalibrating partially parallel acquisitions (GRAPPA) to jointly reconstruct data acquired with different contrast preparations, and show its application in 2D, 3D, and simultaneous multi-slice (SMS) acquisitions. We extend the joint parallel imaging concept to exploit limited support and smooth phase constraints through Joint (J-) LORAKS formulation. J-LORAKS allows joint parallel imaging from limited autocalibration signal region, as well as permitting partial Fourier sampling and calibrationless reconstruction. RESULTS: We demonstrate highly accelerated 2D balanced steady-state free precession with phase cycling, SMS multi-echo spin echo, 3D multi-echo magnetization-prepared rapid gradient echo, and multi-echo gradient recalled echo acquisitions in vivo. Compared to conventional GRAPPA, proposed joint acquisition/reconstruction techniques provide more than 2-fold reduction in reconstruction error. CONCLUSION: JVC-GRAPPA takes advantage of additional spatial encoding from phase information and image similarity, and employs different sampling patterns across acquisitions. J-LORAKS achieves a more parsimonious low-rank representation of local k-space by considering multiple images as additional coils. Both approaches provide dramatic improvement in artifact and noise mitigation over conventional single-contrast parallel imaging reconstruction. Magn Reson Med 80:619-632, 2018. © 2018 International Society for Magnetic Resonance in Medicine.

A theoretical signal processing framework for linear diffusion MRI: Implications for parameter estimation and experiment design.

The data measured in diffusion MRI can be modeled as the Fourier transform of the Ensemble Average Propagator (EAP), a probability distribution that summarizes the molecular diffusion behavior of the spins within each voxel. This Fourier relationship is potentially advantageous because of the extensive theory that has been developed to characterize the sampling requirements, accuracy, and stability of linear Fourier reconstruction methods. However, existing diffusion MRI data sampling and signal estimation methods have largely been developed and tuned without the benefit of such theory, instead relying on approximations, intuition, and extensive empirical evaluation. This paper aims to address this discrepancy by introducing a novel theoretical signal processing framework for diffusion MRI. The new framework can be used to characterize arbitrary linear diffusion estimation methods with arbitrary q-space sampling, and can be used to theoretically evaluate and compare the accuracy, resolution, and noise-resilience of different data acquisition and parameter estimation techniques. The framework is based on the EAP, and makes very limited modeling assumptions. As a result, the approach can even provide new insight into the behavior of model-based linear diffusion estimation methods in contexts where the modeling assumptions are inaccurate. The practical usefulness of the proposed framework is illustrated using both simulated and real diffusion MRI data in applications such as choosing between different parameter estimation methods and choosing between different q-space sampling schemes.

LORAKS makes better SENSE: Phase-constrained partial fourier SENSE reconstruction without phase calibration.

PURPOSE: Parallel imaging and partial Fourier acquisition are two classical approaches for accelerated MRI. Methods that combine these approaches often rely on prior knowledge of the image phase, but the need to obtain this prior information can place practical restrictions on the data acquisition strategy. In this work, we propose and evaluate SENSE-LORAKS, which enables combined parallel imaging and partial Fourier reconstruction without requiring prior phase information. THEORY AND METHODS: The proposed formulation is based on combining the classical SENSE model for parallel imaging data with the more recent LORAKS framework for MR image reconstruction using low-rank matrix modeling. Previous LORAKS-based methods have successfully enabled calibrationless partial Fourier parallel MRI reconstruction, but have been most successful with nonuniform sampling strategies that may be hard to implement for certain applications. By combining LORAKS with SENSE, we enable highly accelerated partial Fourier MRI reconstruction for a broader range of sampling trajectories, including widely used calibrationless uniformly undersampled trajectories. RESULTS: Our empirical results with retrospectively undersampled datasets indicate that when SENSE-LORAKS reconstruction is combined with an appropriate k-space sampling trajectory, it can provide substantially better image quality at high-acceleration rates relative to existing state-of-the-art reconstruction approaches. CONCLUSION: The SENSE-LORAKS framework provides promising new opportunities for highly accelerated MRI. Magn Reson Med 77:1021-1035, 2017. © 2016 International Society for Magnetic Resonance in Medicine.