Improving parallel imaging by jointly reconstructing multi-contrast data

JUST-Net: Jointly unrolled cross-domain optimization based spatio-temporal reconstruction network for accelerated 3D myelin water imaging

PURPOSE

We introduced a novel reconstruction network, jointly unrolled cross-domain optimization-based spatio-temporal reconstruction network (JUST-Net), aimed at accelerating 3D multi-echo gradient-echo (mGRE) data acquisition and improving the quality of resulting myelin water imaging (MWI) maps.

METHOD

An unrolled cross-domain spatio-temporal reconstruction network was designed. The main idea is to combine frequency and spatio-temporal image feature representations and to sequentially implement convolution layers in both domains. The k-space subnetwork utilizes shared information from adjacent frames, whereas the image subnetwork applies separate convolutions in both spatial and temporal dimensions. The proposed reconstruction network was evaluated for both retrospectively and prospectively accelerated acquisition. Furthermore, it was assessed in simulation studies and real-world cases with k-space corruptions to evaluate its potential for motion artifact reduction.

RESULTS

The proposed JUST-Net enabled highly reproducible and accelerated 3D mGRE acquisition for whole-brain MWI, reducing the acquisition time from fully sampled 15:23 to 2:22 min within a 3-min reconstruction time. The normalized root mean squared error of the reconstructed mGRE images increased by less than 4.0%, and the correlation coefficients for MWI showed a value of over 0.68 when compared to the fully sampled reference. Additionally, the proposed method demonstrated a mitigating effect on both simulated and clinical motion-corrupted cases.

CONCLUSION

The proposed JUST-Net has demonstrated the capability to achieve high acceleration factors for 3D mGRE-based MWI, which is expected to facilitate widespread clinical applications of MWI.

Joint MAPLE: Accelerated joint T1 and T 2 * $$ {{\mathrm{T}}_2}^{\ast } $$ mapping with scan-specific self-supervised networks

PURPOSE

Quantitative MRI finds important applications in clinical and research studies. However, it is encoding intensive and may suffer from prohibitively long scan times. Accelerated MR parameter mapping techniques have been developed to help address these challenges. Here, an accelerated joint T1,T 2 * $$ {{\mathrm{T}}_2}^{\ast } $$ , frequency and proton density mapping technique with scan-specific self-supervised network reconstruction is proposed to synergistically combine parallel imaging, model-based, and deep learning approaches to speed up parameter mapping.

METHODS

Proposed framework, Joint MAPLE, includes parallel imaging, signal modeling, and data consistency blocks which are optimized jointly in a combined loss function. A scan-specific self-supervised reconstruction is embedded into the framework, which takes advantage of multi-contrast data from a multi-echo, multi-flip angle, gradient echo acquisition.

RESULTS

In comparison with parallel reconstruction techniques powered by low-rank methods, emerging scan specific networks, and model-basedT 2 * $$ {{\mathrm{T}}_2}^{\ast } $$ estimation approaches, the proposed framework reduces the reconstruction error in parameter maps by approximately two-fold on average at acceleration rates as high as R = 16 with uniform sampling. It can outperform evaluated parallel reconstruction techniques up to four-fold on average in the presence of challenging sub-sampling masks. It is observed that Joint MAPLE performs well at extreme acceleration rates of R = 25 and R = 36 with error values less than 20%.

CONCLUSION

Joint MAPLE enables higher fidelity parameter estimation at high acceleration rates by synergistically combining parallel imaging and model-based parameter mapping and exploiting multi-echo, multi-flip angle datasets. Utilizing a scan-specific self-supervised reconstruction obviates the need for large data sets for training while improving the parameter estimation ability.

Time-efficient, high-resolution 3T whole-brain relaxometry using 3D-QALAS with wave-CAIPI readouts

PURPOSE

Volumetric, high-resolution, quantitative mapping of brain-tissue relaxation properties is hindered by long acquisition times and SNR challenges. This study combines time-efficient wave-controlled aliasing in parallel imaging (wave-CAIPI) readouts with the 3D quantification using an interleaved Look-Locker acquisition sequence with a T2 preparation pulse (3D-QALAS), enabling full-brain quantitative T1 , T2 , and proton density (PD) maps at 1.15-mm3 isotropic voxels in 3 min.

METHODS

Wave-CAIPI readouts were embedded in the standard 3D-QALAS encoding scheme, enabling full-brain quantitative parameter maps (T1 , T2 , and PD) at acceleration factors of R = 3 × 2 with minimum SNR loss due to g-factor penalties. The quantitative parameter maps were estimated using a dictionary-based mapping algorithm incorporating inversion efficiency and B1 -field inhomogeneity effects. The parameter maps using the accelerated protocol were quantitatively compared with those obtained from the conventional 3D-QALAS sequence using GRAPPA acceleration of R = 2 in the ISMRM/NIST phantom, and in 10 healthy volunteers.

RESULTS

When tested in both the ISMRM/NIST phantom and 10 healthy volunteers, the quantitative maps using the accelerated protocol showed excellent agreement against those obtained from conventional 3D-QALAS at RGRAPPA  = 2.

CONCLUSION

Three-dimensional QALAS enhanced with wave-CAIPI readouts enables time-efficient, full-brain quantitative T1 , T2 , and PD mapping at 1.15 mm3 in 3 min at R = 3 × 2 acceleration. The quantitative maps obtained from the accelerated wave-CAIPI 3D-QALAS protocol showed very similar values to those from the standard 3D-QALAS (R = 2) protocol, alluding to the robustness and reliability of the proposed method.

Three-dimensional echo-shifted EPI with simultaneous blip-up and blip-down acquisitions for correcting geometric distortion

PURPOSE

EPI with blip-up/down acquisition (BUDA) can provide high-quality images with minimal distortions by using two readout trains with opposing phase-encoding gradients. Because of the need for two separate acquisitions, BUDA doubles the scan time and degrades the temporal resolution when compared to single-shot EPI, presenting a major challenge for many applications, particularly fMRI. This study aims at overcoming this challenge by developing an echo-shifted EPI BUDA (esEPI-BUDA) technique to acquire both blip-up and blip-down datasets in a single shot.

METHODS

A 3D esEPI-BUDA pulse sequence was designed by using an echo-shifting strategy to produce two EPI readout trains. These readout trains produced a pair of k-space datasets whose k-space trajectories were interleaved with opposite phase-encoding gradient directions. The two k-space datasets were separately reconstructed using a 3D SENSE algorithm, from which time-resolved B0 -field maps were derived using TOPUP in FSL and then input into a forward model of joint parallel imaging reconstruction to correct for geometric distortion. In addition, Hankel structured low-rank constraint was incorporated into the reconstruction framework to improve image quality by mitigating the phase errors between the two interleaved k-space datasets.

RESULTS

The 3D esEPI-BUDA technique was demonstrated in a phantom and an fMRI study on healthy human subjects. Geometric distortions were effectively corrected in both phantom and human brain images. In the fMRI study, the visual activation volumes and their BOLD responses were comparable to those from conventional 3D echo-planar images.

CONCLUSION

The improved imaging efficiency and dynamic distortion correction capability afforded by 3D esEPI-BUDA are expected to benefit many EPI applications.

High-fidelity mesoscale in-vivo diffusion MRI through gSlider-BUDA and circular EPI with S-LORAKS reconstruction

PURPOSE

To develop a high-fidelity diffusion MRI acquisition and reconstruction framework with reduced echo-train-length for less T2* image blurring compared to typical highly accelerated echo-planar imaging (EPI) acquisitions at sub-millimeter isotropic resolution.

METHODS

