Local SAR management strategies to use two-channel RF shimming for fetal MRI at 3 T.

PURPOSE: This study evaluates the imaging performance of two-channel RF-shimming for fetal MRI at 3 T using four different local specific absorption rate (SAR) management strategies. METHODS: Due to the ambiguity of safe local SAR levels for fetal MRI, local SAR limits for RF shimming were determined based on either each individual's own SAR levels in standard imaging mode (CP mode) or the maximum SAR level observed across seven pregnant body models in CP mode. Local SAR was constrained either indirectly by further constraining the whole-body SAR (wbSAR) or directly by using subject-specific local SAR models. Each strategy was evaluated by the improvement of the transmit field efficiency (average |B1 + |) and nonuniformity (|B1 + | variation) inside the fetus compared with CP mode for the same wbSAR. RESULTS: Constraining wbSAR when using RF shimming decreases B1 + efficiency inside the fetus compared with CP mode (by 12%-30% on average), making it inefficient for SAR management. Using subject-specific models with SAR limits based on each individual's own CP mode SAR value, B1 + efficiency and nonuniformity are improved on average by 6% and 13% across seven pregnant models. In contrast, using SAR limits based on maximum CP mode SAR values across seven models, B1 + efficiency and nonuniformity are improved by 13% and 25%, compared with the best achievable improvement without SAR constraints: 15% and 26%. CONCLUSION: Two-channel RF-shimming can safely and significantly improve the transmit field inside the fetus when subject-specific models are used with local SAR limits based on maximum CP mode SAR levels in the pregnant population.

Cross-Sectional Observational Study of Typical in utero Fetal Movements Using Machine Learning.

Early variations of fetal movements are the hallmark of a healthy developing central nervous system. However, there are no automatic methods to quantify the complex 3D motion of the developing fetus in utero. The aim of this prospective study was to use machine learning (ML) on in utero MRI to perform quantitative kinematic analysis of fetal limb movement, assessing the impact of maternal, placental, and fetal factors. In this cross-sectional, observational study, we used 76 sets of fetal (24-40 gestational weeks [GW]) blood oxygenation level-dependent (BOLD) MRI scans of 52 women (18-45 years old) during typical pregnancies. Pregnant women were scanned for 5-10 min while breathing room air (21% O2) and for 5-10 min while breathing 100% FiO2 in supine and/or lateral position. BOLD acquisition time was 20 min in total with effective temporal resolution approximately 3 s. To quantify upper and lower limb kinematics, we used a 3D convolutional neural network previously trained to track fetal key points (wrists, elbows, shoulders, ankles, knees, hips) on similar BOLD time series. Tracking was visually assessed, errors were manually corrected, and the absolute movement time (AMT) for each joint was calculated. To identify variables that had a significant association with AMT, we constructed a mixed-model ANOVA with interaction terms. Fetuses showed significantly longer duration of limb movements during maternal hyperoxia. We also found a significant centrifugal increase of AMT across limbs and significantly longer AMT of upper extremities <31 GW and longer AMT of lower extremities >35 GW. In conclusion, using ML we successfully quantified complex 3D fetal limb motion in utero and across gestation, showing maternal factors (hyperoxia) and fetal factors (gestational age, joint) that impact movement. Quantification of fetal motion on MRI is a potential new biomarker of fetal health and neuromuscular development.

Stochastic-offset-enhanced restricted slice excitation and 180° refocusing designs with spatially non-linear ΔB0 shim array fields.

PURPOSE: Developing a general framework with a novel stochastic offset strategy for the design of optimized RF pulses and time-varying spatially non-linear ΔB0 shim array fields for restricted slice excitation and refocusing with refined magnetization profiles within the intervals of the fixed voxels. METHODS: Our framework uses the decomposition property of the Bloch equations to enable joint design of RF-pulses and shim array fields for restricted slice excitation and refocusing with auto-differentiation optimization. Bloch simulations are performed independently on orthogonal basis vectors, Mx, My, and Mz, which enables designs for arbitrary initial magnetizations. Requirements for refocusing pulse designs are derived from the extended phase graph formalism obviating time-consuming sub-voxel isochromatic simulations to model the effects of crusher gradients. To refine resultant slice-profiles because of voxelwise optimization functions, we propose an algorithm that stochastically offsets spatial points at which loss is computed during optimization. RESULTS: We first applied our proposed design framework to standard slice-selective excitation and refocusing pulses in the absence of non-linear ΔB0 shim array fields and compared them against pulses designed with Shinnar-Le Roux algorithm. Next, we demonstrated our technique in a simulated setup of fetal brain imaging in pregnancy for restricted-slice excitation and refocusing of the fetal brain. CONCLUSIONS: Our proposed framework for optimizing RF pulse and time-varying spatially non-linear ΔB0 shim array fields achieve high fidelity restricted-slice excitation and refocusing for fetal MRI, which could enable zoomed fast-spin-echo-MRI and other applications.

