SMASH, SENSE, PILS, GRAPPA: how to choose the optimal method

High-Frequency Space Diffusion Model for Accelerated MRI

Diffusion models with continuous stochastic differential equations (SDEs) have shown superior performances in image generation. It can serve as a deep generative prior to solving the inverse problem in magnetic resonance (MR) reconstruction. However, low-frequency regions of k -space data are typically fully sampled in fast MR imaging, while existing diffusion models are performed throughout the entire image or k -space, inevitably introducing uncertainty in the reconstruction of low-frequency regions. Additionally, existing diffusion models often demand substantial iterations to converge, resulting in time-consuming reconstructions. To address these challenges, we propose a novel SDE tailored specifically for MR reconstruction with the diffusion process in high-frequency space (referred to as HFS-SDE). This approach ensures determinism in the fully sampled low-frequency regions and accelerates the sampling procedure of reverse diffusion. Experiments conducted on the publicly available fastMRI dataset demonstrate that the proposed HFS-SDE method outperforms traditional parallel imaging methods, supervised deep learning, and existing diffusion models in terms of reconstruction accuracy and stability. The fast convergence properties are also confirmed through theoretical and experimental validation. Our code and weights are available at https://github.com/Aboriginer/HFS-SDE.

Ferumoxytol-Enhanced Cardiac Cine MRI Reconstruction Using a Variable-Splitting Spatiotemporal Network

BACKGROUND

Balanced steady-state free precession (bSSFP) imaging is commonly used in cardiac cine MRI but prone to image artifacts. Ferumoxytol-enhanced (FE) gradient echo (GRE) has been proposed as an alternative. Utilizing the abundance of bSSFP images to develop a computationally efficient network that is applicable to FE GRE cine would benefit future network development.

PURPOSE

To develop a variable-splitting spatiotemporal network (VSNet) for image reconstruction, trained on bSSFP cine images and applicable to FE GRE cine images.

STUDY TYPE

Retrospective and prospective.

SUBJECTS

41 patients (26 female, 53 ± 19 y/o) for network training, 31 patients (19 female, 49 ± 17 y/o) and 5 healthy subjects (5 female, 30 ± 7 y/o) for testing.

FIELD STRENGTH/SEQUENCE

1.5T and 3T, bSSFP and GRE.

ASSESSMENT

VSNet was compared to VSNet with total variation loss, compressed sensing and low rank methods for 14× accelerated data. The GRAPPA×2/×3 images served as the reference. Peak signal-to-noise-ratio (PSNR), structural similarity index (SSIM), left ventricular (LV) and right ventricular (RV) end-diastolic volume (EDV), end-systolic volume (ESV), and ejection fraction (EF) were measured. Qualitative image ranking and scoring were independently performed by three readers. Latent scores were calculated based on scores of each method relative to the reference.

STATISTICS

Linear mixed-effects regression, Tukey method, Fleiss' Kappa, Bland-Altman analysis, and Bayesian categorical cumulative probit model. A P-value <0.05 was considered statistically significant.

RESULTS

VSNet achieved significantly higher PSNR (32.7 ± 0.2), SSIM (0.880 ± 0.004), rank (2.14 ± 0.06), and latent scores (-1.72 ± 0.22) compared to other methods (rank >2.90, latent score < -2.63). Fleiss' Kappa was 0.52 for scoring and 0.61 for ranking. VSNet showed no significantly different LV and RV ESV (P = 0.938) and EF (P = 0.143) measurements, but statistically significant different (2.62 mL) EDV measurements compared to the reference.

CONCLUSION

VSNet produced the highest image quality and the most accurate functional measurements for FE GRE cine images among the tested 14× accelerated reconstruction methods.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 1.

Alternating Look-Locker for quantitative T1 , T1ρ and B1 3D MRI mapping

PURPOSE

To develop a new MRI method, entitled alternating Look-Locker (aLL), for quantitativeT 1 $$ {T}_1 $$ ,T 1 ρ $$ {T}_{1\uprho} $$ , andB 1 $$ {B}_1 $$ 3D mapping.

METHODS

A Look-Locker scheme that alternates magnetization from +Z and -Z axes of the laboratory frame is utilized in combination with a 3D Multi-Band Sweep Imaging with Fourier Transformation (MB-SWIFT) readout. The analytical solution describing the spin evolution during aLL, as well as the correction required for segmented acquisition were derived. The simultaneousB 1 $$ {B}_1 $$ andT 1 $$ {T}_1 $$ mapping are demonstrated on an agar/saline phantom and on an in-vivo rat head.T 1 ρ $$ {T}_{1\uprho} $$ relaxation was achieved by cyclically applying magnetization preparation (MP) modules consisting of two adiabatic pulses.T 1 ρ $$ {T}_{1\uprho} $$ values in the rat brain in-vivo and in a gadobenate dimeglumine (Gd-DTPA) phantom were compared to those obtained with a previously introduced steady-state (SS) method.

RESULTS

The accuracy and precision of the analytical solution was tested by Bloch simulations. With the application of MP modules, the aLL method provides simultaneousT 1 $$ {T}_1 $$ andT 1 ρ $$ {T}_{1\uprho} $$ maps. Conversely, without it, the method can be used for simultaneousT 1 $$ {T}_1 $$ andB 1 $$ {B}_1 $$ mapping.T 1 ρ $$ {T}_{1\uprho} $$ values were similar with both aLL and SS techniques. However, the aLL method resulted in more robust quantitative mapping compared to the SS method. Unlike the SS method, the aLL method does not require additional scans for generatingT 1 $$ {T}_1 $$ maps.

CONCLUSION

The proposed method offers a new flexible tool for quantitative mapping ofT 1 $$ {T}_1 $$ ,T 1 ρ $$ {T}_{1\uprho} $$ , andB 1 $$ {B}_1 $$ . The aLL method can also be used with readout schemes different from MB-SWIFT.

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

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

Compressed SVD-based L + S model to reconstruct undersampled dynamic MRI data using parallel architecture

BACKGROUND

Magnetic Resonance Imaging (MRI) is a highly demanded medical imaging system due to high resolution, large volumetric coverage, and ability to capture the dynamic and functional information of body organs e.g. cardiac MRI is employed to assess cardiac structure and evaluate blood flow dynamics through the cardiac valves. Long scan time is the main drawback of MRI, which makes it difficult for the patients to remain still during the scanning process.

OBJECTIVE

By collecting fewer measurements, MRI scan time can be shortened, but this undersampling causes aliasing artifacts in the reconstructed images. Advanced image reconstruction algorithms have been used in literature to overcome these undersampling artifacts. These algorithms are computationally expensive and require a long time for reconstruction which makes them infeasible for real-time clinical applications e.g. cardiac MRI. However, exploiting the inherent parallelism in these algorithms can help to reduce their computation time.

METHODS

Low-rank plus sparse (L+S) matrix decomposition model is a technique used in literature to reconstruct the highly undersampled dynamic MRI (dMRI) data at the expense of long reconstruction time. In this paper, Compressed Singular Value Decomposition (cSVD) model is used in L+S decomposition model (instead of conventional SVD) to reduce the reconstruction time. The results provide improved quality of the reconstructed images. Furthermore, it has been observed that cSVD and other parts of the L+S model possess highly parallel operations; therefore, a customized GPU based parallel architecture of the modified L+S model has been presented to further reduce the reconstruction time.

RESULTS