We first proposed a circular-EPI trajectory with partial Fourier sampling on both the readout and phase-encoding directions to minimize the echo-train-length and echo time. We then utilized this trajectory in an interleaved two-shot EPI acquisition with reversed phase-encoding polarity, to aid in the correction of off-resonance-induced image distortions and provide complementary k-space coverage in the missing partial Fourier regions. Using model-based reconstruction with structured low-rank constraint and smooth phase prior, we corrected the shot-to-shot phase variations across the two shots and recover the missing k-space data. Finally, we combined the proposed acquisition/reconstruction framework with an SNR-efficient RF-encoded simultaneous multi-slab technique, termed gSlider, to achieve high-fidelity 720 µm and 500 µm isotropic resolution in-vivo diffusion MRI.

RESULTS

Both simulation and in-vivo results demonstrate the effectiveness of the proposed acquisition and reconstruction framework to provide distortion-corrected diffusion imaging at the mesoscale with markedly reduced T2*-blurring. The in-vivo results of 720 µm and 500 µm datasets show high-fidelity diffusion images with reduced image blurring and echo time using the proposed approaches.

CONCLUSIONS

The proposed method provides high-quality distortion-corrected diffusion-weighted images with ∼40% reduction in the echo-train-length and T2* blurring at 500µm-isotropic-resolution compared to standard multi-shot EPI.

Calibrationless reconstruction of uniformly‐undersampled multi‐channel MR data with deep learning estimated ESPIRiT maps

PURPOSE

To develop a truly calibrationless reconstruction method that derives An Eigenvalue Approach to Autocalibrating Parallel MRI (ESPIRiT) maps from uniformly-undersampled multi-channel MR data by deep learning.

METHODS

ESPIRiT, one commonly used parallel imaging reconstruction technique, forms the images from undersampled MR k-space data using ESPIRiT maps that effectively represents coil sensitivity information. Accurate ESPIRiT map estimation requires quality coil sensitivity calibration or autocalibration data. We present a U-Net based deep learning model to estimate the multi-channel ESPIRiT maps directly from uniformly-undersampled multi-channel multi-slice MR data. The model is trained using fully-sampled multi-slice axial brain datasets from the same MR receiving coil system. To utilize subject-coil geometric parameters available for each dataset, the training imposes a hybrid loss on ESPIRiT maps at the original locations as well as their corresponding locations within the standard reference multi-slice axial stack. The performance of the approach was evaluated using publicly available T1-weighed brain and cardiac data.

RESULTS

The proposed model robustly predicted multi-channel ESPIRiT maps from uniformly-undersampled k-space data. They were highly comparable to the reference ESPIRiT maps directly computed from 24 consecutive central k-space lines. Further, they led to excellent ESPIRiT reconstruction performance even at high acceleration, exhibiting a similar level of errors and artifacts to that by using reference ESPIRiT maps.

CONCLUSION

A new deep learning approach is developed to estimate ESPIRiT maps directly from uniformly-undersampled MR data. It presents a general strategy for calibrationless parallel imaging reconstruction through learning from the coil and protocol-specific data.

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.

Iterative training of robust k-space interpolation networks for improved image reconstruction with limited scan specific training samples

PURPOSE

To evaluate an iterative learning approach for enhanced performance of robust artificial-neural-networks for k-space interpolation (RAKI), when only a limited amount of training data (auto-calibration signals [ACS]) are available for accelerated standard 2D imaging.

METHODS

In a first step, the RAKI model was tailored for the case of limited training data amount. In the iterative learning approach (termed iterative RAKI [iRAKI]), the tailored RAKI model is initially trained using original and augmented ACS obtained from a linear parallel imaging reconstruction. Subsequently, the RAKI convolution filters are refined iteratively using original and augmented ACS extracted from the previous RAKI reconstruction. Evaluation was carried out on 200 retrospectively undersampled in vivo datasets from the fastMRI neuro database with different contrast settings.

RESULTS

For limited training data (18 and 22 ACS lines for R = 4 and R = 5, respectively), iRAKI outperforms standard RAKI by reducing residual artifacts and yields better noise suppression when compared to standard parallel imaging, underlined by quantitative reconstruction quality metrics. Additionally, iRAKI shows better performance than both GRAPPA and standard RAKI in case of pre-scan calibration with varying contrast between training- and undersampled data.

CONCLUSION

RAKI benefits from the iterative learning approach, which preserves the noise suppression feature, but requires less original training data for the accurate reconstruction of standard 2D images thereby improving net acceleration.

Virtual coil augmentation for MR coil extrapoltion via deep learning

Magnetic resonance imaging (MRI) is a widely used medical imaging modality. However, due to the limitations in hardware, scan time, and throughput, it is often clinically challenging to obtain high-quality MR images. In this article, we propose a method of using artificial intelligence to expand the coils to achieve the goal of generating the virtual coils. The main characteristic of our work is utilizing dummy variable technology to expand/extrapolate the receive coils in both image and k-space domains. The high-dimensional information formed by coil expansion is used as the prior information to improve the reconstruction performance of parallel imaging. Two main components are incorporated into the network design, namely variable augmentation technology and sum of squares (SOS) objective function. Variable augmentation provides the network with more high-dimensional prior information, which is helpful for the network to extract the deep feature information of the data. The SOS objective function is employed to solve the deficiency of k-space data training while speeding up convergence. Experimental results demonstrated its great potentials in accelerating parallel imaging reconstruction.

A unified model for reconstruction and R2* mapping of accelerated 7T data using the quantitative recurrent inference machine

Quantitative MRI (qMRI) acquired at the ultra-high field of 7 Tesla has been used in visualizing and analyzing subcortical structures. qMRI relies on the acquisition of multiple images with different scan settings, leading to extended scanning times. Data redundancy and prior information from the relaxometry model can be exploited by deep learning to accelerate the imaging process. We propose the quantitative Recurrent Inference Machine (qRIM), with a unified forward model for joint reconstruction and R₂^*-mapping from sparse data, embedded in a Recurrent Inference Machine (RIM), an iterative inverse problem-solving network. To study the dependency of the proposed extension of the unified forward model to network architecture, we implemented and compared a quantitative End-to-End Variational Network (qE2EVN). Experiments were performed with high-resolution multi-echo gradient echo data of the brain at 7T of a cohort study covering the entire adult life span. The error in reconstructed R₂^* from undersampled data relative to reference data significantly decreased for the unified model compared to sequential image reconstruction and parameter fitting using the RIM. With increasing acceleration factor, an increasing reduction in the reconstruction error was observed, pointing to a larger benefit for sparser data. Qualitatively, this was following an observed reduction of image blurriness in R₂^*-maps. In contrast, when using the U-Net as network architecture, a negative bias in R₂^* in selected regions of interest was observed. Compressed Sensing rendered accurate, but less precise estimates of R₂^*. The qE2EVN showed slightly inferior reconstruction quality compared to the qRIM but better quality than the U-Net and Compressed Sensing. Subcortical maturation over age measured by a linearly increasing interquartile range of R₂^* in the striatum was preserved up to an acceleration factor of 9. With the integrated prior of the unified forward model, the proposed qRIM can exploit the redundancy among repeated measurements and shared information between tasks, facilitating relaxometry in accelerated MRI.

Wave-Encoded Model-Based Deep Learning for Highly Accelerated Imaging with Joint Reconstruction

