Fast and silent MRI using nonlinear gradient fields at the ultrasonic gradient switching frequency of 20 kHz with a Point Spread Function framework reconstruction.

PURPOSE: To demonstrate the feasibility of using a nonlinear gradient field for spatial encoding at the ultrasonic switching frequency of 20 kHz and present a framework to reconstruct data acquired in this way. METHODS: Nonlinear encoding at 20 kHz was realized by using a single-axis silent gradient insert for imaging in the periphery, that, is the nonlinear region, of the gradient field. The gradient insert induces a rapidly oscillating gradient field in the phase-encode direction, which enables nonlinear encoding when combined with a Cartesian readout from the linear whole-body gradients. Data from a 2D gradient echo sequence were reconstructed using a point spread function (PSF) framework. Accelerated scans were also simulated via retrospective undersampling (R = 1 to R = 8) to determine the effectiveness of the PSF-framework for accelerated imaging. RESULTS: Using a nonlinear gradient field switched at 20 kHz and the PSF-framework resulted in images of comparable quality to images from conventional Cartesian linear encoding. At increased acceleration factors (R ≤ 8), the PSF-framework outperformed linear SENSE reconstructions by improved controlling of aliasing artifacts. CONCLUSION: Using the PSF-framework, images of comparable quality to conventional SENSE reconstructions are possible via combining traditional linear and ultrasonic oscillating nonlinear encoding fields. Using nonlinear gradient fields relaxes the demand for strictly linear gradient fields, enabling much higher slew rates with a reduced risk of peripheral nerve stimulation or cardiac stimulation, which could aid in extension to ultrasonic whole-body MRI. The lack of aliasing artifacts also highlights the potential of accelerated imaging using the PSF-framework.

Improving the lesion appearance on FLAIR images synthetized from quantitative MRI: a fast, hybrid approach.

OBJECTIVE: The image quality of synthetized FLAIR (fluid attenuated inversion recovery) images is generally inferior to its conventional counterpart, especially regarding the lesion contrast mismatch. This work aimed to improve the lesion appearance through a hybrid methodology. MATERIALS AND METHODS: We combined a full brain 5-min MR-STAT acquisition followed by FLAIR synthetization step with an ultra-under sampled conventional FLAIR sequence and performed the retrospective and prospective analysis of the proposed method on the patient datasets and a healthy volunteer. RESULTS: All performance metrics of the proposed hybrid FLAIR images on patient datasets were significantly higher than those of the physics-based FLAIR images (p < 0.005), and comparable to those of conventional FLAIR images. The small difference between prospective and retrospective analysis on a healthy volunteer demonstrated the validity of the retrospective analysis of the hybrid method as presented for the patient datasets. DISCUSSION: The proposed hybrid FLAIR achieved an improved lesion appearance in the clinical cases with neurological diseases compared to the physics-based FLAIR images, Future prospective work on patient data will address the validation of the method from a diagnostic perspective by radiological inspection of the new images over a larger patient cohort.

Towards retrospective motion correction and reconstruction for clinical 3D brain MRI protocols with a reference contrast.

OBJECT: In a typical MR session, several contrasts are acquired. Due to the sequential nature of the data acquisition process, the patient may experience some discomfort at some point during the session, and start moving. Hence, it is quite common to have MR sessions where some contrasts are well-resolved, while other contrasts exhibit motion artifacts. Instead of repeating the scans that are corrupted by motion, we introduce a reference-guided retrospective motion correction scheme that takes advantage of the motion-free scans, based on a generalized rigid registration routine. MATERIALS AND METHODS: We focus on various existing clinical 3D brain protocols at 1.5 Tesla MRI based on Cartesian sampling. Controlled experiments with three healthy volunteers and three levels of motion are performed. RESULTS: Radiological inspection confirms that the proposed method consistently ameliorates the corrupted scans. Furthermore, for the set of specific motion tests performed in this study, the quality indexes based on PSNR and SSIM shows only a modest decrease in correction quality as a function of motion complexity. DISCUSSION: While the results on controlled experiments are positive, future applications to patient data will ultimately clarify whether the proposed correction scheme satisfies the radiological requirements.