Four cardiac MRI datasets (three different cardiac perfusion acquired from different patients and one cardiac cine data), each with different acceleration factors of 2, 6 and 8 are used for experiments in this paper. Experimental results demonstrate that using the proposed parallel architecture for the reconstruction of cardiac perfusion data provides a speed-up factor up to 19.15× (with memory latency) and 70.55× (without memory latency) in comparison to the conventional CPU reconstruction with no compromise on image quality.

CONCLUSION

The proposed method is well-suited for real-time clinical applications, offering a substantial reduction in reconstruction time.

Simultaneous multislice EPI prospective motion correction by real-time receiver phase correction and coil sensitivity map interpolation

PURPOSE

To improve the image reconstruction for prospective motion correction (PMC) of simultaneous multislice (SMS) EPI of the brain, an update of receiver phase and resampling of coil sensitivities are proposed and evaluated.

METHODS

A camera-based system was used to track head motion (3 translations and 3 rotations) and dynamically update the scan position and orientation. We derived the change in receiver phase associated with a shifted field of view (FOV) and applied it in real-time to each k-space line of the EPI readout trains. Second, for the SMS reconstruction, we adapted resampled coil sensitivity profiles reflecting the movement of slices. Single-shot gradient-echo SMS-EPI scans were performed in phantoms and human subjects for validation.

RESULTS

Brain SMS-EPI scans in the presence of motion with PMC and no phase correction for scan plane shift showed noticeable artifacts. These artifacts were visually and quantitatively attenuated when corrections were enabled. Correcting misaligned coil sensitivity maps improved the temporal SNR (tSNR) of time series by 24% (p = 0.0007) for scans with large movements (up to ˜35 mm and 30°). Correcting the receiver phase improved the tSNR of a scan with minimal head movement by 50% from 50 to 75 for a United Kingdom biobank protocol.

CONCLUSION

Reconstruction-induced motion artifacts in single-shot SMS-EPI scans acquired with PMC can be removed by dynamically adjusting the receiver phase of each line across EPI readout trains and updating coil sensitivity profiles during reconstruction. The method may be a valuable tool for SMS-EPI scans in the presence of subject motion.

Optimizing variable flip angles in magnetization-prepared gradient-echo sequences for efficient 3D-T1ρ mapping

PURPOSE

To optimize the choice of the flip angles of magnetization-prepared gradient-echo sequences for improved accuracy, precision, and speed of 3D-T1ρ mapping.

METHODS

We propose a new optimization approach for finding variable flip-angle values that improve magnetization-prepared gradient-echo sequences used for 3D-T1ρ mapping. This new approach can improve the accuracy and SNR, while reducing filtering effects. We demonstrate the concept in the three different versions of the magnetization-prepared gradient-echo sequences that are typically used for 3D-T1ρ mapping and evaluate their performance in model agarose phantoms (n = 4) and healthy volunteers (n = 5) for knee joint imaging. We also tested the optimization with sequence parameters targeting faster acquisitions.

RESULTS

Our results show that optimized variable flip angle can improve the accuracy and the precision of the sequences, seen as a reduction of the mean of normalized absolute difference from about 5%-6% to 3%-4% in model phantoms and from 15%-16% to 11%-13% in the knee joint, and improving SNR from about 12-28 to 22-32 in agarose phantoms and about 7-14 to 13-17 in healthy volunteers. The optimization can also compensate for the loss in quality caused by making the sequence faster. This results in sequence configurations that acquire more data per unit of time with SNR and mean of normalized absolute difference measurements close to its slower versions.

CONCLUSION

The optimization of the variable flip angle can be used to increase accuracy and precision, and to improve the speed of the typical imaging sequences used for quantitative 3D-T1ρ mapping of the knee joint.

Influence of half Fourier and elliptical scanning (radial scan) on magnetic resonance images

This study aimed to investigate the effect of using slice partial Fourier (SPF), phase partial Fourier (PPF), and radial scan (Elliptical scanning) methods on image quality. Changes in signal-to-noise ratio (SNR), effective slice thickness, and in-plane resolution were measured in 3D-gradient echo when SPF, PPF, and radial scan were used. Effective slice thickness increased and SNR increased when SPF was used; in-plane resolution decreased and SNR decreased when PPF was used; effective slice thickness did not change, in-plane resolution decreased, and SNR increased when the radial scan method was used. The radial scan method reduces image quality and imaging time compared to those in the SPF and PPF methods.

Review of oxygen-enhanced lung mri: Pulse sequences for image acquisition and T1 measurement

Oxygen-enhanced MR imaging (OE-MRI) is a special proton imaging technique that can be performed without modifying the scanner hardware. Many fundamental studies have been conducted following the initial reporting of this technique in 1996, illustrating the high potential for its clinical application. This review aims to summarise and analyse current pulse sequences and T1 measurement methods for OE-MRI, including fundamental theories, existing pulse sequences applied to OE-MRI acquisition and T1 mapping. Wash-in and wash-out time identify lung function and are sensitive to ventilation; thus, dynamic OE-MRI is also discussed in this review. We compare OE-MRI with the primary competitive technique, hyperpolarised gas MRI. Finally, an overview of lower-field applications of OE-MRI is highlighted, as relatively recent publications demonstrated positive results. Lower-field OE-MRI, which is lower than 1.5 T, could be an alternative modality for detecting lung diseases. This educational review is aimed at researchers who want a quick summary of the steps needed to perform pulmonary OE-MRI with a particular focus on sequence design, settings, and quantification methods.

Psilocybin desynchronizes brain networks

1The relationship between the acute effects of psychedelics and their persisting neurobiological and psychological effects is poorly understood. Here, we tracked brain changes with longitudinal precision functional mapping in healthy adults before, during, and for up to 3 weeks after oral psilocybin and methylphenidate (17 MRI visits per participant) and again 6+ months later. Psilocybin disrupted connectivity across cortical networks and subcortical structures, producing more than 3-fold greater acute changes in functional networks than methylphenidate. These changes were driven by desynchronization of brain activity across spatial scales (area, network, whole brain). Psilocybin-driven desynchronization was observed across association cortex but strongest in the default mode network (DMN), which is connected to the anterior hippocampus and thought to create our sense of self. Performing a perceptual task reduced psilocybin-induced network changes, suggesting a neurobiological basis for grounding, connecting with physical reality during psychedelic therapy. The acute brain effects of psilocybin are consistent with distortions of space-time and the self. Psilocybin induced persistent decrease in functional connectivity between the anterior hippocampus and cortex (and DMN in particular), lasting for weeks but normalizing after 6 months. Persistent suppression of hippocampal-DMN connectivity represents a candidate neuroanatomical and mechanistic correlate for psilocybin's pro-plasticity and anti-depressant effects.

Accelerating GRAPPA reconstruction using SoC design for real-time cardiac MRI

