Feasibility of undersampled spiral trajectories in MREPT for fast conductivity imaging.

PURPOSE: To investigate spiral-based imaging including trajectories with undersampling as a fast and robust alternative for phase-based magnetic resonance electrical properties tomography (MREPT) techniques. METHODS: Spiral trajectories with various undersampling ratios were prescribed to acquire images from an experimental phantom and a healthy volunteer at 3T. The non-Cartesian acquisitions were reconstructed using SPIRiT, and conductivity maps were derived using phase-based cr-MREPT. The resulting maps were compared between different sampling trajectories. Additionally, a conductivity map was obtained using a Cartesian balanced SSFP acquisition from the volunteer to comparatively demonstrate the robustness of the proposed method. RESULTS: The phantom and volunteer results illustrate the benefits of the spiral acquisitions. Specifically, undersampled spiral acquisitions display improved robustness against field inhomogeneity artifacts and lowered SD values with shortened readout times. Furthermore, average of conductivity values measured for the cerebrospinal fluid with the spiral acquisitions were 1.703 S/m, indicating a close agreement with the theoretical values of 1.794 S/m. CONCLUSION: A spiral-based acquisition framework for conductivity imaging with and without undersampling is presented. Overall, spiral-based acquisitions improved robustness against field inhomogeneity artifacts, while achieving whole head coverage with multiple averages in less than a minute.

Multidimensional compressed sensing to advance 23 Na multi-quantum coherences MRI.

PURPOSE: Sodium (23 Na) multi-quantum coherences (MQC) MRI was accelerated using three-dimensional (3D) and a dedicated five-dimensional (5D) compressed sensing (CS) framework for simultaneous Cartesian single (SQ) and triple quantum (TQ) sodium imaging of in vivo human brain at 3.0 and 7.0 T. THEORY AND METHODS: 3D 23 Na MQC MRI requires multi-echo paired with phase-cycling and exhibits thus a multidimensional space. A joint reconstruction framework to exploit the sparsity in all imaging dimensions by extending the conventional 3D CS framework to 5D was developed. 3D MQC images of simulated brain, phantom and healthy brain volunteers obtained from 3.0 T and 7.0 T were retrospectively and prospectively undersampled. Performance of the CS models were analyzed by means of structural similarity index (SSIM), root mean squared error (RMSE), signal-to-noise ratio (SNR) and signal quantification of tissue sodium concentration and TQ/SQ ratio. RESULTS: It was shown that an acceleration of three-fold, leading to less than 2 × 10 $$ 2\times 10 $$ min of scan time with a resolution of 8 × 8 × 20 mm 3 $$ 8\times 8\times 20\;{\mathrm{mm}}^3 $$ at 3.0 T, are possible. 5D CS improved SSIM by 3%, 5%, 1% and reduced RMSE by 50%, 30%, 8% for in vivo SQ, TQ, and TQ/SQ ratio maps, respectively. Furthermore, for the first time prospective undersampling enabled unprecedented high resolution from 8 × 8 × 20 mm 3 $$ 8\times 8\times 20\;{\mathrm{mm}}^3 $$ to 6 × 6 × 10 mm 3 $$ 6\times 6\times 10\;{\mathrm{mm}}^3 $$ MQC images of in vivo human brain at 7.0 T without extending acquisition time. CONCLUSION: 5D CS proved to allow up to three-fold acceleration retrospectively on 3.0 T data. 2-fold acceleration was demonstrated prospectively at 7.0 T to reach higher spatial resolution of 23 Na MQC MRI.

Sodium triple quantum MR signal extraction using a single-pulse sequence with single quantum time efficiency.