Latent signal models: Learning compact representations of signal evolution for improved time-resolved, multi-contrast MRI.

PURPOSE: To improve time-resolved reconstructions by training auto-encoders to learn compact representations of Bloch-simulated signal evolution and inserting the decoder into the forward model. METHODS: Building on model-based nonlinear and linear subspace techniques, we train auto-encoders on dictionaries of simulated signal evolution to learn compact, nonlinear, latent representations. The proposed latent signal model framework inserts the decoder portion of the auto-encoder into the forward model and directly reconstructs the latent representation. Latent signal models essentially serve as a proxy for fast and feasible differentiation through the Bloch equations used to simulate signal. This work performs experiments in the context of T2 -shuffling, gradient echo EPTI, and MPRAGE-shuffling. We compare how efficiently auto-encoders represent signal evolution in comparison to linear subspaces. Simulation and in vivo experiments then evaluate if reducing degrees of freedom by incorporating our proxy for the Bloch equations, the decoder portion of the auto-encoder, into the forward model improves reconstructions in comparison to subspace constraints. RESULTS: An auto-encoder with 1 real latent variable represents single-tissue fast spin echo, EPTI, and MPRAGE signal evolution to within 0.15% normalized RMS error, enabling reconstruction problems with 3 degrees of freedom per voxel (real latent variable + complex scaling) in comparison to linear models with 4-8 degrees of freedom per voxel. In simulated/in vivo T2 -shuffling and in vivo EPTI experiments, the proposed framework achieves consistent quantitative normalized RMS error improvement over linear approaches. From qualitative evaluation, the proposed approach yields images with reduced blurring and noise amplification in MPRAGE-shuffling experiments. CONCLUSION: Directly solving for nonlinear latent representations of signal evolution improves time-resolved MRI reconstructions.

Scan-specific artifact reduction in k-space (SPARK) neural networks synergize with physics-based reconstruction to accelerate MRI.

PURPOSE: To develop a scan-specific model that estimates and corrects k-space errors made when reconstructing accelerated MRI data. METHODS: Scan-specific artifact reduction in k-space (SPARK) trains a convolutional-neural-network to estimate and correct k-space errors made by an input reconstruction technique by back-propagating from the mean-squared-error loss between an auto-calibration signal (ACS) and the input technique's reconstructed ACS. First, SPARK is applied to generalized autocalibrating partially parallel acquisitions (GRAPPA) and demonstrates improved robustness over other scan-specific models, such as robust artificial-neural-networks for k-space interpolation (RAKI) and residual-RAKI. Subsequent experiments demonstrate that SPARK synergizes with residual-RAKI to improve reconstruction performance. SPARK also improves reconstruction quality when applied to advanced acquisition and reconstruction techniques like 2D virtual coil (VC-) GRAPPA, 2D LORAKS, 3D GRAPPA without an integrated ACS region, and 2D/3D wave-encoded imaging. RESULTS: SPARK yields SSIM improvement and 1.5 - 2× root mean squared error (RMSE) reduction when applied to GRAPPA and improves robustness to ACS size for various acceleration rates in comparison to other scan-specific techniques. When applied to advanced reconstruction techniques such as residual-RAKI, 2D VC-GRAPPA and LORAKS, SPARK achieves up to 20% RMSE improvement. SPARK with 3D GRAPPA also improves RMSE performance by ~2×, SSIM performance, and perceived image quality without a fully sampled ACS region. Finally, SPARK synergizes with non-Cartesian, 2D and 3D wave-encoding imaging by reducing RMSE between 20% and 25% and providing qualitative improvements. CONCLUSION: SPARK synergizes with physics-based acquisition and reconstruction techniques to improve accelerated MRI by training scan-specific models to estimate and correct reconstruction errors in k-space.