GPU-accelerated Bloch simulations and MR-STAT reconstructions using the Julia programming language.

PURPOSE: MR-STAT is a relatively new multiparametric quantitative MRI technique in which quantitative paramater maps are obtained by solving a large-scale nonlinear optimization problem. Managing reconstruction times is one of the main challenges of MR-STAT. In this work we leverage GPU hardware to reduce MR-STAT reconstruction times. A highly optimized, GPU-compatible Bloch simulation toolbox is developed as part of this work that can be utilized for other quantitative MRI techniques as well. METHODS: The Julia programming language was used to develop a flexible yet highly performant and GPU-compatible Bloch simulation toolbox called BlochSimulators.jl. The runtime performance of the toolbox is benchmarked against other Bloch simulation toolboxes. Furthermore, a (partially matrix-free) modification of a previously presented (matrix-free) MR-STAT reconstruction algorithm is proposed and implemented using the Julia language on GPU hardware. The proposed algorithm is combined with BlochSimulators.jl and the resulting MR-STAT reconstruction times on GPU hardware are compared to previously presented MR-STAT reconstruction times. RESULTS: The BlochSimulators.jl package demonstrates superior runtime performance on both CPU and GPU hardware when compared to other existing Bloch simulation toolboxes. The GPU-accelerated partially matrix-free MR-STAT reconstruction algorithm, which relies on BlochSimulators.jl, allows for reconstructions of 68 seconds per two-dimensional (2D slice). CONCLUSION: By combining the proposed Bloch simulation toolbox and the partially matrix-free reconstruction algorithm, 2D MR-STAT reconstructions can be performed in the order of one minute on a modern GPU card. The Bloch simulation toolbox can be utilized for other quantitative MRI techniques as well, for example for online dictionary generation for MR Fingerprinting.

High SNR full brain relaxometry at 7T by accelerated MR-STAT.

PURPOSE: To demonstrate the feasibility and robustness of the Magnetic Resonance Spin TomogrAphy in Time-domain (MR-STAT) framework for fast, high SNR relaxometry at 7T. METHODS: To deploy MR-STAT on 7T-systems, we designed optimized flip-angles using the BLAKJac-framework that incorporates the SAR-constraints. Transmit RF-inhomogeneities were mitigated by including a measured B 1 + $$ {B}_1^{+} $$ -map in the reconstruction. Experiments were performed on a gel-phantom and on five volunteers to explore the robustness of the sequence and its sensitivity to B 1 + $$ {B}_1^{+} $$ inhomogeneities. The SNR-gain at 7T was explored by comparing phantom and in vivo results to MR-STAT at 3T in terms of SNR-efficiency. RESULTS: The higher SNR at 7T enabled two-fold acceleration with respect to current 2D MR-STAT protocols at lower field strengths. The resulting scan had whole-brain coverage, with 1 x 1 x 3 mm3 resolution (1.5 mm slice-gap) and was acquired within 3 min including the B 1 + $$ {B}_1^{+} $$ -mapping. After B 1 + $$ {B}_1^{+} $$ -correction, the estimated T1 and T2 in a phantom showed a mean relative error of, respectively, 1.7% and 4.4%. In vivo, the estimated T1 and T2 in gray and white matter corresponded to the range of values reported in literature with a variation over the subjects of 1.0%-2.1% (WM-GM) for T1 and 4.3%-5.3% (WM-GM) for T2. We measured a higher SNR-efficiency at 7T (R = 2) than at 3T for both T1 and T2 with, respectively, a 4.1 and 2.3 times increase in SNR-efficiency. CONCLUSION: We presented an accelerated version of MR-STAT tailored to high field (7T) MRI using a low-SAR flip-angle train and showed high quality parameter maps with an increased SNR-efficiency compared to MR-STAT at 3T.