A recently introduced model-based deep learning (MoDL) technique successfully incorporates convolutional neural network (CNN)-based regularizers into physics-based parallel imaging reconstruction using a small number of network parameters. Wave-controlled aliasing in parallel imaging (CAIPI) is an emerging parallel imaging method that accelerates imaging acquisition by employing sinusoidal gradients in the phase- and slice/partition-encoding directions during the readout to take better advantage of 3D coil sensitivity profiles. We propose wave-encoded MoDL (wave-MoDL) combining the wave-encoding strategy with unrolled network constraints for highly accelerated 3D imaging while enforcing data consistency. We extend wave-MoDL to reconstruct multicontrast data with CAIPI sampling patterns to leverage similarity between multiple images to improve the reconstruction quality. We further exploit this to enable rapid quantitative imaging using an interleaved look-locker acquisition sequence with T₂ preparation pulse (3D-QALAS). Wave-MoDL enables a 40 s MPRAGE acquisition at 1 mm resolution at 16-fold acceleration. For quantitative imaging, wave-MoDL permits a 1:50 min acquisition for T₁, T₂, and proton density mapping at 1 mm resolution at 12-fold acceleration, from which contrast-weighted images can be synthesized as well. In conclusion, wave-MoDL allows rapid MR acquisition and high-fidelity image reconstruction and may facilitate clinical and neuroscientific applications by incorporating unrolled neural networks into wave-CAIPI reconstruction.

Virtual Conjugate Coil for Improving KerNL Reconstruction

Ker NL is a general kernel-based framework for auto calibrated reconstruction method, which does not need any explicit formulas of the kernel function for characterizing nonlinear relationships between acquired and unacquired k-space data. It is non-iterative without requiring a large amount of computational costs. Since the limited autocalibration signals (ACS) are acquired to perform KerNL calibration and the calibration suffers from the overfitting problem, more training data can improve the kernel model accuracy. In this work, virtual conjugate coil data are incorporated into the KerNL calibration and estimation process for enhancing reconstruction performance. Experimental results show that the proposed method can further suppress noise and aliasing artifacts with fewer ACS data and higher acceleration factors. Computation efficiency is still retained to keep fast reconstruction with the random projection.

Simultaneously optimizing sampling pattern for joint acceleration of multi-contrast MRI using model-based deep learning

BACKGROUND

Acceleration of MR imaging (MRI) is a popular research area, and usage of deep learning for acceleration has become highly widespread in the MR community. Joint acceleration of multiple-acquisition MRI was proven to be effective over a single-acquisition approach. Also, optimization in the sampling pattern demonstrated its advantage over conventional undersampling pattern. However, optimizing the sampling patterns for joint acceleration of multiple-acquisition MRI has not been investigated well.

PURPOSE

To develop a model-based deep learning scheme to optimize sampling patterns for a joint acceleration of multi-contrast MRI.

METHODS

The proposed scheme combines sampling pattern optimization and multi-contrast MRI reconstruction. It was extended from the physics-guided method of the joint model-based deep learning (J-MoDL) scheme to optimize the separate sampling pattern for each of multiple contrasts simultaneously for their joint reconstruction. Tests were performed with three contrasts of T2-weighted, FLAIR, and T1-weighted images. The proposed multi-contrast method was compared to (i) single-contrast method with sampling optimization (baseline J-MoDL), (ii) multi-contrast method without sampling optimization, and (iii) multi-contrast method with single common sampling optimization for all contrasts. The optimized sampling patterns were analyzed for sampling location overlap across contrasts. The scheme was also tested in a data-driven scenario, where the inversion between input and label was learned from the under-sampled data directly and tested on knee datasets for generalization test.

RESULTS

The proposed scheme demonstrated a quantitative and qualitative advantage over the single-contrast scheme with sampling pattern optimization and the multi-contrast scheme without sampling pattern optimization. Optimizing the separate sampling pattern for each of the multi-contrasts was superior to optimizing only one common sampling pattern for all contrasts. The proposed scheme showed less overlap in sampling locations than the single-contrast scheme. The main hypothesis was also held in the data-driven situation as well. The brain-trained model worked well on the knee images, demonstrating its generalizability.

CONCLUSION

Our study introduced an effective scheme that combines the sampling optimization and the multi-contrast acceleration. The seamless combination resulted in superior performance over the other existing methods.

Accelerated 3D myelin water imaging using joint spatio-temporal reconstruction

PURPOSE

To enable acceleration in 3D multi-echo gradient echo (mGRE) acquisition for myelin water imaging (MWI) by combining joint parallel imaging (JPI) and joint deep learning (JDL).

METHODS

We implemented a multistep reconstruction process using both advanced parallel imaging and deep learning network which can utilize joint spatiotemporal components between the multi-echo images to further accelerate 3D mGRE acquisition for MWI. In the first step, JPI was performed to estimate missing k-space lines. Next, JDL was implemented to reduce residual artifacts and produce high-fidelity reconstruction by using variable splitting optimization consisting of spatiotemporal denoiser block, data consistency block, and weighted average block. The proposed method was evaluated for MWI with 2D Cartesian uniform under-sampling for each echo, enabling scan times of up to approximately 2 min for 2 mm × 2 mm × 2 mm $2\ {\rm mm} \times 2\ {\rm mm} \times 2\ {\rm mm}$ 3D coverage.

RESULTS

The proposed method showed acceptable MWI quality with improved quantitative values compared to both JPI and JDL methods individually. The improved performance of the proposed method was demonstrated by the low normalized mean-square error and high-frequency error norm values of the reconstruction with high similarity to the fully sampled MWI.

CONCLUSION

Joint spatiotemporal reconstruction approach by combining JPI and JDL can achieve high acceleration factors for 3D mGRE-based MWI.

Multi-echo balanced SSFP with a sequential phase-encoding order for functional MR imaging at 7T

PURPOSE

To develop a 2D multi-echo passband balanced SSFP (bSSFP) sequence using an echo-train readout with a sequential phase-encoding order (sequential multi-echo bSSFP), and evaluate its performance in fast functional brain imaging at 7 T.

METHODS

As images of sequential multi-echo bSSFP exhibit multiple ghosts due to periodic k-space modulations, a GRAPPA-based reconstruction method was proposed to eliminate ghosting artifacts. MRI experiments were performed to compare the image quality of multi-echo bSSFP and conventional single-echo bSSFP. Submillimeter-resolution fMRI using a checkerboard visual stimulus was conducted to compare the activation characteristics of multi-echo bSSFP, conventional single-echo bSSFP and standard gradient-echo EPI (GE-EPI).

RESULTS

A higher mean structural similarity index was found between images of single-echo bSSFP and multi-echo bSSFP with a shorter echo train length (ETL). Multi-echo bSSFP (ETL = 3) showed higher temporal SNR (tSNR) values than GRAPPA-accelerated single-echo bSSFP (R = 2). In submillimeter-resolution fMRI experiments, multi-echo bSSFP (ETL = 3) approached the imaging speed of GRAPPA-accelerated single-echo bSSFP (R = 2), but without tSNR penalty and reduced activation due to acceleration. The median t-value and the number of significantly activated voxels were comparable between GE-EPI and multi-echo bSSFP (ETL = 3) that provides virtually distortion-free functional images and inherits the activation patterns of conventional bSSFP.

CONCLUSION

Sequential multi-echo bSSFP (ETL = 3) is suitable for fast fMRI with submillimeter in-plane resolution, and offers an option to accelerate bSSFP imaging without tSNR penalty like parallel imaging.

Varying undersampling directions for accelerating multiple acquisition magnetic resonance imaging