Selective RF excitation designs enabled by time-varying spatially non-linear ΔB0 fields with applications in fetal MRI.

PURPOSE: To demonstrate, through numerical simulations, novel designs of spatially selective radiofrequency (RF) excitations of the fetal brain by both a restricted 2D slice and 3D inner-volume selection. These designs exploit a single-channel RF pulse, conventional gradient fields, and the spatially non-linear ΔB0 fields of a multi-coil shim array, using an auto-differentiation optimization algorithm. METHODS: The design algorithm jointly optimizes the RF pulse and the time-varying ΔB0 fields, which is produced by a 64-channel multi-coil ΔB0 body array to augment the RF and the linear gradient fields, using an auto-differentiation approach. Two design targets were specified, one a 4-mm thick slice with a limited in-slice extent in one dimension ("restricted slice"), and the other a 3D inner-volume selection encompassing the fetal brain ("inner volume"). The RF duration was limited to 2 ms for the restricted slice excitation and 6 ms for the inner-volume excitation. RESULTS: Excitation profiles were achieved for both the restricted slice excitation task (one-minus-minimum magnitude, 8%) within the region of interest (ROI) and (maximum-minus-zero magnitude, 8%) in the suppressed regions and the fetal brain volume excitation task (13% and 9%, respectively). CONCLUSIONS: The proposed joint design of RF and time-varying, spatially non-linear ΔB0 fields achieves the target excitation profiles with short RF pulse durations and demonstrates the potential to enhance fetal MRI with multi-channel body shim arrays.

Automated detection and reacquisition of motion-degraded images in fetal HASTE imaging at 3 T.

PURPOSE: Fetal brain Magnetic Resonance Imaging suffers from unpredictable and unconstrained fetal motion that causes severe image artifacts even with half-Fourier single-shot fast spin echo (HASTE) readouts. This work presents the implementation of a closed-loop pipeline that automatically detects and reacquires HASTE images that were degraded by fetal motion without any human interaction. METHODS: A convolutional neural network that performs automatic image quality assessment (IQA) was run on an external GPU-equipped computer that was connected to the internal network of the MRI scanner. The modified HASTE pulse sequence sent each image to the external computer, where the IQA convolutional neural network evaluated it, and then the IQA score was sent back to the sequence. At the end of the HASTE stack, the IQA scores from all the slices were sorted, and only slices with the lowest scores (corresponding to the slices with worst image quality) were reacquired. RESULTS: The closed-loop HASTE acquisition framework was tested on 10 pregnant mothers, for a total of 73 acquisitions of our modified HASTE sequence. The IQA convolutional neural network, which was successfully employed by our modified sequence in real time, achieved an accuracy of 85.2% and area under the receiver operator characteristic of 0.899. CONCLUSION: The proposed acquisition/reconstruction pipeline was shown to successfully identify and automatically reacquire only the motion degraded fetal brain HASTE slices in the prescribed stack. This minimizes the overall time spent on HASTE acquisitions by avoiding the need to repeat the entire stack if only few slices in the stack are motion-degraded.

A 31-channel integrated "AC/DC" B0 shim and radiofrequency receive array coil for improved 7T MRI.

PURPOSE: To test an integrated "AC/DC" array approach at 7T, where B0 inhomogeneity poses an obstacle for functional imaging, diffusion-weighted MRI, MR spectroscopy, and other applications. METHODS: A close-fitting 7T 31-channel (31-ch) brain array was constructed and tested using combined Rx and ΔB0 shim channels driven by a set of rapidly switchable current amplifiers. The coil was compared to a shape-matched 31-ch reference receive-only array for RF safety, signal-to-noise ratio (SNR), and inter-element noise correlation. We characterize the coil array's ability to provide global and dynamic (slice-optimized) shimming using ΔB0 field maps and echo planar imaging (EPI) acquisitions. RESULTS: The SNR and average noise correlation were similar to the 31-ch reference array. Global and slice-optimized shimming provide 11% and 40% improvements respectively compared to baseline second-order spherical harmonic shimming. Birdcage transmit coil efficiency was similar for the reference and AC/DC array setups. CONCLUSION: Adding ΔB0 shim capability to a 31-ch 7T receive array can significantly boost 7T brain B0 homogeneity without sacrificing the array's rdiofrequency performance, potentially improving ultra-high field neuroimaging applications that are vulnerable to off-resonance effects.

