General phase regularized reconstruction using phase cycling

Structured low-rank reconstruction for navigator-free water/fat separated multi-shot diffusion-weighted EPI

PURPOSE

Multi-shot diffusion-weighted EPI allows an increase in image resolution and reduced geometric distortions and can be combined with chemical-shift encoding (Dixon) to separate water/fat signals. However, such approaches suffer from physiological motion-induced shot-to-shot phase variations. In this work, a structured low-rank-based navigator-free algorithm is proposed to address the challenge of simultaneously separating water/fat signals and correcting for physiological motion-induced shot-to-shot phase variations in multi-shot EPI-based diffusion-weighted MRI.

THEORY AND METHODS

We propose an iterative, model-based reconstruction pipeline that applies structured low-rank regularization to estimate and eliminate the shot-to-shot phase variations in a data-driven way, while separating water/fat images. The algorithm is tested in different anatomies, including head-neck, knee, brain, and prostate. The performance is validated in simulations and in-vivo experiments in comparison to existing approaches.

RESULTS

In-vivo experiments and simulations demonstrated the effectiveness of the proposed algorithm compared to extra-navigated and an alternative self-navigation approach. The proposed algorithm demonstrates the capability to reconstruct in the multi-shot/Dixon hybrid space domain under-sampled datasets, using the same number of acquired EPI shots compared to conventional fat-suppression techniques but eliminating fat signals through chemical-shift encoding. In addition, partial Fourier reconstruction can also be achieved by using the concept of virtual conjugate coils in conjunction with the proposed algorithm.

CONCLUSION

The proposed algorithm effectively eliminates the shot-to-shot phase variations and separates water/fat images, making it a promising solution for future DWI on different anatomies.

Distortion-corrected image reconstruction with deep learning on an MRI-Linac

PURPOSE

MRI is increasingly utilized for image-guided radiotherapy due to its outstanding soft-tissue contrast and lack of ionizing radiation. However, geometric distortions caused by gradient nonlinearities (GNLs) limit anatomical accuracy, potentially compromising the quality of tumor treatments. In addition, slow MR acquisition and reconstruction limit the potential for effective image guidance. Here, we demonstrate a deep learning-based method that rapidly reconstructs distortion-corrected images from raw k-space data for MR-guided radiotherapy applications.

METHODS

We leverage recent advances in interpretable unrolling networks to develop a Distortion-Corrected Reconstruction Network (DCReconNet) that applies convolutional neural networks (CNNs) to learn effective regularizations and nonuniform fast Fourier transforms for GNL-encoding. DCReconNet was trained on a public MR brain dataset from 11 healthy volunteers for fully sampled and accelerated techniques, including parallel imaging (PI) and compressed sensing (CS). The performance of DCReconNet was tested on phantom, brain, pelvis, and lung images acquired on a 1.0T MRI-Linac. The DCReconNet, CS-, PI-and UNet-based reconstructed image quality was measured by structural similarity (SSIM) and RMS error (RMSE) for numerical comparisons. The computation time and residual distortion for each method were also reported.

RESULTS

Imaging results demonstrated that DCReconNet better preserves image structures compared to CS- and PI-based reconstruction methods. DCReconNet resulted in the highest SSIM (0.95 median value) and lowest RMSE (<0.04) on simulated brain images with four times acceleration. DCReconNet is over 10-times faster than iterative, regularized reconstruction methods.

CONCLUSIONS

DCReconNet provides fast and geometrically accurate image reconstruction and has the potential for MRI-guided radiotherapy applications.

Phase-regularized and displacement-regularized compressed sensing for fast magnetic resonance elastography