PURPOSE: Sodium triple quantum (TQ) signal has been shown to be a valuable biomarker for cell viability. Despite its clinical potential, application of Sodium TQ signal is hindered by complex pulse sequences with long scan times. This study proposes a method to approximate the TQ signal using a single excitation pulse without phase cycling. METHODS: The proposed method is based on a single excitation pulse and a comparison of the free induction decay (FID) with the integral of the FID combined with a shifting reconstruction window. The TQ signal is calculated from this FID only. As a proof of concept, the method was also combined with a multi-echo UTE imaging sequence on a 9.4 T preclinical MRI scanner for the possibility of fast TQ MRI. RESULTS: The extracted Sodium TQ signals of single-pulse and spin echo FIDs were in close agreement with theory and TQ measurement by traditional three-pulse sequence (TQ time proportional phase increment [TQTPPI)]. For 2%, 4%, and 6% agar samples, the absolute deviations of the maximum TQ signals between SE and theoretical (time proportional phase increment TQTPPI) TQ signals were less than 1.2% (2.4%), and relative deviations were less than 4.6% (6.8%). The impact of multi-compartment systems and noise on the accuracy of the TQ signal was small for simulated data. The systematic error was <3.4% for a single quantum (SQ) SNR of 5 and at maximum <2.5% for a multi-compartment system. The method also showed the potential of fast in vivo SQ and TQ imaging. CONCLUSION: Simultaneous SQ and TQ MRI using only a single-pulse sequence and SQ time efficiency has been demonstrated. This may leverage the full potential of the Sodium TQ signal in clinical applications.

Low-rank reconstruction for simultaneous double half-echo 23Na and undersampled 23Na multi-quantum coherences MRI.

PURPOSE: To develop a new sequence to simultaneously acquire Cartesian sodium (23Na) MRI and accelerated Cartesian single (SQ) and triple quantum (TQ) sodium MRI of in vivo human brain at 7 T by leveraging two dedicated low-rank reconstruction frameworks. THEORY AND METHODS: The Double Half-Echo technique enables short echo time Cartesian 23Na MRI and acquires two k-space halves, reconstructed by a low-rank coupling constraint. Additionally, three-dimensional (3D) 23Na Multi-Quantum Coherences (MQC) MRI requires multi-echo sampling paired with phase-cycling, exhibiting a redundant multidimensional space. Simultaneous Autocalibrating and k-Space Estimation (SAKE) were used to reconstruct highly undersampled 23Na MQC MRI. Reconstruction performance was assessed against five-dimensional (5D) CS, evaluating structural similarity index (SSIM), root mean squared error (RMSE), signal-to-noise ratio (SNR), and quantification of tissue sodium concentration and TQ/SQ ratio in silico, in vitro, and in vivo. RESULTS: The proposed sequence enabled the simultaneous acquisition of fully sampled 23Na MRI while leveraging prospective undersampling for 23Na MQC MRI. SAKE improved TQ image reconstruction regarding SSIM by 6% and reduced RMSE by 35% compared to 5D CS in vivo. Thanks to prospective undersampling, the spatial resolution of 23Na MQC MRI was enhanced from 8 × 8 × 15 $$ 8\times 8\times 15 $$ mm3 to 8 × 8 × 8 $$ 8\times 8\times 8 $$ mm3 while reducing acquisition time from 2 × 31 $$ 2\times 31 $$ min to 2 × 23 $$ 2\times 23 $$ min. CONCLUSION: The proposed sequence, coupled with low-rank reconstructions, provides an efficient framework for comprehensive whole-brain sodium MRI, combining TSC, T2*, and TQ/SQ ratio estimations. Additionally, low-rank matrix completion enables the reconstruction of highly undersampled 23Na MQC MRI, allowing for accelerated acquisition or enhanced spatial resolution.

Improved reproducibility for myocardial ASL: Impact of physiological and acquisition parameters.

