High-fidelity intravoxel incoherent motion parameter mapping using locally low-rank and subspace modeling.

PURPOSE: Intravoxel incoherent motion (IVIM) is a quantitative magnetic resonance imaging (MRI) method used to quantify perfusion properties of tissue non-invasively without contrast. However, clinical applications are limited by unreliable parameter estimates, particularly for the perfusion fraction (f) and pseudodiffusion coefficient (D*). This study aims to develop a high-fidelity reconstruction for reliable estimation of IVIM parameters. The proposed method is versatile and amenable to various acquisition schemes and fitting methods. METHODS: To address current challenges with IVIM, we adapted several advanced reconstruction techniques. We used a low-rank approximation of IVIM images and temporal subspace modeling to constrain the magnetization dynamics of the bi-exponential diffusion signal decay. In addition, motion-induced phase variations were corrected between diffusion directions and b-values, facilitating the use of high SNR real-valued diffusion data. The proposed method was evaluated in simulations and in vivo brain acquisitions in six healthy subjects and six individuals with a history of SARS-CoV-2 infection and compared with the conventionally reconstructed magnitude data. Following reconstruction, IVIM parameters were estimated voxel-wise. RESULTS: Our proposed method reduced noise contamination in simulations, resulting in a 60%, 58.9%, and 83.9% reduction in the NRMSE for D, f, and D*, respectively, compared to the conventional reconstruction. In vivo, anisotropic properties of D, f, and D* were preserved with the proposed method, highlighting microvascular differences in gray matter between individuals with a history of COVID-19 and those without (p = 0.0210), which wasn't observed with the conventional reconstruction. CONCLUSION: The proposed method yielded a more reliable estimation of IVIM parameters with less noise than the conventional reconstruction. Further, the proposed method preserved anisotropic properties of IVIM parameter estimates and demonstrated differences in microvascular perfusion in COVID-affected subjects, which weren't observed with conventional reconstruction methods.

Multi-band multi-shot diffusion MRI reconstruction with joint usage of structured low-rank constraints and explicit phase mapping.

PURPOSE: To develop a joint reconstruction method for multi-band multi-shot diffusion MRI. THEORY AND METHODS: Multi-band multi-shot EPI acquisition is an effective approach for high-resolution diffusion MRI, but requires specific algorithms to correct the inter-shot phase variations. The phase correction can be done by first estimating the explicit phase map and then feeding it into the k-space signal formulation model. Alternatively, the phase information can be used indirectly as structured low-rank constraints in k-space. The 2 methods differ in reconstruction accuracy and efficiency. We aim to combine the 2 different approaches for improved image quality and reconstruction efficiency simultaneously, termed "joint usage of structured low-rank constraints and explicit phase mapping" (JULEP). The proposed JULEP reconstruction is tested on both single-band and multi-band, multi-shot diffusion data, with different resolutions and b values. The results of JULEP are compared with conventional methods with explicit phase mapping (i.e., multiplexed sensitivity-encoding [MUSE]) and structured low-rank constraints (i.e., MUSSELS), and another joint reconstruction method (i.e., network estimated artifacts for tempered reconstruction [NEATR]). RESULTS: JULEP improves the quality of the navigator and subsequently facilitates the reconstruction of final diffusion images. Compared with all 3 other methods (MUSE, MUSSELS, and NEATR), JULEP mitigates residual structural bias and improves temporal SNRs in the final diffusion image, particularly at high multi-band factors. Compared with MUSSELS, JULEP also improves computational efficiency. CONCLUSION: The proposed JULEP method significantly improves the image quality and reconstruction efficiency of multi-band multi-shot diffusion MRI, which can promote a broader application of high-resolution diffusion MRI.

Multi-band- and in-plane-accelerated diffusion MRI enabled by model-based deep learning in q-space and its extension to learning in the spherical harmonic domain.