In this study, we propose a new sampling strategy for efficiently accelerating multiple acquisition MRI. The new sampling strategy is to obtain data along different phase-encoding directions across multiple acquisitions. The proposed sampling strategy was evaluated in multicontrast MR imaging (T1, T2, proton density) and multiple phase-cycled (PC) balanced steady-state free precession (bSSFP) imaging by using convolutional neural networks with central and random sampling patterns. In vivo MRI acquisitions as well as a public database were used to test the concept. Based on both visual inspection and quantitative analysis, the proposed sampling strategy showed better performance than sampling along the same phase-encoding direction in both multicontrast MR imaging and multiple PC-bSSFP imaging, regardless of sampling pattern (central, random) or datasets (public, retrospective and prospective in vivo). For the prospective in vivo applications, acceleration was performed by sampling along different phase-encoding directions at the time of acquisition with a conventional rectangular field of view, which demonstrated the advantage of the proposed sampling strategy in the real environment. Preliminary trials on compressed sensing (CS) also demonstrated improvement of CS with the proposed idea. Sampling along different phase-encoding directions across multiple acquisitions is advantageous for accelerating multiacquisition MRI, irrespective of sampling pattern or datasets, with further improvement through transfer learning.

Arterial spin labeling MR image denoising and reconstruction using unsupervised deep learning

Arterial spin labeling (ASL) imaging is a powerful magnetic resonance imaging technique that allows to quantitatively measure blood perfusion non-invasively, which has great potential for assessing tissue viability in various clinical settings. However, the clinical applications of ASL are currently limited by its low signal-to-noise ratio (SNR), limited spatial resolution, and long imaging time. In this work, we propose an unsupervised deep learning-based image denoising and reconstruction framework to improve the SNR and accelerate the imaging speed of high resolution ASL imaging. The unique feature of the proposed framework is that it does not require any prior training pairs but only the subject's own anatomical prior, such as T1-weighted images, as network input. The neural network was trained from scratch in the denoising or reconstruction process, with noisy images or sparely sampled k-space data as training labels. Performance of the proposed method was evaluated using in vivo experiment data obtained from 3 healthy subjects on a 3T MR scanner, using ASL images acquired with 44-min acquisition time as the ground truth. Both qualitative and quantitative analyses demonstrate the superior performance of the proposed txtc framework over the reference methods. In summary, our proposed unsupervised deep learning-based denoising and reconstruction framework can improve the image quality and accelerate the imaging speed of ASL imaging.

SNR Enhancement for Multi-TE MRSI Using Joint Low-Dimensional Model and Spatial Constraints

We present a novel method to enhance the SNR for multi-TE MR spectroscopic imaging (MRSI) data by integrating learned nonlinear low-dimensional model and spatial constraints. A deep complex convolutional autoen-coder (DCCAE) was developed to learn a nonlinear low-dimensional representation of the high-dimensional multi-TE 1 H spectroscopy signals. The learned model significantly reduces the data dimension thus serving as an effective constraint for noise reduction. A reconstruction formulation was proposed to integrate the spatiospectral encoding model, the learned model, and a spatial constraint for an SNR-enhancing reconstruction from multi-TE data. The proposed method has been evaluated using both numerical simulations and in vivo brain MRSI experiments. The superior denoising performance of the proposed over alternative methods was demonstrated, both qualitatively and quantitatively. In vivo multi-TE data was used to assess the improved metabolite quantification reproducibility and accuracy achieved by the proposed method. We expect the proposed SNR-enhancing reconstruction to enable faster and/or higher-resolution multi-TE 1 H-MRSI of the brain, potentially useful for various clinical applications.

Automated Parameter Selection for Accelerated MRI Reconstruction via Low-Rank Modeling of Local k-Space Neighborhoods

PURPOSE

Image quality in accelerated MRI rests on careful selection of various reconstruction parameters. A common yet tedious and error-prone practice is to hand-tune each parameter to attain visually appealing reconstructions. Here, we propose a parameter tuning strategy to automate hybrid parallel imaging (PI) - compressed sensing (CS) reconstructions via low-rank modeling of local k-space neighborhoods (LORAKS) supplemented with sparsity regularization in wavelet and total variation (TV) domains.

METHODS

For low-rank regularization, we leverage a soft-thresholding operation based on singular values for matrix rank selection in LORAKS. For sparsity regularization, we employ Stein's unbiased risk estimate criterion to select the wavelet regularization parameter and local standard deviation of reconstructions to select the TV regularization parameter. Comprehensive demonstrations are presented on a numerical brain phantom and in vivo brain and knee acquisitions. Quantitative assessments are performed via PSNR, SSIM and NMSE metrics.

RESULTS

The proposed hybrid PI-CS method improves reconstruction quality compared to PI-only techniques, and it achieves on par image quality to reconstructions with brute-force optimization of reconstruction parameters. These results are prominent across several different datasets and the range of examined acceleration rates.

CONCLUSION

A data-driven parameter tuning strategy to automate hybrid PI-CS reconstructions is presented. The proposed method achieves reliable reconstructions of accelerated multi-coil MRI datasets without the need for exhaustive hand-tuning of reconstruction parameters.

Rapid simultaneous acquisition of macromolecular tissue volume, susceptibility, and relaxometry maps

PURPOSE

A major obstacle to the clinical implementation of quantitative MR is the lengthy acquisition time required to derive multi-contrast parametric maps. We sought to reduce the acquisition time for QSM and macromolecular tissue volume by acquiring both contrasts simultaneously by leveraging their redundancies. The joint virtual coil concept with GRAPPA (JVC-GRAPPA) was applied to reduce acquisition time further.

METHODS

Three adult volunteers were imaged on a 3 Tesla scanner using a multi-echo 3D GRE sequence acquired at 3 head orientations. Macromolecular tissue volume, QSM, R2∗ , T₁ , and proton density maps were reconstructed. The same sequence (GRAPPA R = 4) was performed in subject 1 with a single head orientation for comparison. Fully sampled data was acquired in subject 2, from which retrospective undersampling was performed (R = 6 GRAPPA and R = 9 JVC-GRAPPA). Prospective undersampling was performed in subject 3 (R = 6 GRAPPA and R = 9 JVC-GRAPPA) using gradient blips to shift k-space sampling in later echoes.

RESULTS

Subject 1's multi-orientation and single-orientation macromolecular tissue volume maps were not significantly different based on RMSE. For subject 2, the retrospectively undersampled JVC-GRAPPA and GRAPPA generated similar results as fully sampled data. This approach was validated with the prospectively undersampled images in subject 3. Using QSM, R2∗ , and macromolecular tissue volume, the contributions of myelin and iron content to susceptibility were estimated.

CONCLUSION

We have developed a novel strategy to simultaneously acquire data for the reconstruction of 5 intrinsically coregistered 1-mm isotropic resolution multi-parametric maps, with a scan time of 6 min using JVC-GRAPPA.

Isotropic multichannel total variation framework for joint reconstruction of multicontrast parallel MRI

Purpose: To develop a synergistic image reconstruction framework that exploits multicontrast (MC), multicoil, and compressed sensing (CS) redundancies in magnetic resonance imaging (MRI). Approach: CS, MC acquisition, and parallel imaging (PI) have been individually well developed, but the combination of the three has not been equally well studied, much less the potential benefits of isotropy within such a setting. Inspired by total variation theory, we introduce an isotropic MC image regularizer and attain its full potential by integrating it into compressed MC multicoil MRI. A convex optimization problem is posed to model the new variational framework and a first-order algorithm is developed to solve the problem. Results: It turns out that the proposed isotropic regularizer outperforms many of the state-of-the-art reconstruction methods not only in terms of rotation-invariance preservation of symmetrical features, but also in suppressing noise or streaking artifacts, which are normally encountered in PI methods at aggressive undersampling rates. Moreover, the new framework significantly prevents intercontrast leakage of contrast-specific details, which seems to be a difficult situation to handle for some variational and low-rank MC reconstruction approaches. Conclusions: The new framework is a viable option for image reconstruction in fast protocols of MC parallel MRI, potentially reducing patient discomfort in otherwise long and time-consuming scans.

NeuroMix-A single-scan brain exam

PURPOSE