Liver magnetic resonance elastography (MRE) is a noninvasive stiffness measurement technique that captures the tissue displacement in the phase of the signal. To limit the scanning time to a single breath-hold, liver MRE usually involves advanced readout techniques such as simultaneous multislice (SMS) or multishot methods. Furthermore, all these readout techniques require additional in-plane acceleration using either parallel imaging capabilities, such as sensitivity encoding (SENSE), or k -space undersampling, such as compressed sensing (CS). However, these methods apply a single regularization function on the complex image. This study aims to design and evaluate methods that use separate regularization on the magnitude and phase of MRE to exploit their distinct spatiotemporal characteristics. Specifically, we introduce two compressed sensing methods. The first method, termed phase-regularized compressed sensing (PRCS), applies a two-dimensional total variation (TV) prior to the magnitude and two-dimensional wavelet regularization to the phase. The second method, termed displacement-regularized compressed sensing (DRCS), exploits the spatiotemporal redundancy using 3D total variation on the magnitude. Additionally, DRCS includes a displacement fitting function to apply wavelet regularization to the displacement phasor. Both DRCS and PRCS were evaluated with different levels of compression factors in three datasets: an in silico abdomen dataset, an in vitro tissue-mimicking phantom, and an in vivo liver dataset. The reconstructed images were compared with the full sampled reconstruction, zero-filling reconstruction, wavelet-regularized compressed sensing, and a low rank plus sparse reconstruction. The metrics used for quantitative evaluation were the structural similarity index (SSIM) of magnitude (M-SSIM), displacement (D-SSIM), and shear modulus (S-SSIM), and mean shear modulus. Results from highly undersampled in silico and in vitro datasets demonstrate that the DRCS method provides higher reconstruction quality than the conventional compressed sensing method for a wide range of stiffness values. Notably, DRCS provides 24% and 22% increase in D-SSIM compared with CS for the in silico and in vitro datasets, respectively. Comparison with liver stiffness measured from full sampled data and highly undersampled data (CR=4) demonstrates that the DRCS method provided the strongest correlation ( R 2 =0.95), second-lowest mean bias (-0.18 kPa, lowest for CS with -0.16 kPa), and lowest coefficient of variation (CV=3.6%). Our results demonstrate the potential of using DRCS to improve the reconstruction quality of accelerated MRE.

Unsupervised cycle-consistent network using restricted subspace field map for removing susceptibility artifacts in EPI

PURPOSE

To design an unsupervised deep neural model for correcting susceptibility artifacts in single-shot Echo Planar Imaging (EPI) and evaluate the model for preclinical and clinical applications.

METHODS

This work proposes an unsupervised cycle-consistent model based on the restricted subspace field map to take advantage of both the deep learning (DL) and the reverse polarity-gradient (RPG) method for single-shot EPI. The proposed model consists of three main components: (1) DLRPG neural network (DLRPG-net) to obtain field maps based on a pair of images acquired with reversed phase encoding; (2) spin physical model-based modules to obtain the corrected undistorted images based on the learned field map; and (3) cycle-consistency loss between the input images and back-calculated images from each cycle is explored for network training. In addition, the field maps generated by DLRPG-net belong to a restricted subspace, which is a span of predefined cubic splines to ensure the smoothness of the field maps and avoid blurring in the corrected images. This new method is trained and validated on both preclinical and clinical datasets for diffusion MRI.

RESULTS

The proposed network could effectively generate smooth field maps and correct susceptibility artifacts in single-shot EPI. Simulated and in vivo preclinical/clinical experiments demonstrated that our method outperforms the state-of-the-art susceptibility artifact correction methods. Furthermore, the ablation experiments of the cycle-consistent network and the restricted subspace in generating field maps did show the advantages of DLRPG-net.

CONCLUSION

The proposed method (DLRPG-net) can effectively correct susceptibility artifacts for preclinical and clinical single-shot EPI sequences.

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.

Motion-corrected 4D-Flow MRI for neurovascular applications

Neurovascular 4D-Flow MRI has emerged as a powerful tool for comprehensive cerebrovascular hemodynamic characterization. Clinical studies in at risk populations such as aging adults indicate hemodynamic markers can be confounded by motion-induced bias. This study develops and characterizes a high fidelity 3D self-navigation approach for retrospective rigid motion correction of neurovascular 4D-Flow data. A 3D radial trajectory with pseudorandom ordering was combined with a multi-resolution low rank regularization approach to enable high spatiotemporal resolution self-navigators from extremely undersampled data. Phantom and volunteer experiments were performed at 3.0T to evaluate the ability to correct for different amounts of induced motions. In addition, the approach was applied to clinical-research exams from ongoing aging studies to characterize performance in the clinical setting. Simulations, phantom and volunteer experiments with motion correction produced images with increased vessel conspicuity, reduced image blurring, and decreased variability in quantitative measures. Clinical exams revealed significant changes in hemodynamic parameters including blood flow rates, flow pulsatility index, and lumen areas after motion correction in probed cerebral arteries (Flow: P<0.001 Lt ICA, P=0.002 Rt ICA, P=0.004 Lt MCA, P=0.004 Rt MCA; Area: P<0.001 Lt ICA, P<0.001 Rt ICA, P=0.004 Lt MCA, P=0.004 Rt MCA; flow pulsatility index: P=0.042 Rt ICA, P=0.002 Lt MCA). Motion induced bias can lead to significant overestimation of hemodynamic markers in cerebral arteries. The proposed method reduces measurement bias from rigid motion in neurovascular 4D-Flow MRI in challenging populations such as aging adults.