PURPOSE: To propose a new method for the recovery of combined in-plane- and multi-band (MB)-accelerated diffusion MRI data. METHODS: Combining MB acceleration with in-plane acceleration is crucial to improve the time efficiency of high (angular and spatial) resolution diffusion scans. However, as the MB factor and in-plane acceleration increase, the reconstruction becomes challenging due to the heavy aliasing. The new reconstruction utilizes an additional q-space prior to constrain the recovery, which is derived from the previously proposed qModeL framework. Specifically, the qModeL prior provides a pre-learned representation of the diffusion signal space to which the measured data belongs. We show that the pre-learned q-space prior along with a model-based iterative reconstruction that accommodate multi-band unaliasing, can efficiently reconstruct the in-plane- and MB-accelerated data. The power of joint reconstruction is maximally utilized by using an incoherent under-sampling pattern in the k-q domain. We tested the proposed method on single- and multi-shell data, with high/low angular resolution, high/low spatial resolution, healthy/abnormal tissues, and 3T/7T field strengths. Furthermore, the learning is extended to the spherical harmonic basis, to provide a rotational invariant learning framework. RESULTS: The qModeL joint reconstruction is shown to simultaneously unalias and jointly recover DWIs with reasonable accuracy in all the cases studied. The reconstruction error from 18-fold accelerated multi-shell datasets was <3%. The microstructural maps derived from the accelerated acquisitions also exhibit reasonable accuracy for both healthy and abnormal tissues. The deep learning (DL)-enabled reconstructions are comparable to those derived using traditional methods. CONCLUSION: qModeL enables the joint recovery of combined in-plane- and MB-accelerated dMRI utilizing DL.

qModeL: A plug-and-play model-based reconstruction for highly accelerated multi-shot diffusion MRI using learned priors.

PURPOSE: To introduce a joint reconstruction method for highly undersampled multi-shot diffusion weighted (msDW) scans. METHODS: Multi-shot EPI methods enable higher spatial resolution for diffusion MRI, but at the expense of long scan-time. Highly accelerated msDW scans are needed to enable their utilization in advanced microstructure studies, which require high q-space coverage. Previously, joint k-q undersampling methods coupled with compressed sensing were shown to enable very high acceleration factors. However, the reconstruction of this data using sparsity priors is challenging and is not suited for multi-shell data. We propose a new reconstruction that recovers images from the combined k-q data jointly. The proposed qModeL reconstruction brings together the advantages of model-based iterative reconstruction and machine learning, extending the idea of plug-and-play algorithms. Specifically, qModeL works by prelearning the signal manifold corresponding to the diffusion measurement space using deep learning. The prelearned manifold prior is incorporated into a model-based reconstruction to provide a voxel-wise regularization along the q-dimension during the joint recovery. Notably, the learning does not require in vivo training data and is derived exclusively from biophysical modeling. Additionally, a plug-and-play total variation denoising provides regularization along the spatial dimension. The proposed framework is tested on k-q undersampled single-shell and multi-shell msDW acquisition at various acceleration factors. RESULTS: The qModeL joint reconstruction is shown to recover DWIs from 8-fold accelerated msDW acquisitions with error less than 5% for both single-shell and multi-shell data. Advanced microstructural analysis performed using the undersampled reconstruction also report reasonable accuracy. CONCLUSION: qModeL enables the joint recovery of highly accelerated multi-shot dMRI utilizing learning-based priors. The bio-physically driven approach enables the use of accelerated multi-shot imaging for multi-shell sampling and advanced microstructure studies.

Improved MUSSELS reconstruction for high-resolution multi-shot diffusion weighted imaging.