Real-time cardiac MRI is a rapidly developing area of research that has the potential to improve the diagnosis and treatment of cardiovascular diseases. However, the acquisition of high-quality real-time cardiac MR (CMR) images is challenging as it requires a high frame rate and temporal resolution. To overcome this challenge, there have been recent efforts on several approaches including hardware-based improvements and image reconstruction techniques such as compressed sensing and parallel MRI. The use of parallel MRI techniques such as GRAPPA (Generalized Autocalibrating Partial Parallel Acquisition) is a promising approach for improving the temporal resolution of MRI and expanding its applications in clinical practice. However, the GRAPPA algorithm involves a significant amount of computation, particularly for high acceleration factors and large datasets. This can result in long reconstruction times, which can limit the ability to achieve real-time imaging or high frame rates. One solution to this challenge is to use specialized hardware i.e. field-programmable gate arrays (FPGAs). In this work, a novel 32-bit floating-point FPGA-based GRAPPA accelerator is proposed with an aim to reconstruct high-quality cardiac MR images at higher frame rates, making it well suited for real-time clinical applications. The proposed FPGA-based accelerator consists of custom-designed data processing units named as dedicated computational engines (DCEs) that allow for a continuous flow of data between the calibration and synthesis stages of GRAPPA reconstruction process. This greatly increases the throughput and reduces the latency of the overall proposed system. Moreover, a high-speed memory module (DDR4-SDRAM) is integrated with the proposed architecture to store the multi-coil MR data. An on-chip quad-core ARM Cortex-A53 processor is used to manage access control information required for data transfer between the DCEs and DDR4-SDRAM. The proposed accelerator is implemented on Xilinx Zynq UltraScale + MPSoC using high-level synthesis (HLS) and hardware descriptive language (HDL) with an aim to explore the trade-offs between the reconstruction time, resource utilization and design effort. Several experiments have been performed using in-vivo cardiac datasets i.e. 18-receiver coil and 30-receiver coil to evaluate the performance of the proposed accelerator. A comparison is performed with the contemporary CPU and GPU-based GRAPPA reconstruction methods in terms of reconstruction time, frames-per-second and reconstruction accuracy (RMSE and SNR). The results show that the proposed accelerator achieves speed-up factors up to 121× and 9× as compared to the contemporary CPU-based and GPU-based GRAPPA reconstruction methods, respectively. Moreover, it has been demonstrated that the proposed accelerator can achieve reconstruction rates of up to ∼27 frames-per-second while maintaining the visual quality of the reconstructed images.

A Deep Learning Framework for Cardiac MR Under-Sampled Image Reconstruction with a Hybrid Spatial and k-Space Loss Function

Magnetic resonance imaging (MRI) is an efficient, non-invasive diagnostic imaging tool for a variety of disorders. In modern MRI systems, the scanning procedure is time-consuming, which leads to problems with patient comfort and causes motion artifacts. Accelerated or parallel MRI has the potential to minimize patient stress as well as reduce scanning time and medical costs. In this paper, a new deep learning MR image reconstruction framework is proposed to provide more accurate reconstructed MR images when under-sampled or aliased images are generated. The proposed reconstruction model is designed based on the conditional generative adversarial networks (CGANs) where the generator network is designed in a form of an encoder–decoder U-Net network. A hybrid spatial and k-space loss function is also proposed to improve the reconstructed image quality by minimizing the L1-distance considering both spatial and frequency domains simultaneously. The proposed reconstruction framework is directly compared when CGAN and U-Net are adopted and used individually based on the proposed hybrid loss function against the conventional L1-norm. Finally, the proposed reconstruction framework with the extended loss function is evaluated and compared against the traditional SENSE reconstruction technique using the evaluation metrics of structural similarity (SSIM) and peak signal to noise ratio (PSNR). To fine-tune and evaluate the proposed methodology, the public Multi-Coil k-Space OCMR dataset for cardiovascular MR imaging is used. The proposed framework achieves a better image reconstruction quality compared to SENSE in terms of PSNR by 6.84 and 9.57 when U-Net and CGAN are used, respectively. Similarly, it demonstrates SSIM of the reconstructed MR images comparable to the one provided by the SENSE algorithm when U-Net and CGAN are used. Comparing cases where the proposed hybrid loss function is used against the cases with the simple L1-norm, the reconstruction performance can be noticed to improve by 6.84 and 9.57 for U-Net and CGAN, respectively. To conclude this, the proposed framework using CGAN provides the best reconstruction performance compared with U-Net or the conventional SENSE reconstruction techniques. The proposed framework seems to be useful for the practical reconstruction of cardiac images since it can provide better image quality in terms of SSIM and PSNR.

BigBrain-MR: a new digital phantom with anatomically-realistic magnetic resonance properties at 100-µm resolution for magnetic resonance methods development

The benefits, opportunities and growing availability of ultra-high field magnetic resonance imaging (MRI) for humans have prompted an expansion in research and development efforts towards increasingly more advanced high-resolution imaging techniques. To maximize their effectiveness, these efforts need to be supported by powerful computational simulation platforms that can adequately reproduce the biophysical characteristics of MRI, with high spatial resolution. In this work, we have sought to address this need by developing a novel digital phantom with realistic anatomical detail up to 100-µm resolution, including multiple MRI properties that affect image generation. This phantom, termed BigBrain-MR, was generated from the publicly available BigBrain histological dataset and lower-resolution in-vivo 7T-MRI data, using a newly-developed image processing framework that allows mapping the general properties of the latter into the fine anatomical scale of the former. Overall, the mapping framework was found to be effective and robust, yielding a diverse range of realistic "in-vivo-like" MRI contrasts and maps at 100-µm resolution. BigBrain-MR was then tested in three imaging applications (motion effects and interpolation, super-resolution imaging, and parallel imaging reconstruction) to investigate its properties, value and validity as a simulation platform. The results consistently showed that BigBrain-MR can closely approximate the behavior of real in-vivo data, more realistically and with more extensive features than a more classic option such as the Shepp-Logan phantom. Its flexibility in simulating different contrast mechanisms and artifacts may also prove valuable for educational applications. BigBrain-MR is therefore deemed a favorable choice to support methodological development and demonstration in brain MRI, and has been made freely available to the community.

Joint K-space and Image-space Parallel Imaging (KIPI) for accelerated chemical exchange saturation transfer acquisition

PURPOSE

To develop an auto-calibrated technique by joint K-space and Image-space Parallel Imaging (KIPI) for accelerated CEST acquisition.

THEORY AND METHODS

The KIPI method selects a calibration frame with a low acceleration factor (AF) and auto-calibration signals (ACS) acquired, from which the coil sensitivity profiles and artifact correction maps are calculated after restoring the k-space by GRAPPA. Then the other frames with high AF and without ACS can be reconstructed by SENSE and artifact suppression. The signal leakage due to the T₂ -decay filtering in k-space compromises the SENSE reconstruction, which can be corrected by the artifact suppression algorithm of KIPI. The 2D and 3D imaging experiments were done on the phantom, healthy volunteer, and brain tumor patient with a 3T scanner.

RESULTS

The proposed KIPI method was evaluated by retrospectively undersampled data with variable AFs and compared against existing parallel imaging methods (SENSE/auto, GRAPPA, and ESPIRiT). KIPI enabled CEST frames with random AFs to achieve similar image quality, eliminated the strong aliasing artifacts, and generated significantly smaller errors than the other methods (p < 0.01). The KIPI method permitted an AF up to 12-fold in both phase-encoding and slice-encoding directions for 3D CEST source images, achieving an overall 8.2-fold speedup in scan time.

CONCLUSION

KIPI is a novel auto-calibrated parallel imaging method that enables variable AFs for different CEST frames, achieves a significant reduction in scan time, and does not compromise the accuracy of CEST maps.

Multi-channel GAN-based calibration-free diffusion-weighted liver imaging with simultaneous coil sensitivity estimation and reconstruction

Introduction

Diffusion-weighted imaging (DWI) with parallel reconstruction may suffer from a mismatch between the coil calibration scan and imaging scan due to motions, especially for abdominal imaging.

Methods

This study aimed to construct an iterative multichannel generative adversarial network (iMCGAN)-based framework for simultaneous sensitivity map estimation and calibration-free image reconstruction. The study included 106 healthy volunteers and 10 patients with tumors.

Results