A three-dimensional Magnetic Resonance Spin Tomography in Time-domain protocol for high-resolution multiparametric quantitative magnetic resonance imaging.

Magnetic Resonance Spin TomogrAphy in Time-domain (MR-STAT) is a multiparametric quantitative MR framework, which allows for simultaneously acquiring quantitative tissue parameters such as T1, T2, and proton density from one single short scan. A typical two-dimensional (2D) MR-STAT acquisition uses a gradient-spoiled, gradient-echo sequence with a slowly varying RF flip-angle train and Cartesian readouts, and the quantitative tissue maps are reconstructed by an iterative, model-based optimization algorithm. In this work, we design a three-dimensional (3D) MR-STAT framework based on previous 2D work, in order to achieve better image signal-to-noise ratio, higher though-plane resolution, and better tissue characterization. Specifically, we design a 7-min, high-resolution 3D MR-STAT sequence, and the corresponding two-step reconstruction algorithm for the large-scale dataset. To reduce the long acquisition time, Cartesian undersampling strategies such as SENSE are adopted in our transient-state quantitative framework. To reduce the computational burden, a data-splitting scheme is designed for decoupling the 3D reconstruction problem into independent 2D reconstructions. The proposed 3D framework is validated by numerical simulations, phantom experiments, and in vivo experiments. High-quality knee quantitative maps with 0.8 × 0.8 × 1.5 mm3 resolution and bilateral lower leg maps with 1.6 mm isotropic resolution can be acquired using the proposed 7-min acquisition sequence and the 3-min-per-slice decoupled reconstruction algorithm. The proposed 3D MR-STAT framework could have wide clinical applications in the future.

Water diffusion and T2 quantification in transient-state MRI: the effect of RF pulse sequence.

In quantitative measurement of the T 2 value of tissues, the diffusion of water molecules has been recognized as a confounder. This is most notably so for transient-state quantitative mapping techniques, which allow simultaneous estimation of T 1 and T 2 . In prior work, apparently conflicting conclusions are presented on the level of diffusion-induced bias on the T2 estimate. So far there is a lack of studies on the effect of the RF pulse angle sequence on the level of diffusion-induced bias. In this work, we show that the specific transient-state RF pulse sequence has a large effect on this level of bias. In particular, the bias level is strongly influenced by the mean value of the RF pulse angles. Also, for realistic values of the spoiling gradient area, we infer that the diffusion-induced bias is negligible for non-liquid human tissues; yet, for phantoms, the effect can be substantial (15% of the true T 2 value) for some RF pulse sequences. This should be taken into account in validation procedures.

Generalizable synthetic MRI with physics-informed convolutional networks.

BACKGROUND: Magnetic resonance imaging (MRI) provides state-of-the-art image quality for neuroimaging, consisting of multiple separately acquired contrasts. Synthetic MRI aims to accelerate examinations by synthesizing any desirable contrast from a single acquisition. PURPOSE: We developed a physics-informed deep learning-based method to synthesize multiple brain MRI contrasts from a single 5-min acquisition and investigate its ability to generalize to arbitrary contrasts. METHODS: A dataset of 55 subjects acquired with a clinical MRI protocol and a 5-min transient-state sequence was used. The model, based on a generative adversarial network, maps data acquired from the five-minute scan to "effective" quantitative parameter maps (q*-maps), feeding the generated PD, T1 , and T2 maps into a signal model to synthesize four clinical contrasts (proton density-weighted, T1 -weighted, T2 -weighted, and T2 -weighted fluid-attenuated inversion recovery), from which losses are computed. The synthetic contrasts are compared to an end-to-end deep learning-based method proposed by literature. The generalizability of the proposed method is investigated for five volunteers by synthesizing three contrasts unseen during training and comparing these to ground truth acquisitions via qualitative assessment and contrast-to-noise ratio (CNR) assessment. RESULTS: The physics-informed method matched the quality of the end-to-end method for the four standard contrasts, with structural similarity metrics above 0.75 ± 0.08 (±std), peak signal-to-noise ratios above 22.4 ± 1.9, representing a portion of compact lesions comparable to standard MRI. Additionally, the physics-informed method enabled contrast adjustment, and similar signal contrast and comparable CNRs to the ground truth acquisitions for three sequences unseen during model training. CONCLUSIONS: The study demonstrated the feasibility of physics-informed, deep learning-based synthetic MRI to generate high-quality contrasts and generalize to contrasts beyond the training data. This technology has the potential to accelerate neuroimaging protocols.