PURPOSE: MUSSELS is a one-step iterative reconstruction method for multishot diffusion weighted (msDW) imaging. The current work presents an efficient implementation, termed IRLS MUSSELS, that enables faster reconstruction to enhance its utility for high-resolution diffusion MRI studies. METHODS: The recently proposed MUSSELS reconstruction belongs to a new class of parallel imaging-based methods that recover artifact-free DWIs from msDW data without needing phase compensation. The reconstruction is achieved via structured low-rank matrix completion algorithms, which are computationally demanding due to the large size of the Hankel matrices and their associated computations involving singular value decompositions. Because of this, computational demands of the MUSSELS reconstruction scales as the matrix size and the number of shots increases, which hinders its practical utility for high-resolution applications. In this work, we derive a computationally efficient MUSSELS formulation by modifying the iterative reweighted least squares (IRLS) method that were proposed earlier to solve such problems. Using whole-brain in vivo data, we show the utility of the IRLS MUSSELS for routine high-resolution studies with reduced computational burden. RESULTS: IRLS MUSSELS provides about five times faster reconstruction for matrix sizes 192 × 192 and 256 × 256 compared to the earlier MUSSELS implementation. The widely employed conjugate symmetry priors can also be incorporated into IRLS MUSSELS to reduce blurring of the partial Fourier acquisitions, without incurring much computational burden. CONCLUSIONS: The proposed method is observed to be computationally efficient to enable routine high-resolution studies. The computational complexity matches the traditional msDWI reconstruction methods and provides improved reconstruction results with the additional constraints.

SMS MUSSELS: A navigator-free reconstruction for simultaneous multi-slice-accelerated multi-shot diffusion weighted imaging.

PURPOSE: To introduce a novel reconstruction method for simultaneous multi-slice (SMS)-accelerated multi-shot diffusion weighted imaging (ms-DWI). METHODS: SMS acceleration using blipped-CAIPI schemes have been proposed to speed up the acquisition of ms-DWIs. The reconstruction of the data requires (a) phase compensation to combine data from different shots and (b) slice unfolding to separate the data of different slices. The traditional approaches first estimate the phase maps corresponding to each shot and slice which are then employed to iteratively recover the slice unfolded DWIs without phase artifacts. In contrast, the proposed reconstruction directly recovers the slice-unfolded k-space data of the multiple shots for each slice in a single-step recovery scheme. The proposed method is enabled by the low-rank property inherent in the k-space samples of ms-DW acquisition. This enabled to formulate a joint recovery scheme that simultaneously (a) unfolds the k-space data of each slice using a SENSE-based scheme and (b) recover the missing k-space samples in each slice of the multi-shot acquisition employing a structured low-rank matrix completion. Additional smoothness regularization is also utilized for higher acceleration factors. The proposed joint recovery is tested on simulated and in vivo data and compared to similar un-navigated methods. RESULTS: Our experiments show effective slice unfolding and successful recovery of DWIs with minimal phase artifacts using the proposed method. The performance is comparable to existing methods at low acceleration factors and better than existing methods for higher acceleration factors. CONCLUSIONS: For the slice accelerations considered in this study, the proposed method can successfully recover DWIs from SMS-accelerated ms-DWI acquisitions.

Structured Low-Rank Algorithms: Theory, Magnetic Resonance Applications, and Links to Machine Learning.

In this survey, we provide a detailed review of recent advances in the recovery of continuous domain multidimensional signals from their few non-uniform (multichannel) measurements using structured low-rank matrix completion formulation. This framework is centered on the fundamental duality between the compactness (e.g., sparsity) of the continuous signal and the rank of a structured matrix, whose entries are functions of the signal. This property enables the reformulation of the signal recovery as a low-rank structured matrix completion, which comes with performance guarantees. We will also review fast algorithms that are comparable in complexity to current compressed sensing methods, which enables the application of the framework to large-scale magnetic resonance (MR) recovery problems. The remarkable flexibility of the formulation can be used to exploit signal properties that are difficult to capture by current sparse and low-rank optimization strategies. We demonstrate the utility of the framework in a wide range of MR imaging (MRI) applications, including highly accelerated imaging, calibration-free acquisition, MR artifact correction, and ungated dynamic MRI.

A general algorithm for compensation of trajectory errors: Application to radial imaging.