Implement a fast, motion-robust pulse sequence that acquires T₁ -weighted, T₂ -weighted, T₂ ^* -weighted, T₂ fluid-attenuated inversion recovery, and DWI data in one run with only one prescription and one prescan.

METHODS

A software framework was developed that configures and runs several sequences in one main sequence. Based on that framework, the NeuroMix sequence was implemented, containing motion robust single-shot sequences using EPI and fast spin echo (FSE) readouts (without EPI distortions). Optional multi-shot sequences that provide better contrast, higher resolution, or isotropic resolution could also be run within the NeuroMix sequence. An optimized acquisition order was implemented that minimizes times where no data is acquired.

RESULTS

NeuroMix is customizable and takes between 1:20 and 4 min for a full brain scan. A comparison with the predecessor EPIMix revealed significant improvements for T₂ -weighted and T₂ fluid-attenuated inversion recovery, while taking only 8 s longer for a similar configuration. The optional contrasts were less motion robust but offered a significant increase in quality, detail, and contrast. Initial clinical scans on 1 pediatric and 1 adult patient showed encouraging image quality.

CONCLUSION

The single-shot FSE readouts for T₂ -weighted and T₂ fluid-attenuated inversion recovery and the optional multishot FSE and 3D-EPI contrasts significantly increased diagnostic value compared with EPIMix, allowing NeuroMix to be considered as a standalone brain MRI application.

Quad-Contrast Imaging: Simultaneous Acquisition of Four Contrast-Weighted Images (PD-Weighted, T₂-Weighted, PD-FLAIR and T₂-FLAIR Images) With Synthetic T₁-Weighted Image, T₁- and T₂-Maps

Magnetic resonance imaging (MRI) can provide multiple contrast-weighted images using different pulse sequences and protocols. However, a long acquisition time of the images is a major challenge. To address this limitation, a new pulse sequence referred to as quad-contrast imaging is presented. The quad-contrast sequence enables the simultaneous acquisition of four contrast-weighted images (proton density (PD)-weighted, T₂-weighted, PD-fluid attenuated inversion recovery (FLAIR), and T₂-FLAIR), and the synthesis of T₁-weighted images and T₁- and T₂-maps in a single scan. The scan time is less than 6 min and is further reduced to 2 min 50 s using a deep learning-based parallel imaging reconstruction. The natively acquired quad contrasts demonstrate high quality images, comparable to those from the conventional scans. The deep learning-based reconstruction successfully reconstructed highly accelerated data (acceleration factor 6), reporting smaller normalized root mean squared errors (NRMSEs) and higher structural similarities (SSIMs) than those from conventional generalized autocalibrating partially parallel acquisitions (GRAPPA)-reconstruction (mean NRMSE of 4.36% vs. 10.54% and mean SSIM of 0.990 vs. 0.953). In particular, the FLAIR contrast is natively acquired and does not suffer from lesion-like artifacts at the boundary of tissue and cerebrospinal fluid, differentiating the proposed method from synthetic imaging methods. The quad-contrast imaging method may have the potentials to be used in a clinical routine as a rapid diagnostic tool.

Accelerating whole-heart 3D T2 mapping: Impact of undersampling strategies and reconstruction techniques

PURPOSE

We aim to determine an advantageous approach for the acceleration of high spatial resolution 3D cardiac T2 relaxometry data by comparing the performance of different undersampling patterns and reconstruction methods over a range of acceleration rates.

METHODS

Multi-volume 3D high-resolution cardiac images were acquired fully and undersampled retrospectively using 1) optimal CAIPIRINHA and 2) a variable density random (VDR) sampling. Data were reconstructed using 1) multi-volume sensitivity encoding (SENSE), 2) joint-sparsity SENSE and 3) model-based SENSE. Four metrics were calculated on 3 naïve swine and 8 normal human subjects over a whole left-ventricular region of interest: root-mean-square error (RMSE) of image signal intensity, RMSE of T2, the bias of mean T2, and standard deviation (SD) of T2. Fully sampled data and volume-by-volume SENSE with standard equally spaced undersampling were used as references. The Jaccard index calculated from one swine with acute myocardial infarction (MI) was used to demonstrate preservation of segmentation of edematous tissues with elevated T2.

RESULTS

In naïve swine and normal human subjects, all methods had similar performance when the net reduction factor (Rnet) <2.5. VDR sampling with model-based SENSE showed the lowest RMSEs (10.5%-14.2%) and SDs (+1.7-2.4 ms) of T2 when Rnet>2.5, while VDR sampling with the joint-sparsity SENSE had the lowest bias of mean T2 (0.0-1.1ms) when Rnet>3. The RMSEs of parametric T2 values (9.2%-24.6%) were larger than for image signal intensities (5.2%-18.4%). In the swine with MI, VDR sampling with either joint-sparsity or model-based SENSE showed consistently higher Jaccard index for all Rnet (0.71-0.50) than volume-by-volume SENSE (0.68-0.30).

CONCLUSIONS

Retrospective exploration of undersampling and reconstruction in 3D whole-heart T2 parametric mapping revealed that maps were more sensitive to undersampling than images, presenting a more stringent limiting factor on Rnet. The combination of VDR sampling patterns with model-based or joint-sparsity SENSE reconstructions were more robust for Rnet>3.

Accelerated calibrationless parallel transmit mapping using joint transmit and receive low-rank tensor completion

Purpose

To evaluate an algorithm for calibrationless parallel imaging to reconstruct undersampled parallel transmit field maps for the body and brain.

Methods

Using a combination of synthetic data and in vivo measurements from brain and body, 3 different approaches to a joint transmit and receive low-rank tensor completion algorithm are evaluated. These methods included: 1) virtual coils using the product of receive and transmit sensitivities, 2) joint-receiver coils that enforces a low rank structure across receive coils of all transmit modes, and 3) transmit low rank that uses a low rank structure for both receive and transmit modes simultaneously. The performance of each is investigated for different noise levels and different acceleration rates on an 8-channel parallel transmit 7 Tesla system.

Results

The virtual coils method broke down with increasing noise levels or acceleration rates greater than 2, producing normalized RMS error greater than 0.1. The joint receiver coils method worked well up to acceleration factors of 4, beyond which the normalized RMS error exceeded 0.1. Transmit low rank enabled an eightfold acceleration, with most normalized RMS errors remaining below 0.1.

Conclusion

This work demonstrates that undersampling factors of up to eightfold are feasible for transmit array mapping and can be reconstructed using calibrationless parallel imaging methods.

Efficient T2 mapping with blip-up/down EPI and gSlider-SMS (T2 -BUDA-gSlider)

PURPOSE

To rapidly obtain high isotropic-resolution T₂ maps with whole-brain coverage and high geometric fidelity.

METHODS

A T₂ blip-up/down EPI acquisition with generalized slice-dithered enhanced resolution (T₂ -BUDA-gSlider) is proposed. A RF-encoded multi-slab spin-echo (SE) EPI acquisition with multiple TEs was developed to obtain high SNR efficiency with reduced TR. This was combined with an interleaved 2-shot EPI acquisition using blip-up/down phase encoding. An estimated field map was incorporated into the joint multi-shot EPI reconstruction with a structured low rank constraint to achieve distortion-free and robust reconstruction for each slab without navigation. A Bloch simulated subspace model was integrated into gSlider reconstruction and used for T₂ quantification.

RESULTS

In vivo results demonstrated that the T₂ values estimated by the proposed method were consistent with gold standard spin-echo acquisition. Compared to the reference 3D fast spin echo (FSE) images, distortion caused by off-resonance and eddy current effects were effectively mitigated.

CONCLUSION