Data‐driven motion‐corrected brain MRI incorporating pose‐dependent B 0 fields

Purpose

To develop a fully data‐driven retrospective intrascan motion‐correction framework for volumetric brain MRI at ultrahigh field (7 Tesla) that includes modeling of pose‐dependent changes in polarizing magnetic (B₀) fields.

Theory and Methods

Tissue susceptibility induces spatially varying B₀ distributions in the head, which change with pose. A physics‐inspired B₀ model has been deployed to model the B₀ variations in the head and was validated in vivo. This model is integrated into a forward parallel imaging model for imaging in the presence of motion. Our proposal minimizes the number of added parameters, enabling the developed framework to estimate dynamic B₀ variations from appropriately acquired data without requiring navigators. The effect on data‐driven motion correction is validated in simulations and in vivo.

Results

The applicability of the physics‐inspired B₀ model was confirmed in vivo. Simulations show the need to include the pose‐dependent B₀ fields in the reconstruction to improve motion‐correction performance and the feasibility of estimating B₀ evolution from the acquired data. The proposed motion and B₀ correction showed improved image quality for strongly corrupted data at 7 Tesla in simulations and in vivo.

Conclusion

We have developed a motion‐correction framework that accounts for and estimates pose‐dependent B₀ fields. The method improves current state‐of‐the‐art data‐driven motion‐correction techniques when B₀ dependencies cannot be neglected. The use of a compact physics‐inspired B₀ model together with leveraging the parallel imaging encoding redundancy and previously proposed optimized sampling patterns enables a purely data‐driven approach.

MRI Reconstruction with Separate Magnitude and Phase Priors Based on Dual-Tree Complex Wavelet Transform

The methods of compressed sensing magnetic resonance imaging (CS-MRI) can be divided into two categories roughly based on the number of target variables. One group devotes to estimating the complex-valued MRI image. And the other calculates the magnitude and phase parts of the complex-valued MRI image, respectively, by enforcing separate penalties on them. We propose a new CS-based method based on dual-tree complex wavelet (DT CWT) sparsity, which is under the frame of the second class of CS-MRI. Owing to the separate regularization frame, this method reduces the impact of the phase jumps (that means the jumps or discontinuities of phase values) on magnitude reconstruction. Moreover, by virtue of the excellent features of DT CWT, such as nonoscillating envelope of coefficients and multidirectional selectivity, the proposed method is capable of capturing more details in the magnitude and phase images. The experimental results show that the proposed method recovers the image contour and edges information well and can eliminate the artifacts in magnitude results caused by phase jumps.

Velocity vector reconstruction for real-time phase-contrast MRI with radial Maxwell correction

PURPOSE

To develop an auto-calibrated image reconstruction for highly accelerated multi-directional phase-contrast (PC) MRI that compensates for (1) reconstruction instabilities occurring for phase differences near ± π and (2) phase errors by concomitant magnetic fields that differ for individual radial spokes.

THEORY AND METHODS

A model-based image reconstruction for real-time PC MRI based on nonlinear inversion is extended to multi-directional flow by exploiting multiple flow-encodings for the estimation of velocity vectors. An initial smoothing constraint during iterative optimization is introduced to resolve the ambiguity of the solution space by penalizing phase wraps. Maxwell terms are considered as part of the signal model on a line-by-line basis to address phase errors by concomitant magnetic fields. The reconstruction methods are evaluated using simulated data and cross-sectional imaging of a rotating-disc, as well as in vivo for the aortic arch and cervical spinal canal at 3T.

RESULTS