PURPOSE: To investigate and mitigate the influence of physiological and acquisition-related parameters on myocardial blood flow (MBF) measurements obtained with myocardial Arterial Spin Labeling (myoASL). METHODS: A Flow-sensitive Alternating Inversion Recovery (FAIR) myoASL sequence with bSSFP and spoiled GRE (spGRE) readout is investigated for MBF quantification. Bloch-equation simulations and phantom experiments were performed to evaluate how variations in acquisition flip angle (FA), acquisition matrix size (AMS), heart rate (HR) and blood T 1 $$ {\mathrm{T}}_1 $$ relaxation time ( T 1 , B $$ {\mathrm{T}}_{1,B} $$ ) affect quantification of myoASL-MBF. In vivo myoASL-images were acquired in nine healthy subjects. A corrected MBF quantification approach was proposed based on subject-specific T 1 , B $$ {\mathrm{T}}_{1,B} $$ values and, for spGRE imaging, subtracting an additional saturation-prepared baseline from the original baseline signal. RESULTS: Simulated and phantom experiments showed a strong dependence on AMS and FA ( R 2 $$ {R}^2 $$ >0.73), which was eliminated in simulations and alleviated in phantom experiments using the proposed saturation-baseline correction in spGRE. Only a very mild HR dependence ( R 2 $$ {R}^2 $$ >0.59) was observed which was reduced when calculating MBF with individual T 1 , B $$ {\mathrm{T}}_{1,B} $$ . For corrected spGRE, in vivo mean global spGRE-MBF ranged from 0.54 to 2.59 mL/g/min and was in agreement with previously reported values. Compared to uncorrected spGRE, the intra-subject variability within a measurement (0.60 mL/g/min), between measurements (0.45 mL/g/min), as well as the inter-subject variability (1.29 mL/g/min) were improved by up to 40% and were comparable with conventional bSSFP. CONCLUSION: Our results show that physiological and acquisition-related factors can lead to spurious changes in myoASL-MBF if not accounted for. Using individual T 1 , B $$ {\mathrm{T}}_{1,B} $$ and a saturation-baseline can reduce these variations in spGRE and improve reproducibility of FAIR-myoASL against acquisition parameters.

Comparison of triple quantum (TQ) TPPI and inversion recovery TQ TPPI pulse sequences at 9.4 and 21.1 T.

PURPOSE: Both sodium T1 triple quantum (TQ) signal and T1 relaxation pathways have a unique sensitivity to the sodium molecular environment. In this study an inversion recovery time proportional phase increment (IRTQTPPI) pulse sequence was investigated for simultaneous and reliable quantification of sodium TQ signal and bi-exponential T1 relaxation times. METHODS: The IRTQTPPI sequence combines inversion recovery TQ filtering and time proportional phase increment. The reliable and reproducible results were achieved by the pulse sequence optimized in three ways: (1) optimization of the nonlinear fit for the determination of both T1-TQ signal and T1 relaxation times; (2) suppression of unwanted signals by assessment of four different phase cycles; (3) nonlinear sampling during evolution time for optimal scan time without any compromises in fit accuracy. The relaxation times T1 and T2 and the TQ signals from IRTQTPPI and TQTPPI were compared between 9.4 and 21.1 T. The motional environment of the sodium nuclei was evaluated by calculation of correlation times and nuclear quadrupole interaction strengths. RESULTS: Reliable measurements of the T1-TQ signals and T1 bi-exponential relaxation times were demonstrated. The fit parameters for all four phase cycles were in good agreement with one another, with a negligible influence of unwanted signals. The agar samples yielded normalized T1-TQ signals from 3% to 16% relative to single quantum (SQ) signals at magnetic fields of both 9.4 and 21.1 T. In comparison, the normalized T2-TQ signal was in the range 15%-35%. The TQ/SQ signal ratio was decreased at 21.1 T as compared with 9.4 T for both T1 and T2 relaxation pathways. The bi-exponential T1 relaxation time separation ranged from 15 to 18 ms at 9.4 T and 15 to 21 ms at 21.1 T. The T2 relaxation time separation was larger, ranging from 28 to 35 ms at 9.4 T and 37 to 40 ms at 21.1 T. CONCLUSION: The IRTQTPPI sequence, while providing a less intensive TQ signal than TQTPPI, allows a simultaneous and reliable quantification of both the T1-TQ signal and T1 relaxation times. The unique sensitivities of the T1 and T2 relaxation pathways to different types of molecular motion provide a deeper understanding of the sodium MR environment.

Linear multi-scale modeling of diffusion MRI data: A framework for characterization of oriented structures across length scales.