BUDA-gSlider SE-EPI acquisition and gSlider-subspace joint reconstruction enabled distortion-free whole-brain T₂ mapping in 2 min at ~1 mm³ isotropic resolution, which could bring significant benefits to related clinical and neuroscience applications.

Accelerated multicontrast reconstruction for synthetic MRI using joint parallel imaging and variable splitting networks

PURPOSE

Synthetic magnetic resonance imaging (MRI) requires the acquisition of multicontrast images to estimate quantitative parameter maps, such as T₁ , T₂ , and proton density (PD). The study aims to develop a multicontrast reconstruction method based on joint parallel imaging (JPI) and joint deep learning (JDL) to enable further acceleration of synthetic MRI.

METHODS

The JPI and JDL methods are extended and combined to improve reconstruction for better-quality, synthesized images. JPI is performed as a first step to estimate the missing k-space lines, and JDL is then performed to correct and refine the previous estimate with a trained neural network. For the JDL architecture, the original variable splitting network (VS-Net) is modified and extended to form a joint variable splitting network (JVS-Net) to apply to multicontrast reconstructions. The proposed method is designed and tested for multidynamic multiecho (MDME) images with Cartesian uniform under-sampling using acceleration factors between 4 and 8.

RESULTS

It is demonstrated that the normalized root-mean-square error (nRMSE) is lower and the structural similarity index measure (SSIM) values are higher with the proposed method compared to both the JPI and JDL methods individually. The method also demonstrates the potential to produce a set of synthesized contrast-weighted images that closely resemble those from the fully sampled acquisition without erroneous artifacts.

CONCLUSION

Combining JPI and JDL enables the reconstruction of highly accelerated synthetic MRIs.

Minimizing echo and repetition times in magnetic resonance imaging using a double half-echo k-space acquisition and low-rank reconstruction

Sampling k-space asymmetrically (ie, partial Fourier sampling) in the readout direction is a common way to reduce the echo time (TE) during magnetic resonance image acquisitions. This technique requires overlap around the center of k-space to provide a calibration region for reconstruction, which limits the minimum fractional echo to ~60% before artifacts are observed. The present study describes a method for reconstructing images from exact half echoes using two separate acquisitions with reversed readout polarity, effectively providing a full line of k-space without additional data around central k-space. This approach can benefit sequences or applications that prioritize short TE, short inter-echo spacing or short repetition time. An example of the latter is demonstrated to reduce banding artifacts in balanced steady-state free precession.

Distortion-free, high-isotropic-resolution diffusion MRI with gSlider BUDA-EPI and multicoil dynamic B0 shimming

PURPOSE

We combine SNR-efficient acquisition and model-based reconstruction strategies with newly available hardware instrumentation to achieve distortion-free in vivo diffusion MRI of the brain at submillimeter-isotropic resolution with high fidelity and sensitivity on a clinical 3T scanner.

METHODS

We propose blip-up/down acquisition (BUDA) for multishot EPI using interleaved blip-up/blip-down phase encoding and incorporate B₀ forward-modeling into structured low-rank reconstruction to enable distortion-free and navigator-free diffusion MRI. We further combine BUDA-EPI with an SNR-efficient simultaneous multislab acquisition (generalized slice-dithered enhanced resolution ["gSlider"]), to achieve high-isotropic-resolution diffusion MRI. To validate gSlider BUDA-EPI, whole-brain diffusion data at 860-μm and 780-μm data sets were acquired. Finally, to improve the conditioning and minimize noise penalty in BUDA reconstruction at very high resolutions where B₀ inhomogeneity can have a detrimental effect, the level of B₀ inhomogeneity was reduced by incorporating slab-by-slab dynamic shimming with a 32-channel AC/DC coil into the acquisition. Whole-brain 600-μm diffusion data were then acquired with this combined approach of gSlider BUDA-EPI with dynamic shimming.

RESULTS

The results of 860-μm and 780-μm datasets show high geometry fidelity with gSlider BUDA-EPI. With dynamic shimming, the BUDA reconstruction's noise penalty was further alleviated. This enables whole-brain 600-μm isotropic resolution diffusion imaging with high image quality.

CONCLUSIONS

The gSlider BUDA-EPI method enables high-quality, distortion-free diffusion imaging across the whole brain at submillimeter resolution, where the use of multicoil dynamic B₀ shimming further improves reconstruction performance, which can be particularly useful at very high resolutions.

APIR4EMC: Autocalibrated parallel imaging reconstruction for extended multi-contrast imaging

PURPOSE

To improve image quality of multi-contrast imaging with the proposed Autocalibrated Parallel Imaging Reconstruction for Extended Multi-Contrast Imaging (APIR4EMC).

METHODS

APIR4EMC reconstructs multi-contrast images in an autocalibrated parallel imaging reconstruction framework by adding contrasts as virtual coils. Compensation of signal evolution along the echo train of different contrasts is performed to improve signal prediction for missing samples. As a proof of concept, we performed prospectively accelerated phantom and in-vivo brain acquisitions with T1, T1-fat saturated (Fatsat), T2, PD, and FLAIR contrasts. The k-space sampling patterns of these acquisitions were jointly optimized. Images were jointly reconstructed with the proposed APIR4EMC method as well as individually with GRAPPA. Root mean square error (RMSE) to fully sampled reference images and g-factor maps were computed for both methods in the phantom experiment. Visual evaluation was performed in the in-vivo experiment.

RESULTS

Compared to GRAPPA, APIR4EMC reduced artifacts and improved SNR of the reconstructed images in the phantom acquisitions. Quantitatively, APIR4EMC substantially reduced noise amplification (g-factor) as well as RMSE compared to GRAPPA. Signal evolution compensation reduced artifacts. In the in-vivo experiments, 1 mm³ isotropic 3D images with contrasts of T1, T1-Fatsat, T2, PD, and FLAIR were acquired in as little as 7.5 min with the acceleration factor of 9. Reconstruction quality was consistent with the phantom results.

CONCLUSION

Compared to single contrast reconstruction with GRAPPA, APIR4EMC reduces artifacts and noise amplification in accelerated multi-contrast imaging.

Joint calibrationless reconstruction of highly undersampled multicontrast MR datasets using a low-rank Hankel tensor completion framework

PURPOSE

To jointly reconstruct highly undersampled multicontrast two-dimensional (2D) datasets through a low-rank Hankel tensor completion framework.

METHODS

A multicontrast Hankel tensor completion (MC-HTC) framework is proposed to exploit the shareable information in multicontrast datasets with respect to their highly correlated image structure, common spatial support, and shared coil sensitivity for joint reconstruction. This is achieved by first organizing multicontrast k-space datasets into a single block-wise Hankel tensor. Subsequent low-rank tensor approximation via higher-order singular value decomposition (HOSVD) uses the image structural correlation by considering different contrasts as virtual channels. Meanwhile, the HOSVD imposes common spatial support and shared coil sensitivity by treating data from different contrasts as from additional k-space kernels. The missing k-space data are then recovered by iteratively performing such low-rank approximation and enforcing data consistency. This joint reconstruction framework was evaluated using multicontrast multichannel 2D human brain datasets (T₁ -weighted, T₂ -weighted, fluid-attenuated inversion recovery, and T₁ -weighted-inversion recovery) of identical image geometry with random and uniform undersampling schemes.

RESULTS

The proposed method offered high acceleration, exhibiting significantly less residual errors when compared with both single-contrast SAKE (simultaneous autocalibrating and k-space estimation) and multicontrast J-LORAKS (joint parallel-imaging-low-rank matrix modeling of local k-space neighborhoods) low-rank reconstruction. Furthermore, the MC-HTC framework was applied uniquely to Cartesian uniform undersampling by incorporating a novel complementary k-space sampling strategy where the phase-encoding direction among different contrasts is orthogonally alternated.