Real-time myocardial landmark tracking for MRI-guided cardiac radio-ablation using Gaussian Processes.

Objective.The high speed of cardiorespiratory motion introduces a unique challenge for cardiac stereotactic radio-ablation (STAR) treatments with the MR-linac. Such treatments require tracking myocardial landmarks with a maximum latency of 100 ms, which includes the acquisition of the required data. The aim of this study is to present a new method that allows to track myocardial landmarks from few readouts of MRI data, thereby achieving a latency sufficient for STAR treatments.Approach.We present a tracking framework that requires only few readouts of k-space data as input, which can be acquired at least an order of magnitude faster than MR-images. Combined with the real-time tracking speed of a probabilistic machine learning framework called Gaussian Processes, this allows to track myocardial landmarks with a sufficiently low latency for cardiac STAR guidance, including both the acquisition of required data, and the tracking inference.Main results.The framework is demonstrated in 2D on a motion phantom, andin vivoon volunteers and a ventricular tachycardia (arrhythmia) patient. Moreover, the feasibility of an extension to 3D was demonstrated byin silico3D experiments with a digital motion phantom. The framework was compared with template matching-a reference, image-based, method-and linear regression methods. Results indicate an order of magnitude lower total latency (<10 ms) for the proposed framework in comparison with alternative methods. The root-mean-square-distances and mean end-point-distance with the reference tracking method was less than 0.8 mm for all experiments, showing excellent (sub-voxel) agreement.Significance.The high accuracy in combination with a total latency of less than 10 ms-including data acquisition and processing-make the proposed method a suitable candidate for tracking during STAR treatments. Additionally, the probabilistic nature of the Gaussian Processes also gives access to real-time prediction uncertainties, which could prove useful for real-time quality assurance during treatments.

Gaussian Processes for real-time 3D motion and uncertainty estimation during MR-guided radiotherapy.

Respiratory motion during radiotherapy causes uncertainty in the tumor's location, which is typically addressed by an increased radiation area and a decreased dose. As a result, the treatments' efficacy is reduced. The recently proposed hybrid MR-linac scanner holds the promise to efficiently deal with such respiratory motion through real-time adaptive MR-guided radiotherapy (MRgRT). For MRgRT, motion-fields should be estimated from MR-data and the radiotherapy plan should be adapted in real-time according to the estimated motion-fields. All of this should be performed with a total latency of maximally 200 ms, including data acquisition and reconstruction. A measure of confidence in such estimated motion-fields is highly desirable, for instance to ensure the patient's safety in case of unexpected and undesirable motion. In this work, we propose a framework based on Gaussian Processes to infer 3D motion-fields and uncertainty maps in real-time from only three readouts of MR-data. We demonstrated an inference frame rate up to 69 Hz including data acquisition and reconstruction, thereby exploiting the limited amount of required MR-data. Additionally, we designed a rejection criterion based on the motion-field uncertainty maps to demonstrate the framework's potential for quality assurance. The framework was validated in silico and in vivo on healthy volunteer data (n=5) acquired using an MR-linac, thereby taking into account different breathing patterns and controlled bulk motion. Results indicate end-point-errors with a 75th percentile below 1 mm in silico, and a correct detection of erroneous motion estimates with the rejection criterion. Altogether, the results show the potential of the framework for application in real-time MR-guided radiotherapy with an MR-linac.