Diffusion-weighted magnetic resonance imaging (DW-MRI) has evolved to provide increasingly sophisticated investigations of the human brain's structural connectome in vivo. Restriction spectrum imaging (RSI) is a method that reconstructs the orientation distribution of diffusion within tissues over a range of length scales. In its original formulation, RSI represented the signal as consisting of a spectrum of Gaussian diffusion response functions. Recent technological advances have enabled the use of ultra-high b-values on human MRI scanners, providing higher sensitivity to intracellular water diffusion in the living human brain. To capture the complex diffusion time dependence of the signal within restricted water compartments, we expand upon the RSI approach to represent restricted water compartments with non-Gaussian response functions, in an extended analysis framework called linear multi-scale modeling (LMM). The LMM approach is designed to resolve length scale and orientation-specific information with greater specificity to tissue microstructure in the restricted and hindered compartments, while retaining the advantages of the RSI approach in its implementation as a linear inverse problem. Using multi-shell, multi-diffusion time DW-MRI data acquired with a state-of-the-art 3 T MRI scanner equipped with 300 mT/m gradients, we demonstrate the ability of the LMM approach to distinguish different anatomical structures in the human brain and the potential to advance mapping of the human connectome through joint estimation of the fiber orientation distributions and compartment size characteristics.

Dynamic mode decomposition of dynamic MRI for assessment of pulmonary ventilation and perfusion.

PURPOSE: To introduce dynamic mode decomposition (DMD) as a robust alternative for the assessment of pulmonary functional information from dynamic non-contrast-enhanced acquisitions. METHODS: Pulmonary fractional ventilation and normalized perfusion maps were obtained using DMD from simulated phantoms as well as in vivo dynamic acquisitions of healthy volunteers at 1.5T. The performance of DMD was compared with conventional Fourier decomposition (FD) and matrix pencil (MP) methods in estimating functional map values. The proposed method was evaluated based on estimated signal amplitude in functional maps across varying number of measurements. RESULTS: Quantitative assessments performed on phantoms and in vivo measurements indicate that DMD is capable of successfully obtaining pulmonary functional maps. Specifically, compared to FD and MP methods, DMD is able to reduce variations in estimated amplitudes across different number of measurements. This improvement is evident in the fractional ventilation and normalized perfusion maps obtain from phantom simulations with frequency variations and noise, as well as in the maps obtained from in vivo measurements. CONCLUSIONS: A robust method for accurately estimating pulmonary ventilation and perfusion related signal changes in dynamic acquisitions is presented. The proposed method uses DMD to obtain functional maps reliably, while reducing amplitude variations caused by differences in number of measurements.

Prediction of experimental cardiac magnetostimulation thresholds using pig-specific body models.

PURPOSE: Modern high-amplitude gradient systems can be limited by the International Electrotechnical Commission 60601-2-33 cardiac stimulation (CS) limit, which was set in a conservative manner based on electrode experiments and E-field simulations in uniform ellipsoidal body models. Here, we show that coupled electromagnetic-electrophysiological modeling in detailed body and heart models can predict CS thresholds, suggesting that such modeling might lead to more detailed threshold estimates in humans. Specifically, we compare measured and predicted CS thresholds in eight pigs. METHODS: We created individualized porcine body models using MRI (Dixon for the whole body, CINE for the heart) that replicate the anatomy and posture of the animals used in our previous experimental CS study. We model the electric fields induced along cardiac Purkinje and ventricular muscle fibers and predict the electrophysiological response of these fibers, yielding CS threshold predictions in absolute units for each animal. Additionally, we assess the total modeling uncertainty through a variability analysis of the 25 main model parameters. RESULTS: Predicted and experimental CS thresholds agree within 19% on average (normalized RMS error), which is smaller than the 27% modeling uncertainty. No significant difference was found between the modeling predictions and experiments (p < 0.05, paired t-test). CONCLUSION: Predicted thresholds matched the experimental data within the modeling uncertainty, supporting the model validity. We believe that our modeling approach can be applied to study CS thresholds in humans for various gradient coils, body shapes/postures, and waveforms, which is difficult to do experimentally.

Functional lung imaging of 2-year-old children after congenital diaphragmatic hernia repair using dynamic mode decomposition MRI.