CONCLUSION

The proposed MC-HTC approach presents an effective tensor completion framework to jointly reconstruct highly undersampled multicontrast 2D datasets without coil-sensitivity calibration.

On the regularization of feature fusion and mapping for fast MR multi-contrast imaging via iterative networks

Multi-contrast (MC) Magnetic Resonance Imaging (MRI) of the same patient usually requires long scanning times, despite the images sharing redundant information. In this work, we propose a new iterative network that utilizes the sharable information among MC images for MRI acceleration. The proposed network has reinforced data fidelity control and anatomy guidance through an iterative optimization procedure of Gradient Descent, leading to reduced uncertainties and improved reconstruction results. Through a convolutional network, the new method incorporates a learnable regularization unit that is capable of extracting, fusing, and mapping shareable information among different contrasts. Specifically, a dilated inception block is proposed to promote multi-scale feature extractions and increase the receptive field diversity for contextual information incorporation. Lastly, an optimal MC information feeding protocol is built through the design of a complementary feature extractor block. Comprehensive experiments demonstrated the superiority of the proposed network, both qualitatively and quantitatively.

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.

Fast multicomponent 3D-T1ρ relaxometry

NMR relaxometry can provide information about the relaxation of the magnetization in different tissues, increasing our understanding of molecular dynamics and biochemical composition in biological systems. In general, tissues have complex and heterogeneous structures composed of multiple pools. As a result, bulk magnetization returns to its original state with different relaxation times, in a multicomponent relaxation. Recovering the distribution of relaxation times in each voxel is a difficult inverse problem; it is usually unstable and requires long acquisition time, especially on clinical scanners. MRI can also be viewed as an inverse problem, especially when compressed sensing (CS) is used. The solution of these two inverse problems, CS and relaxometry, can be obtained very efficiently in a synergistically combined manner, leading to a more stable multicomponent relaxometry obtained with short scan times. In this paper, we will discuss the details of this technique from the viewpoint of inverse problems.

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.

Joint multi-contrast variational network reconstruction (jVN) with application to rapid 2D and 3D imaging

PURPOSE

To improve the image quality of highly accelerated multi-channel MRI data by learning a joint variational network that reconstructs multiple clinical contrasts jointly.

METHODS

Data from our multi-contrast acquisition were embedded into the variational network architecture where shared anatomical information is exchanged by mixing the input contrasts. Complementary k-space sampling across imaging contrasts and Bunch-Phase/Wave-Encoding were used for data acquisition to improve the reconstruction at high accelerations. At 3T, our joint variational network approach across T1w, T2w and T2-FLAIR-weighted brain scans was tested for retrospective under-sampling at R = 6 (2D) and R = 4 × 4 (3D) acceleration. Prospective acceleration was also performed for 3D data where the combined acquisition time for whole brain coverage at 1 mm isotropic resolution across three contrasts was less than 3 min.

RESULTS

Across all test datasets, our joint multi-contrast network better preserved fine anatomical details with reduced image-blurring when compared to the corresponding single-contrast reconstructions. Improvement in image quality was also obtained through complementary k-space sampling and Bunch-Phase/Wave-Encoding where the synergistic combination yielded the overall best performance as evidenced by exemplary slices and quantitative error metrics.

CONCLUSION

By leveraging shared anatomical structures across the jointly reconstructed scans, our joint multi-contrast approach learnt more efficient regularizers, which helped to retain natural image appearance and avoid over-smoothing. When synergistically combined with advanced encoding techniques, the performance was further improved, enabling up to R = 16-fold acceleration with good image quality. This should help pave the way to very rapid high-resolution brain exams.

Fast and high-resolution myelin water imaging: Accelerating multi-echo GRASE with CAIPIRINHA

PURPOSE

Although several MRI methods have been explored to achieve in vivo myelin quantification, imaging the whole brain in clinically acceptable times and sufficiently high resolution remains challenging. To address this problem, this work investigates the acceleration of multi-echo T₂ acquisitions based on the multi-echo gradient and spin echo (GRASE) sequence using CAIPIRINHA undersampling and adapted k-space reordering patterns.

METHODS

A prototype multi-echo GRASE sequence supporting CAIPIRINHA parallel imaging was implemented. Multi-echo T₂ data were acquired from 12 volunteers using the implemented sequence (1.6 × 1.6 × 1.6 mm³ , 84 slices, acquisition time [TA] = 10:30 min) and a multi-echo spin echo (MESE) sequence as reference (1.6 × 1.6 × 3.2 mm³ , single-slice, TA = 5:41 min). Myelin water fraction (MWF) maps derived from both acquisitions were compared via correlation and Bland-Altman analyses. In addition, scan-rescan datasets were acquired to evaluate the repeatability of the derived maps.

RESULTS

Resulting maps from the MESE and multi-echo GRASE sequences were found to be correlated (r = 0.83). The Bland-Altman analysis revealed a mean bias of -0.2% (P = .24) with the limits of agreement ranging from -3.7% to 3.3%. The Pearson's correlation coefficient among MWF values obtained from the scan-rescan datasets was found to be 0.95 and the mean bias equal to 0.11% (P = .32), indicating good repeatability of the retrieved maps.

CONCLUSION

By combining a 3D multi-echo GRASE sequence with CAIPIRINHA sampling, whole-brain MWF maps were obtained in 10:30 min with 1.6 mm isotropic resolution. The good correlation with conventional MESE-based maps demonstrates that the implemented sequence may be a promising alternative to time-consuming MESE acquisitions.

Navigator-based reacquisition and estimation of motion-corrupted data: Application to multi-echo spin echo for carotid wall MRI

PURPOSE

To assess whether artifacts in multi-slice multi-echo spin echo neck imaging, thought to be caused by brief motion events such as swallowing, can be corrected by reacquiring corrupted central k-space data and estimating the remainder with parallel imaging.

METHODS

A single phase-encode line (ky = 0, phase-encode direction anteroposterior) navigator echo was used to identify motion-corrupted data and guide the online reacquisition. If motion corruption was detected in the 7 central k-space lines, they were replaced with reacquired data. Subsequently, GRAPPA reconstruction was trained on the updated central portion of k-space and then used to estimate the remaining motion-corrupted k-space data from surrounding uncorrupted data. Similar compressed sensing-based approaches have been used previously to compensate for respiration in cardiac imaging. The g-factor noise amplification was calculated for the parallel imaging reconstruction of data acquired with a 10-channel neck coil. The method was assessed in scans with 9 volunteers and 12 patients.

RESULTS

The g-factor analysis showed that GRAPPA reconstruction of 2 adjacent motion-corrupted lines causes high noise amplification; therefore, the number of 2-line estimations should be limited. In volunteer scans, median ghosting reduction of 24% was achieved with 2 adjacent motion-corrupted lines correction, and image quality was improved in 2 patient scans that had motion corruption close to the center of k-space.

CONCLUSION

Motion-corrupted echo-trains can be identified with a navigator echo. Combined reacquisition and parallel imaging estimation reduced motion artifacts in multi-slice MESE when there were brief motion events, especially when motion corruption was close to the center of k-space.

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.

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.

Ultimate MRI

The basic principles of Magnetic Resonance have been understood for over 70 years and a mainstay of medical imaging for over 40. At this point, it's no longer about simply porting these principles to medical imaging. But we are by no means confined to simply polishing either. Significant innovation and even revolution can come to old technologies. The recent revolution in optical microscopy shattered the resolution constraint imposed by a seemingly fundamental physical law (the diffraction limit) and reinvigorated a 500-year-old modality. Progress comes from re-examining old-ways and sidestepping underlying assumptions. This is already underway for MRI; and is fueled by advances in image reconstruction. Reconstruction increasingly employs sophisticated general models often using subtle and hopefully innocuous prior knowledge about the object. This allows a careful re-examination of some basic prerequisites for MRI such as uniform static fields, linear encoding fields, full Nyquist sampling, or even a stationary object. These powerful reconstruction tools are driving changes in acquisition strategy and basic hardware. The scanner of the future will know more about itself and its patient and his/her biology than ever before. This strategy emboldens relaxed hardware constraints and more specialized scanners, hopefully expanding the reach and value offered by MR imaging.