Real-time three-directional velocity mapping in the aortic arch is achieved at 1.8 × 1.8 × 6 mm³ spatial and 60 ms temporal resolution. Artificial phase wraps are avoided in all cases using the smoothness constraint. Inter-spoke differences of concomitant magnetic fields are effectively compensated for by the model-based image reconstruction with integrated radial Maxwell correction.

CONCLUSION

Velocity vector reconstructions based on nonlinear inversion allow for high degrees of radial data undersampling paving the way for multi-directional PC MRI in real time. Whether a spoke-wise treatment of Maxwell terms is required or a computationally cheaper frame-wise approach depends on the individual application.

Accelerating quantitative susceptibility and R2* mapping using incoherent undersampling and deep neural network reconstruction

Quantitative susceptibility mapping (QSM) and R2* mapping are MRI post-processing methods that quantify tissue magnetic susceptibility and transverse relaxation rate distributions. However, QSM and R2* acquisitions are relatively slow, even with parallel imaging. Incoherent undersampling and compressed sensing reconstruction techniques have been used to accelerate traditional magnitude-based MRI acquisitions; however, most do not recover the full phase signal, as required by QSM, due to its non-convex nature. In this study, a learning-based Deep Complex Residual Network (DCRNet) is proposed to recover both the magnitude and phase images from incoherently undersampled data, enabling high acceleration of QSM and R2* acquisition. Magnitude, phase, R2*, and QSM results from DCRNet were compared with two iterative and one deep learning methods on retrospectively undersampled acquisitions from six healthy volunteers, one intracranial hemorrhage and one multiple sclerosis patients, as well as one prospectively undersampled healthy subject using a 7T scanner. Peak signal to noise ratio (PSNR), structural similarity (SSIM), root-mean-squared error (RMSE), and region-of-interest susceptibility and R2* measurements are reported for numerical comparisons. The proposed DCRNet method substantially reduced artifacts and blurring compared to the other methods and resulted in the highest PSNR, SSIM, and RMSE on the magnitude, R2*, local field, and susceptibility maps. Compared to two iterative and one deep learning methods, the DCRNet method demonstrated a 3.2% to 9.1% accuracy improvement in deep grey matter susceptibility when accelerated by a factor of four. The DCRNet also dramatically shortened the reconstruction time of single 2D brain images from 36-140 seconds using conventional approaches to only 15-70 milliseconds.

Automatic determination of the regularization weighting for wavelet-based compressed sensing MRI reconstructions

PURPOSE

To present a method that automatically, rapidly, and in a noniterative manner determines the regularization weighting for wavelet-based compressed sensing reconstructions. This method determines level-specific regularization weighting factors from the wavelet transform of the image obtained from zero-filling in k-space.

METHODS

We compare reconstruction results obtained by our method, λ auto , to the ones obtained by the L-curve, λ Lcurve , and the minimum NMSE, λ NMSE . The comparisons are done using in vivo data; then, simulations are used to analyze the impact of undersampling and noise. We use NMSE, Pearson's correlation coefficient, high-frequency error norm, and structural similarity as reconstruction quality indices.

RESULTS

Our method, λ auto , provides improved reconstructed image quality to that obtained by λ Lcurve regardless of undersampling or SNR and comparable quality to λ NMSE at high SNR. The method determines the regularization weighting prospectively with negligible computational time.

CONCLUSION

Our main finding is an automatic, fast, noniterative, and robust procedure to determine the regularization weighting. The impact of this method is to enable prospective and tuning-free wavelet-based compressed sensing reconstructions.

PALMNUT: An Enhanced Proximal Alternating Linearized Minimization Algorithm with Application to Separate Regularization of Magnitude and Phase

We introduce a new algorithm for complex image reconstruction with separate regularization of the image magnitude and phase. This optimization problem is interesting in many different image reconstruction contexts, although is nonconvex and can be difficult to solve. In this work, we first describe a novel implementation of the previous proximal alternating linearized minimization (PALM) algorithm to solve this optimization problem. We then make enhancements to PALM, leading to a new algorithm named PALMNUT that combines the PALM together with Nesterov's momentum and a novel approach that relies on uncoupled coordinatewise step sizes derived from coordinatewise Lipschitz-like bounds. Theoretically, we establish that a version of PALMNUT (without Nesterov's momentum) monotonically decreases the objective function, guaranteeing convergence of the cost function value. Empirical results obtained in the context of magnetic resonance imaging demonstrate that PALMNUT has computational advantages over common existing approaches like alternating minimization. Although our focus is on the application to separate magnitude and phase regularization, we expect that the same approach may also be useful in other nonconvex optimization problems with similar objective function structure.