OBJECTIVES: To investigate the feasibility of non-contrast-enhanced functional lung imaging in 2-year-old children after congenital diaphragmatic hernia (CDH) repair. METHODS: Fifteen patients after CDH repair were examined using non-contrast-enhanced dynamic magnetic resonance imaging (MRI). For imaging two protocols were used during free-breathing: Protocol A with high temporal resolution and Protocol B with high spatial resolution. The dynamic images were then analysed through a recently developed post-processing method called dynamic mode decomposition (DMD) to obtain ventilation and perfusion maps. The ventilation ratios (VRatio) and perfusion ratios (QRatio) of ipsilateral to contralateral lung were compared to evaluate functional differences. Lastly, DMD MRI-based perfusion results were compared with perfusion parameters obtained using dynamic contrast-enhanced (DCE) MRI to assess agreement between methods. RESULTS: Both imaging protocols successfully generated pulmonary ventilation (V) and perfusion (Q) maps in all patients. Overall, the VRatio and QRatio values were 0.84 ± 0.19 and 0.70 ± 0.24 for Protocol A, and 0.88 ± 0.18 and 0.72 ± 0.23 for Protocol B, indicating reduced ventilation ([Formula: see text]) and perfusion ([Formula: see text]) on the ipsilateral side. Moreover, there is a very strong positive correlation ([Formula: see text]) and close agreement between DMD MRI-based perfusion values and DCE MRI-based perfusion parameters. CONCLUSIONS: DMD MRI can obtain pulmonary functional information in 2-year-old CDH patients. The results obtained with DMD MRI correlate with DCE MRI, without the need for ionising radiation or exposure to contrast agents. While further studies with larger cohorts are warranted, DMD MRI is a promising option for functional lung imaging in CDH patients. CLINICAL RELEVANCE STATEMENT: We demonstrate that pulmonary ventilation and perfusion information can be obtained in 2-year-old patients after CDH repair, without the need for ionising radiation or contrast agents by utilising non-contrast-enhanced MRI acquisitions together with dynamic mode decomposition analysis. KEY POINTS: • Non-contrast-enhanced functional MR imaging is a promising option for functional lung imaging in 2-year-old children after congenital diaphragmatic hernia. • DMD MRI can generate pulmonary ventilation and perfusion maps from free-breathing dynamic acquisitions without the need for ionising radiation or contrast agents. • Lung perfusion parameters obtained with DMD MRI correlate with perfusion parameters obtained using dynamic contrast-enhanced MRI.

VirtMRI: A Tool for Teaching MRI.

Magnetic resonance image formation is not trivial and remains a difficult subject for teaching. Therefore, we saw an urgent need to facilitate teaching by developing a practical and easily accessible MR image generator. Due to the increasing interest in X-nuclei MRI, sodium image generation is also offered. The tool is implemented as a web application that is compatible with all standard desktop browsers and is open source. The user interface focuses on the parameters needed for the creation and display of the resulting images. Available MR sequences range from the standard Spin Echo and Inversion Recovery over steady-state to conventional sodium and more advanced single and triple quantum sequences. Additionally, the user interface has parameters to alter the resolution, the noise, and the k-space sampling. Our software is free to use and specifically suited for teaching purposes.

Phase-cycled balanced SSFP imaging for non-contrast-enhanced functional lung imaging.

PURPOSE: To introduce phase-cycled balanced SSFP (bSSFP) acquisition as an alternative in Fourier decomposition MRI for improved robustness against field inhomogeneities. METHODS: Series 2D dynamic lung images were acquired in 5 healthy volunteers at 1.5 T and 3 T using bSSFP sequence with multiple RF phase increments and compared with conventional single RF phase increment acquisitions. The approach was evaluated based on functional map homogeneity analysis, while ensuring image and functional map quality by means of SNR and contrast-to-noise ratio analyses. RESULTS: At both field strengths, functional maps obtained with phase-cycled acquisitions displayed improved robustness against local signal losses compared with single-phase acquisitions. The coefficient of variation (mean ± SD, across volunteers) measured in the ventilation maps resulted in 29.7 ± 2.6 at 1.5 T and 37.5 ± 3.1 at 3 T for phase-cycled acquisitions, compared with 39.9 ± 5.2 at 1.5 T and 49.5 ± 3.7 at 3 T for single-phase acquisitions, indicating a significant improvement ( p<0.05$$ p<0.05 $$ ) in ventilation map homogeneity. CONCLUSIONS: Phase-cycled bSSFP acquisitions improve robustness against field inhomogeneity artifacts and significantly improve ventilation map homogeneity at both field strengths. As such, phase-cycled bSSFP may serve as a robust alternative in lung function assessments.