Cartesian vs radial MR-STAT: An efficiency and robustness study.

MR Spin TomogrAphy in Time-domain ("MR-STAT") is quantitative MR technique in which multiple quantitative parameters are estimated from a single short scan by solving a large-scale non-linear optimization problem. In this work we extended the MR-STAT framework to non-Cartesian gradient trajectories. Cartesian MR-STAT and radial MR-STAT were compared in terms of time-efficiency and robustness in simulations, gel phantom measurements and in vivo measurements. In simulations, we observed that both Cartesian and radial MR-STAT are highly robust against undersampling. Radial MR-STAT does have a lower spatial encoding power because the outer corners of k-space are never sampled. However, especially in T2, this is compensated by a higher dynamic encoding power that comes from sampling the k-space center with each readout. In gel phantom measurements, Cartesian MR-STAT was observed to be robust against overfitting whereas radial MR-STAT suffered from high-frequency artefacts in the parameter maps at later iterations. These artefacts are hypothesized to be related to hardware imperfections and were (partially) suppressed with image filters. The time-efficiencies were higher for Cartesian MR-STAT in all vials. In-vivo, the radial reconstruction again suffered from overfitting artefacts. The robustness of Cartesian MR-STAT over the entire range of experiments may make it preferable in a clinical setting, despite radial MR-STAT resulting in a higher T1 time-efficiency in white matter.

Synthetic MRI with Magnetic Resonance Spin TomogrAphy in Time-Domain (MR-STAT): Results from a Prospective Cross-Sectional Clinical Trial.

BACKGROUND: Magnetic Resonance Spin TomogrAphy in Time-domain (MR-STAT) can reconstruct whole-brain multi-parametric quantitative maps (eg, T1 , T2 ) from a 5-minute MR acquisition. These quantitative maps can be leveraged for synthetization of clinical image contrasts. PURPOSE: The objective was to assess image quality and overall diagnostic accuracy of synthetic MR-STAT contrasts compared to conventional contrast-weighted images. STUDY TYPE: Prospective cross-sectional clinical trial. POPULATION: Fifty participants with a median age of 45 years (range: 21-79 years) consisting of 10 healthy participants and 40 patients with neurological diseases (brain tumor, epilepsy, multiple sclerosis or stroke). FIELD STRENGTH/SEQUENCE: 3T/Conventional contrast-weighted imaging (T1 /T2 weighted, proton density [PD] weighted, and fluid-attenuated inversion recovery [FLAIR]) and a MR-STAT acquisition (2D Cartesian spoiled gradient echo with varying flip angle preceded by a non-selective inversion pulse). ASSESSMENT: Quantitative T1 , T2 , and PD maps were computed from the MR-STAT acquisition, from which synthetic contrasts were generated. Three neuroradiologists blinded for image type and disease randomly and independently evaluated synthetic and conventional datasets for image quality and diagnostic accuracy, which was assessed by comparison with the clinically confirmed diagnosis. STATISTICAL TESTS: Image quality and consequent acceptability for diagnostic use was assessed with a McNemar's test (one-sided α = 0.025). Wilcoxon signed rank test with a one-sided α = 0.025 and a margin of Δ = 0.5 on the 5-level Likert scale was used to assess non-inferiority. RESULTS: All data sets were similar in acceptability for diagnostic use (≥3 Likert-scale) between techniques (T1 w:P = 0.105, PDw:P = 1.000, FLAIR:P = 0.564). However, only the synthetic MR-STAT T2 weighted images were significantly non-inferior to their conventional counterpart; all other synthetic datasets were inferior (T1 w:P = 0.260, PDw:P = 1.000, FLAIR:P = 1.000). Moreover, true positive/negative rates were similar between techniques (conventional: 88%, MR-STAT: 84%). DATA CONCLUSION: MR-STAT is a quantitative technique that may provide radiologists with clinically useful synthetic contrast images within substantially reduced scan time. EVIDENCE LEVEL: 1 Technical Efficacy: Stage 2.

Efficient performance analysis and optimization of transient-state sequences for multiparametric magnetic resonance imaging.

In transient-state multiparametric MRI sequences such as Magnetic Resonance Spin TomogrAphy in Time-domain (MR-STAT), MR fingerprinting, or hybrid-state imaging, the flip angle pattern of the RF excitation varies over the sequence. This gives considerable freedom to choose an optimal pattern of flip angles. For pragmatic reasons, most optimization methodologies choose for a single-voxel approach (i.e., without taking the spatial encoding scheme into account). Particularly in MR-STAT, the context of spatial encoding is important. In the current study, we present a methodology, called BLock Analysis of a K-space-domain Jacobian (BLAKJac), which is sufficiently fast to optimize a sequence in the context of a predetermined phase-encoding pattern. Based on MR-STAT acquisitions and reconstructions, we show that sequences optimized using BLAKJac are more reliable in terms of actually achieved precision than conventional single-voxel-optimized sequences. In addition, BLAKJac provides analytical tools that give insights into the performance of the sequence in a very limited computation time. Our experiments are based on MR-STAT, but the theory is equally valid for other transient-state multiparametric methods.

⊥-loss: A symmetric loss function for magnetic resonance imaging reconstruction and image registration with deep learning.

Convolutional neural networks (CNNs) are increasingly adopted in medical imaging, e.g., to reconstruct high-quality images from undersampled magnetic resonance imaging (MRI) acquisitions or estimate subject motion during an examination. MRI is naturally acquired in the complex domain C, obtaining magnitude and phase information in k-space. However, CNNs in complex regression tasks are almost exclusively trained to minimize the L2 loss or maximizing the magnitude structural similarity (SSIM), which are possibly not optimal as they do not take full advantage of the magnitude and phase information present in the complex domain. This work identifies that minimizing the L2 loss in the complex field has an asymmetry in the magnitude/phase loss landscape and is biased, underestimating the reconstructed magnitude. To resolve this, we propose a new loss function for regression in the complex domain called ⊥-loss, which adds a novel phase term to established magnitude loss functions, e.g., L2 or SSIM. We show ⊥-loss is symmetric in the magnitude/phase domain and has favourable properties when applied to regression in the complex domain. Specifically, we evaluate the ⊥+ℓ2-loss and ⊥+SSIM-loss for complex undersampled MR image reconstruction tasks and MR image registration tasks. We show that training a model to minimize the ⊥+ℓ2-loss outperforms models trained to minimize the L2 loss and results in similar performance compared to models trained to maximize the magnitude SSIM while offering high-quality phase reconstruction. Moreover, ⊥-loss is defined in Rn, and we apply the loss function to the R2 domain by learning 2D deformation vector fields for image registration. We show that a model trained to minimize the ⊥+ℓ2-loss outperforms models trained to minimize the end-point error loss.

Acceleration Strategies for MR-STAT: Achieving High-Resolution Reconstructions on a Desktop PC Within 3 Minutes.