The performance of iMCGAN was evaluated in healthy participants and patients and compared with the SAKE, ALOHA-net, and DeepcomplexMRI reconstructions. The peak signal-to-noise ratio (PSNR), structural similarity index measure (SSIM), root mean squared error (RMSE), and histograms of apparent diffusion coefficient (ADC) maps were calculated for assessing image qualities. The proposed iMCGAN outperformed the other methods in terms of the PSNR (iMCGAN: 41.82 ± 2.14; SAKE: 17.38 ± 1.78; ALOHA-net: 20.43 ± 2.11 and DeepcomplexMRI: 39.78 ± 2.78) for b = 800 DWI with an acceleration factor of 4. Besides, the ghosting artifacts in the SENSE due to the mismatch between the DW image and the sensitivity maps were avoided using the iMCGAN model.

Discussion

The current model iteratively refined the sensitivity maps and the reconstructed images without additional acquisitions. Thus, the quality of the reconstructed image was improved, and the aliasing artifact was alleviated when motions occurred during the imaging procedure.

A 2D-GRAPPA Algorithm with a Boomerang Kernel for 3D MRI Data Accelerated along Two Phase-Encoding Directions

For the reconstruction of 3D MRI data that are accelerated along the two phase-encoding directions, the 2D-generalized autocalibrating partially parallel acquisitions (GRAPPA) algorithm can be used to estimate the missing data in the k-space. We propose a new boomerang-shaped kernel based on theoretic and systemic analyses of the shape and dimensions of the kernel. The reconstruction efficiency of the 2D-GRAPPA algorithm with the proposed boomerang-shaped kernel (i.e., boomerang kernel (BK)-2D-GRAPPA) was compared with other 2D-GRAPPA algorithms that utilize different types of kernels (i.e., EX-2D-GRAPPA and SK-2D-GRAPPA) based on computer simulation, phantom and in vivo experiments. The proposed method was validated for different sets of ACS lines with acceleration factors from four to eight and various sizes of the kernels. A quantitative analysis was also performed by comparing the normalized root mean squared error (nRMSE) in the images and the undersampled edges. Computer simulation, in vivo and phantom experiments, and the quantitative analysis, showed that the proposed method could reduce aliasing artifacts without reducing the SNRs of the reconstructed images.

[Applications of Artificial Intelligence in MR Image Acquisition and Reconstruction]

Recently, artificial intelligence (AI) technology has shown potential clinical utility in a wide range of MRI fields. In particular, AI models for improving the efficiency of the image acquisition process and the quality of reconstructed images are being actively developed by the MR research community. AI is expected to further reduce acquisition times in various MRI protocols used in clinical practice when compared to current parallel imaging techniques. Additionally, AI can help with tasks such as planning, parameter optimization, artifact reduction, and quality assessment. Furthermore, AI is being actively applied to automate MR image analysis such as image registration, segmentation, and object detection. For this reason, it is important to consider the effects of protocols or devices in MR image analysis. In this review article, we briefly introduced issues related to AI application of MR image acquisition and reconstruction.

Utility of Wavelet Denoising with Geometry Factor Weighting for Gadoxetic Acid-enhanced Hepatobiliary-phase MR Imaging

PURPOSE

The wavelet denoising with geometry factor weighting (g-denoising) method can reduce the image noise by adapting to spatially varying noise levels induced by parallel imaging. The aim of this study was to investigate the clinical applicability of g-denoising on hepatobiliary-phase (HBP) images with gadoxetic acid.

METHODS

We subjected 53 patients suspected of harboring hepatic neoplastic lesions to gadoxetic acid-enhanced HBP imaging with and without g-denoising (g⁺HBP and g^-HBP). The matrix size was reduced for g⁺HBP images to avoid prolonging the scanning time. Two radiologists calculated the SNR, the portal vein-, and paraspinal muscle contrast-to-noise ratio (CNR) relative to the hepatic parenchyma (liver-to-portal vein- and liver-to-muscle CNR). Two other radiologists independently graded the sharpness of the liver edge, the visibility of intrahepatic vessels, the image noise, the homogeneity of liver parenchyma, and the overall image quality using a 5-point scale. Differences between g^-HBP and g⁺HBP images were determined with the two-sided Wilcoxon signed-rank test.

RESULTS

The liver-to-portal- and liver-to-muscle CNR and the SNR were significantly higher on g⁺HBP- than g^-HBP images (P < 0.01), as was the qualitative score for the image noise, homogeneity of liver parenchyma, and overall image quality (P < 0.01). Although there were no significant differences in the scores for the sharpness of the liver edge or the score assigned for the visibility of intrahepatic vessels (P = 0.05, 0.43), with g⁺HBP the score was lower in three patients for the sharpness of the liver edge and in six patients for the visibility of intrahepatic vessels.

CONCLUSION

At gadoxetic acid-enhanced HBP imaging, g-denoising yielded a better image quality than conventional HBP imaging although the anatomic details may be degraded.

Multi-level pooling encoder-decoder convolution neural network for MRI reconstruction

MRI reconstruction is one of the critical processes of MRI machines, along with the acquisition. Due to a slow processing time of signal acquiring, parallel imaging and reconstruction techniques are applied for acceleration. To accelerate the acquisition process, fewer raw data are sampled simultaneously with all RF coils acquisition. Then, the reconstruction uses under-sampled data from all RF coils to restore the final MR image that resembles the fully sampled MR image. These processes have been a traditional procedure inside the MRI system since the invention of the multi-coils MRI machine. This paper proposes the deep learning technique with a lightweight network. The deep neural network is capable of generating the high-quality reconstructed MR image with a high peak signal-to-noise ratio (PSNR). This also opens a high acceleration factor for MR data acquisition. The lightweight network is called Multi-Level Pooling Encoder-Decoder Net (MLPED Net). The proposed network outperforms the traditional encoder-decoder networks on 4-fold acceleration with a significant margin on every evaluation metric. The network can be trained end-to-end, and it is a lightweight structure that can reduce training time significantly. Experimental results are based on a publicly available MRI Knee dataset from the fastMRI competition.

SLOW: A novel spectral editing method for whole-brain MRSI at ultra high magnetic field

PURPOSE

At ultra-high field (UHF), B₁ ⁺ -inhomogeneities and high specific absorption rate (SAR) of adiabatic slice-selective RF-pulses make spatial resolved spectral-editing extremely challenging with the conventional MEGA-approach. The purpose of the study was to develop a whole-brain resolved spectral-editing MRSI at UHF (UHF, B₀  ≥ 7T) within clinical acceptable measurement-time and minimal chemical-shift-displacement-artifacts (CSDA) allowing for simultaneous GABA/Glx-, 2HG-, and PE-editing on a clinical approved 7T-scanner.

METHODS

Slice-selective adiabatic refocusing RF-pulses (2π-SSAP) dominate the SAR to the patient in (semi)LASER based MEGA-editing sequences, causing large CSDA and long measurement times to fulfill SAR requirements, even using SAR-minimized GOIA-pulses. Therefore, a novel type of spectral-editing, called SLOW-editing, using two different pairs of phase-compensated chemical-shift selective adiabatic refocusing-pulses (2π-CSAP) with different refocusing bandwidths were investigated to overcome these problems.

RESULTS

Compared to conventional echo-planar spectroscopic imaging (EPSI) and MEGA-editing, SLOW-editing shows robust refocusing and editing performance despite to B₁ ⁺ -inhomogeneity, and robustness to B₀ -inhomogeneities (0.2 ppm ≥ ΔB₀  ≥ -0.2 ppm). The narrow bandwidth (∼0.6-0.8 kHz) CSAP reduces the SAR by 92%, RF peak power by 84%, in-excitation slab CSDA by 77%, and has no in-plane CSDA. Furthermore, the CSAP implicitly dephases water, lipid and all the other signals outside of range (≥ 4.6 ppm and ≤1.4 ppm), resulting in additional water and lipid suppression (factors ≥ 1000s) at zero SAR-cost, and no spectral aliasing artifacts.

CONCLUSION

A new spectral-editing has been developed that is especially suitable for UHF, and was successfully applied for 2HG, GABA+, PE, and Glx-editing within 10 min clinical acceptable measurement time.

MULTI-parametric MR imaging with fLEXible design (MULTIPLEX)

PURPOSE

To introduce a gradient echo (GRE) -based method, namely MULTIPLEX, for single-scan 3D multi-parametric MRI with high resolution, signal-to-noise ratio (SNR), accuracy, efficiency, and acquisition flexibility.

THEORY

With a comprehensive design with dual-repetition time (TR), dual flip angle (FA), multi-echo, and optional flow modulation features, the MULTIPLEX signals contain information on radiofrequency (RF) B_1t fields, proton density, T₁ , susceptibility and blood flows, facilitating multiple qualitative images and parametric maps.

METHODS

MULTIPLEX was evaluated on system phantom and human brains, via visual inspection for image contrasts and quality or quantitative evaluation via simulation, phantom scans and literature comparison. Region-of-interest (ROI) analysis was performed on parametric maps of the system phantom and brain scans, extracting the mean and SD of the T₁ , T2∗ , proton density (PD), and/or quantitative susceptibility mapping (QSM) values for comparison with reference values or literature.

RESULTS

One MULTIPLEX scan offers multiple sets of images, including but not limited to: composited PDW/T₁ W/ T2∗ W, aT₁ W, SWI, MRA (optional), B_1t map, T₁ map, T2∗ / R2∗ maps, PD map, and QSM. The quantitative error of phantom T₁ , T2∗ and PD mapping were <5%, and those in brain scans were in good agreement with literature. MULTIPLEX scan times for high resolution (0.68 × 0.68 × 2 mm³ ) whole brain coverage were about 7.5 min, while processing times were <1 min. With flow modulation, additional MRA images can be obtained without affecting the quality or accuracy of other images.

CONCLUSION

The proposed MUTLIPLEX method possesses great potential for multi-parametric MR imaging.

Alternating Learning Approach for Variational Networks and Undersampling Pattern in Parallel MRI Applications

This work proposes an alternating learning approach to learn the sampling pattern (SP) and the parameters of variational networks (VN) in accelerated parallel magnetic resonance imaging (MRI). We investigate four variations of the learning approach, that alternates between improving the SP, using bias-accelerated subset selection, and improving parameters of the VN, using ADAM. The variations include the use of monotone or non-monotone alternating steps and systematic reduction of learning rates. The algorithms learn an effective pair to be used in future scans, including an SP that captures fewer k-space samples in which the generated undersampling artifacts are removed by the VN reconstruction. The quality of the VNs and SPs obtained by the proposed approaches is compared against different methods, including other kinds of joint learning methods and state-of-art reconstructions, on two different datasets at various acceleration factors (AF). We observed improvements visually and in three different figures of merit commonly used in deep learning (RMSE, SSIM, and HFEN) on AFs from 2 to 20 with brain and knee joint datasets when compared to the other approaches. The improvements ranged from 1% to 62% over the next best approach tested with VNs. The proposed approach has shown stable performance, obtaining similar learned SPs under different initial training conditions. We observe that the improvement is not only due to the learned sampling density, it is also due to the learned position of samples in k-space. The proposed approach was able to learn effective pairs of SPs and reconstruction VNs, improving 3D Cartesian accelerated parallel MRI applications.

Image quality and acquisition time assessments for phase oversampling in compressed sensing sensitivity encoding: Comparison with conventional SENSE

This study compared sensitivity encoding (SENSE) and compressed sensing sensitivity encoding (CS-SENSE) for phase oversampling distance and assessed its impact on image quality and image acquisition time. The experiment was performed with a large diameter phantom using 16-channel anterior body coils. All imaging data were divided into three groups according to the parallel imaging technique and oversampling distances: groups A (SENSE with phase oversampling distance of 150 mm), B (CS-SENSE with phase oversampling distance of 100 mm), and C (CS-SENSE with phase oversampling distance of 75 mm). No statistically significant differences were observed among groups A, B, and C regarding both T2 and T1 turbo spin-echo (TSE) sequences using an acceleration factor (AF) of 2 (p = 0.301 and 0.289, respectively). In comparison with AF 2 of group A, the scan time of AF 2 of groups B and C was reduced by 11.2% and 23.5% (T2 TSE) and 15.8% and 22.7% (T1 TSE), respectively, while providing comparable image quality. Significant image noise and aliasing artifact were more evident at AF ≥ 2 in group A compared with groups B and C. CS-SENSE with a less phase oversampling distance can reduce image acquisition time without image quality degradation compared with that of SENSE, despite the increase in aliasing artifact as the AF increased in both CS-SENSE and SENSE.

Fetal Neuroimaging Updates

MR imaging is used in conjunction with ultrasound screening for fetal brain abnormalities because it offers better contrast, higher resolution, and has multiplanar capabilities that increase the accuracy and confidence of diagnosis. Fetal motion still severely limits the MR imaging sequences that can be acquired. We outline the current acquisition strategies for fetal brain MR imaging and discuss the near term advances that will improve its reliability. Prospective and retrospective motion correction aim to make the complement of MR neuroimaging modalities available for fetal diagnosis, improve the performance of existing modalities, and open new horizons to understanding in utero brain development.

Fast data-driven learning of parallel MRI sampling patterns for large scale problems

In this study, a fast data-driven optimization approach, named bias-accelerated subset selection (BASS), is proposed for learning efficacious sampling patterns (SPs) with the purpose of reducing scan time in large-dimensional parallel MRI. BASS is applicable when Cartesian fully-sampled k-space measurements of specific anatomy are available for training and the reconstruction method for undersampled measurements is specified; such information is used to define the efficacy of any SP for recovering the values at the non-sampled k-space points. BASS produces a sequence of SPs with the aim of finding one of a specified size with (near) optimal efficacy. BASS was tested with five reconstruction methods for parallel MRI based on low-rankness and sparsity that allow a free choice of the SP. Three datasets were used for testing, two of high-resolution brain images ([Formula: see text]-weighted images and, respectively, [Formula: see text]-weighted images) and another of knee images for quantitative mapping of the cartilage. The proposed approach has low computational cost and fast convergence; in the tested cases it obtained SPs up to 50 times faster than the currently best greedy approach. Reconstruction quality increased by up to 45% over that provided by variable density and Poisson disk SPs, for the same scan time. Optionally, the scan time can be nearly halved without loss of reconstruction quality. Quantitative MRI and prospective accelerated MRI results show improvements. Compared with greedy approaches, BASS rapidly learns effective SPs for various reconstruction methods, using larger SPs and larger datasets; enabling better selection of sampling-reconstruction pairs for specific MRI problems.

A multi-dimensional integration (MDI) strategy for MR T2 * mapping