Variable flip angle echo planar time-resolved imaging (vFA-EPTI) for fast high-resolution gradient echo myelin water imaging

Myelin water imaging techniques based on multi-compartment relaxometry have been developed as an important tool to measure myelin concentration in vivo, but are limited by the long scan time of multi-contrast multi-echo acquisition. In this work, a fast imaging technique, termed variable flip angle Echo Planar Time-Resolved Imaging (vFA-EPTI), is developed to acquire multi-echo and multi-flip-angle gradient-echo data with significantly reduced acquisition time, providing rich information for multi-compartment analysis of gradient-echo myelin water imaging (GRE-MWI). The proposed vFA-EPTI method achieved 26 folds acceleration with good accuracy by utilizing an efficient continuous readout, optimized spatiotemporal encoding across echoes and flip angles, as well as a joint subspace reconstruction. An approach to estimate off-resonance field changes between different flip-angle acquisitions was also developed to ensure high-quality joint reconstruction across flip angles. The accuracy of myelin water fraction (MWF) estimate under high acceleration was first validated by a retrospective undersampling experiment using a lengthy fully-sampled data as reference. Prospective experiments were then performed where whole-brain MWF and multi-compartment quantitative maps were obtained in 5 min at 1.5 mm isotropic resolution and 24 min at 1 mm isotropic resolution at 3T. Additionally, ultra-high resolution data at 600 μm isotropic resolution were acquired at 7T, which show detailed structures within the cortex such as the line of Gennari, demonstrating the ability of the proposed method for submillimeter GRE-MWI that can be used to study cortical myeloarchitecture in vivo.

A Modified Phase Cycling Method for Complex-Valued MRI Reconstruction

The phase cycling method is a state-of-the-art method to reconstruct complex-valued MR image. However, when it follows practical two-dimensional (2D) subsampling Cartesian acquisition which is only enforcing random sampling in the phase-encoding direction, a number of artifacts in magnitude appear. A modified approach is proposed to remove these artifacts under practical MRI subsampling, by adding one-dimensional total variation (TV) regularization into the phase cycling method to "pre-process" the magnitude component before its update. Furthermore, an operation used in SFISTA is employed to update the magnitude and phase images for better solutions. The results of the experiments show the ability of the proposed method to eliminate the ring artifacts and improve the magnitude reconstruction.

Echo planar time-resolved imaging with subspace reconstruction and optimized spatiotemporal encoding

PURPOSE

To develop new encoding and reconstruction techniques for fast multi-contrast/quantitative imaging.

METHODS

The recently proposed Echo Planar Time-resolved Imaging (EPTI) technique can achieve fast distortion- and blurring-free multi-contrast/quantitative imaging. In this work, a subspace reconstruction framework is developed to improve the reconstruction accuracy of EPTI at high encoding accelerations. The number of unknowns in the reconstruction is significantly reduced by modeling the temporal signal evolutions using low-rank subspace. As part of the proposed reconstruction approach, a B₀ -update algorithm and a shot-to-shot B₀ variation correction method are developed to enable the reconstruction of high-resolution tissue phase images and to mitigate artifacts from shot-to-shot phase variations. Moreover, the EPTI concept is extended to 3D k-space for 3D GE-EPTI, where a new "temporal-variant" of CAIPI encoding is proposed to further improve performance.

RESULTS

The effectiveness of the proposed subspace reconstruction was demonstrated first in 2D GESE EPTI, where the reconstruction achieved higher accuracy when compared to conventional B₀ -informed GRAPPA. For 3D GE-EPTI, a retrospective undersampling experiment demonstrates that the new temporal-variant CAIPI encoding can achieve up to 72× acceleration with close to 2× reduction in reconstruction error when compared to conventional spatiotemporal-CAIPI encoding. In a prospective undersampling experiment, high-quality whole-brain T 2 ∗ and tissue phase maps at 1 mm isotropic resolution were acquired in 52 seconds at 3T using 3D GE-EPTI with temporal-variant CAIPI encoding.

CONCLUSION

The proposed subspace reconstruction and optimized temporal-variant CAIPI encoding can further improve the performance of EPTI for fast quantitative mapping.

Constraints in estimating the proton density fat fraction

The study evaluates four physically motivated constraints in the estimation of the proton density fat fraction (PDFF). Least squares approaches were developed for constraining the parameters in PDFF quantification based on the physics of magnetic resonance imaging. These were smooth fieldmap, smooth initial phase, nonnegative proton density and moderate R₂^∗ values. The constraints were evaluated in terms of their influence on the bias and standard deviation of the estimated parameters using numerical simulations and in vivo data acquired at 0.35 T. Results show that unconstrained least squares estimation is noisy and biased and that constraints can be effective at reducing both the standard deviation and bias.

Multi-shot diffusion-weighted MRI reconstruction with magnitude-based spatial-angular locally low-rank regularization (SPA-LLR)

PURPOSE

To resolve the motion-induced phase variations in multi-shot multi-direction diffusion-weighted imaging (DWI) by applying regularization to magnitude images.

THEORY AND METHODS

A nonlinear model was developed to estimate phase and magnitude images separately. A locally low-rank regularization (LLR) term was applied to the magnitude images from all diffusion-encoding directions to exploit the spatial and angular correlation. In vivo experiments with different resolutions and b-values were performed to validate the proposed method.

RESULTS

The proposed method significantly reduces the noise level compared to the conventional reconstruction method and achieves submillimeter (0.8mm and 0.9mm isotropic resolutions) DWI with a b-value of 1,000  s / mm 2 and 1-mm isotropic DWI with a b-value of 2,000  s / mm 2 without modification of the sequence.

CONCLUSIONS

A joint reconstruction method with spatial-angular LLR regularization on magnitude images substantially improves multi-direction DWI reconstruction, simultaneously removes motion-induced phase artifacts, and denoises images.

Highly accelerated multishot echo planar imaging through synergistic machine learning and joint reconstruction

PURPOSE

To introduce a combined machine learning (ML)- and physics-based image reconstruction framework that enables navigator-free, highly accelerated multishot echo planar imaging (msEPI) and demonstrate its application in high-resolution structural and diffusion imaging.

METHODS

Single-shot EPI is an efficient encoding technique, but does not lend itself well to high-resolution imaging because of severe distortion artifacts and blurring. Although msEPI can mitigate these artifacts, high-quality msEPI has been elusive because of phase mismatch arising from shot-to-shot variations which preclude the combination of the multiple-shot data into a single image. We utilize deep learning to obtain an interim image with minimal artifacts, which permits estimation of image phase variations attributed to shot-to-shot changes. These variations are then included in a joint virtual coil sensitivity encoding (JVC-SENSE) reconstruction to utilize data from all shots and improve upon the ML solution.

RESULTS

Our combined ML + physics approach enabled Rinplane × multiband (MB) = 8- × 2-fold acceleration using 2 EPI shots for multiecho imaging, so that whole-brain T2 and T2 * parameter maps could be derived from an 8.3-second acquisition at 1 × 1 × 3-mm3 resolution. This has also allowed high-resolution diffusion imaging with high geometrical fidelity using 5 shots at Rinplane × MB = 9- × 2-fold acceleration. To make these possible, we extended the state-of-the-art MUSSELS reconstruction technique to simultaneous multislice encoding and used it as an input to our ML network.

CONCLUSION

Combination of ML and JVC-SENSE enabled navigator-free msEPI at higher accelerations than previously possible while using fewer shots, with reduced vulnerability to poor generalizability and poor acceptance of end-to-end ML approaches.

Accelerating cardiac diffusion tensor imaging combining local low-rank and 3D TV constraint

OBJECTIVE

Diffusion tensor magnetic resonance imaging (DT-MRI, or DTI) is a promising technique for invasively probing biological tissue structures. However, DTI is known to suffer from much longer acquisition time with respect to conventional MRI and the problem is worsened when dealing with in vivo acquisitions. Therefore, faster DTI for both ex vivo and in vivo scans is highly desired.

MATERIALS AND METHODS