Improving D-2-hydroxyglutarate MR spectroscopic imaging in mutant isocitrate dehydrogenase glioma patients with multiplexed RF-receive/B0 -shim array coils at 3 T.

MR spectroscopic imaging (MRSI) noninvasively maps the metabolism of human brains. In particular, the imaging of D-2-hydroxyglutarate (2HG) produced by glioma isocitrate dehydrogenase (IDH) mutations has become a key application in neuro-oncology. However, the performance of full field-of-view MRSI is limited by B0 spatial nonuniformity and lipid artifacts from tissues surrounding the brain. Array coils that multiplex RF-receive and B0 -shim electrical currents (AC/DC mixing) over the same conductive loops provide many degrees of freedom to improve B0 uniformity and reduce lipid artifacts. AC/DC coils are highly efficient due to compact design, requiring low shim currents (<2 A) that can be switched fast (0.5 ms) with high interscan reproducibility (10% coefficient of variation for repeat measurements). We measured four tumor patients and five volunteers at 3 T and show that using AC/DC coils in addition to the vendor-provided second-order spherical harmonics shim provides 19% narrower spectral linewidth, 6% higher SNR, and 23% less lipid content for unrestricted field-of-view MRSI, compared with the vendor-provided shim alone. We demonstrate that improvement in MRSI data quality led to 2HG maps with higher contrast-to-noise ratio for tumors that coincide better with the FLAIR-enhancing lesions in mutant IDH glioma patients. Smaller Cramér-Rao lower bounds for 2HG quantification are obtained in tumors by AC/DC shim, corroborating with simulations that predicted improved accuracy and precision for narrower linewidths. AC/DC coils can be used synergistically with optimized acquisition schemes to improve metabolic imaging for precision oncology of glioma patients. Furthermore, this methodology has broad applicability to other neurological disorders and neuroscience.

Safety and imaging performance of two-channel RF shimming for fetal MRI at 3T.

PURPOSE: This study investigates whether two-channel radiofrequency (RF) shimming can improve imaging without increasing specific absorption rate (SAR) for fetal MRI at 3T. METHODS: Transmit field ( B1+ ) average and variation in the fetus was simulated in seven numerical pregnant body models. Safety was quantified by maternal and fetal peak local SAR and fetal average SAR. The shim parameter space was divided into improved B1+ (magnitude and homogeneity) and improved SAR regions, and an overlap where RF shimming improved both classes of metrics compared with birdcage mode was assessed. Additionally, the effect of fetal position, tissue detail, and dielectric properties on transmit field and SAR was studied. RESULTS: A region of subject-specific RF shim parameter space improving both B1+ and SAR metrics was found for five of the seven models. Optimizing only B1+ metrics improved B1+ efficiency across models by 15% on average and 28% for the best-case model. B1+ variation improved by 26% on average and 49% for the best case. However, for these shim settings, fetal SAR increased by up to 106%. The overlap region, where both B1+ and SAR metrics improve, showed an average B1+ efficiency improvement of 6% on average across models and 19% for the best-case model. B1+ variation improved by 13% on average and 40% for the best case. RFS could also decrease maternal/fetal SAR by up to 49%/58%. CONCLUSION: RF shimming can improve imaging compared with birdcage mode without increasing fetal and maternal SAR when a patient-specific SAR model is incorporated into the shimming procedure.

Individualized SAR calculations using computer vision-based MR segmentation and a fast electromagnetic solver.

PURPOSE: We propose a fast, patient-specific workflow for on-line specific absorption rate (SAR) supervision. An individualized electromagnetic model is created while the subject is on the table, followed by rapid SAR estimates for that individual. Our goal is an improved correspondence between the patient and model, reducing reliance on general anatomical body models. METHODS: A 3D fat-water 3T acquisition (~2 minutes) is automatically segmented using a computer vision algorithm (~1 minute) into what we found to be the most important electromagnetic tissue classes: air, bone, fat, and soft tissues. We then compute the individual's EM field exposure and global and local SAR matrices using a fast electromagnetic integral equation solver. We assess the approach in 10 volunteers and compare to the SAR seen in a standard generic body model (Duke). RESULTS: The on-the-table workflow averaged 7'44″. Simulation of the simplified Duke models confirmed that only air, bone, fat, and soft tissue classes are needed to estimate global and local SAR with an error of 6.7% and 2.7%, respectively, compared to the full model. In contrast, our volunteers showed a 16.0% and 20.3% population variability in global and local SAR, respectively, which was mostly underestimated by the Duke model. CONCLUSION: Timely construction and deployment of a patient-specific model is computationally feasible. The benefit of resolving the population heterogeneity compared favorably to the modest modeling error incurred. This suggests that individualized SAR estimates can improve electromagnetic safety in MRI and possibly reduce conservative safety margins that account for patient-model mismatch, especially in non-standard patients.