An anthropomorphic pelvis phantom for MR-guided prostate interventions.

PURPOSE: To design and manufacture a pelvis phantom for magnetic resonance (MR)-guided prostate interventions, such as MRGB (MR-guided biopsy) or brachytherapy seed placement. METHODS: The phantom was designed to mimic the human pelvis incorporating bones, bladder, prostate with four lesions, urethra, arteries, veins, and six lymph nodes embedded in ballistic gelatin. A hollow rectum enables transrectal access to the prostate. To demonstrate the feasibility of the phantom for minimal invasive MRI-guided interventions, a targeted inbore MRGB was performed. The needle probe was rectally inserted and guided using an MRI-compatible remote controlled manipulator (RCM). RESULTS: The presented pelvis phantom has realistic imaging properties for MR imaging (MRI), computed tomography (CT) and ultrasound (US). In the targeted inbore MRGB, a prostate lesion was successfully hit with an accuracy of 3.5 mm. The experiment demonstrates that the limited size of the rectum represents a realistic impairment for needle placements. CONCLUSION: The phantom provides a valuable platform for evaluating the performance of MRGB systems. Interventionalists can use the phantom to learn how to deal with challenging situations, without risking harm to patients.

Measurement of magnetostimulation thresholds in the porcine heart.

PURPOSE: Powerful MRI gradient systems can surpass the International Electrotechnical Commission (IEC) 60601-2-33 limit for cardiac stimulation (CS), which was determined by simple electromagnetic simulations and electrode stimulation experiments. Only a few canine studies measured magnetically induced CS thresholds in vivo and extrapolating them to human safety limits can be challenging. METHODS: We measured cardiac magnetostimulation thresholds in 10 healthy, anesthetized pigs using capacitors discharged into a flat spiral coil to produce damped sinusoidal waveforms with effective stimulus duration ts,eff  = 0.45 ms. Electrocardiography (ECG), blood pressure, and peripheral oximetry signals were recorded to determine threshold coil currents yielding cardiac capture. Dixon and CINE MR volumes from each animal were segmented to generate porcine-specific electromagnetic models to calculate dB/dt and E-field values in the porcine heart at threshold. For comparison, we also simulated maximum dB/dt and E-field values created by three MRI gradient systems in the heart of a human body model. RESULTS: The average dB/dt threshold estimated in the porcine heart was 1.66 ± 0.23 kT/s, which is 11-fold greater than the IEC dB/dt limit at ts,eff  = 0.45 ms, and 31-fold greater than the maximum value created by the investigated MRI gradients in the human heart. The average E-field threshold estimated in the porcine heart was 92.9 ± 13.5 V/m, which is 6-fold greater than the IEC E-field limit at ts,eff  = 0.45 ms and 37-fold greater than the maximum gradient-induced E-field in the human heart. CONCLUSION: This first measurement of cardiac magnetostimulation thresholds in pigs indicates that the IEC cardiac safety limit is conservative for the investigated stimulus duration (ts,eff  = 0.45 ms).

Acceleration of Magnetic Resonance Fingerprinting Reconstruction Using Denoising and Self-Attention Pyramidal Convolutional Neural Network.

Magnetic resonance fingerprinting (MRF) based on echo-planar imaging (EPI) enables whole-brain imaging to rapidly obtain T1 and T2* relaxation time maps. Reconstructing parametric maps from the MRF scanned baselines by the inner-product method is computationally expensive. We aimed to accelerate the reconstruction of parametric maps for MRF-EPI by using a deep learning model. The proposed approach uses a two-stage model that first eliminates noise and then regresses the parametric maps. Parametric maps obtained by dictionary matching were used as a reference and compared with the prediction results of the two-stage model. MRF-EPI scans were collected from 32 subjects. The signal-to-noise ratio increased significantly after the noise removal by the denoising model. For prediction with scans in the testing dataset, the mean absolute percentage errors between the standard and the final two-stage model were 3.1%, 3.2%, and 1.9% for T1, and 2.6%, 2.3%, and 2.8% for T2* in gray matter, white matter, and lesion locations, respectively. Our proposed two-stage deep learning model can effectively remove noise and accurately reconstruct MRF-EPI parametric maps, increasing the speed of reconstruction and reducing the storage space required by dictionaries.

Investigating cardiac stimulation limits of MRI gradient coils using electromagnetic and electrophysiological simulations in human and canine body models.

PURPOSE: Cardiac stimulation (CS) limits to gradient coil switching speed are difficult to measure in humans; instead, current regulatory guidelines (IEC 60601-2-33) are based on animal experiments and electric field-to-dB/dt conversion factors computed for a simple, homogeneous body model. We propose improvement to this methodology by using more detailed CS modeling based on realistic body models and electrophysiological models of excitable cardiac fibers. METHODS: We compute electric fields induced by a solenoid, coplanar loops, and a commercial gradient coil in two human body models and a canine model. The canine simulations mimic previously published experiments. We generate realistic fiber topologies for the cardiac Purkinje and ventricular muscle fiber networks using rule-based algorithms, and evaluate CS thresholds using validated electrodynamic models of these fibers. RESULTS: We were able to reproduce the average measured canine CS thresholds within 5%. In all simulations, the Purkinje fibers were stimulated before the ventricular fibers, and therefore set the effective CS threshold. For the investigated gradient coil, simulated CS thresholds for the x-, y-, and z-axis were at least one order of magnitude greater than the International Electrotechnical Commission limit. CONCLUSION: We demonstrate an approach to simulate gradient-induced CS using a combination of electromagnetic and electrophysiological modeling. Pending additional validation, these simulations could guide the assessment of CS limits to MRI gradient coil switching speed. Such an approach may lead to less conservative, but still safe, operation limits, enabling the use of the maximum gradient amplitude versus slew rate parameter space of recent, powerful gradient systems.

Free-breathing simultaneous T1 , T2 , and T2∗ quantification in the myocardium.

PURPOSE: To implement a free-breathing sequence for simultaneous quantification of T1 , T2 , and T2∗ for comprehensive tissue characterization of the myocardium in a single scan using a multi-gradient-echo readout with saturation and T2 preparation pulses. METHODS: In the proposed Saturation And T2 -prepared Relaxometry with Navigator-gating (SATURN) technique, a series of multi-gradient-echo (GRE) images with different magnetization preparations was acquired during free breathing. A total of 35 images were acquired in 26.5 ± 14.9 seconds using multiple saturation times and T2 preparation durations and with imaging at 5 echo times. Bloch simulations and phantom experiments were used to validate a 5-parameter fit model for accurate relaxometry. Free-breathing simultaneous T1 , T2 , and T2∗ measurements were performed in 10 healthy volunteers and 2 patients using SATURN at 3T and quantitatively compared to conventional single-parameter methods such as SASHA for T1 , T2 -prepared bSSFP, and multi-GRE for T2∗ . RESULTS: Simulations confirmed accurate fitting with the 5-parameter model. Phantom measurements showed good agreement with the reference methods in the relevant range for in vivo measurements. Compared to single-parameter methods comparable accuracy was achieved. SATURN produced in vivo parameter maps that were visually comparable to single-parameter methods. No significant difference between T1 , T2 , and T2∗ times acquired with SATURN and single-parameter methods was shown in quantitative measurements (SATURN T1=1573±86ms , T2=33.2±3.6ms , T2∗=25.3±6.1ms ; conventional methods: T1=1544±107ms , T2=33.2±3.6ms , T2∗=23.8±5.5ms ; P>.2 ) CONCLUSION: SATURN enables simultaneous quantification of T1 , T2 , and T2∗ in the myocardium for comprehensive tissue characterization with co-registered maps, in a single scan with good agreement to single-parameter methods.

Accelerated white matter lesion analysis based on simultaneous T1 and T2∗ quantification using magnetic resonance fingerprinting and deep learning.

PURPOSE: To develop an accelerated postprocessing pipeline for reproducible and efficient assessment of white matter lesions using quantitative magnetic resonance fingerprinting (MRF) and deep learning. METHODS: MRF using echo-planar imaging (EPI) scans with varying repetition and echo times were acquired for whole brain quantification of T1 and T2∗ in 50 subjects with multiple sclerosis (MS) and 10 healthy volunteers along 2 centers. MRF T1 and T2∗ parametric maps were distortion corrected and denoised. A CNN was trained to reconstruct the T1 and T2∗ parametric maps, and the WM and GM probability maps. RESULTS: Deep learning-based postprocessing reduced reconstruction and image processing times from hours to a few seconds while maintaining high accuracy, reliability, and precision. Mean absolute error performed the best for T1 (deviations 5.6%) and the logarithmic hyperbolic cosinus loss the best for T2∗ (deviations 6.0%). CONCLUSIONS: MRF is a fast and robust tool for quantitative T1 and T2∗ mapping. Its long reconstruction and several postprocessing steps can be facilitated and accelerated using deep learning.

23 Na MRI in ischemic stroke: Acquisition time reduction using postprocessing with convolutional neural networks.

Quantitative 23 Na magnetic resonance imaging (MRI) provides tissue sodium concentration (TSC), which is connected to cell viability and vitality. Long acquisition times are one of the most challenging aspects for its clinical establishment. K-space undersampling is an approach for acquisition time reduction, but generates noise and artifacts. The use of convolutional neural networks (CNNs) is increasing in medical imaging and they are a useful tool for MRI postprocessing. The aim of this study is 23 Na MRI acquisition time reduction by k-space undersampling. CNNs were applied to reduce the resulting noise and artifacts. A retrospective analysis from a prospective study was conducted including image datasets from 46 patients (aged 72 ± 13 years; 25 women, 21 men) with ischemic stroke; the 23 Na MRI acquisition time was 10 min. The reconstructions were performed with full dataset (FI) and with a simulated dataset an image that was acquired in 2.5 min (RI). Eight different CNNs with either U-Net-based or ResNet-based architectures were implemented with RI as input and FI as label, using batch normalization and the number of filters as varying parameters. Training was performed with 9500 samples and testing included 400 samples. CNN outputs were evaluated based on signal-to-noise ratio (SNR) and structural similarity (SSIM). After quantification, TSC error was calculated. The image quality was subjectively rated by three neuroradiologists. Statistical significance was evaluated by Student's t-test. The average SNR was 21.72 ± 2.75 (FI) and 10.16 ± 0.96 (RI). U-Nets increased the SNR of RI to 43.99 and therefore performed better than ResNet. SSIM of RI to FI was improved by three CNNs to 0.91 ± 0.03. CNNs reduced TSC error by up to 15%. The subjective rating of CNN-generated images showed significantly better results than the subjective image rating of RI. The acquisition time of 23 Na MRI can be reduced by 75% due to postprocessing with a CNN on highly undersampled data.

Characterization of chronic active multiple sclerosis lesions with sodium (23 Na) magnetic resonance imaging-preliminary observations.

BACKGROUND AND PURPOSE: There has been an increasing interest in chronic active multiple sclerosis (MS) lesions as a new magnetic resonance imaging (MRI) marker of disease progression. Chronic active lesions are characterized by progressive tissue matrix damage, axonal loss and chronic inflammation. Sodium (23 Na) MRI provides a biochemical marker of cell integrity and tissue viability in a quantitative manner. The aim of this study was to investigate with 23 Na MRI tissue abnormalities in chronic active lesions as indicators of tissue destruction. METHODS: To identify chronic active lesions, two 3D magnetization-prepared rapid acquisition gradient-echo datasets obtained 12 months apart were processed using the voxel-guided morphometry algorithm. Cross-sectional 23 Na MRI was performed during the 12-month follow-up period. Total sodium concentration was calculated in chronic active lesions compared to shrinking, chronic stable and acute contrast-enhancing lesions. RESULTS: Overall, 70 MS lesions (21 chronic active, 10 shrinking, 29 chronic stable lesions, 10 acute contrast-enhancing lesions) in 12 patients were included. Total sodium concentration in chronic active lesions (49.57 ± 8.47 mM) was significantly higher than in shrinking (42.16 ± 3.9 mM; p = 0.03) and chronic stable lesions (39.92 ± 4.82 mM; p < 0.001). Chronic active lesions showed similar sodium values compared to acute contrast-enhancing lesions (48.06 ± 6.65 mM; p = 0.97). No differences between shrinking and chronic stable lesions were observed (p = 0.89). CONCLUSION: High sodium values in chronic active MS lesions may be an indicator of ongoing inflammation and tissue damage.