This paper proposes a new compressed sensing (CS) reconstruction method that employs local low-rank (LLR) model and three-dimensional (3D) total variation (TV) constraint to reconstruct cardiac diffusion-weighted (DW) images from highly undersampled k-space data. The LLR model takes the set of DW images corresponding to different diffusion gradient directions as a 3D image volume and decomposes the latter into overlapping 3D blocks. Then, the 3D blocks are stacked as two-dimensional (2D) matrix. Finally, low-rank property is applied to each block matrix and the 3D TV constraint to the 3D image volume. The underlying constrained optimization problem is finally solved using the first-order fast method. The proposed method is evaluated on real ex vivo cardiac DTI data as a prerequisite to in vivo cardiac DTI applications.

RESULTS

The results on real human ex vivo cardiac DTI images demonstrate that the proposed method exhibits lower reconstruction errors for DTI indices, including fractional anisotropy (FA), mean diffusivities (MD), transverse angle (TA), and helix angle (HA), compared to existing CS-based DTI image reconstruction techniques.

CONCLUSION

The proposed method provides better reconstruction quality and more accurate DTI indices in comparison with the state-of-the-art CS-based DW image reconstruction methods.

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.

Dynamic water/fat separation and B 0 inhomogeneity mapping-joint estimation using undersampled triple-echo multi-spoke radial FLASH

PURPOSE

To achieve dynamic water/fat separation and B 0 field inhomogeneity mapping via model-based reconstructions of undersampled triple-echo multi-spoke radial FLASH acquisitions.

METHODS

This work introduces an undersampled triple-echo multi-spoke radial FLASH sequence, which uses (i) complementary radial spokes per echo train for faster spatial encoding, (ii) asymmetric echoes for flexible and nonuniform echo spacing, and (iii) a golden angle increment across frames for optimal k-space coverage. Joint estimation of water, fat, B 0 inhomogeneity, and coil sensitivity maps from undersampled triple-echo data poses a nonlinear and non-convex inverse problem which is solved by a model-based reconstruction with suitable regularization. The developed methods are validated using phantom experiments with different degrees of undersampling. Real-time MRI studies of the knee, liver, and heart are conducted without prospective gating or retrospective data sorting at temporal resolutions of 70, 158, and 40 ms, respectively.

RESULTS

Up to 18-fold undersampling is achieved in this work. Even in the presence of rapid physiological motion, large B 0 field inhomogeneities, and phase wrapping, the model-based reconstruction yields reliably separated water/fat maps in conjunction with spatially smooth inhomogeneity maps.

CONCLUSIONS

The combination of a triple-echo acquisition and joint reconstruction technique provides a practical solution to time-resolved and motion robust water/fat separation at high spatial and temporal resolution.

On the limitations of partial Fourier acquisition in phase-contrast MRI of turbulent kinetic energy

PURPOSE

To investigate limitations of partial Fourier acquisition in phase-contrast MRI of turbulent kinetic energy (TKE).

METHODS

To assess the validity of partial Fourier reconstruction of TKE and phase images, computational fluid dynamics data of mean and turbulent velocities in a stenotic U-bend phantom was used. Partial Fourier acquisition with 75% k-space coverage was simulated and TKE data were reconstructed using zero-filling, homodyne reconstruction, and the method of projections onto convex sets (POCS). Results were compared to data from fully sampled k-space and 75% symmetric sampling. In addition, compressed sensing (CS) reconstruction was compared for a standard variable density sampling pattern and a variable density sampling pattern combined with 75% partial Fourier. For illustration purposes, in vivo examples of velocity magnitude and TKE maps of aortic flow reconstructed with the different methods are provided.

RESULTS

In accordance with theory, partial Fourier reconstruction of TKE maps from phase-contrast data results in artifacts relative to fully sampled data. It is demonstrated that neither homodyne reconstruction nor POCS can improve reconstruction of TKE data with respect to zero-filling reconstruction when compared to ground-truth (RMS error: 4.70%, 4.34%, and 2.45% for homodyne, POCS, and zero-filling reconstruction of in vivo data, respectively). CS reconstruction from data acquired with partial Fourier did not recover the resolution loss incurred by partial Fourier sampling.

CONCLUSION

Partial Fourier reconstruction of TKE maps from phase-contrast data does not yield a benefit over zero-filling reconstruction. In consequence, symmetric sampling is preferred over partial Fourier acquisition for a given number of phase-encodes in phase-contrast MRI.