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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

MRI-based transfer function determination through the transfer matrix by jointly fitting the incident and scattered B1+ field.

PURPOSE: A purely experimental method for MRI-based transfer function (TF) determination is presented. A TF characterizes the potential for radiofrequency heating of a linear implant by relating the incident tangential electric field to a scattered electric field at its tip. We utilize the previously introduced transfer matrix (TM) to determine transfer functions solely from the MR measurable quantities, that is, the B1+ and transceive phase distributions. This technique can extend the current practice of phantom-based TF assessment with dedicated experimental setup toward MR-based methods that have the potential to assess the TF in more realistic situations. THEORY AND METHODS: An analytical description of the B1+ magnitude and transceive phase distribution around a wire-like implant was derived based on the TM. In this model, the background field is described using a superposition of spherical and cylindrical harmonics while the transfer matrix is parameterized using a previously introduced attenuated wave model. This analytical description can be used to estimate the transfer matrix and transfer function based on the measured B1+ distribution. RESULTS: The TF was successfully determined for 2 mock-up implants: a 20-cm bare copper wire and a 20-cm insulated copper wire with 10 mm of insulation stripped at both endings in respectively 4 and 3 different trajectories. The measured TFs show a strong correlation with a reference determined from simulations and between the separate experiments with correlation coefficients above 0.96 between all TFs. Compared to the simulated TF, the maximum deviation in the estimated tip field is 9.4% and 12.2% for the bare and insulated wire, respectively. CONCLUSIONS: A method has been developed to measure the TF of medical implants using MRI experiments. Jointly fitting the incident and scattered B1+ distributions with an analytical description based on the transfer matrix enables accurate determination of the TF of 2 test implants. The presented method no longer needs input from simulated data and can therefore, in principle, be used to measure TF's in test animals or corpses.

Conditional safety margins for less conservative peak local SAR assessment: A probabilistic approach.

PURPOSE: The introduction of a linear safety factor to address peak local specific absorption rate (pSAR10g ) uncertainties (eg, intersubject variation, modeling inaccuracies) bears one considerable drawback: It often results in over-conservative scanning constraints. We present a more efficient approach to define a variable safety margin based on the conditional probability density function of the effectively obtained pSAR10g value, given the estimated pSAR10g value. METHODS: The conditional probability density function can be estimated from previously simulated data. A representative set of true and estimated pSAR10g samples was generated by means of our database of 23 subject-specific models with an 8-fractionated dipole array for prostate imaging at 7 T. The conditional probability density function was calculated for each possible estimated pSAR10g value and used to determine the corresponding safety margin with an arbitrary low probability of underestimation. This approach was applied to five state-of-the-art local SAR estimation methods, namely: (1) using just the generic body model "Duke"; (2) using our model library to assess the maximum pSAR10g value over all models; (3) using the most representative "local SAR model"; (4) using the five most representative local SAR models; and (5) using a recently developed deep learning-based method. RESULTS: Compared with the more conventional safety factor, the conditional safety-margin approach results in lower (up to 30%) mean overestimation for all investigated local SAR estimation methods. CONCLUSION: The proposed probabilistic approach for pSAR10g correction allows more accurate local SAR assessment with much lower overestimation, while a predefined level of underestimation is accepted (eg, 0.1%).

High-resolution in vivo MR-STAT using a matrix-free and parallelized reconstruction algorithm.

MR-STAT is a recently proposed framework that allows the reconstruction of multiple quantitative parameter maps from a single short scan by performing spatial localisation and parameter estimation on the time-domain data simultaneously, without relying on the fast Fourier transform (FFT). To do this at high resolution, specialized algorithms are required to solve the underlying large-scale nonlinear optimisation problem. We propose a matrix-free and parallelized inexact Gauss-Newton based reconstruction algorithm for this purpose. The proposed algorithm is implemented on a high-performance computing cluster and is demonstrated to be able to generate high-resolution (1 mm × 1 mm in-plane resolution) quantitative parameter maps in simulation, phantom, and in vivo brain experiments. Reconstructed T1 and T2 values for the gel phantoms are in agreement with results from gold standard measurements and, for the in vivo experiments, the quantitative values show good agreement with literature values. In all experiments, short pulse sequences with robust Cartesian sampling are used, for which MR fingerprinting reconstructions are shown to fail.

Accuracy and precision of electrical permittivity mapping at 3T: the impact of three B1+ mapping techniques.

PURPOSE: To investigate the sequence-specific impact of B1+ amplitude mapping on the accuracy and precision of permittivity reconstruction at 3T in the pelvic region. METHODS: B1+ maps obtained with actual flip angle imaging (AFI), Bloch-Siegert (BS), and dual refocusing echo acquisition mode (DREAM) sequences, set to a clinically feasible scan time of 5 minutes, were compared in terms of accuracy and precision with electromagnetic and Bloch simulations and MR measurements. Permittivity maps were reconstructed based on these B1+ maps with Helmholtz-based electrical properties tomography. Accuracy and precision in permittivity were assessed. A 2-compartment phantom with properties and size similar to the human pelvis was used for both simulations and measurements. Measurements were also performed on a female volunteer's pelvis. RESULTS: Accuracy was evaluated with noiseless simulations on the phantom. The maximum B1+ bias relative to the true B1+ distribution was 1% for AFI and BS and 6% to 15% for DREAM. This caused an average permittivity bias relative to the true permittivity of 7% to 20% for AFI and BS and 12% to 35% for DREAM. Precision was assessed in MR experiments. The lowest standard deviation in permittivity, found in the phantom for BS, measured 22.4 relative units and corresponded to a standard deviation in B1+ of 0.2% of the B1+ average value. As regards B1+ precision, in vivo and phantom measurements were comparable. CONCLUSIONS: Our simulation framework quantitatively predicts the different impact of B1+ mapping techniques on permittivity reconstruction and shows high sensitivity of permittivity reconstructions to sequence-specific bias and noise perturbation in the B1+ map. These findings are supported by the experimental results.

