Large-scale 3D non-Cartesian coronary MRI reconstruction using distributed memory-efficient physics-guided deep learning with limited training data.

OBJECT: To enable high-quality physics-guided deep learning (PG-DL) reconstruction of large-scale 3D non-Cartesian coronary MRI by overcoming challenges of hardware limitations and limited training data availability. MATERIALS AND METHODS: While PG-DL has emerged as a powerful image reconstruction method, its application to large-scale 3D non-Cartesian MRI is hindered by hardware limitations and limited availability of training data. We combine several recent advances in deep learning and MRI reconstruction to tackle the former challenge, and we further propose a 2.5D reconstruction using 2D convolutional neural networks, which treat 3D volumes as batches of 2D images to train the network with a limited amount of training data. Both 3D and 2.5D variants of the PG-DL networks were compared to conventional methods for high-resolution 3D kooshball coronary MRI. RESULTS: Proposed PG-DL reconstructions of 3D non-Cartesian coronary MRI with 3D and 2.5D processing outperformed all conventional methods both quantitatively and qualitatively in terms of image assessment by an experienced cardiologist. The 2.5D variant further improved vessel sharpness compared to 3D processing, and scored higher in terms of qualitative image quality. DISCUSSION: PG-DL reconstruction of large-scale 3D non-Cartesian MRI without compromising image size or network complexity is achieved, and the proposed 2.5D processing enables high-quality reconstruction with limited training data.

Combined thalamic and pallidal deep brain stimulation in diabetic hemiballism/hemichorea.

Hemiballism/hemichorea (HH) is a hyperkinetic movement disorder observed mostly in older adults with cerebrovascular diseases. Although the symptoms improve without any treatment, lesioning or DBS (deep brain stimulation) may be rarely required to provide symptomatic relief for patients with severe involuntary movements. HH is a rare complication of uncontrolled diabetes. There are only a few reported cases of diabetic HH that have been surgically treated. Thus, herein, we have reported the case of a 75-year-old female with type-II diabetes mellitus that presented with disabling involuntary limb movements of the left side, despite being treated conservatively for six months. DBS targeting the globus pallidus internus (GPi) and ventral intermediate (Vim) thalamic nucleus was performed. Complete resolution of symptoms was achieved with a combined stimulation of the thalamic Vim nucleus (at 1.7 mA) and GPi (at 2.4 mA). The combined stimulation of the Vim nucleus and GPi effectively resolved the diabetes-induced HH symptoms in our patient. Thus, although certain conclusions cannot be drawn due to the rarity of the surgically treated patients with HH, the combined stimulation is a novel treatment option for resistant HH.

Gamma Knife radiosurgery for the treatment of trigeminal neuralgia: A single center-experience.

INTRODUCTION AND OBJECTIVES: We aimed to assess the outcomes of patients with trigeminal neuralgia (TGN) who underwent Gamma Knife radiosurgery (GKRS). MATERIALS AND METHODS: Fifty-three patients with typical TGN underwent GKRS from May 2012 until December 2022. Among these patients, 45 patients who were follow-up for at least 12 months were included in the study. A mean dose of 87.5 Gy (range, 80-90) was administered to the trigeminal nerve. Postoperatively, outcome was considered excellent if the patient was pain- and medication-free. RESULTS: The mean symtpom duration was 9.53 years, and the mean patient age was 59.8 years (range, 34-85). The mean follow-up period was 46.8 months (range, 12-127 months). 46.7% of patients had a history of previous surgical interventions. A single nerve division was affected in 14 patients (31.1%), and multiple divisions were affected in 31 patients (68.9%). The rate of initial pain relief was 80%. Hypoesthesia in the area of trigeminal nerve developed in 30 (66.7%). Twenty patients (44.4%) exhibited excellent results within 72.4 months. Recurrence occurred in 11 patients (24.4%) with 27.6 months. CONCLUSIONS: Our results suggest that GKRS is a safe and effective procedure. Thus, it is an attractive first- and second-line treatment choice for TGN.

Unsupervised deep learning with convolutional neural networks for static parallel transmit design: A retrospective study.

PURPOSE: To mitigate B 1 + $$ {B}_1^{+} $$ inhomogeneity at 7T for multi-channel transmit arrays using unsupervised deep learning with convolutional neural networks (CNNs). METHODS: Deep learning parallel transmit (pTx) pulse design has received attention, but such methods have relied on supervised training and did not use CNNs for multi-channel B 1 + $$ {B}_1^{+} $$ maps. In this work, we introduce an alternative approach that facilitates the use of CNNs with multi-channel B 1 + $$ {B}_1^{+} $$ maps while performing unsupervised training. The multi-channel B 1 + $$ {B}_1^{+} $$ maps are concatenated along the spatial dimension to enable shift-equivariant processing amenable to CNNs. Training is performed in an unsupervised manner using a physics-driven loss function that minimizes the discrepancy of the Bloch simulation with the target magnetization, which eliminates the calculation of reference transmit RF weights. The training database comprises 3824 2D sagittal, multi-channel B 1 + $$ {B}_1^{+} $$ maps of the healthy human brain from 143 subjects. B 1 + $$ {B}_1^{+} $$ data were acquired at 7T using an 8Tx/32Rx head coil. The proposed method is compared to the unregularized magnitude least-squares (MLS) solution for the target magnetization in static pTx design. RESULTS: The proposed method outperformed the unregularized MLS solution for RMS error and coefficient-of-variation and had comparable energy consumption. Additionally, the proposed method did not show local phase singularities leading to distinct holes in the resulting magnetization unlike the unregularized MLS solution. CONCLUSION: Proposed unsupervised deep learning with CNNs performs better than unregularized MLS in static pTx for speed and robustness.

Locally low-rank denoising in transform domains.

Purpose: To develop an extension to locally low rank (LLR) denoising techniques based on transform domain processing that reduces the number of images required in the MR image series for high-quality denoising. Theory and Methods: LLR methods with random matrix theory-based thresholds are successfully used in the denoising of MR image series in a number of applications. The performance of these methods depend on how well the LLR assumption is satisfied, which deteriorates with few numbers of images, as is commonly encountered in quantitative MRI applications. We propose a transform-domain approach for denoising of MR image series to represent the underlying signal with higher fidelity when using a locally low rank approximation. The efficacy of the method is demonstrated for fully-sampled k-space, undersampled k-space, DICOM images, and complex-valued SENSE-1 images in quantitative MRI applications with as few as 4 images. Results: For both MSK and brain applications, the transform domain denoising preserves local subtle variability, whereas the quantitative maps based on image domain LLR methods tend to be locally more homogeneous. Conclusion: A transform domain extension to LLR denoising produces high quality images and is compatible with both raw k-space data and vendor reconstructed data. This allows for improved imaging and more accurate quantitative analyses and parameters obtained therefrom.

High-fidelity Database-free Deep Learning Reconstruction for Real-time Cine Cardiac MRI.

Real-time cine cardiac MRI provides an ECG-free free-breathing alternative to clinical gold-standard ECG-gated breath-hold segmented cine MRI for evaluation of heart function. Real-time cine MRI data acquisition during free breathing snapshot imaging enables imaging of patient cohorts that cannot be imaged with segmented or breath-hold acquisitions, but requires rapid imaging to achieve sufficient spatial-temporal resolutions. However, at high acceleration rates, conventional reconstruction techniques suffer from residual aliasing and temporal blurring, including advanced methods such as compressed sensing with radial trajectories. Recently, deep learning (DL) reconstruction has emerged as a powerful tool in MRI. However, its utility for free-breathing real-time cine MRI has been limited, as database-learning of spatio-temporal correlations with varying breathing and cardiac motion patterns across subjects has been challenging. Zero-shot self-supervised physics-guided deep learning (PG-DL) reconstruction has been proposed to overcome such challenges of database training by enabling subject-specific training. In this work, we adapt zero-shot PG-DL for real-time cine MRI with a spatio-temporal regularization. We compare our method to TGRAPPA, locally low-rank (LLR) regularized reconstruction and database-trained PG-DL reconstruction, both for retrospectively and prospectively accelerated datasets. Results on highly accelerated real-time Cartesian cine MRI show that the proposed method outperforms other reconstruction methods, both visibly in terms of noise and aliasing, and quantitatively.

Improved Simultaneous Multi-slice imaging with Composition of k-space Interpolations (SMS-COOKIE) for myocardial T1 mapping.

The aim of this study is to develop and evaluate a regularized Simultaneous Multi-Slice (SMS) reconstruction method for improved Cardiac Magnetic Resonance Imaging (CMR). The proposed reconstruction method, SMS with COmpOsition of k-space IntErpolations (SMS-COOKIE) combines the advantages of Iterative Self-consistent Parallel Imaging Reconstruction (SPIRiT) and split slice-Generalized Autocalibrating Partially Parallel Acquisitions (GRAPPA), while allowing regularization for further noise reduction. The proposed SMS-COOKIE was implemented with and without regularization, and validated using a Saturation Pulse-Prepared Heart rate Independent inversion REcovery (SAPPHIRE) myocardial T1 mapping sequence. The performance of the proposed reconstruction method was compared to ReadOut (RO)-SENSE-GRAPPA and split slice-GRAPPA, on both retrospectively and prospectively three-fold SMS-accelerated data with an additional two-fold in-plane acceleration. All SMS reconstruction methods yielded similar T1 values compared to single band imaging. SMS-COOKIE showed lower spatial variability in myocardial T1 with significant improvement over RO-SENSE-GRAPPA and split slice-GRAPPA (P < 10-4). The proposed method with additional locally low rank (LLR) regularization reduced the spatial variability, again with significant improvement over RO-SENSE-GRAPPA and split slice-GRAPPA (P < 10-4). In conclusion, improved reconstruction quality was achieved with the proposed SMS-COOKIE, which also provided lower spatial variability with significant improvement over split slice-GRAPPA.

Physics-Driven Deep Learning for Computational Magnetic Resonance Imaging: Combining physics and machine learning for improved medical imaging.

Physics-driven deep learning methods have emerged as a powerful tool for computational magnetic resonance imaging (MRI) problems, pushing reconstruction performance to new limits. This article provides an overview of the recent developments in incorporating physics information into learning-based MRI reconstruction. We consider inverse problems with both linear and non-linear forward models for computational MRI, and review the classical approaches for solving these. We then focus on physics-driven deep learning approaches, covering physics-driven loss functions, plug-and-play methods, generative models, and unrolled networks. We highlight domain-specific challenges such as real- and complex-valued building blocks of neural networks, and translational applications in MRI with linear and non-linear forward models. Finally, we discuss common issues and open challenges, and draw connections to the importance of physics-driven learning when combined with other downstream tasks in the medical imaging pipeline.

Robust cardiac T 1 ρ $$ {\mathrm{T}}_{1_{\boldsymbol{\rho}}} $$ mapping at 3T using adiabatic spin-lock preparations.

PURPOSE: The aim of this study is to develop and optimize an adiabatic T 1 ρ $$ {\mathrm{T}}_{1\uprho} $$ ( T 1 ρ , adiab $$ {\mathrm{T}}_{1\uprho, \mathrm{adiab}} $$ ) mapping method for robust quantification of spin-lock (SL) relaxation in the myocardium at 3T. METHODS: Adiabatic SL (aSL) preparations were optimized for resilience against B 0 $$ {\mathrm{B}}_0 $$ and B 1 + $$ {\mathrm{B}}_1^{+} $$ inhomogeneities using Bloch simulations. Optimized B 0 $$ {\mathrm{B}}_0 $$ -aSL, Bal-aSL and B 1 $$ {\mathrm{B}}_1 $$ -aSL modules, each compensating for different inhomogeneities, were first validated in phantom and human calf. Myocardial T 1 ρ $$ {\mathrm{T}}_{1\uprho} $$ mapping was performed using a single breath-hold cardiac-triggered bSSFP-based sequence. Then, optimized T 1 ρ , adiab $$ {\mathrm{T}}_{1\uprho, \mathrm{adiab}} $$ preparations were compared to each other and to conventional SL-prepared T 1 ρ $$ {\mathrm{T}}_{1\uprho} $$ maps (RefSL) in phantoms to assess repeatability, and in 13 healthy subjects to investigate image quality, precision, reproducibility and intersubject variability. Finally, aSL and RefSL sequences were tested on six patients with known or suspected cardiovascular disease and compared with LGE, T 1 $$ {\mathrm{T}}_1 $$ , and ECV mapping. RESULTS: The highest T 1 ρ , adiab $$ {\mathrm{T}}_{1\uprho, \mathrm{adiab}} $$ preparation efficiency was obtained in simulations for modules comprising 2 HS pulses of 30 ms each. In vivo T 1 ρ , adiab $$ {\mathrm{T}}_{1\uprho, \mathrm{adiab}} $$ maps yielded significantly higher quality than RefSL maps. Average myocardial T 1 ρ , adiab $$ {\mathrm{T}}_{1\uprho, \mathrm{adiab}} $$ values were 183.28 ± $$ \pm $$ 25.53 ms, compared with 38.21 ± $$ \pm $$ 14.37 ms RefSL-prepared T 1 ρ $$ {\mathrm{T}}_{1\uprho} $$ . T 1 ρ , adiab $$ {\mathrm{T}}_{1\uprho, \mathrm{adiab}} $$ maps showed a significant improvement in precision (avg. 14.47 ± $$ \pm $$ 3.71% aSL, 37.61 ± $$ \pm $$ 19.42% RefSL, p < 0.01) and reproducibility (avg. 4.64 ± $$ \pm $$ 2.18% aSL, 47.39 ± $$ \pm $$ 12.06% RefSL, p < 0.0001), with decreased inter-subject variability (avg. 8.76 ± $$ \pm $$ 3.65% aSL, 51.90 ± $$ \pm $$ 15.27% RefSL, p < 0.0001). Among aSL preparations, B 0 $$ {\mathrm{B}}_0 $$ -aSL achieved the better inter-subject variability. In patients, B 1 $$ {\mathrm{B}}_1 $$ -aSL preparations showed the best artifact resilience among the adiabatic preparations. T 1 ρ , adiab $$ {\mathrm{T}}_{1\uprho, \mathrm{adiab}} $$ times show focal alteration colocalized with areas of hyper-enhancement in the LGE images. CONCLUSION: Adiabatic preparations enable robust in vivo quantification of myocardial SL relaxation times at 3T.

High-fidelity Database-free Deep Learning Reconstruction for Real-time Cine Cardiac MRI.

Real-time cine cardiac MRI provides an ECG-free free-breathing alternative to clinical gold-standard ECG-gated breath-hold segmented cine MRI for evaluation of heart function. Real-time cine MRI data acquisition during free breathing snapshot imaging enables imaging of patient cohorts that cannot be imaged with segmented or breath-hold acquisitions, but requires rapid imaging to achieve sufficient spatial-temporal resolutions. However, at high acceleration rates, conventional reconstruction techniques suffer from residual aliasing and temporal blurring, including advanced methods such as compressed sensing with radial trajectories. Recently, deep learning (DL) reconstruction has emerged as a powerful tool in MRI. However, its utility for free-breathing real-time cine MRI has been limited, as database-learning of spatio-temporal correlations with varying breathing and cardiac motion patterns across subjects has been challenging. Zero-shot self-supervised physics-guided deep learning (PG-DL) reconstruction has been proposed to overcome such challenges of database training by enabling subject-specific training. In this work, we adapt zero-shot PG-DL for real-time cine MRI with a spatio-temporal regularization. We compare our method to TGRAPPA, locally low-rank (LLR) regularized reconstruction and database-trained PG-DL reconstruction, both for retrospectively and prospectively accelerated datasets. Results on highly accelerated real-time Cartesian cine MRI show that the proposed method outperforms other reconstruction methods, both visibly in terms of noise and aliasing, and quantitatively.

Evaluating increases in sensitivity from NORDIC for diverse fMRI acquisition strategies.

As the neuroimaging field moves towards detecting smaller effects at higher spatial resolutions, and faster sampling rates, there is increased attention given to the deleterious contribution of unstructured, thermal noise. Here, we critically evaluate the performance of a recently developed reconstruction method, termed NORDIC, for suppressing thermal noise using datasets acquired with various field strengths, voxel sizes, sampling rates, and task designs. Following minimal preprocessing, statistical activation (t-values) of NORDIC processed data was compared to the results obtained with alternative denoising methods. Additionally, we examined the consistency of the estimates of task responses at the single-voxel, single run level, using a finite impulse response (FIR) model. To examine the potential impact on effective image resolution, the overall smoothness of the data processed with different methods was estimated. Finally, to determine if NORDIC alters or removes temporal information important for modeling responses, we employed an exhaustive leave-p-out cross validation approach, using FIR task responses to predict held out timeseries, quantified using R2. After NORDIC, the t-values are increased, an improvement comparable to what could be achieved by 1.5 voxels smoothing, and task events are clearly visible and have less cross-run error. These advantages are achieved with smoothness estimates increasing by less than 4%, while 1.5 voxel smoothing is associated with increases of over 140%. Cross-validated R2s based on the FIR models show that NORDIC is not measurably distorting the temporal structure of the data under this approach and is the best predictor of non-denoised time courses. The results demonstrate that analyzing 1 run of data after NORDIC produces results equivalent to using 2 to 3 original runs and that NORDIC performs equally well across a diverse array of functional imaging protocols. Significance Statement: For functional neuroimaging, the increasing availability of higher field strengths and ever higher spatiotemporal resolutions has led to concomitant increase in concerns about the deleterious effects of thermal noise. Historically this noise source was suppressed using methods that reduce spatial precision such as image blurring or averaging over a large number of trials or sessions, which necessitates large data collection efforts. Here, we critically evaluate the performance of a recently developed reconstruction method, termed NORDIC, which suppresses thermal noise. Across datasets varying in field strength, voxel sizes, sampling rates, and task designs, NORDIC produces substantial gains in data quality. Both conventional t-statistics derived from general linear models and coefficients of determination for predicting unseen data are improved. These gains match or even exceed those associated with 1 voxel Full Width Half Max image smoothing, however, even such small amounts of smoothing are associated with a 52% reduction in estimates of spatial precision, whereas the measurable difference in spatial precision is less than 4% following NORDIC.

Improved Quantitative Parameter Estimation for Prostate T2 Relaxometry using Convolutional Neural Networks.

This work seeks to evaluate multiple methods for quantitative parameter estimation from standard T2 mapping acquisitions in the prostate. The T2 estimation performance of methods based on neural networks (NN) was quantitatively compared to that of conventional curve fitting techniques. Large physics-based synthetic datasets simulating T2 mapping acquisitions were generated for training NNs and for quantitative performance comparisons. Ten combinations of different NN architectures, training strategies, and training corpora were implemented and compared with four different curve fitting strategies. All methods were compared quantitatively using synthetic data with known ground truth, and further compared on in vivo test data, with and without noise augmentation, to evaluate feasibility and noise robustness. In the evaluation on synthetic data, a convolutional neural network (CNN), trained in a supervised fashion using synthetic data generated from naturalistic images, showed the highest overall accuracy and precision amongst all the methods. On in vivo data, this best-performing method produced low-noise T2 maps and showed the least deterioration with increasing input noise levels. This study showed that a CNN, trained with synthetic data in a supervised manner, may provide superior T2 estimation performance compared to conventional curve fitting, especially in low signal-to-noise regions.

Signal intensity informed multi-coil encoding operator for physics-guided deep learning reconstruction of highly accelerated myocardial perfusion CMR.

PURPOSE: To develop a physics-guided deep learning (PG-DL) reconstruction strategy based on a signal intensity informed multi-coil (SIIM) encoding operator for highly-accelerated simultaneous multislice (SMS) myocardial perfusion cardiac MRI (CMR). METHODS: First-pass perfusion CMR acquires highly-accelerated images with dynamically varying signal intensity/SNR following the administration of a gadolinium-based contrast agent. Thus, using PG-DL reconstruction with a conventional multi-coil encoding operator leads to analogous signal intensity variations across different time-frames at the network output, creating difficulties in generalization for varying SNR levels. We propose to use a SIIM encoding operator to capture the signal intensity/SNR variations across time-frames in a reformulated encoding operator. This leads to a more uniform/flat contrast at the output of the PG-DL network, facilitating generalizability across time-frames. PG-DL reconstruction with the proposed SIIM encoding operator is compared to PG-DL with conventional encoding operator, split slice-GRAPPA, locally low-rank (LLR) regularized reconstruction, low-rank plus sparse (L + S) reconstruction, and regularized ROCK-SPIRiT. RESULTS: Results on highly accelerated free-breathing first pass myocardial perfusion CMR at three-fold SMS and four-fold in-plane acceleration show that the proposed method improves upon the reconstruction methods use for comparison. Substantial noise reduction is achieved compared to split slice-GRAPPA, and aliasing artifacts reduction compared to LLR regularized reconstruction, L + S reconstruction and PG-DL with conventional encoding. Furthermore, a qualitative reader study indicated that proposed method outperformed all methods. CONCLUSION: PG-DL reconstruction with the proposed SIIM encoding operator improves generalization across different time-frames /SNRs in highly accelerated perfusion CMR.

Cardiac phase-resolved late gadolinium enhancement imaging.

Late gadolinium enhancement (LGE) with cardiac magnetic resonance (CMR) imaging is the clinical reference for assessment of myocardial scar and focal fibrosis. However, current LGE techniques are confined to imaging of a single cardiac phase, which hampers assessment of scar motility and does not allow cross-comparison between multiple phases. In this work, we investigate a three step approach to obtain cardiac phase-resolved LGE images: (1) Acquisition of cardiac phase-resolved imaging data with varying T 1 weighting. (2) Generation of semi-quantitative T 1 * maps for each cardiac phase. (3) Synthetization of LGE contrast to obtain functional LGE images. The proposed method is evaluated in phantom imaging, six healthy subjects at 3T and 20 patients at 1.5T. Phantom imaging at 3T demonstrates consistent contrast throughout the cardiac cycle with a coefficient of variation of 2.55 ± 0.42%. In-vivo results show reliable LGE contrast with thorough suppression of the myocardial tissue is healthy subjects. The contrast between blood and myocardium showed moderate variation throughout the cardiac cycle in healthy subjects (coefficient of variation 18.2 ± 3.51%). Images were acquired at 40-60 ms and 80 ms temporal resolution, at 3T and 1.5, respectively. Functional LGE images acquired in patients with myocardial scar visualized scar tissue throughout the cardiac cycle, albeit at noticeably lower imaging resolution and noise resilience than the reference technique. The proposed technique bears the promise of integrating the advantages of phase-resolved CMR with LGE imaging, but further improvements in the acquisition quality are warranted for clinical use.

Unsupervised Deep Learning Methods for Biological Image Reconstruction and Enhancement: An overview from a signal processing perspective.

Recently, deep learning approaches have become the main research frontier for biological image reconstruction and enhancement problems thanks to their high performance, along with their ultra-fast inference times. However, due to the difficulty of obtaining matched reference data for supervised learning, there has been increasing interest in unsupervised learning approaches that do not need paired reference data. In particular, self-supervised learning and generative models have been successfully used for various biological imaging applications. In this paper, we overview these approaches from a coherent perspective in the context of classical inverse problems, and discuss their applications to biological imaging, including electron, fluorescence and deconvolution microscopy, optical diffraction tomography and functional neuroimaging.

Adiabatic spin-lock preparations enable robust in vivo cardiac T1ρ-mapping at 3T.

Magnetic Resonance Imaging (MRI) is the clinical gold standard for the assessment of myocardial viability but requires injection of exogenous gadolinium-based contrast agents. Recently, T1ρ-mapping has been proposed as a fully non-invasive alternative for imaging myocardial fibrosis without the need for contrast agent injection. However, its applicability at high fields is hindered by susceptibility to MRI system imperfections, such as inhomogeneities in the B0 and B1+ fields. In this work we propose a single breath-hold ECG-triggered single-shot bSSFP sequence to enable T1ρ-mapping in vivo at 3T. Adiabatic T1ρ preparations are evaluated to reduce B0 and B1+ sensitivity in comparison with conventional spin-lock (SL) modules. Numerical Bloch simulations were performed to identify optimal parameters for the adiabatic pulses. Experiments yield T1ρ values in the myocardium equal to 48.13±54.08 ms for the best adiabatic preparation and 16.01±20.75 ms for the reference non-adiabatic SL, with 26.91% against 89.74% relative difference in T1ρ values across two shimming conditions. Both phantom and in vivo measurements show increased myocardium/blood contrast and improved resilience against system imperfections compared to non-adiabatic T1ρ preparations, enabling the use at 3T. Clinical relevance- Adiabatically-prepared T1ρ-mapping sequences form a promising candidate for non-contrast evaluation of ischemic and non-ischemic cardiomyopathies at 3T.

Signal-Intensity Informed Multi-Coil MRI Encoding Operator for Improved Physics-Guided Deep Learning Reconstruction of Dynamic Contrast-Enhanced MRI.

Dynamic contrast enhanced (DCE) MRI acquires a series of images following the administration of a contrast agent, and plays an important clinical role in diagnosing various diseases. DCE MRI typically necessitates rapid imaging to provide sufficient spatio-temporal resolution and coverage. Conventional MRI acceleration techniques exhibit limited image quality at such high acceleration rates. Recently, deep learning (DL) methods have gained interest for improving highly-accelerated MRI. However, DCE MRI series show substantial variations in SNR and contrast across images. This hinders the quality and generalizability of DL methods, when applied across time frames. In this study, we propose signal intensity informed multi-coil MRI encoding operator for improved DL reconstruction of DCE MRI. The output of the corresponding inverse problem for this forward operator leads to more uniform contrast across time frames, since the proposed operator captures signal intensity variations across time frames while not altering the coil sensitivities. Our results in perfusion cardiac MRI show that high-quality images are reconstructed at very high acceleration rates, with substantial improvement over existing methods.

Myocardial Approximate Spin-lock Dispersion Mapping using a Simultaneous T2 and TRAFF2 Mapping at 3T MRI.

Ischemic heart disease (IHD) is one of the leading causes of death worldwide. Myocardial infarction (MI) represents a third of all IHD cases, and cardiac magnetic resonance imaging (MRI) is often used to assess its damage to myocardial viability. Late gadolinium enhancement (LGE) is the current gold standard, but the use of gadolinium-based agents limits the clinical applicability in some patients. Spin-lock (SL) dispersion has recently been proposed as a promising non-contrast biomarker for the assessment of MI. However, at 3T, the required range of SL preparations acquired at different amplitudes suffers from specific absorption rate (SAR) limitations and off-resonance artifacts. Relaxation Along a Fictitious Field (RAFF) is an alternative to SL preparations with lower SAR requirements, while still sampling relaxation in the rotating frame. In this study, a single breath-hold simultaneous TRAFF2 and T2 mapping sequence is proposed for SL dispersion mapping at 3T. Excellent reproducibility (coefficient of variations lower than 10%) was achieved in phantom experiments, indicating good intrascan repeatability. The average myocardial TRAFF2, T2, and SL dispersion obtained with the proposed sequence (68.0±10.7 ms, 44.0±4.0 ms, and 0.4±0.2 ×10-4 s2, respectively) were comparable to the reference methods (62.7±11.7 ms, 41.2±2.4 ms, and 0.3±0.2x 10-4s2, respectively). High visual map quality, free of B0 and B1+ related artifacts, for T2, TRAFF2, and SL dispersion maps were obtained in phantoms and in vivo, suggesting promise in clinical use at 3T. Clinical relevance - and imaging promises non-contrast assessment of scar and focal fibrosis in a single breath-hold using approximate spin-lock dispersion mapping.

Revisiting [Formula: see text]-wavelet compressed-sensing MRI in the era of deep learning.

Following their success in numerous imaging and computer vision applications, deep-learning (DL) techniques have emerged as one of the most prominent strategies for accelerated MRI reconstruction. These methods have been shown to outperform conventional regularized methods based on compressed sensing (CS). However, in most comparisons, CS is implemented with two or three hand-tuned parameters, while DL methods enjoy a plethora of advanced data science tools. In this work, we revisit [Formula: see text]-wavelet CS reconstruction using these modern tools. Using ideas such as algorithm unrolling and advanced optimization methods over large databases that DL algorithms utilize, along with conventional insights from wavelet representations and CS theory, we show that [Formula: see text]-wavelet CS can be fine-tuned to a level close to DL reconstruction for accelerated MRI. The optimized [Formula: see text]-wavelet CS method uses only 128 parameters compared to >500,000 for DL, employs a convex reconstruction at inference time, and performs within <1% of a DL approach that has been used in multiple studies in terms of quantitative quality metrics.

Multi-mask self-supervised learning for physics-guided neural networks in highly accelerated magnetic resonance imaging.

Self-supervised learning has shown great promise because of its ability to train deep learning (DL) magnetic resonance imaging (MRI) reconstruction methods without fully sampled data. Current self-supervised learning methods for physics-guided reconstruction networks split acquired undersampled data into two disjoint sets, where one is used for data consistency (DC) in the unrolled network, while the other is used to define the training loss. In this study, we propose an improved self-supervised learning strategy that more efficiently uses the acquired data to train a physics-guided reconstruction network without a database of fully sampled data. The proposed multi-mask self-supervised learning via data undersampling (SSDU) applies a holdout masking operation on the acquired measurements to split them into multiple pairs of disjoint sets for each training sample, while using one of these pairs for DC units and the other for defining loss, thereby more efficiently using the undersampled data. Multi-mask SSDU is applied on fully sampled 3D knee and prospectively undersampled 3D brain MRI datasets, for various acceleration rates and patterns, and compared with the parallel imaging method, CG-SENSE, and single-mask SSDU DL-MRI, as well as supervised DL-MRI when fully sampled data are available. The results on knee MRI show that the proposed multi-mask SSDU outperforms SSDU and performs as well as supervised DL-MRI. A clinical reader study further ranks the multi-mask SSDU higher than supervised DL-MRI in terms of signal-to-noise ratio and aliasing artifacts. Results on brain MRI show that multi-mask SSDU achieves better reconstruction quality compared with SSDU. The reader study demonstrates that multi-mask SSDU at R = 8 significantly improves reconstruction compared with single-mask SSDU at R = 8, as well as CG-SENSE at R = 2.