Joint multi-contrast variational network reconstruction (jVN) with application to rapid 2D and 3D imaging.

PURPOSE: To improve the image quality of highly accelerated multi-channel MRI data by learning a joint variational network that reconstructs multiple clinical contrasts jointly. METHODS: Data from our multi-contrast acquisition were embedded into the variational network architecture where shared anatomical information is exchanged by mixing the input contrasts. Complementary k-space sampling across imaging contrasts and Bunch-Phase/Wave-Encoding were used for data acquisition to improve the reconstruction at high accelerations. At 3T, our joint variational network approach across T1w, T2w and T2-FLAIR-weighted brain scans was tested for retrospective under-sampling at R = 6 (2D) and R = 4 × 4 (3D) acceleration. Prospective acceleration was also performed for 3D data where the combined acquisition time for whole brain coverage at 1 mm isotropic resolution across three contrasts was less than 3 min. RESULTS: Across all test datasets, our joint multi-contrast network better preserved fine anatomical details with reduced image-blurring when compared to the corresponding single-contrast reconstructions. Improvement in image quality was also obtained through complementary k-space sampling and Bunch-Phase/Wave-Encoding where the synergistic combination yielded the overall best performance as evidenced by exemplary slices and quantitative error metrics. CONCLUSION: By leveraging shared anatomical structures across the jointly reconstructed scans, our joint multi-contrast approach learnt more efficient regularizers, which helped to retain natural image appearance and avoid over-smoothing. When synergistically combined with advanced encoding techniques, the performance was further improved, enabling up to R = 16-fold acceleration with good image quality. This should help pave the way to very rapid high-resolution brain exams.

Individual variation in simulated fetal SAR assessed in multiple body models.

PURPOSE: We generate 12 models from 4 pregnant individuals to evaluate individual differences in local specific absorption rate (SAR) for differing body habitus and fetal and maternal positions. METHODS: Structural MR images from 4 pregnant subjects (including supine and left-lateral maternal positions) were manually segmented to create 12 body models by rotating the fetus, modifying the fat content, and altering the maternal arm position in 1 of the subjects. Electromagnetic simulations modeled at 3 Tesla determined the average and peak local SAR in the maternal trunk, fetus, fetal brain, and amniotic fluid. RESULTS: We observed a significant range of fetal and maternal peak local SAR across the models (maternal trunk: 19.14-44.03 watts/kg, fetus: 9.93-18.79 watts/kg, fetal brain 3.36-10.3 watts/kg). We found that maternal body habitus changes introduced a significant variation in the maternal peak local SAR but not the fetal local SAR. However, the maternal position (either rotating the mother to left-lateral position or altering the arm position) introduced changes in fetal peak local SAR (range: 11.9-17.9 watts/kg). Rotating the fetus also introduced variation in the fetal and fetal brain peak local SAR. CONCLUSION: The observed variation in SAR emphasizes the need for more anatomical models to enable better safety management of individuals during fetal MRI, including a wider range of gestational ages.

Nonlinear dipole inversion (NDI) enables robust quantitative susceptibility mapping (QSM).

High-quality Quantitative Susceptibility Mapping (QSM) with Nonlinear Dipole Inversion (NDI) is developed with pre-determined regularization while matching the image quality of state-of-the-art reconstruction techniques and avoiding over-smoothing that these techniques often suffer from. NDI is flexible enough to allow for reconstruction from an arbitrary number of head orientations and outperforms COSMOS even when using as few as 1-direction data. This is made possible by a nonlinear forward-model that uses the magnitude as an effective prior, for which we derived a simple gradient descent update rule. We synergistically combine this physics-model with a Variational Network (VN) to leverage the power of deep learning in the VaNDI algorithm. This technique adopts the simple gradient descent rule from NDI and learns the network parameters during training, hence requires no additional parameter tuning. Further, we evaluate NDI at 7 T using highly accelerated Wave-CAIPI acquisitions at 0.5 mm isotropic resolution and demonstrate high-quality QSM from as few as 2-direction data.

Functional oxygen extraction fraction (OEF) imaging with turbo gradient spin echo QUIXOTIC (Turbo QUIXOTIC).

PURPOSE: QUantitative Imaging of eXtraction of Oxygen and TIssue Consumption (QUIXOTIC) is a recent technique that measures voxel-wise oxygen extraction fraction (OEF) but suffers from long scan times, limiting its application. We implemented multiecho QUIXOTIC dubbed turbo QUIXOTIC (tQUIXOTIC) that reduces scan time eightfold and then applied it in functional MRI. METHODS: tQUIXOTIC utilizes a novel turbo gradient spin echo readout enabling measurement of venular blood transverse relaxation rate in a single tag-control acquisition. Using tQUIXOTIC, we estimated cortical gray matter (GM) OEF, created voxel-by-voxel GM OEF maps, and quantified changes in visual cortex OEF during a blocked design flashing checkerboard visual stimulus. Contamination from cerebrospinal fluid partial volume averaging was estimated and corrected. RESULTS: The average cortical GM OEF was estimated as 0.38 ± 0.06 (n = 8) using a 3.4-min acquisition. The average OEF in the visual cortex was estimated as 0.43 ± 0.04 at baseline and 0.35 ± 0.05 during activation, with an average %ΔOEF of -20%. These values are consistent with those of past studies. CONCLUSION: tQUIXOTIC successfully estimated cortical GM OEF in clinical scan times and detected changes in OEF during blocked design visual stimulation. tQUIXOTIC will be useful to monitor regional OEF clinically and in blocked design or event-related functional MRI experiments. Magn Reson Med 79:2713-2723, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Improved magnetic resonance fingerprinting reconstruction with low-rank and subspace modeling.

PURPOSE: This article introduces a constrained imaging method based on low-rank and subspace modeling to improve the accuracy and speed of MR fingerprinting (MRF). THEORY AND METHODS: A new model-based imaging method is developed for MRF to reconstruct high-quality time-series images and accurate tissue parameter maps (e.g., T1 , T2 , and spin density maps). Specifically, the proposed method exploits low-rank approximations of MRF time-series images, and further enforces temporal subspace constraints to capture magnetization dynamics. This allows the time-series image reconstruction problem to be formulated as a simple linear least-squares problem, which enables efficient computation. After image reconstruction, tissue parameter maps are estimated via dictionary-based pattern matching, as in the conventional approach. RESULTS: The effectiveness of the proposed method was evaluated with in vivo experiments. Compared with the conventional MRF reconstruction, the proposed method reconstructs time-series images with significantly reduced aliasing artifacts and noise contamination. Although the conventional approach exhibits some robustness to these corruptions, the improved time-series image reconstruction in turn provides more accurate tissue parameter maps. The improvement is pronounced especially when the acquisition time becomes short. CONCLUSIONS: The proposed method significantly improves the accuracy of MRF, and also reduces data acquisition time. Magn Reson Med 79:933-942, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Motion-robust sub-millimeter isotropic diffusion imaging through motion corrected generalized slice dithered enhanced resolution (MC-gSlider) acquisition.

PURPOSE: To develop an efficient MR technique for ultra-high resolution diffusion MRI (dMRI) in the presence of motion. METHODS: gSlider is an SNR-efficient high-resolution dMRI acquisition technique. However, subject motion is inevitable during a prolonged scan for high spatial resolution, leading to potential image artifacts and blurring. In this study, an integrated technique termed Motion Corrected gSlider (MC-gSlider) is proposed to obtain high-quality, high-resolution dMRI in the presence of large in-plane and through-plane motion. A motion-aware reconstruction with spatially adaptive regularization is developed to optimize the conditioning of the image reconstruction under difficult through-plane motion cases. In addition, an approach for intra-volume motion estimation and correction is proposed to achieve motion correction at high temporal resolution. RESULTS: Theoretical SNR and resolution analysis validated the efficiency of MC-gSlider with regularization, and aided in selection of reconstruction parameters. Simulations and in vivo experiments further demonstrated the ability of MC-gSlider to mitigate motion artifacts and recover detailed brain structures for dMRI at 860 µm isotropic resolution in the presence of motion with various ranges. CONCLUSION: MC-gSlider provides motion-robust, high-resolution dMRI with a temporal motion correction sensitivity of 2 s, allowing for the recovery of fine detailed brain structures in the presence of large subject movements.