MRI signals are intrinsically multi-dimensional, and signal behavior may be orthogonal among different dimensions. Such dimensional orthogonality can be utilized to eliminate unwanted effects and facilitate mathematical simplicity during image processing for improved outcomes. In this work, we will demonstrate and analyze the principles and performance of a newly developed multi-dimensional integration (MDI) strategy in MR T₂ * mapping. By constructing a complex signal function to extract the inter-echo signal changes, MDI solves an optimization problem by processing all signal dimensions (eg echoes, flip angles and coil channels) in one integrative step. MDI was compared with routine curve fitting methods on noise behavior, quantification accuracy and computational efficiency. All methods were tested and compared on simulation, phantom and knee data. Monte Carlo simulations were performed on simulation and all MRI data to investigate noise propagation from k space to T₂ * maps. For phantom tests, T₂ * values in regions of interest were extracted on a voxel-wise basis and analyzed using a paired t-test between scanning parameters and mapping methods, with p < 0.05 being significantly different. MDI facilitated a straightforward processing procedure, yielding homogeneous, high-signal-to-noise-ratio (SNR) and artifact-free T₂ * maps without explicit coil combination or additional measures. Compared with routine fitting methods, MDI offered significantly (p < 0.05) improved SNR and quantitative accuracy/robustness, with two to three orders higher computational efficiency. MDI also represented low-SNR signals with low T₂ * values, avoiding misinterpretation with long-T₂ * species.

Diagnostic Accuracy of Quantitative Multicontrast 5-Minute Knee MRI Using Prospective Artificial Intelligence Image Quality Enhancement

BACKGROUND. Potential approaches for abbreviated knee MRI, including prospective acceleration with deep learning, have achieved limited clinical implementation. OBJECTIVE. The objective of this study was to evaluate the interreader agreement between conventional knee MRI and a 5-minute 3D quantitative double-echo steady-state (qDESS) sequence with automatic T2 mapping and deep learning super-resolutionaugmentation and to compare the diagnostic performance of the two methods regarding findings from arthroscopic surgery. METHODS. Fifty-one patients with knee pain underwent knee MRI that included an additional 3D qDESS sequence with automatic T2 mapping. Fourier interpolation was followed by prospective deep learning super resolution to enhance qDESS slice resolution twofold. A musculoskeletal radiologist and a radiology resident performed independent retrospective evaluations of articular cartilage, menisci, ligaments, bones, extensor mechanism, and synovium using conventional MRI. Following a 2-month washout period, readers reviewed qDESS images alone followed by qDESS with the automatic T2 maps. Interreader agreement between conventional MRI and qDESS was computed using percentage agreement and Cohen kappa. The sensitivity and specificity of conventional MRI, qDESS alone, and qDESS plus T2 mapping were compared with arthroscopic findings using exact McNemar tests. RESULTS. Conventional MRI and qDESS showed 92% agreement in evaluating all tissues. Kappa was 0.79 (95% CI, 0.76-0.81) across all imaging findings. In 43 patients who underwent arthroscopy, sensitivity and specificity were not significantly different (p = .23 to > .99) between conventional MRI (sensitivity, 58-93%; specificity, 27-87%) and qDESS alone (sensitivity, 54-90%; specificity, 23-91%) for cartilage, menisci, ligaments, and synovium. For grade 1 cartilage lesions, sensitivity and specificity were 33% and 56%, respectively, for conventional MRI; 23% and 53% for qDESS (p = .81); and 46% and 39% for qDESS with T2 mapping (p = .80). For grade 2A lesions, values were 27% and 53% for conventional MRI, 26% and 52% for qDESS (p = .02), and 58% and 40% for qDESS with T2 mapping (p < .001). CONCLUSION. The qDESS method prospectively augmented with deep learning showed strong interreader agreement with conventional knee MRI and near-equivalent diagnostic performance regarding arthroscopy. The ability of qDESS to automatically generate T2 maps increases sensitivity for cartilage abnormalities. CLINICAL IMPACT. Using prospective artificial intelligence to enhance qDESS image quality may facilitate an abbreviated knee MRI protocol while generating quantitative T2 maps.

Motion correction methods for MRS: experts' consensus recommendations

Long acquisition times due to intrinsically low signal-to-noise ratio and the need for highly homogeneous B₀ field make MRS particularly susceptible to motion or scanner instability compared with MRI. Motion-induced changes in both localization and shimming (ie B₀ homogeneity) degrade MRS data quality. To mitigate the effects of motion three approaches can be employed: (1) subject immobilization, (2) retrospective correction, and (3) prospective real-time correction using internal and/or external tracking methods. Prospective real-time correction methods can simultaneously update localization and the B₀ field to improve MRS data quality. While localization errors can be corrected with both internal (navigators) and external (optical camera, NMR probes) tracking methods, the B₀ field correction requires internal navigator methods to measure the B₀ field inside the imaged volume and the possibility to update the scanner shim hardware in real time. Internal and external tracking can rapidly update the MRS localization with submillimeter and subdegree precision, while scanner frequency and first-order shims of scanner hardware can be updated by internal methods every sequence repetition. These approaches are most well developed for neuroimaging, for which rigid transformation is primarily applicable. Real-time correction greatly improves the stability of MRS acquisition and quantification, as shown in clinical studies on subjects prone to motion, including children and patients with movement disorders, enabling robust measurement of metabolite signals including those with low concentrations, such as gamma-aminobutyric acid and glutathione. Thus, motion correction is recommended for MRS users and calls for tighter integration and wider availability of such methods by MR scanner manufacturers.

Susceptibility Weighted Imaging for evaluation of musculoskeletal lesions

The presence of blood or calcium in the musculoskeletal (MSK) system may be linked to specific pathological conditions. The ability of MRI for calcium detection is usually limited compared with other techniques such as CT. In a similar manner, the accuracy of MRI for detection and evaluation of hemorrhage in soft tissues is closely linked to the degree of degradation of blood products. Blood and calcium are substances that cause local inhomogeneity of the magnetic field resulting in susceptibility artifacts. To try to evaluate these substances, specific MRI sequences which are highly sensitive to these local magnetic field inhomogeneities such as Susceptibility Weighted Imaging (SWI) have been developed and successfully applied in the Central Nervous System, but scarcely used in MSK. SWI may increase the overall sensitivity of MRI to detect blood and calcium in several clinical scenarios such as degenerative joint disease or bone and soft tissue lesion assessment and discriminate between both compounds, something which is not always possible with conventional MRI approaches. In this paper, physical basis and technical adjustment for SWI acquisition at MSK are detailed reviewing the potential application of SWI in different MSK clinical scenarios.

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.

Mapping the rest of the human connectome: Atlasing the spinal cord and peripheral nervous system

The emergence of diffusion, structural, and functional neuroimaging methods has enabled major multi-site efforts to map the human connectome, which has heretofore been defined as containing all neural connections in the central nervous system (CNS). However, these efforts are not structured to examine the richness and complexity of the peripheral nervous system (PNS), which arguably forms the (neglected) rest of the connectome. Despite increasing interest in an atlas of the spinal cord (SC) and PNS which is simultaneously stereotactic, interactive, electronically dissectible, scalable, population-based and deformable, little attention has thus far been devoted to this task of critical importance. Nevertheless, the atlasing of these complete neural structures is essential for neurosurgical planning, neurological localization, and for mapping those components of the human connectome located outside of the CNS. Here we recommend a modification to the definition of the human connectome to include the SC and PNS, and argue for the creation of an inclusive atlas to complement current efforts to map the brain's human connectome, to enhance clinical education, and to assist progress in neuroscience research. In addition to providing a critical overview of existing neuroimaging techniques, image processing methodologies and algorithmic advances which can be combined for the creation of a full connectome atlas, we outline a blueprint for ultimately mapping the entire human nervous system and, thereby, for filling a critical gap in our scientific knowledge of neural connectivity.

Accelerating cardiac cine MRI using a deep learning-based ESPIRiT reconstruction

PURPOSE

To propose a novel combined parallel imaging and deep learning-based reconstruction framework for robust reconstruction of highly accelerated 2D cardiac cine MRI data.