Accelerated whole-brain perfusion imaging using a simultaneous multislice spin-echo and gradient-echo sequence with joint virtual coil reconstruction

PURPOSE

Dynamic susceptibility contrast imaging requires high temporal sampling, which poses limits on achievable spatial coverage and resolution. Additionally, more encoding-intensive multi-echo acquisitions for quantitative imaging are desired to mitigate contrast leakage effects, which further limits spatial encoding. We present an accelerated sequence that provides whole-brain coverage at an improved spatio-temporal resolution, to allow for dynamic quantitative R2 and R2 * mapping during contrast-enhanced imaging.

METHODS

A multi-echo spin and gradient-echo sequence was implemented with simultaneous multislice acquisition. Complementary k-space sampling between repetitions and joint virtual coil reconstruction were used along with a dynamic phase-matching technique to achieve high-quality reconstruction at 9-fold acceleration, which enabled 2 × 2 × 5 mm whole-brain imaging at TR of 1.5 to 1.7 seconds. The multi-echo images from this sequence were fit to achieve quantitative R2 and R2 * maps for each repetition, and subsequently used to find perfusion measures including cerebral blood flow and cerebral blood volume.

RESULTS

Images reconstructed using joint virtual coil show improved image quality and g-factor compared with conventional reconstruction methods, resulting in improved quantitative maps with a 9-fold acceleration factor and whole-brain coverage during the dynamic perfusion acquisition.

CONCLUSION

The method presented shows the advantage of using a joint virtual coil-GRAPPA reconstruction to allow for high acceleration factors while maintaining reliable image quality for quantitative perfusion mapping, with the potential to improve tumor diagnostics and monitoring.

High-fidelity, high-isotropic-resolution diffusion imaging through gSlider acquisition with B 1 + and T1 corrections and integrated ΔB0 /Rx shim array

PURPOSE

B 1 + and T₁ corrections and dynamic multicoil shimming approaches were proposed to improve the fidelity of high-isotropic-resolution generalized slice-dithered enhanced resolution (gSlider) diffusion imaging.

METHODS

An extended reconstruction incorporating B 1 + inhomogeneity and T₁ recovery information was developed to mitigate slab-boundary artifacts in short-repetition time (TR) gSlider acquisitions. Slab-by-slab dynamic B₀ shimming using a multicoil integrated ΔB₀ /Rx shim array and high in-plane acceleration (R_inplane = 4) achieved with virtual-coil GRAPPA were also incorporated into a 1-mm isotropic resolution gSlider acquisition/reconstruction framework to achieve a significant reduction in geometric distortion compared to single-shot echo planar imaging (EPI).

RESULTS

The slab-boundary artifacts were alleviated by the proposed B 1 + and T₁ corrections compared to the standard gSlider reconstruction pipeline for short-TR acquisitions. Dynamic shimming provided >50% reduction in geometric distortion compared to conventional global second-order shimming. One-millimeter isotropic resolution diffusion data show that the typically problematic temporal and frontal lobes of the brain can be imaged with high geometric fidelity using dynamic shimming.

CONCLUSIONS

The proposed B 1 + and T₁ corrections and local-field control substantially improved the fidelity of high-isotropic-resolution diffusion imaging, with reduced slab-boundary artifacts and geometric distortion compared to conventional gSlider acquisition and reconstruction. This enabled high-fidelity whole-brain 1-mm isotropic diffusion imaging with 64 diffusion directions in 20 min using a 3T clinical scanner.

Improving GRAPPA reconstruction using joint nonlinear kernel mapped and phase conjugated virtual coils

To improve the reconstruction condition and alleviate the noise amplification of GRAPPA reconstruction by aggregating the phase conjugated and nonlinear kernel mapped coils with the original physical coil. Nonlinear GRAPPA (NL-GRAPPA) is a kernel-based non-iterative approach which can reduce noise-induced error in GRAPPA reconstruction. And virtual conjugate coil (VCC) embeds the conjugate symmetric property of k-space into GRAPPA data synthesis to improve reconstruction condition. This work proposed NL-VCC-GRAPPA to jointly utilize the nonlinear mapped virtual coil and phase conjugated virtual coil to further reduce noise amplification in parallel imaging. In vivo static and dynamic 2D imaging accelerated by uniform undersampling schemes were performed to evaluate the proposed method in terms of visual image quality, root-mean-square-error (RMSE), and geometry factor (g-factor). The effects of acceleration factors, calibration data size and kernel shape on the proposed model were also separately analyzed and discussed. The proposed method illustrated improved visual image quality evidenced by reduced retrospective RMSE and prospective g-factor comparing with conventional GRAPPA and the recently proposed iterative SENSE-LORAKS reconstructions. Although a larger amount of calibration data and smaller kernel size were required to stabilize the calibration of fourfold extended kernel for the proposed method, it was non-iterative and relatively insensitive to parameter adjustment in the applications. The proposed NL-VCC-extension to conventional GRAPPA brings visible improvements for imaging scenarios accelerated by the widely available uniform undersampling schemes in a practically efficient manner without iteration.

Projection onto Epigraph Sets for Rapid Self-Tuning Compressed Sensing MRI

The compressed sensing (CS) framework leverages the sparsity of MR images to reconstruct from the undersampled acquisitions. CS reconstructions involve one or more regularization parameters that weigh sparsity in transform domains against fidelity to acquired data. While parameter selection is critical for reconstruction quality, the optimal parameters are subject and dataset specific. Thus, commonly practiced heuristic parameter selection generalizes poorly to independent datasets. Recent studies have proposed to tune parameters by estimating the risk of removing significant image coefficients. Line searches are performed across the parameter space to identify the parameter value that minimizes this risk. Although effective, these line searches yield prolonged reconstruction times. Here, we propose a new self-tuning CS method that uses computationally efficient projections onto epigraph sets of the l₁ and total-variation norms to simultaneously achieve parameter selection and regularization. In vivo demonstrations are provided for balanced steady-state free precession, time-of-flight, and T1-weighted imaging. The proposed method achieves an order of magnitude improvement in computational efficiency over line-search methods while maintaining near-optimal parameter selection.

Statistically Segregated k-Space Sampling for Accelerating Multiple-Acquisition MRI

A central limitation of multiple-acquisition magnetic resonance imaging (MRI) is the degradation in scan efficiency as the number of distinct datasets grows. Sparse recovery techniques can alleviate this limitation via randomly undersampled acquisitions. A frequent sampling strategy is to prescribe for each acquisition a different random pattern drawn from a common sampling density. However, naive random patterns often contain gaps or clusters across the acquisition dimension that, in turn, can degrade reconstruction quality or reduce scan efficiency. To address this problem, a statistically segregated sampling method is proposed for multiple-acquisition MRI. This method generates multiple patterns sequentially while adaptively modifying the sampling density to minimize k-space overlap across patterns. As a result, it improves incoherence across acquisitions while still maintaining similar sampling density across the radial dimension of k-space. Comprehensive simulations and in vivo results are presented for phase-cycled balanced steady-state free precession and multi-echo [Formula: see text]-weighted imaging. Segregated sampling achieves significantly improved quality in both Fourier and compressed-sensing reconstructions of multiple-acquisition datasets.