MR-STAT is an emerging quantitative magnetic resonance imaging technique which aims at obtaining multi-parametric tissue parameter maps from single short scans. It describes the relationship between the spatial-domain tissue parameters and the time-domain measured signal by using a comprehensive, volumetric forward model. The MR-STAT reconstruction solves a large-scale nonlinear problem, thus is very computationally challenging. In previous work, MR-STAT reconstruction using Cartesian readout data was accelerated by approximating the Hessian matrix with sparse, banded blocks, and can be done on high performance CPU clusters with tens of minutes. In the current work, we propose an accelerated Cartesian MR-STAT algorithm incorporating two different strategies: firstly, a neural network is trained as a fast surrogate to learn the magnetization signal not only in the full time-domain but also in the compressed low-rank domain; secondly, based on the surrogate model, the Cartesian MR-STAT problem is re-formulated and split into smaller sub-problems by the alternating direction method of multipliers. The proposed method substantially reduces the computational requirements for runtime and memory. Simulated and in-vivo balanced MR-STAT experiments show similar reconstruction results using the proposed algorithm compared to the previous sparse Hessian method, and the reconstruction times are at least 40 times shorter. Incorporating sensitivity encoding and regularization terms is straightforward, and allows for better image quality with a negligible increase in reconstruction time. The proposed algorithm could reconstruct both balanced and gradient-spoiled in-vivo data within 3 minutes on a desktop PC, and could thereby facilitate the translation of MR-STAT in clinical settings.

Real-Time Non-Rigid 3D Respiratory Motion Estimation for MR-Guided Radiotherapy Using MR-MOTUS.

The MR-Linac is a combination of an MR-scanner and radiotherapy linear accelerator (Linac) which holds the promise to increase the precision of radiotherapy treatments with MR-guided radiotherapy by monitoring motion during radiotherapy with MRI, and adjusting the radiotherapy plan accordingly. Optimal MR-guidance for respiratory motion during radiotherapy requires MR-based 3D motion estimation with a latency of 200-500 ms. Currently this is still challenging since typical methods rely on MR-images, and are therefore limited by the 3D MR-imaging latency. In this work, we present a method to perform non-rigid 3D respiratory motion estimation with 170 ms latency, including both acquisition and reconstruction. The proposed method called real-time low-rank MR-MOTUS reconstructs motion-fields directly from k -space data, and leverages an explicit low-rank decomposition of motion-fields to split the large scale 3D+t motion-field reconstruction problem posed in our previous work into two parts: (I) a medium-scale offline preparation phase and (II) a small-scale online inference phase which exploits the results of the offline phase for real-time computations. The method was validated on free-breathing data of five volunteers, acquired with a 1.5T Elekta Unity MR-Linac. Results show that the reconstructed 3D motion-field are anatomically plausible, highly correlated with a self-navigation motion surrogate ( R=0.975 ±0.0110 ), and can be reconstructed with a total latency of 170 ms that is sufficient for real-time MR-guided abdominal radiotherapy.

Fast and accurate modeling of transient-state, gradient-spoiled sequences by recurrent neural networks.

Fast and accurate modeling of MR signal responses are typically required for various quantitative MRI applications, such as MR fingerprinting. This work uses a new extended phase graph (EPG)-Bloch model for accurate simulation of transient-state, gradient-spoiled MR sequences, and proposes a recurrent neural network (RNN) as a fast surrogate of the EPG-Bloch model for computing large-scale MR signals and derivatives. The computational efficiency of the RNN model is demonstrated by comparisons with other existing models, showing one to three orders of acceleration compared with the latest GPU-accelerated, open-source EPG package. By using numerical and in vivo brain data, two used cases, namely, MRF dictionary generation and optimal experimental design, are also provided. Results show that the RNN surrogate model can be efficiently used for computing large-scale dictionaries of transient-state signals and derivatives within tens of seconds, resulting in several orders of magnitude acceleration with respect to state-of-the-art implementations. The practical application of transient-state quantitative techniques can therefore be substantially facilitated.

Nonrigid 3D motion estimation at high temporal resolution from prospectively undersampled k-space data using low-rank MR-MOTUS.

PURPOSE: With the recent introduction of the MR-LINAC, an MR-scanner combined with a radiotherapy LINAC, MR-based motion estimation has become of increasing interest to (retrospectively) characterize tumor and organs-at-risk motion during radiotherapy. To this extent, we introduce low-rank MR-MOTUS, a framework to retrospectively reconstruct time-resolved nonrigid 3D+t motion fields from a single low-resolution reference image and prospectively undersampled k-space data acquired during motion. THEORY: Low-rank MR-MOTUS exploits spatiotemporal correlations in internal body motion with a low-rank motion model, and inverts a signal model that relates motion fields directly to a reference image and k-space data. The low-rank model reduces the degrees-of-freedom, memory consumption, and reconstruction times by assuming a factorization of space-time motion fields in spatial and temporal components. METHODS: Low-rank MR-MOTUS was employed to estimate motion in 2D/3D abdominothoracic scans and 3D head scans. Data were acquired using golden-ratio radial readouts. Reconstructed 2D and 3D respiratory motion fields were, respectively, validated against time-resolved and respiratory-resolved image reconstructions, and the head motion against static image reconstructions from fully sampled data acquired right before and right after the motion. RESULTS: Results show that 2D+t respiratory motion can be estimated retrospectively at 40.8 motion fields per second, 3D+t respiratory motion at 7.6 motion fields per second and 3D+t head-neck motion at 9.3 motion fields per second. The validations show good consistency with image reconstructions. CONCLUSIONS: The proposed framework can estimate time-resolved nonrigid 3D motion fields, which allows to characterize drifts and intra and inter-cycle patterns in breathing motion during radiotherapy, and could form the basis for real-time MR-guided radiotherapy.

Real-time assessment of potential peak local specific absorption rate value without phase monitoring: Trigonometric maximization method for worst-case local specific absorption rate determination.

PURPOSE: Multi-transmit MRI systems are typically equipped with dedicated hardware to sample the reflected/lost power in the transmit channels. After extensive calibration, the amplitude and phase of the signal at the feed of each array element can be accurately determined. However, determining the phase is more difficult and monitoring errors can lead to a hazardous peak local specific absorption rate (pSAR10g ) underestimation. For this purpose, methods were published for online maximum potential pSAR10g estimation without relying on phase monitoring, but these methods produce considerable overestimation. We present a trigonometric maximization method to determine the actual worst-case pSAR10g without any overestimation. THEORY AND METHOD: The proposed method takes advantage of the sinusoidal relation between the SAR10g in each voxel and the phases of input signals, to return the maximum achievable SAR10g in a few iterations. The method is applied to determine the worst-case pSAR10g for three multi-transmit array configurations at 7T: (1) body array with eight fractionated dipoles; (2) head array with eight fractionated dipoles; (3) head array with eight rectangular loops. The obtained worst-case pSAR10g values are compared with the pSAR10g values determined with a commonly used method and with a more efficient method based on reference-phases. RESULTS: For each voxel, the maximum achievable SAR10g is determined in less than 0.1 ms. Compared to the reference-phases-based method, the proposed method reduces the mean overestimation of the actual pSAR10g up to 52%, while never underestimating the true pSAR10g . CONCLUSION: The proposed method can widely improve the performance of parallel transmission MRI systems without phase monitoring.

Accelerated MR-STAT Reconstructions Using Sparse Hessian Approximations.

MR-STAT is a quantitative magnetic resonance imaging framework for obtaining multi-parametric quantitative tissue parameter maps using data from single short scans. A large-scale optimization problem is solved in which spatial localization of signal and estimation of tissue parameters are performed simultaneously by directly fitting a Bloch-based volumetric signal model to measured time-domain data. In previous work, a highly parallelized, matrix-free Gauss-Newton reconstruction algorithm was presented that can solve the large-scale optimization problem for high-resolution scans. The main computational bottleneck in this matrix-free method is solving a linear system involving (an approximation to) the Hessian matrix at each iteration. In the current work, we analyze the structure of the Hessian matrix in relation to the dynamics of the spin system and derive conditions under which the (approximate) Hessian admits a sparse structure. In the case of Cartesian sampling patterns with smooth RF trains we demonstrate how exploiting this sparsity can reduce MR-STAT reconstruction times by approximately an order of magnitude.