METHODS

We propose DL-ESPIRiT, an unrolled neural network architecture that utilizes an extended coil sensitivity model to address SENSE-related field-of-view (FOV) limitations in previously proposed deep learning-based reconstruction frameworks. Additionally, we propose a novel neural network design based on (2+1)D spatiotemporal convolutions to produce more accurate dynamic MRI reconstructions than conventional 3D convolutions. The network is trained on fully sampled 2D cardiac cine datasets collected from 11 healthy volunteers with IRB approval. DL-ESPIRiT is compared against a state-of-the-art parallel imaging and compressed sensing method known as l 1 -ESPIRiT. The reconstruction accuracy of both methods is evaluated on retrospectively undersampled datasets (R = 12) with respect to standard image quality metrics as well as automatic deep learning-based segmentations of left ventricular volumes. Feasibility of DL-ESPIRiT is demonstrated on two prospectively undersampled datasets acquired in a single heartbeat per slice.

RESULTS

The (2+1)D DL-ESPIRiT method produces higher fidelity image reconstructions when compared to l 1 -ESPIRiT reconstructions with respect to standard image quality metrics (P < .001). As a result of improved image quality, segmentations made from (2+1)D DL-ESPIRiT images are also more accurate than segmentations from l 1 -ESPIRiT images.

CONCLUSIONS

DL-ESPIRiT synergistically combines a robust parallel imaging model and deep learning-based priors to produce high-fidelity reconstructions of retrospectively undersampled 2D cardiac cine data acquired with reduced FOV. Although a proof-of-concept is shown, further experiments are necessary to determine the efficacy of DL-ESPIRiT in prospectively undersampled 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.

Rapid Knee MRI Acquisition and Analysis Techniques for Imaging Osteoarthritis

Osteoarthritis (OA) of the knee is a major source of disability that has no known treatment or cure. Morphological and compositional MRI is commonly used for assessing the bone and soft tissues in the knee to enhance the understanding of OA pathophysiology. However, it is challenging to extend these imaging methods and their subsequent analysis techniques to study large population cohorts due to slow and inefficient imaging acquisition and postprocessing tools. This can create a bottleneck in assessing early OA changes and evaluating the responses of novel therapeutics. The purpose of this review article is to highlight recent developments in tools for enhancing the efficiency of knee MRI methods useful to study OA. Advances in efficient MRI data acquisition and reconstruction tools for morphological and compositional imaging, efficient automated image analysis tools, and hardware improvements to further drive efficient imaging are discussed in this review. For each topic, we discuss the current challenges as well as potential future opportunities to alleviate these challenges. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY STAGE: 3.

Methodological Considerations for Neuroimaging in Deep Brain Stimulation of the Subthalamic Nucleus in Parkinson's Disease Patients

Deep brain stimulation (DBS) of the subthalamic nucleus is a neurosurgical intervention for Parkinson's disease patients who no longer appropriately respond to drug treatments. A small fraction of patients will fail to respond to DBS, develop psychiatric and cognitive side-effects, or incur surgery-related complications such as infections and hemorrhagic events. In these cases, DBS may require recalibration, reimplantation, or removal. These negative responses to treatment can partly be attributed to suboptimal pre-operative planning procedures via direct targeting through low-field and low-resolution magnetic resonance imaging (MRI). One solution for increasing the success and efficacy of DBS is to optimize preoperative planning procedures via sophisticated neuroimaging techniques such as high-resolution MRI and higher field strengths to improve visualization of DBS targets and vasculature. We discuss targeting approaches, MRI acquisition, parameters, and post-acquisition analyses. Additionally, we highlight a number of approaches including the use of ultra-high field (UHF) MRI to overcome limitations of standard settings. There is a trade-off between spatial resolution, motion artifacts, and acquisition time, which could potentially be dissolved through the use of UHF-MRI. Image registration, correction, and post-processing techniques may require combined expertise of traditional radiologists, clinicians, and fundamental researchers. The optimization of pre-operative planning with MRI can therefore be best achieved through direct collaboration between researchers and clinicians.

Calibrationless parallel imaging reconstruction for multislice MR data using low-rank tensor completion

PURPOSE

To provide joint calibrationless parallel imaging reconstruction of highly accelerated multislice 2D MR k-space data.

METHODS

Adjacent image slices in multislice MR data have similar coil sensitivity maps, spatial support, and image content. Such similarities can be utilized to improve image quality by reconstructing multiple slices jointly with low-rank tensor completion. Specifically, the multichannel k-space data from multiple slices are constructed into a block-wise Hankel tensor and iteratively updated by promoting tensor low-rankness through higher-order SVD. This multislice block-wise Hankel tensor completion was implemented for 2D spiral and Cartesian k-space undersampling where sampling patterns vary between adjacent slices. The approach was evaluated with human brain MR data and compared to the traditional single-slice simultaneous autocalibrating and k-space estimation reconstruction.

RESULTS

The proposed multislice block-wise Hankel tensor completion approach robustly reconstructed highly undersampled multislice 2D spiral and Cartesian data. It produced substantially lower level of artifacts compared to the traditional single-slice simultaneous autocalibrating and k-space estimation reconstruction. Quantitative evaluation using error maps and root mean square error demonstrated its significantly improved performance in terms of residual artifacts and root mean square error.

CONCLUSION

Our proposed multislice block-wise Hankel tensor completion method exploits the similar coil sensitivity and image content within multislice MR data through a tensor completion framework. It offers a new and effective approach to acquire and reconstruct highly undersampled multislice MR data in a calibrationless manner.

Simultaneous use of individual and joint regularization terms in compressive sensing: Joint reconstruction of multi-channel multi-contrast MRI acquisitions

Multi-contrast images are commonly acquired together to maximize complementary diagnostic information, albeit at the expense of longer scan times. A time-efficient strategy to acquire high-quality multi-contrast images is to accelerate individual sequences and then reconstruct undersampled data with joint regularization terms that leverage common information across contrasts. However, these terms can cause features that are unique to a subset of contrasts to leak into the other contrasts. Such leakage-of-features may appear as artificial tissues, thereby misleading diagnosis. The goal of this study is to develop a compressive sensing method for multi-channel multi-contrast magnetic resonance imaging (MRI) that optimally utilizes shared information while preventing feature leakage. Joint regularization terms group sparsity and colour total variation are used to exploit common features across images while individual sparsity and total variation are also used to prevent leakage of distinct features across contrasts. The multi-channel multi-contrast reconstruction problem is solved via a fast algorithm based on Alternating Direction Method of Multipliers. The proposed method is compared against using only individual and only joint regularization terms in reconstruction. Comparisons were performed on single-channel simulated and multi-channel in-vivo datasets in terms of reconstruction quality and neuroradiologist reader scores. The proposed method demonstrates rapid convergence and improved image quality for both simulated and in-vivo datasets. Furthermore, while reconstructions that solely use joint regularization terms are prone to leakage-of-features, the proposed method reliably avoids leakage via simultaneous use of joint and individual terms, thereby holding great promise for clinical use.

Does higher sampling rate (multiband + SENSE) improve group statistics - An example from social neuroscience block design at 3T

Multiband (MB) or Simultaneous multi-slice (SMS) acquisition schemes allow the acquisition of MRI signals from more than one spatial coordinate at a time. Commercial availability has brought this technique within the reach of many neuroscientists and psychologists. Most early evaluation of the performance of MB acquisition employed resting state fMRI or the most basic tasks. In this study, we tested whether the advantages of using MB acquisition schemes generalize to group analyses using a cognitive task more representative of typical cognitive neuroscience applications. Twenty-three subjects were scanned on a Philips 3 ​T scanner using five sequences, up to eight-fold acceleration with MB-factors 1 to 4, SENSE factors up to 2 and corresponding TRs of 2.45s down to 0.63s, while they viewed (i) movie blocks showing complex actions with hand object interactions and (ii) control movie blocks without hand object interaction. Data were processed using a widely used analysis pipeline implemented in SPM12 including the unified segmentation and canonical HRF modelling. Using random effects group-level, voxel-wise analysis we found that all sequences were able to detect the basic action observation network known to be recruited by our task. The highest t-values were found for sequences with MB4 acceleration. For the MB1 sequence, a 50% bigger voxel volume was needed to reach comparable t-statistics. The group-level t-values for resting state networks (RSNs) were also highest for MB4 sequences. Here the MB1 sequence with larger voxel size did not perform comparable to the MB4 sequence. Altogether, we can thus recommend the use of MB4 (and SENSE 1.5 or 2) on a Philips scanner when aiming to perform group-level analyses using cognitive block design fMRI tasks and voxel sizes in the range of cortical thickness (e.g. 2.7 ​mm isotropic). While results will not be dramatically changed by the use of multiband, our results suggest that MB will bring a moderate but significant benefit.

Utility of deep learning super-resolution in the context of osteoarthritis MRI biomarkers

BACKGROUND

Super-resolution is an emerging method for enhancing MRI resolution; however, its impact on image quality is still unknown.

PURPOSE

To evaluate MRI super-resolution using quantitative and qualitative metrics of cartilage morphometry, osteophyte detection, and global image blurring.

STUDY TYPE

Retrospective.

POPULATION

In all, 176 MRI studies of subjects at varying stages of osteoarthritis.

FIELD STRENGTH/SEQUENCE

Original-resolution 3D double-echo steady-state (DESS) and DESS with 3× thicker slices retrospectively enhanced using super-resolution and tricubic interpolation (TCI) at 3T.

ASSESSMENT

A quantitative comparison of femoral cartilage morphometry was performed for the original-resolution DESS, the super-resolution, and the TCI scans in 17 subjects. A reader study by three musculoskeletal radiologists assessed cartilage image quality, overall image sharpness, and osteophytes incidence in all three sets of scans. A referenceless blurring metric evaluated blurring in all three image dimensions for the three sets of scans.

STATISTICAL TESTS

Mann-Whitney U-tests compared Dice coefficients (DC) of segmentation accuracy for the DESS, super-resolution, and TCI images, along with the image quality readings and blurring metrics. Sensitivity, specificity, and diagnostic odds ratio (DOR) with 95% confidence intervals compared osteophyte detection for the super-resolution and TCI images, with the original-resolution as a reference.

RESULTS

DC for the original-resolution (90.2 ± 1.7%) and super-resolution (89.6 ± 2.0%) were significantly higher (P < 0.001) than TCI (86.3 ± 5.6%). Segmentation overlap of super-resolution with the original-resolution (DC = 97.6 ± 0.7%) was significantly higher (P < 0.0001) than TCI overlap (DC = 95.0 ± 1.1%). Cartilage image quality for sharpness and contrast levels, and the through-plane quantitative blur factor for super-resolution images, was significantly (P < 0.001) better than TCI. Super-resolution osteophyte detection sensitivity of 80% (76-82%), specificity of 93% (92-94%), and DOR of 32 (22-46) was significantly higher (P < 0.001) than TCI sensitivity of 73% (69-76%), specificity of 90% (89-91%), and DOR of 17 (13-22).

DATA CONCLUSION

Super-resolution appears to consistently outperform naïve interpolation and may improve image quality without biasing quantitative biomarkers.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2020;51:768-779.

FPGA-based hardware accelerator for SENSE (a parallel MR image reconstruction method)

SENSE (Sensitivity Encoding) is a parallel MRI (pMRI) technique that allows accelerated data acquisition using multiple receiver coils and reconstructs the artifact-free images from the acquired under-sampled data. However, an increasing number of receiver coils has raised the computational demands of pMRI techniques to an extent where the reconstruction time on general purpose computers becomes impractically long for real-time MRI. Field Programmable Gate Arrays (FPGAs) have recently emerged as a viable hardware platform for accelerating pMRI algorithms (e.g. SENSE). However, recent efforts to accelerate SENSE using FPGAs have been focused on a fixed number of receiver coils (L=8) and acceleration factor (A_f=2). This paper presents a novel 32-bit floating-point FPGA-based hardware accelerator for SENSE (HW-ACC-SENSE); having an ability to work in coordination with an on-chip ARM processor performing reconstructions for different values of L and A_f. Moreover, the proposed design provides flexibility to integrate multiple units of HW-ACC-SENSE with an on-chip ARM processor, for low-latency image reconstruction. The VIVADO High-Level-Synthesis (HLS) tool has been used to design and implement the HW-ACC-SENSE on the Xilinx FPGA development board (ZCU102). A series of experiments has been performed on in-vivo datasets acquired using 8, 12 and 30 receiver coil elements. The performance of the proposed architecture is compared with the single thread and multi-thread CPU-based implementations of SENSE. The results show that the proposed design withstands the reconstruction quality of the SENSE algorithm while demonstrating a maximum speed-gain up to 298× over the CPU counterparts in our experiments.

Parallel magnetic resonance image reconstruction from a single-element parametric amplifier

In magnetic resonance imaging (MRI), acquisition speed is always an important issue. In this paper, we propose a promising technique to achieve parallel MRI (pMRI) on a single-channel spectrometer, using a novel Wireless Amplified Nuclear MR Detector (WAND) for spatial encoding in image reconstruction. For this, a planar structure double frequency WAND is designed and fabricated, where two of its frequencies - 'signal', ω1 and 'idler', ω2 are effectively utilized as two separate "channels" for accelerated acquisition. We provided a thorough background needed for the method and subsequently parallel imaging algorithms. Sum-of-Squares (SoS) reconstruction and GeneRalized Autocalibrating Partially Parallel Acquisition (GRAPPA) reconstruction are used to reconstruct as well as to analyze the SNR in the resulting images and validate our hypothesis. Experimental results using phantom datasets demonstrate that the proposed method of parallel imaging yield a better sensitivity for the combined images ('idler' + 'signal') than the sensitivity acquired for each individual image and thus significantly improving the reconstruction quality with optimal signal-to-noise ratio. We also demonstrated the achievable acceleration factor of this approach.

Time optimal control-based RF pulse design under gradient imperfections

Purpose

This study incorporates a gradient system imperfection model into an optimal control framework for radio frequency (RF) pulse design.

Theory and Methods

The joint design of minimum‐time RF and slice selective gradient shapes is posed as an optimal control problem. Hardware limitations such as maximal amplitudes for RF and slice selective gradient or its slew rate are included as hard constraints to assure practical applicability of the optimized waveforms. In order to guarantee the performance of the optimized waveform with possible gradient system disturbances such as limited system bandwidth and eddy currents, a measured gradient impulse response function (GIRF) for a specific system is integrated into the optimization.

Results

The method generates optimized RF and pre‐distorted slice selective gradient shapes for refocusing that are able to fully compensate the modeled imperfections of the gradient system under investigation. The results nearly regenerate the optimal results of an idealized gradient system. The numerical Bloch simulations are validated by phantom and in‐vivo experiments on 2 3T scanners.

Conclusions

The presented design approach demonstrates the successful correction of gradient system imperfections within an optimal control framework for RF pulse design.