PURPOSE: To reconstruct artifact-free images from measured k-space data, when the actual k-space trajectory deviates from the nominal trajectory due to gradient imperfections. METHODS: Trajectory errors arising from eddy currents and gradient delays introduce phase inconsistencies in several fast scanning MR pulse sequences, resulting in image artifacts. The proposed algorithm provides a novel framework to compensate for this phase distortion. The algorithm relies on the construction of a multi-block Hankel matrix, where each block is constructed from k-space segments with the same phase distortion. In the presence of spatially smooth phase distortions between the segments, the complete block-Hankel matrix is known to be highly low-rank. Since each k-space segment is only acquiring part of the k-space data, the reconstruction of the phase compensated image from their partially parallel measurements is posed as a structured low-rank matrix optimization problem, assuming the coil sensitivities to be known. RESULTS: The proposed formulation is tested on radial acquisitions in several settings including partial Fourier and golden-angle acquisitions. The experiments demonstrate the ability of the algorithm to successfully remove the artifacts arising from the trajectory errors, without the need for trajectory or phase calibration. The quality of the reconstruction was comparable to corrections achieved using the Trajectory Auto-Corrected Image Reconstruction (TrACR) for radial acquisitions. CONCLUSION: The proposed method provides a general framework for the recovery of artifact-free images from radial trajectories without the need for trajectory calibration.

Alterations of the cerebellum and basal ganglia in bipolar disorder mood states detected by quantitative T1ρ mapping.

OBJECTIVES: Quantitative mapping of T1 relaxation in the rotating frame (T1ρ) is a magnetic resonance imaging technique sensitive to pH and other cellular and microstructural factors, and is a potentially valuable tool for identifying brain alterations in bipolar disorder. Recently, this technique identified differences in the cerebellum and cerebral white matter of euthymic patients vs healthy controls that were consistent with reduced pH in these regions, suggesting an underlying metabolic abnormality. The current study built upon this prior work to investigate brain T1ρ differences across euthymic, depressed, and manic mood states of bipolar disorder. METHODS: Forty participants with bipolar I disorder and 29 healthy control participants matched for age and gender were enrolled. Participants with bipolar disorder were imaged in one or more mood states, yielding 27, 12, and 13 imaging sessions in euthymic, depressed, and manic mood states, respectively. Three-dimensional, whole-brain anatomical images and T1ρ maps were acquired for all participants, enabling voxel-wise evaluation of T1ρ differences between bipolar mood state and healthy control groups. RESULTS: All three mood state groups had increased T1ρ relaxation times in the cerebellum compared to the healthy control group. Additionally, the depressed and manic groups had reduced T1ρ relaxation times in and around the basal ganglia compared to the control and euthymic groups. CONCLUSIONS: The study implicated the cerebellum and basal ganglia in the pathophysiology of bipolar disorder and its mood states, the roles of which are relatively unexplored. These findings motivate further investigation of the underlying cause of the abnormalities, and the potential role of altered metabolic activity in these regions.

Multi-shot sensitivity-encoded diffusion data recovery using structured low-rank matrix completion (MUSSELS).

PURPOSE: To introduce a novel method for the recovery of multi-shot diffusion weighted (MS-DW) images from echo-planar imaging (EPI) acquisitions. METHODS: Current EPI-based MS-DW reconstruction methods rely on the explicit estimation of the motion-induced phase maps to recover artifact-free images. In the new formulation, the k-space data of the artifact-free DWI is recovered using a structured low-rank matrix completion scheme, which does not require explicit estimation of the phase maps. The structured matrix is obtained as the lifting of the multi-shot data. The smooth phase-modulations between shots manifest as null-space vectors of this matrix, which implies that the structured matrix is low-rank. The missing entries of the structured matrix are filled in using a nuclear-norm minimization algorithm subject to the data-consistency. The formulation enables the natural introduction of smoothness regularization, thus enabling implicit motion-compensated recovery of the MS-DW data. RESULTS: Our experiments on in-vivo data show effective removal of artifacts arising from inter-shot motion using the proposed method. The method is shown to achieve better reconstruction than the conventional phase-based methods. CONCLUSION: We demonstrate the utility of the proposed method to effectively recover artifact-free images from Cartesian fully/under-sampled and partial Fourier acquired data without the use of explicit phase estimates. Magn Reson Med 78:494-507, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Efficient parallel reconstruction for high resolution multishot spiral diffusion data with low rank constraint.

PURPOSE: To propose a novel reconstruction method using parallel imaging with low rank constraint to accelerate high resolution multishot spiral diffusion imaging. THEORY AND METHODS: The undersampled high resolution diffusion data were reconstructed based on a low rank (LR) constraint using similarities between the data of different interleaves from a multishot spiral acquisition. The self-navigated phase compensation using the low resolution phase data in the center of k-space was applied to correct shot-to-shot phase variations induced by motion artifacts. The low rank reconstruction was combined with sensitivity encoding (SENSE) for further acceleration. The efficiency of the proposed joint reconstruction framework, dubbed LR-SENSE, was evaluated through error quantifications and compared with ℓ1 regularized compressed sensing method and conventional iterative SENSE method using the same datasets. RESULTS: It was shown that with a same acceleration factor, the proposed LR-SENSE method had the smallest normalized sum-of-squares errors among all the compared methods in all diffusion weighted images and DTI-derived index maps, when evaluated with different acceleration factors (R = 2, 3, 4) and for all the acquired diffusion directions. CONCLUSION: Robust high resolution diffusion weighted image can be efficiently reconstructed from highly undersampled multishot spiral data with the proposed LR-SENSE method. Magn Reson Med 77:1359-1366, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Fast iterative algorithm for the reconstruction of multishot non-cartesian diffusion data.

PURPOSE: To accelerate the motion-compensated iterative reconstruction of multishot non-Cartesian diffusion data. METHOD: The motion-compensated recovery of multishot non-Cartesian diffusion data is often performed using a modified iterative sensitivity-encoded algorithm. Specifically, the encoding matrix is replaced with a combination of nonuniform Fourier transforms and composite sensitivity functions, which account for the motion-induced phase errors. The main challenge with this scheme is the significantly increased computational complexity, which is directly related to the total number of composite sensitivity functions (number of shots × number of coils). The dimensionality of the composite sensitivity functions and hence the number of Fourier transforms within each iteration is reduced using a principal component analysis-based scheme. Using a Toeplitz-based conjugate gradient approach in combination with an augmented Lagrangian optimization scheme, a fast algorithm is proposed for the sparse recovery of diffusion data. RESULTS: The proposed simplifications considerably reduce the computation time, especially in the recovery of diffusion data from under-sampled reconstructions using sparse optimization. By choosing appropriate number of basis functions to approximate the composite sensitivities, faster reconstruction (close to 9 times) with effective motion compensation is achieved. CONCLUSION: The proposed enhancements can offer fast motion-compensated reconstruction of multishot diffusion data for arbitrary k-space trajectories.

Acceleration of high angular and spatial resolution diffusion imaging using compressed sensing with multichannel spiral data.

PURPOSE: To accelerate the acquisition of simultaneously high spatial and angular resolution diffusion imaging. METHODS: Accelerated imaging is achieved by recovering the diffusion signal at all voxels simultaneously from under-sampled k-q space data using a compressed sensing algorithm. The diffusion signal at each voxel is modeled as a sparse complex Gaussian mixture model. The joint recovery scheme enables incoherent under-sampling of the 5-D k-q space, obtained by randomly skipping interleaves of a multishot variable density spiral trajectory. This sampling and reconstruction strategy is observed to provide considerably improved reconstructions than classical k-q under-sampling and reconstruction schemes. The complex model enables to account for the noise statistics without compromising the computational efficiency and theoretical convergence guarantees. The reconstruction framework also incorporates compensation of motion induced phase errors that result from the multishot acquisition. RESULTS: Reconstructions of the diffusion signal from under-sampled data using the proposed method yields accurate results with errors less that 5% for different accelerations and b-values. The proposed method is also shown to perform better than standard k-q acceleration schemes. CONCLUSIONS: The proposed scheme can significantly accelerate the acquisition of high spatial and angular resolution diffusion imaging by accurately reconstructing crossing fiber architectures from under-sampled data.

2D CAIPI accelerated 3D multi-slab diffusion weighted EPI combined with qModeL reconstruction for fast high resolution microstructure imaging.

PURPOSE: To develop acceleration strategies for 3D multi-slab diffusion weighted imaging (3D ms-DWI) for enabling applications that require simultaneously high spatial (1 mm isotropic) and angular (> 30 directions) resolutions. METHODS: 3D ms-DWI offers high SNR-efficiency, with the ability to achieve high isotropic spatial resolution, yet suffers from long scan-times for studies requiring high angular resolutions. We develop 6D k-q space acceleration strategies to reduce the scan-time. Specifically, we develop non-uniform 3D ky-kz under-sampling employing a shot-selective 2D CAIPI sampling approach. To achieve inter-shot phase-compensation, 2D navigators were employed that utilize the same CAIPI trajectory. An iterative model-based 3D multi-shot reconstruction was designed by incorporating phase into the forward encoding process. Additionally, the shot-selective non-uniform ky-kz CAIPI acceleration was randomized along the q-dimension. The 3D model-based multi-shot reconstruction is then extended to a joint reconstruction that simultaneously reconstructs all the q-space points, with the help of a spatial total variation and deep-learned q-space regularization. RESULTS: The proposed reconstruction is shown to achieve adequate phase-compensation in both 2D CAIPI accelerated and additional ky-kz under-sampled cases. Using retrospective under-sampling experiments, we show that k-q accelerations close a factor of 12 can be achieved with a reconstruction error < 3% for both single and multi-shell data. This translates to a scan-time reduction by 3-fold for experiments with simultaneously high spatial and angular resolutions. CONCLUSION: The proposed method facilitates the utilization of 3D ms-DWI for simultaneously high k-q resolution applications with close to 3× reduced scan-time.

Validation of deep learning techniques for quality augmentation in diffusion MRI for clinical studies.

The objective of this study is to evaluate the efficacy of deep learning (DL) techniques in improving the quality of diffusion MRI (dMRI) data in clinical applications. The study aims to determine whether the use of artificial intelligence (AI) methods in medical images may result in the loss of critical clinical information and/or the appearance of false information. To assess this, the focus was on the angular resolution of dMRI and a clinical trial was conducted on migraine, specifically between episodic and chronic migraine patients. The number of gradient directions had an impact on white matter analysis results, with statistically significant differences between groups being drastically reduced when using 21 gradient directions instead of the original 61. Fourteen teams from different institutions were tasked to use DL to enhance three diffusion metrics (FA, AD and MD) calculated from data acquired with 21 gradient directions and a b-value of 1000 s/mm2. The goal was to produce results that were comparable to those calculated from 61 gradient directions. The results were evaluated using both standard image quality metrics and Tract-Based Spatial Statistics (TBSS) to compare episodic and chronic migraine patients. The study results suggest that while most DL techniques improved the ability to detect statistical differences between groups, they also led to an increase in false positive. The results showed that there was a constant growth rate of false positives linearly proportional to the new true positives, which highlights the risk of generalization of AI-based tasks when assessing diverse clinical cohorts and training using data from a single group. The methods also showed divergent performance when replicating the original distribution of the data and some exhibited significant bias. In conclusion, extreme caution should be exercised when using AI methods for harmonization or synthesis in clinical studies when processing heterogeneous data in clinical studies, as important information may be altered, even when global metrics such as structural similarity or peak signal-to-noise ratio appear to suggest otherwise.

MoDL-MUSSELS: Model-Based Deep Learning for Multishot Sensitivity-Encoded Diffusion MRI.

We introduce a model-based deep learning architecture termed MoDL-MUSSELS for the correction of phase errors in multishot diffusion-weighted echo-planar MR images. The proposed algorithm is a generalization of the existing MUSSELS algorithm with similar performance but significantly reduced computational complexity. In this work, we show that an iterative re-weighted least-squares implementation of MUSSELS alternates between a multichannel filter bank and the enforcement of data consistency. The multichannel filter bank projects the data to the signal subspace, thus exploiting the annihilation relations between shots. Due to the high computational complexity of the self-learned filter bank, we propose replacing it with a convolutional neural network (CNN) whose parameters are learned from exemplary data. The proposed CNN is a hybrid model involving a multichannel CNN in the k-space and another CNN in the image space. The k-space CNN exploits the annihilation relations between the shot images, while the image domain network is used to project the data to an image manifold. The experiments show that the proposed scheme can yield reconstructions that are comparable to state-of-the-art methods while offering several orders of magnitude reduction in run-time.

MULTI-ECHO RECOVERY WITH FIELD INHOMOGENEITY COMPENSATION USING STRUCTURED LOW-RANK MATRIX COMPLETION.

Echo-planar imaging (EPI), which is the main workhorse of functional MRI, suffers from field inhomogeneity-induced geometric distortions. The amount of distortion is proportional to the readout duration, which restricts the maximum achievable spatial resolution. The spatially varying nature of the T 2 * decay makes it challenging for EPI schemes with a single echo time to obtain good sensitivity to functional activations in different brain regions. Despite the use of parallel MRI and multislice acceleration, the number of different echo times that can be acquired in a reasonable TR is limited. The main focus of this work is to introduce a rosette-based acquisition scheme and a structured low-rank reconstruction algorithm to overcome the above challenges. The proposed scheme exploits the exponential structure of the time series to recover distortion-free images from several echoes simultaneously.

Model-Based Deep Learning for Reconstruction of Joint k-q Under-sampled High Resolution Diffusion MRI.

We propose a model-based deep learning architecture for the reconstruction of highly accelerated diffusion magnetic resonance imaging (MRI) that enables high resolution imaging. The proposed reconstruction jointly recovers all the diffusion weighted images in a single step from a joint k-q under-sampled acquisition in a parallel MRI setting. We propose the novel use of a pre-trained denoiser as a regularizer in a model-based reconstruction for the recovery of highly under-sampled data. Specifically, we designed the denoiser based on a general diffusion MRI tissue microstructure model for multi-compartmental modeling. By using a wide range of biologically plausible parameter values for the multi-compartmental microstructure model, we simulated diffusion signal that spans the entire microstructure parameter space. A neural network was trained in an unsupervised manner using an autoencoder to learn the diffusion MRI signal subspace. We employed the autoencoder in a model-based reconstruction and show that the autoencoder provides a strong denoising prior to recover the q-space signal. We show reconstruction results on a simulated brain dataset that shows high acceleration capabilities of the proposed method.

Proton Exchange Magnetic Resonance Imaging: Current and Future Applications in Psychiatric Research.

Proton exchange provides a powerful contrast mechanism for magnetic resonance imaging (MRI). MRI techniques sensitive to proton exchange provide new opportunities to map, with high spatial and temporal resolution, compounds important for brain metabolism and function. Two such techniques, chemical exchange saturation transfer (CEST) and T1 relaxation in the rotating frame (T1ρ), are emerging as promising tools in the study of neurological and psychiatric illnesses to study brain metabolism. This review describes proton exchange for non-experts, highlights the current status of proton-exchange MRI, and presents advantages and drawbacks of these techniques compared to more traditional methods of imaging brain metabolism, including positron emission tomography (PET) and MR spectroscopy (MRS). Finally, this review highlights new frontiers for the use of CEST and T1ρ in brain research.

Calibration-Free B0 Correction of EPI Data Using Structured Low Rank Matrix Recovery.

We introduce a structured low rank algorithm for the calibration-free compensation of field inhomogeneity artifacts in echo planar imaging (EPI) MRI data. We acquire the data using two EPI readouts that differ in echo-time. Using time segmentation, we reformulate the field inhomogeneity compensation problem as the recovery of an image time series from highly undersampled Fourier measurements. The temporal profile at each pixel is modeled as a single exponential, which is exploited to fill in the missing entries. We show that the exponential behavior at each pixel, along with the spatial smoothness of the exponential parameters, can be exploited to derive a 3-D annihilation relation in the Fourier domain. This relation translates to a low rank property on a structured multi-fold Toeplitz matrix, whose entries correspond to the measured k-space samples. We introduce a fast two-step algorithm for the completion of the Toeplitz matrix from the available samples. In the first step, we estimate the null space vectors of the Toeplitz matrix using only its fully sampled rows. The null space is then used to estimate the signal subspace, which facilitates the efficient recovery of the time series of images. We finally demonstrate the proposed approach on spherical MR phantom data and human data and show that the artifacts are significantly reduced.