Computation of ultimate SAR amplification factors for radiofrequency hyperthermia in non-uniform body models: impact of frequency and tumour location.

PURPOSE: We introduce a method for calculation of the ultimate specific absorption rate (SAR) amplification factors (uSAF) in non-uniform body models. The uSAF is the greatest possible SAF achievable by any hyperthermia (HT) phased array for a given frequency, body model and target heating volume. METHODS: First, we generate a basis-set of solutions to Maxwell's equations inside the body model. We place a large number of electric and magnetic dipoles around the body model and excite them with random amplitudes and phases. We then compute the electric fields created in the body model by these excitations using an ultra-fast volume integral solver called MARIE. We express the field pattern that maximises the SAF in the target tumour as a linear combination of these basis fields and optimise the combination weights so as to maximise SAF (concave problem). We compute the uSAFs in the Duke body models at 10 frequencies in the 20-900 MHz range and for twelve 3 cm-diameter tumours located at various depths in the head and neck. RESULTS: For both shallow and deep tumours, the frequency yielding the greatest uSAF was ∼900 MHz. Since this is the greatest frequency that we simulated, we hypothesise that the globally optimal frequency is actually greater. CONCLUSIONS: The uSAFs computed in this work are very large (40-100 for shallow tumours and 4-17 for deep tumours), indicating that there is a large room for improvement of the current state-of-the-art head and neck HT devices.

Construction and modeling of a reconfigurable MRI coil for lowering SAR in patients with deep brain stimulation implants.

Post-operative MRI of patients with deep brain simulation (DBS) implants is useful to assess complications and diagnose comorbidities, however more than one third of medical centers do not perform MRIs on this patient population due to stringent safety restrictions and liability risks. A new system of reconfigurable magnetic resonance imaging head coil composed of a rotatable linearly-polarized birdcage transmitter and a close-fitting 32-channel receive array is presented for low-SAR imaging of patients with DBS implants. The novel system works by generating a region with low electric field magnitude and steering it to coincide with the DBS lead trajectory. We demonstrate that the new coil system substantially reduces the SAR amplification around DBS electrodes compared to commercially available circularly polarized coils in a cohort of 9 patient-derived realistic DBS lead trajectories. We also show that the optimal coil configuration can be reliably identified from the image artifact on B1+ field maps. Our preliminary results suggest that such a system may provide a viable solution for high-resolution imaging of DBS patients in the future. More data is needed to quantify safety limits and recommend imaging protocols before the novel coil system can be used on patients with DBS implants.

Mitigation of partial volume effects in susceptibility-based oxygenation measurements by joint utilization of magnitude and phase (JUMP).

PURPOSE: Susceptibility-based blood oxygenation measurements in small vessels of the brain derive from gradient echo (GRE) phase and can provide localized assessment of brain function and pathology. However, when vessel diameter becomes smaller than the acquisition voxel size, partial volume effects compromise these measurements. The purpose of this study was to develop a technique to improve the reliability of vessel oxygenation estimates in the presence of partial volume effects. METHODS: Intravoxel susceptibility variations are present when a vessel and parenchyma experience partial volume effects, modifying the voxel's GRE phase signal and attenuating the GRE magnitude signal. Using joint utilization of magnitude and phase (JUMP), both vessel susceptibility and voxel partial volume fraction can be estimated, providing measurements of venous oxygen saturation ( Yv) in straight, nearly vertical vessels that have improved robustness to partial volume effects. RESULTS: JUMP was demonstrated by estimating vessel Yv in numerical and in vivo experiments. Deviations from ground truth of Yv measurements in vessels tilted up to 30° from B0 were reduced by over 50% when using JUMP compared with phase-only techniques. CONCLUSION: JUMP exploits both magnitude and phase data in GRE imaging to mitigate partial volume effects in estimation of vessel oxygenation. Magn Reson Med 77:1713-1727, 2017. © 2016 International Society for Magnetic Resonance in Medicine.