Whole-brain 3D FLAIR at 7T using direct signal control.

PURPOSE: Image quality obtained for brain imaging at 7T can be hampered by inhomogeneities in the static magnetic field, B0 , and the RF electromagnetic field, B1 . In imaging sequences such as fluid-attenuated inversion recovery (FLAIR), which is used to assess neurological disorders, these inhomogeneities cause spatial variations in signal that can reduce clinical efficacy. In this work, we aim to correct for signal inhomogeneities to ensure whole-brain coverage with 3D FLAIR at 7T. METHODS: The direct signal control (DSC) framework was used to optimize channel weightings applied to the 8 transmit channels used in this work on a pulse-by-pulse basis through the echo train in the FLAIR sequences. 3D FLAIR brain images were acquired on 5 different subjects and compared with imaging using a quadrature-like mode of the transmit array. Precomputed "universal" DSC solutions calculated from a separate set of 5 subjects were also explored. RESULTS: DSC consistently enabled improved imaging across all subjects, with no dropouts in signal seen over the entire brain volume, which contrasted with imaging in quadrature mode. Further, the universal DSC solutions also consistently improved imaging despite not being optimized specifically for the subject being imaged. CONCLUSION: 3D FLAIR brain imaging at 7T is substantially improved using DSC and is able to recover regions of low signal without increasing imaging time or interecho spacing.

Error analysis of helmholtz-based MR-electrical properties tomography.

PURPOSE: MR electrical properties tomography (MR-EPT) aims to measure tissue electrical properties by computing spatial derivatives of measured B1+ data. This computation is very sensitive to spatial fluctuations caused, for example, by noise and Gibbs ringing. In this work, the error arising from the computation of spatial derivatives using finite difference kernels (FD error) has been investigated. In relation to this FD error, it has also been investigated whether mitigation strategies such as Gibbs ringing correction and Gaussian apodization can be beneficial for conductivity reconstructions. METHODS: Conductivity reconstructions were performed on a phantom (by means of simulations and MR measurements at 3T) and on a human brain model. The accuracy was evaluated as a function of image resolution, FD kernel size, k-space windowing, and signal-to-noise ratio. The impact of mitigation strategies was also investigated. RESULTS: The adopted small FD kernel is highly sensitive to spatial fluctuations, whereas the large FD kernel is more noise-robust. However, large FD kernels lead to extended numerical boundary error propagation, which severely hampers the MR-EPT reconstruction accuracy for highly spatially convoluted tissue structures such as the human brain. Mitigation strategies slightly improve the accuracy of conductivity reconstructions. For the adopted derivative kernels and the investigated scenario, MR-EPT conductivity reconstructions show low accuracy: less than 37% of the voxels have a relative error lower than 30%. CONCLUSION: The numerical error introduced by the computation of spatial derivatives using FD kernels is one of the major causes of limited accuracy in Helmholtz-based MR-EPT reconstructions. Magn Reson Med 80:90-100, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Combining a reduced field of excitation with SENSE-based parallel imaging for maximum imaging efficiency.

PURPOSE: To show that a combination of parallel imaging using sensitivity encoding (SENSE) and inner volume imaging (IVI) combines the known benefits of both techniques. SENSE with a reduced field of excitation (rFOX) is termed rSENSE. THEORY AND METHODS: The noise level in SENSE reconstructions is reduced by removing voxels from the unfolding process that are rendered silent by using rFOX. The silent voxels need to be identified beforehand, this is done by using rFOX in the coil sensitivity maps. In vivo experiments were performed at 7 Tesla using a 32-channel receive coil. RESULTS: Good image quality could be obtained in vivo with rSENSE at acceleration factors that are higher than could be obtained using SENSE or IVI alone. With rSENSE we were also able to accelerate scans using an rFOX that was purposely designed to be imperfect or incompatible at all with IVI. CONCLUSION: rSENSE has been demonstrated in vivo with two-dimensionally selective radiofrequency pulses. Besides allowing additional scan acceleration, it offers a greater robustness and flexibility than IVI. The proposed method can be used with other field strengths, anatomies and other rFOX techniques. Magn Reson Med 78:88-96, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution Non Commercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.

Optimal control design of turbo spin-echo sequences with applications to parallel-transmit systems.

PURPOSE: The design of turbo spin-echo sequences is modeled as a dynamic optimization problem which includes the case of inhomogeneous transmit radiofrequency fields. This problem is efficiently solved by optimal control techniques making it possible to design patient-specific sequences online. THEORY AND METHODS: The extended phase graph formalism is employed to model the signal evolution. The design problem is cast as an optimal control problem and an efficient numerical procedure for its solution is given. The numerical and experimental tests address standard multiecho sequences and pTx configurations. RESULTS: Standard, analytically derived flip angle trains are recovered by the numerical optimal control approach. New sequences are designed where constraints on radiofrequency total and peak power are included. In the case of parallel transmit application, the method is able to calculate the optimal echo train for two-dimensional and three-dimensional turbo spin echo sequences in the order of 10 s with a single central processing unit (CPU) implementation. The image contrast is maintained through the whole field of view despite inhomogeneities of the radiofrequency fields. CONCLUSION: The optimal control design sheds new light on the sequence design process and makes it possible to design sequences in an online, patient-specific fashion. Magn Reson Med 77:361-373, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine.